题目描述:戳这里
题意:求
∑
1
b
∑
1
d
G
C
D
(
x
,
y
)
=
k
\displaystyle\sum_1^b\sum_1^dGCD(x,y)=k
1∑b1∑dGCD(x,y)=k
PS:a,c在题面中已经说明了是没用的[汗]
题解:
这题看似好像可以将x和y同时除以一个k,然后用欧拉函数直接搞。但是,因为是莫比乌斯反演的模板题,我们还是思考一下反演的做法吧。
作为一个初学者,表示一脸懵逼[钨丝反演]。。。
我们先化简一下等式:
∑
x
=
1
B
∑
y
=
1
D
G
C
D
(
x
,
y
)
=
k
\sum_{x=1}^B\sum_{y=1}^DGCD(x,y)=k
x=1∑By=1∑DGCD(x,y)=k
∑
x
=
1
B
k
∑
y
=
1
D
k
G
C
D
(
x
,
y
)
=
1
\sum_{x=1}^\frac{B}{k}\sum_{y=1}^\frac{D}{k}GCD(x,y)=1
x=1∑kBy=1∑kDGCD(x,y)=1
关于莫比乌斯函数,我们知道一个性质:
∑
d
∣
n
μ
(
d
)
=
[
n
=
1
]
\sum_{d\mid n}\mu(d)=[n=1]
d∣n∑μ(d)=[n=1]
那么:
∑
d
∣
g
c
d
(
x
,
y
)
μ
(
d
)
=
1
\sum_{d\mid gcd(x,y)}\mu(d)=1
d∣gcd(x,y)∑μ(d)=1
我们把它代到原式:
∑
x
=
1
B
k
∑
y
=
1
D
k
∑
d
∣
g
c
d
(
x
,
y
)
μ
(
d
)
\sum_{x=1}^\frac{B}{k}\sum_{y=1}^\frac{D}{k}\sum_{d\mid gcd(x,y)}\mu(d)
x=1∑kBy=1∑kDd∣gcd(x,y)∑μ(d)
可以“分解”一下d|gcd(x,y)这一条件:
∑
x
=
1
B
k
∑
y
=
1
D
k
∑
d
∣
x
d
∣
y
μ
(
d
)
\sum_{x=1}^\frac{B}{k}\sum_{y=1}^\frac{D}{k}\sum_{d\mid x\\ d\mid y}\mu(d)
x=1∑kBy=1∑kDd∣xd∣y∑μ(d)
首先枚举d,那么可以写成:
∑
d
=
1
m
i
n
(
B
k
,
D
k
)
μ
(
d
)
⌊
B
d
k
⌋
⌊
D
d
k
⌋
\sum_{d=1}^{min(\frac{B}{k},\frac{D}{k})}\mu(d)\lfloor\frac{B}{dk}\rfloor\lfloor\frac{D}{dk}\rfloor
d=1∑min(kB,kD)μ(d)⌊dkB⌋⌊dkD⌋
这样,我们就化简好了。那么我们只要构造出
μ
(
n
)
\mu(n)
μ(n)这个函数,然后暴力枚举一下d就好了。
但是要注意这样其实是会多枚举的,因为题目中认为GCD(5,7)和GCD(7,5)只算一种。
那么我们只要容斥一下就好了。
不难看出,多算的量就是getmu(
B
k
\frac{B}{k}
kB,
B
k
\frac{B}{k}
kB)/2(假设B小于D)。
所以答案就是:getmu(
B
k
\frac{B}{k}
kB,
D
k
\frac{D}{k}
kD) - getmu(
B
k
\frac{B}{k}
kB,
B
k
\frac{B}{k}
kB)/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;
}