欧拉函数总结学习

你可以看这几篇欧拉函数了解一下，也可以看我的博客
继莫比乌斯反演以后 我又观摩了几个重要的数学知识，虽然自己面对大厦一般的数学还只是浅尝一勺而已。自己也发现不能被一些东西误导，更不能总用多疑的心去面对生活。扯远了，说一下这几个东西吧。

欧拉函数

欧拉函数：指小于等于n的正整数中与n互质的数的个数。

用$\phi (n)$表示欧拉函数。
性质
1 积性函数，这里不给证明。。
2 $\phi (1)=1$；
3 通式：$\phi (x)=x(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_n})$,其中$p_1,p_2$等为x的质因数。
简证：（参考了刘汝佳的紫书）
对于小于n的与n互质的整数个数，可以用容斥原理，为

$$n-\frac{n}{p_1} - \frac{n}{p_2}...+\frac{n}{p_1p_2}+...$$
综述为 $$\phi (n)=\sum_{s\subseteq{\{p_1,p_2,...p_n}\}} (-1)^s \frac{n}{\prod_{p_i\in s}p_i}$$
而对 $\phi (x)=x(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_n})$式子展开有从每一项括号中各取一项加起来的和，就有偶数个元素为正，奇数个为负的法则出现。
所以它与 $\phi (n)=\sum_{s\subseteq{\{p_1,p_2,...p_n}\}} (-1)^s \frac{n}{\prod_{p_i\in s}p_i}$等价。
代码实现：

int phi(int i)
{
    int ans=i;
    int m=(int )sqrt(i+0.5);
    for(int j=2;j<=m;j++)
    {
        if(i%j==0)
        {
            ans=ans/j*(j-1);
            while(i%j==0)
            {
                i/=j;
            }
        }
    }
    if(i>1)
    {
        ans=ans/i*(i-1);//注意ans=ans/i*(i-1)与ans*(1-1/i)等价
    }
    return ans;
}

如果写成数表，则有（不懂的我在注释里标了）

void phi_(int n)
{
    memset(phi,0,sizeof(phi));
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!phi[i])
        //如果phi[i]仍为0 证明i为一质数（因为前面没有数字能为它的约数）
        {
            for(int j=i;j<=n;j+=i)
            //枚举i的倍数并更新它们的值（i一定是j的约数）
            {
                if(!phi[j])
                {
                    phi[j]=j;
                }
                phi[j]=phi[j]/i*(i-1);
            }
        }
    }
}

这样我们可以求出各数的欧拉函数了。
感谢看到这QAQ