UVa 11426 - GCD - Extreme (II)

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2421

代码及其注释：

#include <iostream>

#include <cstdio>

#include <string>

#include <cstring>

#include <cmath>

#include <algorithm>

#include <queue>



#define ll long long

#define lint long long

using namespace std;



const int N=4000005;

int phi[N];//phi[k] 表示从1到k 和k互质的数有几个

ll f[N],sum[N];//f[k] 表示从1到k-1 依次求和k的最大公约数 所以最大公约数的总和

void phin(int n=N-1)

{

    memset(phi,0,sizeof(phi));

    phi[1]=1;

    for(int i=2;i<=n;++i)

    if(!phi[i])

    {

        for(int j=i;j<=n;j=j+i)

        {

            if(!phi[j]) phi[j]=j;

            phi[j]=phi[j]/i*(i-1);

        }

    }

}

int main()

{

    //freopen("data.in","r",stdin);

    phin();

    memset(f,0,sizeof(f));

    memset(sum,0,sizeof(sum));

    for(int i=1;i<N;++i)

    for(int j=i+i,l=2;j<N;j=j+i,++l)

    {

        f[j]+=(ll)phi[l]*(ll)i;//这里的phi[l]表示小于j的数中与j的最大公约数是i的数的个数

    }

    for(int i=1;i<N;++i)

    sum[i]=f[i]+sum[i-1];

    int k;

    while(cin>>k)

    {

        if(!k) break;

        cout<<sum[k]<<endl;

    }

    return 0;

}