/// Sàng Phi hàm Euler:
#define maxa 1000000
#define long long long
long sum[maxa+7];
long phi[maxa+7];

void sieve(){
    forinc(i,1,maxa) phi[i]=i;

    forinc(i,2,maxa) if(phi[i]==i) for(int j=i;j<=maxa;j+=i) phi[j]-=phi[j]/i;

    forinc(i,1,maxa) for(int j=i;j<=maxa;j+=i) sum[j]+=phi[i]*(j/i);

    forinc(i,1,maxa) sum[i]=sum[i]-i+sum[i-1];
}