P3455 [POI2007]ZAP-Queries（莫比乌斯反演）

题目传送门

前置： 对于函数$f(x)$和$F(x)$。

如果有$F(n)=\sum _{d|n}f(d)$，那么有$f(n)=\sum _{d|n}\mu (d)F(\frac{n}{d})$
如果有$F(n)=\sum _{n|d}f(d)$，那么有$f(n)=\sum _{n|d}\mu (\frac{d}{n})F(d)$

题意： $T$组询问，给你 $a,b,d$，求出$1\le x\le a,1\le y\le b,gcd(x,y)=d$ 的数量。

思路： 以下推导中，除法均为整数。

$f(d)=\sum _{x=1}^a\sum _{y=1}^b[gcd(x,y)=d]$
$\begin{array}{rl}F(n)& =\sum _{n|k}f(k)\\ & =\sum _{n|k}\sum _{x=1}^a\sum _{y=1}^b[gcd(x,y)=k]\\ & =\sum _{x=1}^a\sum _{y=1}^b[n|gcd(x,y)]\\ & =\sum _{x=1}^a[n|x]\sum _{y=1}^b[n|y]\\ & =\frac{a}{n}\ast \frac{b}{n}\end{array}$

由于有第一步的 $F(n)=\sum _{n|k}f(k)$ ，所以有$f(d)=\sum _{d|k}\mu (\frac{k}{d})F(k)$。

$\begin{array}{rl}f(d)& =\sum _{d|k}\mu (\frac{k}{d})F(k)\\ & =\sum _{d|k}\mu (\frac{k}{d})\ast \frac{a}{k}\ast \frac{b}{k}\end{array}$

设 $t=\frac{k}{d}$，

$\begin{array}{rl}f(d)& =\sum _{t=1}\mu (t)\frac{a}{t\ast d}\ast \frac{b}{t\ast d}\\ & =\sum _{t=1}\mu (t)\frac{\frac{a}{d}}{t}\ast \frac{\frac{b}{d}}{t}\end{array}$

莫比乌斯函数用欧拉筛筛出来，做一下前缀和，接下来就可以用除法分块计算答案了，时间复杂度为$O(T\sqrt{n})$。

代码：

bool ok[N];
int p[N],mu[N];
void euler()
{
	mu[1] = 1;
	for(int i=2; i<N; i++)
	{
		if(!ok[i])
		{
			p[++p[0]] = i;
			mu[i] = -1;
		}
		for(int j=1; j<=p[0]&&p[j]*i<N; j++)
		{
			ok[p[j]*i] = 1;
			if(i%p[j]==0)
				mu[i*p[j]] = 0;
			else mu[i*p[j]] = -mu[i];
		}
	}
	for(int i=2; i<N; i++)
		mu[i] += mu[i-1];
}


signed main()
{
	int t;
	cin>>t;
	euler();
	while(t--)
	{
		int a,b,d;
		cin>>a>>b>>d;
		a/=d;
		b/=d;
		int mi = min(a,b);
		int res = 0;
		for(int i=1,j; i<=mi; i=j+1)
		{
			j = min(a/(a/i),b/(b/i));
			res += (a/i) * (b/i) * (mu[j]-mu[i-1]);
		}
		cout<<res<<endl;
	}



	return 0;
}