SDOI 2017 数字表格 题解

题目传送门

题目大意： 求 ${\prod }_{i=1}^a{\prod }_{j=1}^{b}p(\gcd (i,j))$，其中 $p$ 是斐波那契数列。

题解

设 $s(d)$ 表示满足 $\gcd (i,j)=d$ 的 $i,j$ 的数量，$F(n)$ 表示满足 $n|\gcd (i,j)$ 的 $i,j$ 的数量，那么有：
$$F(n)=\sum _{n|d}f(d)$$

根据莫比乌斯反演，有：
$$f(d)=\sum _{d|n}F(n)\mu (\frac{n}{d})$$

大力推柿子：
$$\begin{gathered} \prod _{i=1}^a\prod _{j=1}^{b}p(\gcd (i,j)) \\ =\prod _{i=1}^{\min (a,b)}p(i{)}^{f(i)} \\ =\prod _{i=1}^{\min (a,b)}p(i{)}^{{\sum }_{i|d}F(d)\mu (\frac{d}{i})} \end{gathered}$$

显然有$F(n)=\lfloor \frac{a}{n}\rfloor \lfloor \frac{b}{n}\rfloor$，带进去得：
$$=\prod _{i=1}^{\min (a,b)}p(i{)}^{{\sum }_{i|d}\lfloor \frac{a}{d}\rfloor \lfloor \frac{b}{d}\rfloor \mu (\frac{d}{i})}$$

设 $j=\frac{d}{i}$，更换枚举对象，于是有：
$$=\prod _{i=1}^{\min (a,b)}p(i{)}^{{\sum }_{j=1}^{\frac{\min (a,b)}{i}}\mu (j)\lfloor \frac{a}{ij}\rfloor \lfloor \frac{b}{ij}\rfloor }$$

那个 $\sum$ 是可以弄出来的，变成一个 $\prod$
$$=\prod _{i=1}^{\min (a,b)}\prod _{j=1}^{\frac{\min (a,b)}{i}}p(i{)}^{\mu (j)\lfloor \frac{a}{ij}\rfloor \lfloor \frac{b}{ij}\rfloor }$$

根据莫比乌斯反演的一般思想： 两个变量乘在一起是很难看的，比如说里面的 $ij$，于是用一个 $k$ 代替掉，更换枚举对象，有
$$\begin{gathered} =\prod _{k=1}^{\min (a,b)}\prod _{j|k}p(\frac{k}{j}{)}^{\mu (j)\lfloor \frac{a}{k}\rfloor \lfloor \frac{b}{k}\rfloor } \\ =\prod _{k=1}^{\min (a,b)}{\left(\prod _{j|k}p(\frac{k}{j}{)}^{\mu (j)}\right)}^{\lfloor \frac{a}{k}\rfloor \lfloor \frac{b}{k}\rfloor } \end{gathered}$$

设 $S(n)={\prod }_{d|n}p(\frac{n}{d}{)}^{\mu (d)}$，那么柿子变成：
$$=\prod _{k=1}^{\min (a,b)}S(k{)}^{\lfloor \frac{a}{k}\rfloor \lfloor \frac{b}{k}\rfloor }$$

因为 $\mu (d)$ 只有三种取值：$-1,0,1$ ，所以 $S$ 是可以快速地预处理的。

然后除法分块即可。

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 1000010
#define ll long long
#define mod 1000000007

int T,n,m;
int prime[maxn],t=0;
bool v[maxn];
int mu[maxn];
ll f[maxn],S[maxn];
ll ksm(ll x,ll y)
{
	ll re=1,tot=x;
	while(y>0)
	{
		if(y%2==1)re=re*tot%mod;
		tot=tot*tot%mod;
		y/=2;
	}
	return re;
}
ll inv(ll x){return ksm(x,mod-2);}
void work()
{
	mu[1]=1;
	for(int i=2;i<=maxn-10;i++)//线性筛mu
	{
		if(!v[i])prime[++t]=i,mu[i]=-1;
		for(int j=1;j<=t&&i*prime[j]<=maxn-10;j++)
		{
			v[i*prime[j]]=true;
			if(i%prime[j]==0)
			{
				mu[i*prime[j]]=0;
				break;
			}
			mu[i*prime[j]]=-mu[i];
		}
	}
	f[1]=f[2]=1;
	S[0]=S[1]=S[2]=1;
	for(int i=3;i<=maxn-10;i++)
	f[i]=(f[i-1]+f[i-2])%mod,S[i]=1;//预处理斐波那契数列
	
	for(int i=1;i<=maxn-10;i++)
	for(int j=i;j<=maxn-10;j+=i)
	if(mu[j/i]==-1)S[j]=S[j]*inv(f[i])%mod;
	else if(mu[j/i]==1)S[j]=S[j]*f[i]%mod;//处理S
	for(int i=1;i<=maxn-10;i++)//求S的前缀积
	S[i]=S[i]*S[i-1]%mod;
}

int main()
{
	work();
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d %d",&n,&m);
		int l=1,r;
		ll ans=1;
		while(l<=min(n,m))
		{
			r=min(n/(n/l),m/(m/l));
			ans=ans*ksm(S[r]*inv(S[l-1])%mod,(ll)(n/l)*(m/l))%mod;
			l=r+1;
		}
		printf("%lld\n",ans);
	}
}