poj 2478 尤拉函式

尤拉函式基礎之基礎題

 

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iomanip>
using namespace std;
#define Max 1000005
int p[Max];
int ph[Max]; //記錄《=i的且與i互質的個數
void prime() //篩選法求素數
{
    int i,j;
    memset(p,0,sizeof(p));
    p[1]=1;
    for(i=2;i<Max;i++)
    {
        if(p[i]==0)
        {
            for(j=i*2;j<Max;j+=i)
                p[j]=1;
        }
    }
}
void euler()
{
    int i,j;
    for(i=0;i<Max;i++)
    {
        ph[i]=i;
    }
    for(i=2;i<Max;i++)
    {
        if(p[i]==0)
        {
            for(j=i;j<Max;j+=i)
                ph[j]=ph[j]/i*(i-1);
        }
    }
}
int main()
{
    prime();
    euler();
    int n;
    while(cin>>n)
    {
        if(n==0)
            break;
        long long ans=0;
        for(int i=2; i<=n; i++)
            ans+=(ph[i]);
        printf("%lld\n",ans);
    }
    return 0;
}