数学：欧拉函数

一个例题是POJ2478，给定一个数n，求在[1,n]这个范围内两两互质的数的个数

对于这个范围内的每一个数，我们只要找到不超过这个数且与这个数互质的数的个数就可以了

其实就是求欧拉函数

一个一个求的方法在欧拉定理部分有介绍，这里我们是要快速求这个欧拉函数

我在这里整理了两份，一份是利用性质求的，一份是用筛法筛的

先看看有啥性质：

phi(p^k)=(p-1)*p^(k-1)

phi(a*b)=phi(a)*phi(b)（a,b互质）

phi(p)=p-1（p为质数）

很神奇，和质数打交道

那么第一种方法可以顺便把质数表打出来？

第二种就是纯筛了，事实证明第一种方法快一丢丢

#include<cstdio>
const int maxn=1000005;
int tot,n;
bool mark[maxn];
int phi[maxn],pri[maxn];
long long ans;
void getphi()
{
    phi[1]=1;
    for(int i=2;i<=1000000;i++)
    {
        if(!mark[i]) {pri[++tot]=i;phi[i]=i-1;}  //如果是质数 
        for(int j=1;j<=tot;j++)
        {
            int x=pri[j];
            if(i*x>1000000) break;
            mark[i*x]=1;
            if(i%x==0) {phi[i*x]=phi[i]*x;break;}
            else phi[i*x]=phi[i]*phi[x]; 
        } 
    }
}
void getphi2()
{
    for(int i=1;i<=1000000;i++) phi[i]=i;
    for(int i=2;i<=1000000;i++)
    {
        if(phi[i]==i)
            for(int j=i;j<=1000000;j+=i)
                phi[j]=phi[j]/i*(i-1);
    }
}
int main()
{
    getphi2();
    while(scanf("%d",&n))
    {
        if(n==0) break;
        ans=0;
        for(int i=2;i<=n;i++) ans+=phi[i];
        printf("%lld\n",ans);
    }
    return 0;
}