bzoj2226 [Spoj 5971] LCMSum（欧拉函数）

$$\sum\limits_{i=1}^nlcm(i,n)$$
$$=\sum\limits_{i=1}^n\frac{in}{gcd(i,n)}$$
$$=\sum\limits_{d|n}f(n/d)*d*n/d$$
$$=\sum\limits_{d|n}f(d)*n$$
其中f(n)表示1~n中与n互质的数的和，我们发现这个东西就是 $n*\phi(n)/2$，因为对于n>=3,如果(x,n)=1,则一定有(n-x,n)=1.所以是若干对和为n的数。n=1的话特判是1.
所以我们线性素数筛预处理欧拉函数，然后nlogn预处理答案即可。
复杂度 $O(nlogn+T)$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 1000010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int prime[N],tot=0,phi[N];ll sum[N],ans[N];
bool notprime[N];
inline void init(){
    int n=1000000;notprime[1]=1;phi[1]=1;sum[1]=1;
    for(int i=2;i<=n;++i){
        if(!notprime[i]) prime[++tot]=i,phi[i]=i-1;
        for(int j=1;prime[j]*i<=n;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0){phi[prime[j]*i]=phi[i]*prime[j];break;}
            phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }sum[i]=(ll)i*phi[i]>>1;
    }for(int i=1;i<=n;++i)
        for(int j=i;j<=n;j+=i) ans[j]+=sum[i];
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();init();
    while(tst--){
        int n=read();
        printf("%lld\n",ans[n]*n);
    }return 0;
}