UESTC - 621
f[n]:sigma(i,n),i<n
g[n]:sigmaf[i],i<=n
#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e6+11;
const int Limit = 1e6;
typedef unsigned long long ll;
ll phi[maxn];
void euler(int n){
for(int i = 1; i <= n; i++){
phi[i]=i;
}
for(int i = 2; i <= n; i++){
if(phi[i]==i){
for(int j = i; j <= n; j+=i){
phi[j]=phi[j]/i*(i-1);
}
}
}
}
int n;
ll f[maxn],g[maxn];
bool vis[maxn];
int main(){
euler(maxn-11);
ll ans=0;
for(int i = 1; i+i <= Limit; i++){
for(int j = i<<1; j <= Limit; j+=i){
f[j]+=phi[j/i]*i;
}
}
for(int i = 1; i <= Limit; i++) g[i]=g[i-1]+f[i];
int T; scanf("%d",&T);
while(T--){
scanf("%d",&n);
printf("%llu\n",g[n]);
}
return 0;
}

浙公网安备 33010602011771号