欧拉函数详解

欧拉函数$\phi (n)$指是小于或等于n的正整数中与n互质的数的个数$(\phi (1)=1)$。
欧拉函数性质是在整数n上的值等于对n进行素因子分解后，所有的素数幂上的欧拉函数之积。

积性函数

如果一个数论函数$f$有当$gcd(n,m)==1$时

$f(n\ast m)=f(n)f(m)$

就称$f$为积性函数。

当$gcd(n,m)\ne 1$时，也有$f(n\ast m)=f(n)f(m))$时，$f$为完全积性函数。
常见的积性函数有欧拉函数$\phi (n)$，莫比乌斯函数$\mu (n)$，因子和函数$\sigma (n)$等。

欧拉函数

$n={\prod }_{i=1}^r{P}_i^{a_i}（P为质因数）$，$\phi (n)$为积性函数，显然：

$\phi (n)={\prod }_{i=1}^{t}f(P_i^{a_i})$

欧拉函数具有以下性质：

1. 如果$p$为一个素数，$\phi (p)=p-1，$因为$[1，p-1]$之间的素肯定都与互质。反之如果$\phi (p)=p-1，$那么p就为素数。
2. 如果$p$为素数，$a$是一个正整数，那么$\phi (p^a)=p^a-p^{a-1}$
3. 如果$n,m$互质的正整数，那么$\phi (n\ast m)=\phi (n)\phi (m)$
4. 如果$(n,m)=m$，那么$\phi (n\ast m)=m\ast \phi (n)$
5. 如果$n={\prod }_{i=1}^r{p}_i^{a_i}（p为质因数）$，那么：
   $\phi (n)=n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})...(1-\frac{1}{p_r})$
6. 如果$n$为奇数时，$\phi (2n)=\phi (n)$
7. 如果n是一个大于2的正整数，那么$\phi (n)$为偶数。
8. 如果n为正整数，那么:
   $\sum _{d|n}\phi (d)=n$

接下来介绍两个常用的定理：
欧拉定理：对于任何两个互质的正整数$a,m(m>=2)有a^{\phi }(m)\equiv 1(modm)$
费马小定理：当$m$是质数时，$a^m-1\equiv 1(modm)$

根据性质4，我们可以对n进行素因子$o(\sqrt{n})$时间复杂度分解求出$\phi (n)$

int phi(int n)
{
    int rea=n;
    for(int i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            rea=rea-rea/i;
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        rea=rea-rea/n;
    return rea;
}

如果需要用到大量欧拉函数，可以用线性筛$o(n)$求出$n$以内的所有欧拉函数值。

long long prime[maxx+5],phi[maxx+5];
void getphi(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            prime[++prime[0]]=i; phi[i]=i-1;
        }
        for(int j=1;j<=prime[0]&&i*prime[j]<=n;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
}

另外根据性质2，我们也可以直接$o(n)$递推出欧拉函数值。

void getphi(int n)
{
    for(int i=1;i<=maxn;i++) phi[i]=i;
    for(int i=2;i<=maxn;i+=2) phi[i]/=2;
    for(int i=3;i<=maxn;i+=2) 
        if(phi[i]==i)
            for(int j=i;j<=maxn;j+=i)
                phi[j]=phi[j]/i*(i-1);
}