uva 11426

#include<iostream>  
#include<cstdio>  
#include<cstring>  
using namespace std;  
const long long maxn=4100000;  
long long phi[maxn*2];  
long long s[maxn],f[maxn];  
void phi_table(int n)  
{  
    for(int i=2;i<=n;i++)phi[i]=0;  
    phi[1]=1;  
    for(int i=2;i<=n;i++)  
    if(!phi[i])  
    for(int j=i;j<=n;j+=i)  
    {  
        if(!phi[j])phi[j]=j;  
        phi[j]=phi[j]/i*(i-1);  
    }  
}  
int main()  
{  
    phi_table(maxn);  
    memset(f,0,sizeof(f));  
    for(int i=1;i<=maxn;i++)  
     for(int j=i+i;j<=maxn;j+=i)//最大公约数不能为本身  
     f[j]+=i*phi[j/i];  
  
    memset(s,0,sizeof(s));  
    s[1]=0;  
    for(int i=2;i<=maxn;i++)  
    s[i]=s[i-1]+f[i];  
  
    int n;  
    while(scanf("%d",&n),n)  
    {  
        printf("%lld\n",s[n]);  
    }  
    return 0;  
}  

 假设a、b(a<b)互质，那么gcd(a,b)=1，这样当i循环到a、j循环到b时就会向结果中+1，而i循环到2*a、j循环到2*b时就会向结果中+2（gcd(2*a,2*b)=2）...循环到k*a和k*b时就会向结果中+k。这样实际上引起结果变化的根源就在于各对互质的数，当i、j循环到他们自身或者自身的倍数时就会引起结果的改变，那么我们不妨先将每对互质的数对结果的贡献值算出来，最后将各对互质的数对结果的贡献累加起来就可以了。


    假设和b互质的数有n个，也就是n对(?,b)（?和b互质），那么在i、j循环到?、b时结果会增加n，循环到(2*?,2*b)时结果就会增加2*n...当i、j循环到k*?、k*b时结果就会增加k*n。那么我们不妨用a[i]记录各种k、b在满足k*b=i时会增加多少结果，那么最后我们要输出的就是a[2]+a[3]+...+a[N]。


    至于找和b互质的数，就是计算b的欧拉函数的值，然后暴力循环k，并修改对应的a[k*b]即可，整体的复杂度是O(N*logN)的。

题解来自于：http://www.cnblogs.com/staginner/archive/2012/10/29/2745135.html