BZOJ#4816. [Sdoi2017]数字表格

4816: [Sdoi2017]数字表格

Time Limit: 50 Sec Memory Limit: 128 MB
Submit: 1302 Solved: 652

Description

Doris刚刚学习了fibonacci数列。用f[i]表示数列的第i项，那么

f[0]=0

f[1]=1

f[n]=f[n-1]+f[n-2],n>=2

Doris用老师的超级计算机生成了一个n×m的表格，第i行第j列的格子中的数是f[gcd(i,j)]，其中gcd(i,j)表示i,

j的最大公约数。Doris的表格中共有n×m个数，她想知道这些数的乘积是多少。答案对10^9+7取模。

Input

有多组测试数据。

第一个一个数T，表示数据组数。

接下来T行，每行两个数n,m

T<=1000,1<=n,m<=10^6

Output

输出T行，第i行的数是第i组数据的结果

Sample Input

3
2 3
4 5
6 7

Sample Output

1
6
960

Problem：

Solution：

推导：

套路枚举gcd：

根据公式二：

枚举dx：

暴力预处理出：　　

void getmu()
{mu[1]=1;for(int i=2;i<N;i++){if(!vis[i]) prime[++cnt]=i,mu[i]=-1;for(int j=1;j<=cnt;j++){if(i*prime[j]>=N) break;vis[i*prime[j]]=1;if(i%prime[j]==0) {mu[i*prime[j]]=0;break;}else mu[i*prime[j]]=-mu[i];}}f[1]=1;for(int i=2;i<N;i++) f[i]=(f[i-1]+f[i-2])%mod;for(int i=1;i<N;i++) inver[i]=qpow(mod-2,f[i]);fill(g,g+N,1);for(int i=1;i<N;i++)for(int j=1;i*j<N;j++)if(mu[j]) g[i*j]=g[i*j]*(mu[j]==1?f[i]:inver[i])%mod;for(int i=1;i<N;i++) g[i]=g[i]*g[i-1]%mod;
}

附上代码：

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+12;
const int mod=1e9+7;
int vis[N],mu[N],prime[N],cnt;
long long f[N],g[N],inver[N];long long qpow(long long a,long long b)
{long long ans=1LL;while(a){if(a&1) ans=ans*b%mod;b=b*b%mod;a>>=1;    }return ans;
}
void getmu()
{mu[1]=1;for(int i=2;i<N;i++) {if(!vis[i]) prime[++cnt]=i,mu[i]=-1;for(int j=1;j<=cnt;j++) {if(i*prime[j]>=N) break;vis[i*prime[j]]=1;if(i%prime[j]==0) {mu[i*prime[j]]=0;break;}else mu[i*prime[j]]=-mu[i];}}f[1]=1;for(int i=2;i<N;i++) f[i]=(f[i-1]+f[i-2])%mod;for(int i=1;i<N;i++) inver[i]=qpow(mod-2,f[i]);fill(g,g+N,1);for(int i=1;i<N;i++)for(int j=1;i*j<N;j++) if(mu[j]) g[i*j]=g[i*j]*(mu[j]==1?f[i]:inver[i])%mod;for(int i=1;i<N;i++) g[i]=g[i]*g[i-1]%mod;
}int main()
{freopen("a.in","r",stdin);getmu();int T;scanf("%d",&T);while(T--) {int n,m;scanf("%d%d",&n,&m);if(n>m) swap(n,m);long long ans=1;int pos;for(int i=1;i<=n;i=pos+1){pos=min(n/(n/i),m/(m/i));ans=ans*qpow((long long)(n/i)*(m/i),g[pos]*qpow(mod-2,g[i-1])%mod)%mod;}printf("%lld\n",ans);}    return 0;
}

转载于:https://www.cnblogs.com/Heey/p/9100705.html

BZOJ#4816. [Sdoi2017]数字表格相关推荐

BZOJ 4816 [Sdoi2017]数字表格
题目链接 https://lydsy.com/JudgeOnline/problem.php?id=4816 题解反演 ∏T=1min⁡(n,m)(∏d∣Tfib(d)μ(d))⌊n/d⌋⌊m/d⌋ ...
BZOJ 4816[SDOI2017]数字表格（莫比乌斯反演）
题目链接 $Description$ 用$f_i$表示$fibonacci$数列第$i$项,求$\prod_{i=1}^{n}\prod_{j=1}^{m}f[gcd(i,j)]$ ...
SDOI2017 数字表格
SDOI2017 数字表格题意: 题目传送门题解: 答案的式子大致是这样的: \[\prod_{i = 1} ^ n \prod_{j = 1} ^ m f_{gcd(i, j)}\] 然后大力反 ...
[SDOI2017]数字表格
题目描述 Doris刚刚学习了fibonacci数列.用$f[i]$ 表示数列的第$i$ 项,那么 $f[0]=0$ ,$f[1]=1$ , $f[n]=f[n-1]+f[n-2],n\geq 2$ ...
【bzoj4816】[Sdoi2017]数字表格
题目链接 Description Doris刚刚学习了fibonacci数列.用f[i]表示数列的第i项,那么 f[0]=0 f[1]=1 f[n]=f[n-1]+f[n-2],n>=2 Dor ...
「luogu3704」[SDOI2017]数字表格
莫比乌斯反演,其中一些处理比较套路 1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 co ...
并不对劲的bzoj4816:loj2000:p3704[SDOI2017]数字表格
题目大意有函数$f(x)$,$f(0)=0,f(1)=1,f(x)=f(x-1)+f(x-2)$ $t$($t\leq1000$)组询问,每次给定$n,m$(\(n,m\leq1 ...
[bzoj4816][Sdoi2017]数字表格（反演+逆元）
(真不想做莫比乌斯了) 首先根据题意写出式子 ∏(i=1~n)∏(j=1~m)f[gcd(i,j)] 很明显的f可以预处理出来,解决根据套路分析,我们可以先枚举gcd(i,j)==d ∏(d=1~n ...
【BZOJ】【2154】Crash的数字表格
莫比乌斯反演 PoPoQQQ讲义第4题题解:http://www.cnblogs.com/jianglangcaijin/archive/2013/11/27/3446169.html 感觉两次sq ...
BZOJ 2154 [国家集训队]Crash的数字表格 / JZPTAB（莫比乌斯反演，经典好题）（Luogu P1829）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P1829 [国家集训队]Crash的数字表格 / JZPTAB(反演,经典好题) Problem S ...