莫比乌斯反演初步：HDU 1695 GCD

题目描述：戳这里
题意：求$\sum _1^b\sum _1^{d}GCD(x,y)=k$
PS：a,c在题面中已经说明了是没用的[汗]

题解：
这题看似好像可以将x和y同时除以一个k，然后用欧拉函数直接搞。但是，因为是莫比乌斯反演的模板题，我们还是思考一下反演的做法吧。
作为一个初学者，表示一脸懵逼[钨丝反演]。。。
我们先化简一下等式:
$$\sum _{x=1}^B\sum _{y=1}^{D}GCD(x,y)=k$$
$$\sum _{x=1}^{\frac{B}{k}}\sum _{y=1}^{\frac{D}{k}}GCD(x,y)=1$$
关于莫比乌斯函数，我们知道一个性质：

$$\sum _{d\mid n}\mu (d)=[n=1]$$
那么：
$$\sum _{d\mid gcd(x,y)}\mu (d)=1$$

我们把它代到原式：
$$\sum _{x=1}^{\frac{B}{k}}\sum _{y=1}^{\frac{D}{k}}\sum _{d\mid gcd(x,y)}\mu (d)$$
可以“分解”一下d|gcd(x,y)这一条件：
$$\begin{gathered} \sum _{x=1}^{\frac{B}{k}}\sum _{y=1}^{\frac{D}{k}}\sum _{d\mid x \\ d\mid y}\mu (d) \end{gathered}$$
首先枚举d，那么可以写成：
$$\sum _{d=1}^{min(\frac{B}{k},\frac{D}{k})}\mu (d)\lfloor \frac{B}{dk}\rfloor \lfloor \frac{D}{dk}\rfloor$$

这样，我们就化简好了。那么我们只要构造出$\mu (n)$这个函数，然后暴力枚举一下d就好了。
但是要注意这样其实是会多枚举的，因为题目中认为GCD(5,7)和GCD（7，5）只算一种。
那么我们只要容斥一下就好了。
不难看出，多算的量就是getmu($\frac{B}{k}$,$\frac{B}{k}$)/2（假设B小于D）。
所以答案就是：getmu($\frac{B}{k}$,$\frac{D}{k}$) - getmu($\frac{B}{k}$,$\frac{B}{k}$)/2。

ps：因为这题数据范围比较小，直接对于每个询问暴力枚举d,如果数据大的话，考虑到d的取值在某一些区间里是一样的，就可以用除法分块来解决。

代码如下：

#include<cstdio>
#include<string>
using namespace std;
const int maxn=1e5+5;
int n,m,ans,fa[maxn],a[maxn],siz[maxn];
int getfa(int x){
	if (fa[x]!=x) fa[x]=getfa(fa[x]); return fa[x];
}
int main(){
	scanf("%d %d",&n,&m);
	for (int i=1;i<=n;i++) scanf("%d",&a[i]);
	for (int i=1;i<=n;i++) fa[i]=i;
	for (int i=1;i<=m;i++) {
		int x,y;
		scanf("%d %d",&x,&y);
		int fax=getfa(x),fay=getfa(y);
		if (fax!=fay) fa[fax]=fay;
	}
	for (int i=1;i<=n;i++) siz[getfa(i)]++;
	for (int i=1;i<=n;i++)
	if ((siz[fa[i]]!=1&&fa[i]==fa[a[i]])||(siz[fa[i]]==1&&i==a[i])) ans++;
	printf("%d\n",ans);
	return 0;
}