欧拉函数

//欧拉函数：求不超过n且与n互素的正整数的个数
long long phi(long long n)
{
    long long rea = n;
    for (int i=2; i*i<=n; i++)
    {
        if (n%i==0)
        {
            rea-=rea/i;
            do
            n/=i;
            while (n%i==0);
        }
    }
    if (n>1)
        rea-=rea/n;
    return rea;
}

先把素数筛选出来，方便欧拉函数使用

bool boo[50000];
long long p[20000];
void prim()//线性筛素数，在main函数里面先调用函数prim
{
    memset(boo, 0, sizeof (boo));
    boo[0]=boo[1]=1;
    int k=0;
    for (int i=2; i<50000; i++)
    {
        if (!boo[i]) p[k++]=i;
        for (int j=0; j<k&&i*p[j]<50000; j++)
        {
            boo[i*p[j]]=1;
            if (!(i%p[j])) break;
        }
    }
}
long long phi(long long n)
{
    long long rea = n;
    for (int i=0; p[i]*p[i]<=n; i++)//对于一些不是素数的可以不用遍历
    {
        if (n%p[i]==0)
        {
            rea-=rea/p[i];
            do
            n/=p[i];
            while (n%p[i]==0);
        }
    }
    if (n>1)
        rea-=rea/n;
    return rea;
}