poj 2478

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

suressay

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

阅读量687

点赞数

分类专栏： POJ

本文链接：https://blog.youkuaiyun.com/suressay/article/details/8670366

版权

POJ 专栏收录该内容

67 篇文章

订阅专栏

欧拉函数简单应用。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;
const int nMax=1000005;
int tag[nMax],p[100000],N;
long long F[nMax];
void get_prime()
{
    int cnt = 0; tag[1] = 1;
    for (int i = 2; i<nMax; i++)
    {
        if(!tag[i]) p[cnt++]=i;
        for(int j = 0;j<cnt&&p[j]*i<nMax;j++)
        {
            tag[i*p[j]] = 1;
            if (i % p[j] == 0)
            break;
        }
    }
}
int euler(int n)
{
    int res=n;
    for(int i=0;;i++)
    {
        if(n==1) return res;
        if(!tag[n])
        {
            res=res/n*(n-1);
            return res;
        }
        if(n%p[i]==0)
        {
            while(n%p[i]==0)
                n/=p[i];
            res=res/p[i]*(p[i]-1);
        }
    }
}
int main()
{
//    freopen("test.txt","r",stdin);
    get_prime(); F[2]=1;
    for(int i=3;i<=nMax;i++)
        F[i]=F[i-1]+euler(i);

    while(scanf("%d",&N) && N)
        printf("%lld\n",F[N]);
    return 0;

}