400题留念
推一推答案就是
∑T=1n⌊nT⌋2f(T)
f(n)=∑d|nϕ(d)μ(nd)
f(n) 是个两个积性函数的卷积,同样也是积性函数,那么可以考虑用线性筛求
- 若n为素数 f(n)=n−2
- 若n是素数p的幂, f(n)=ϕ(n)−ϕ(n/p)
- 否则记录一下n的最小质因数p及次数c,则 f(n)=f(npc)∗f(pnow∗pc)
线性筛一下 再分块
#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 10000010
using namespace std;
typedef long long ll;
ll f[N];
int p[N],low[N];
int n,t;
inline void rea(int &x){
char c=getchar(); x=0;
for(;c>'9'||c<'0';c=getchar());
for(;c>='0'&&c<='9';x=x*10+c-'0',c=getchar());
}
inline void First(){
f[1]=1;
for(int i=2;i<=N-10;i++){
if(!p[i]) p[++*p]=i,f[i]=i-2,low[i]=i;
for(int j=1;j<=*p&&(ll)i*p[j]<=N-10;j++)
if(p[i*p[j]]=1,i%p[j]) f[i*p[j]]=f[i]*f[p[j]],low[i*p[j]]=p[j];
else{
low[i*p[j]]=low[i]*p[j];
if(low[i]==i) f[i*p[j]]=(i*p[j]-i)-(i*p[j]-i)/p[j];
else f[i*p[j]]=f[i/low[i]]*f[low[i]*p[j]];
break;
}
}
for(int i=2;i<=N-10;i++) f[i]+=f[i-1];
}
int main(){
rea(t); First();
while(t--){
rea(n); ll Ans=0;
for(int i=1,j=1;i<=n;i=j+1){
j=n/(n/i);
Ans+=(f[j]-f[i-1])*(n/i)*(n/i);
}
printf("%lld\n",Ans);
}
}