Description
求解
Solution
易得原式即
根据 时一定有 ,可将原式化为
上述式子中括号内的两个 对应的项相等,故又可以化为
可以将相同的 合并在一起计算,故只需要统计 的个数。当 时, ,所以 的个数有 个。
故答案为
变换求和顺序,设 ,式子化为
设 ,已知 为积性函数,于是可以 预处理。最后枚举 ,统计贡献即可。
时间复杂度
- 欧拉筛预处理:
- 枚举因数并计算贡献:
代码
#include <cstdio>
const int N=1000000;
int tot,p[N+5],phi[N+5];
long long ans[N+5];
bool flg[N+5];
void solve() {
phi[1]=1;
for(int i=2;i<=N;++i) {
if(!flg[i]) p[++tot]=i,phi[i]=i-1;
for(int j=1;j<=tot&&i*p[j]<=N;++j) {
flg[i*p[j]]=1;
if(i%p[j]==0) {
phi[i*p[j]]=phi[i]*p[j];
break;
}
phi[i*p[j]]=phi[i]*(p[j]-1);
}
}
for(int i=1;i<=N;++i) {
for(int j=1;i*j<=N;++j) {
ans[i*j]+=1LL*j*phi[j]/2;
}
}
for(int i=1;i<=N;++i) ans[i]=1LL*i*ans[i]+i;
}
int main() {
int T,n;
solve();
for(scanf("%d",&T);T;--T) {
scanf("%d",&n);
printf("%lld\n",ans[n]);
}
return 0;
}