洛谷P3704 [SDOI2017]数字表格
题目描述
Doris刚刚学习了fibonacci数列。用f[i]f[i] 表示数列的第ii 项,那么
f[0]=0f[0]=0 ,f[1]=1f[1]=1 ,
f[n]=f[n-1]+f[n-2],ngeq 2f[n]=f[n−1]+f[n−2],n≥2
Doris用老师的超级计算机生成了一个n×mn×m 的表格,
第ii 行第jj 列的格子中的数是f[gcd(i,j)]f[gcd(i,j)] ,其中gcd(i,j)gcd(i,j) 表示i,ji,j 的最大公约数。
Doris的表格中共有n×mn×m 个数,她想知道这些数的乘积是多少。
答案对10^9+7109+7 取模。
输入输出格式
输入格式:
有多组测试数据。
第一个一个数TT ,表示数据组数。
接下来TT 行,每行两个数n,mn,m
输出格式:
输出TT 行,第ii 行的数是第ii 组数据的结果
输入输出样例
输入样例#1:
3
2 3
4 5
6 7
输出样例#1:
1
6
960
说明
对10%10% 的数据,1leq n,mleq 1001≤n,m≤100
对30%30% 的数据,1leq n,mleq 10001≤n,m≤1000
另外存在30%30% 的数据,Tleq 3T≤3
对100%100% 的数据,Tleq1000,1leq n,mleq 10^6T≤1000,1≤n,m≤106
时间限制:5s
内存限制:128MB
一道明年也做不出来的反演题,。
参考了一下洛谷题解
$$prod_{d=1}^nprod_{i=1}^nprod_{j=1}^mif(gcd(i,j)==d)f[gcd(i,j)]$$ $$prod_{d=1}^{n}f[d]^{sum_{i=1}^{n/d}sum_{j=1}^{m/d}[gcd(i,j)==1]}$$
观察上面的$$sum_{i=1}^{n/d}sum_{j=1}^{m/d}[gcd(i,j)==1]$$
这是一个经典反演问题,它等价于
$$sum_{i=1}^{n/d}mu(i)[frac{n}{id}][frac{m}{id}]$$
令$T=id$
$$prod_{T=1}^{n}prod_{d|T}f[d]^{[n/T][m/T]mu(T/d)}$$
$$prod_{T=1}^{n}(prod_{d|T}f[d]^{mu(T/d)})^{[n/T][m/T]}$$
然后对于$[n/T]$和$[m/T]$分块
里面的那一个直接暴力,
不明白为什么最后要%mod-1QWQ...
// luogu-judger-enable-o2
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
const int mod=1e9+7;
const int MAXN=1e6;
inline int read()
{
char c=getchar();int x=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
int fastpow(int a,int p)
{
int base=1;
while(p)
{
if(p&1) base=1LL*base%mod*a%mod;
a=1LL*a%mod*a%mod%mod;
p>>=1;
}
return base%mod;
}
int mu[MAXN+10],sum[MAXN+10],fib[MAXN+10],nfib[MAXN+10],vis[MAXN+10],prime[MAXN],tot;
void Pre()
{
fib[1]=1;
for(int i=2;i<=MAXN;i++) fib[i]=(fib[i-1]%mod+fib[i-2]%mod)%mod;
for(int i=1;i<=MAXN;i++)
nfib[i]=(fastpow(fib[i],mod-2)+mod)%mod;
vis[1]=1;mu[1]=1;
for(int i=2;i<=MAXN;i++)
{
if(!vis[i]) prime[++tot]=i,mu[i]=-1;
for(int j=1;j<=tot&&i*prime[j]<=MAXN;j++)
{
vis[i*prime[j]]=1;
if(i%prime[j]==0) {mu[i*prime[j]]=0;break;}
else mu[i*prime[j]]=-mu[i];
}
}
fill(sum,sum+MAXN,1);
for(int i=1;i<=MAXN;i++)
{
if(mu[i]==0) continue;
for(int j=i;j<=MAXN;j+=i)
{
if(mu[i]==-1) sum[j]=1ll*sum[j]%mod*nfib[j/i]%mod;
if(mu[i]==1) sum[j]=1ll*sum[j]%mod*fib[j/i]%mod;
}
}
for(int i=1;i<=MAXN;i++) sum[i]=1ll*sum[i-1]*sum[i]%mod;
}
int main()
{
#ifdef WIN32
freopen("a.in","r",stdin);
#endif
Pre();
int QWQ=read();
while(QWQ--)
{
int N=read(),M=read();
if(N>M) swap(N,M);
long long int ans=1;
for(int i=1,nxt;i<=N;i=nxt+1)
{
nxt=min(N/(N/i),M/(M/i));
long long int pw=1ll*sum[nxt]*fastpow(sum[i-1],mod-2)%mod;
ans=(1ll*ans*(fastpow(pw,1ll*(N/i)*(M/i)%(mod-1))))%mod;
}
printf("%lldn",(ans+mod)%mod);
}
return 0;
}
- 数据库18456错误怎么解决
- JavaWeb(一)Servlet中的ServletConfig与ServletContext
- 语义化HTML:p、h1-6、q、blockquote、hr、address、code、pre、var、cite、dfn和samp
- Win7系统 IIS 调试ASP(aspcmsgbk25) 错误号:3706 提示 “未找到提供程序 该程序可能未正确安装”解决办法
- 普通文件和数据库存储的对比
- JavaWeb(一)之细说Servlet
- jQuery Tools Scrollable使用的限制
- OOAD-设计模式(一)概述
- 导出表结构(数据字典)
- 简单的Ajax封装
- 概率论10 方差与标准差
- MySQL优化原理
- Python深入06 Python的内存管理
- 微信登陆代码
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 通过Session记录数据
- Layui表格单元个编辑事件
- 程序员画图两大神器,贴心老哥推荐
- 被面试官虐了,索引为何使用B+树,你知道吗
- Supermap GIS地图知识点
- 最全一篇Lombok使用讲解,及原理,真香啊
- 过来人对迷茫的程序员一点建议,3种学习方式的优缺点
- Java实现分页模糊查询
- VB.NET 仿spy++拖拽图标获取窗体句柄
- 线程生命周期,五大状态转换分析
- 推荐 2 款超牛逼、炫酷、实用的Docker管理工具!
- [业界方案]用Jaeger来学习分布式追踪系统Opentracing
- 记一次FullGC的排查经历--从日志到业务代码
- [源码分析] OpenTracing之跟踪Redis
- Python 3.8 新功能大揭秘