题目分析
首先是一些莫比乌斯反演的基本套路:
设
n<m
n
<
m
,首先看要求的式子:
∑i=1n∑j=1mijgcd(i,j)
∑
i
=
1
n
∑
j
=
1
m
i
j
g
c
d
(
i
,
j
)
套路1:优先枚举gcd
∑d=1n∑i=1n∑j=1mijd[gcd(i,j)==d]
∑
d
=
1
n
∑
i
=
1
n
∑
j
=
1
m
i
j
d
[
g
c
d
(
i
,
j
)
==
d
]
套路2: gcd(i,j)=d g c d ( i , j ) = d 等价于 gcd(i/d,j/d)=1 g c d ( i / d , j / d ) = 1
∑d=1nd∑i=1⌊nd⌋∑j=1⌊md⌋ij[gcd(i,j)==1]
∑
d
=
1
n
d
∑
i
=
1
⌊
n
d
⌋
∑
j
=
1
⌊
m
d
⌋
i
j
[
g
c
d
(
i
,
j
)
==
1
]
套路3:莫比乌斯函数的性质
若n不为1,则 ∑d|n=1 ∑ d | n = 1 ,否则原式为0.
∑d=1nd∑i=1⌊nd⌋∑j=1⌊md⌋ij∑t|gcd(i,j)μ(t)
∑
d
=
1
n
d
∑
i
=
1
⌊
n
d
⌋
∑
j
=
1
⌊
m
d
⌋
i
j
∑
t
|
g
c
d
(
i
,
j
)
μ
(
t
)
套路4:优先枚举 t t (同时设)
∑d=1nd∑t=1⌊nd⌋t2μ(t)sum(⌊ndt⌋)sum(⌊mdt⌋)
∑
d
=
1
n
d
∑
t
=
1
⌊
n
d
⌋
t
2
μ
(
t
)
s
u
m
(
⌊
n
d
t
⌋
)
s
u
m
(
⌊
m
d
t
⌋
)
套路5:设 dt=T d t = T ,优先枚举 sum(⌊nT⌋)sum(⌊mT⌋) s u m ( ⌊ n T ⌋ ) s u m ( ⌊ m T ⌋ ) (因为这一块可以数论分块根号时间内求,具体方法语言讲不清,看代码吧)
∑T=1nsum(⌊nT⌋)sum(⌊mT⌋)T∑d|Tdμ(d)
∑
T
=
1
n
s
u
m
(
⌊
n
T
⌋
)
s
u
m
(
⌊
m
T
⌋
)
T
∑
d
|
T
d
μ
(
d
)
套路结束。
现在我们要考虑迅速搞出 T∑d|Tdμ(d) T ∑ d | T d μ ( d ) ,那么我们设 f(T)=∑d|Tdμ(d) f ( T ) = ∑ d | T d μ ( d )
这是一个积性函数,可以线性筛。考虑莫比乌斯反演的线性筛,我们为一个数 T T 加入一个新的质因子,如果 T T 有这个质因子的话,那么对于原来的约数 d d ,这个数如果有质因子,加上这个 p p 之后,并不会产生新的质因子。否则,,所以此时 f(Tp)=f(T) f ( T p ) = f ( T ) 。如果 T T 没有这个因子,那么其每一个约数乘上这个因子后,由于,所以 f(Tp)=f(T)−pf(T)=(1−p)f(T)=f(p)f(T) f ( T p ) = f ( T ) − p f ( T ) = ( 1 − p ) f ( T ) = f ( p ) f ( T )
代码
#include<bits/stdc++.h>
using namespace std;
const int mod=20101009,N=10000005;
int f[N],pri[N],is[N],s[N];
int n,m,tot,ans,div2;
int qm(int x) {return x>=mod?x-mod:x;}
void init() {
f[1]=1;
for(int i=2;i<=n;++i) {
if(!is[i]) pri[++tot]=i,f[i]=qm(mod+1-i);
for(int j=1;j<=tot&&i*pri[j]<=n;++j) {
int k=i*pri[j]; is[k]=1;
if(i%pri[j]) f[k]=1LL*f[i]*f[pri[j]]%mod;
else {f[k]=f[i];break;}
}
}
for(int i=1;i<=n;++i) f[i]=qm(f[i-1]+1LL*f[i]*i%mod);
for(int i=1;i<=m;++i) s[i]=qm(s[i-1]+i);
}
int main()
{
scanf("%d%d",&n,&m);
if(n>m) swap(n,m);
init();
for(int i=1,j;i<=n;i=j+1) {
j=min(n/(n/i),m/(m/i));
ans=qm(ans+1LL*qm(f[j]-f[i-1]+mod)*s[n/i]%mod*s[m/i]%mod);
}
printf("%d\n",ans);
return 0;
}