POJ - 2478 Farey Sequence （欧拉函数）

最新推荐文章于 2020-07-17 07:59:37 发布

zhuxiyulu

最新推荐文章于 2020-07-17 07:59:37 发布

阅读量182

点赞数

分类专栏：数论基础

本文链接：https://blog.csdn.net/zhuxiyulu/article/details/76301900

版权

数论基础专栏收录该内容

50 篇文章 0 订阅

订阅专栏

#include<cstdio>
using namespace std;
/*
欧拉函数

n>=2  a/b(a<=a<b<=n&&gcd(a,b)==1)
F2 = {1/2}
F3 = {1/3, 1/2, 2/3}
F4 = {1/4, 1/3, 1/2, 2/3, 3/4}
F5 = {1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5}
求Fn集合中元素个数

*/
typedef long long LL;
const int maxn=1e6+5;
int n;
int phi[maxn];
LL sum[maxn];

void phi_table(int n)
{
    for(int i=0;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);
            }
        }
    }
}
void init()
{
    sum[1]=0;
    for(int i=2;i<=maxn;i++)
    {
        sum[i]=sum[i-1]+phi[i];
    }
}
int main()
{
    phi_table(maxn);
    init();

    while(~scanf("%d",&n))
    {
       if(n==0) break;
       printf("%lld\n",sum[n]);
    }
    return 0;
}