闲话
一定要掌握各种函数暴力求的方法,我就是不会暴力求 μ \mu μ 然后用杜教筛 MLE 了,只能去网上搜……
分析
首先,转移成方程公式,可以得出:
∑
i
=
1
n
(
a
i
×
t
i
)
=
t
n
+
1
×
m
−
1
\sum_{i=1}^{n}(a_i\times t_i)=t_{n+1}\times m-1
i=1∑n(ai×ti)=tn+1×m−1
通过裴蜀定理,可以得出要满足:
gcd
i
=
1
n
(
a
i
)
=
1
\gcd_{i=1}^{n}(a_i)=1
i=1gcdn(ai)=1
这很满足莫比乌斯反演!
但是
n
n
n 个很难处理,所以先来看两个:
∑
a
1
=
1
m
∑
a
2
=
1
m
[
gcd
(
i
1
,
i
2
)
]
=
1
\sum_{a_1=1}^{m}\sum_{a_2=1}^{m}[\gcd(i1,i2)]=1
a1=1∑ma2=1∑m[gcd(i1,i2)]=1
很明显的变换:
∑
d
∣
m
μ
(
d
)
×
⌊
m
d
⌋
2
\sum_{d\mid m}\mu(d)\times\lfloor\frac{m}{d}\rfloor^2
d∣m∑μ(d)×⌊dm⌋2
那么,之后也会是
d
∣
i
1
,
d
∣
i
2
,
d
∣
i
3
…
d\mid i1,d\mid i2,d\mid i3\dots
d∣i1,d∣i2,d∣i3…,那么,可以大胆的猜出公式:
∑
d
∣
m
μ
(
d
)
×
⌊
m
d
⌋
n
\sum_{d\mid m}\mu(d)\times\lfloor\frac{m}{d}\rfloor^n
d∣m∑μ(d)×⌊dm⌋n
其实也很好理解,推出的公式后面全都是:
∑
a
1
=
1
⌊
m
d
⌋
∑
a
2
=
1
⌊
m
d
⌋
∑
a
3
=
1
⌊
m
d
⌋
∑
a
4
=
1
⌊
m
d
⌋
⋯
∑
a
n
=
1
⌊
m
d
⌋
1
=
⌊
m
d
⌋
n
\sum_{a_1=1}^{\lfloor\frac{m}{d}\rfloor}\sum_{a_2=1}^{\lfloor\frac{m}{d}\rfloor}\sum_{a_3=1}^{\lfloor\frac{m}{d}\rfloor}\sum_{a_4=1}^{\lfloor\frac{m}{d}\rfloor}\dots\sum_{a_n=1}^{\lfloor\frac{m}{d}\rfloor}1=\lfloor\frac{m}{d}\rfloor^n
a1=1∑⌊dm⌋a2=1∑⌊dm⌋a3=1∑⌊dm⌋a4=1∑⌊dm⌋⋯an=1∑⌊dm⌋1=⌊dm⌋n
那么,就到这里了,正解。
代码
#include<cstdio>
#include<vector>
#define MAXN 5000005
using namespace std;
typedef long long ll;
int mo[MAXN];
bool flag[MAXN];
vector<int> prim;
inline ll power(ll x,ll y){
ll res=1;
while(y){
if(y&1){
res*=x;
}
x*=x;
y>>=1;
}
return res;
}
inline void prework(){
mo[1]=1;
for(int i=2;i<MAXN;++i){
if(!flag[i]){
mo[i]=-1;
prim.push_back(i);
}
for(int j=0;j<prim.size()&&1ll*i*prim[j]<MAXN;++j){
flag[i*prim[j]]=true;
if(i%prim[j]){
mo[i*prim[j]]=-mo[i];
}else{
mo[i*prim[j]]=0;
break;
}
}
}
}
inline int mobius(int n){
if(n<MAXN){
return mo[n];
}
int ans=1,mul=1;
for(int i=2;i*i<=n;++i){
if(n%i==0){
ans=-ans;
mul*=i;
}
if(n%(i*i)==0){
return 0;
}
}
if(mul!=n){
ans=-ans;
}
return ans;
}
int main(){
prework();
int n,m;
scanf("%d %d",&n,&m);
ll ans=0;
for(int i=1;i<=m;++i){
if(m%i==0){
ans+=mobius(i)*power(m/i,n);
}
}
printf("%lld",ans);
return 0;
}