20160321 POJ2478 Farey Sequence(欧拉函数)

来源:互联网 发布:手机arp扫描软件 编辑:程序博客网 时间:2024/05/22 17:02

题意:

求分母小于等于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;}


0 0
原创粉丝点击