20160321 POJ2478 Farey Sequence（欧拉函数）

题意：

求分母小于等于n的所有无法约分的分数的个数。

思路：

显然有F(2)=phi(2)，F(n)=F(n-1)+phi(n)。

即求2-n欧拉函数的和，模板走你。

#include <iostream>

using namespace std;

const int maxn=1111111;
int minDiv[maxn],phi[maxn];

void genPhi()
{
    for(int i=1;i<maxn;i++)
        minDiv[i]=i;
    for(int i=2;i*i<maxn;i++)
        if(minDiv[i]==i)
           for(int j=i*i;j<maxn;j+=i)
             minDiv[j]=i;
    phi[1]=1;
    for(int i=2;i<maxn;i++)
    {
        phi[i]=phi[i/minDiv[i]];
        if((i/minDiv[i])%minDiv[i]==0)
            phi[i]*=minDiv[i];
        else
            phi[i]*=minDiv[i]-1;
    }
}

long long sum[maxn];

void solve()
{
    sum[2]=phi[2];
    for(int i=3;i<maxn;i++)
        sum[i]=sum[i-1]+phi[i];
}
int main()
{
    genPhi();
    solve();
    int n;
    while(cin>>n)
    {
        if(n==0)break;
        cout<<sum[n]<<endl;
    }
    return 0;
}