欧拉函数小结

遇到了一个欧拉函数的题，于是就四处搜刮资料，拼凑出这个博客。

---

1.函数内容

定义：在数论，对正整数n，欧拉函数是小于或等于n的数中与n互质的数的个数，记为$\phi (n)$。特别的，$\phi (1)=1$。

通式：$\phi (n)=n{\prod }_{i=1}^k(1-\frac{1}{p_i}),p_1,p_2,...,p_k是n的所有质因数$

2.编程实现

①求$\phi (x)$

int euler(int n)
{
    int t=n, sn=sqrt(n+1);
    for(int i=2; i<=sn; ++i)
    {
        if(n%i==0)
        {
            t=t/i*(i-1);
            while(n%i==0) n/=i;
        }
        if(n==1) break;
    }
    if(n>1) t=t/n*(n-1);
    return t;
}

②求$\phi (1)$~$\phi (n)$

void euler(int n)
{
    e[1]=1;
    for(int i=2; i<=n; ++i)
    {
        if(!e[i])///i is a prime
        {
            for(int j=i; j<=n; j+=i)
            {
                if(!e[j]) e[j]=j;
                e[j]=e[j]/i*(i-1);
            }
        }
    }
    return;
}

③线性求1~n的所有素数和$\phi (1)$~$\phi (n)$
下面用到的一些公式会在“相关性质”中说明。

void Euler(int n)
{
    notprime[1]=true;
    e[1]=1;

    for(int i=2; i<=n; ++i)
    {
        if(!notprime[i])
        {
            prime[pn++]=i;
            e[i]=i-1;
        }
        for(int j=0; j<pn&&prime[j]*1ll*i<=n; ++j)
        {
            notprime[prime[j]*i]=true;
            if(i%prime[j]==0)
            {
                e[prime[j]*i]=e[i]*prime[j];
                break;
            }else
            {
                e[prime[j]*i]=e[i]*(prime[j]-1);
            }
        }
    }
    return;
}

以上代码可用百度百科“欧拉函数”词条中的欧拉函数表验证。

3.相关性质

• $m、n互质\to \phi (m\times n)=\phi (m)\times \phi (n)$
• $imodp=0\to \phi (i\times p)=\phi (i)\times p$
• $imodp=0且\frac{i}{p}modp!=0且p为质数\to \phi (i)=\phi (\frac{i}{p})\times (p-1)$
• 欧拉定理：$x与p互质\to x^{\phi (p)}=1(modp)$
  费马小定理：$p为质数\to x^{p-2}=1(modp)$
• $n为奇数\to \phi (n\times 2)=\phi (n)\times 2$
• 小于p且与p互质的数之和为$\frac{\phi (p)\ast p}{2}$，因为若x与p互质，则p-x也与p互质

2020.1.31

---

资料来源：
[1]百度百科“欧拉函数”
[2]数论基础——欧拉函数.csdn博客