欧拉函数公式推导及筛法实现

$x={\sum }_{i=1}{p}_i^{k_i}$，$p_i$ 为 $x$ 的第 $i$ 个质因子，$k_i>0$ 。
则 欧拉函数$\phi (x)=x\ast {\prod }_{i=1}(1-\frac{1}{p_i})$
证明：
欧拉函数即小于 $n$ 且与 $n$ 互质的数的个数，由容斥原理得
$\phi (x)=n-\sum \frac{n}{p_i}+\sum \frac{n}{p_i{p}_j}-\sum \frac{n}{p_i{p}_j{p}_k}+...$，化简即得上式。

筛法求欧拉函数：
线性筛中用每个合数的最小质因子筛掉它，可以通过这点去求出每个数的欧拉函数，设当前合数为 $x$，则
$x=p_1^{k_1}{p}_2^{k_2}{p}_3^{k_3}...$，$p_1$ 为最小质因子，$k>0$
$\phi (x)=x(1-\frac{1}{p_1})(1-\frac{1}{p_2})(1-\frac{1}{p_3})...$
设 $t=x/p_1$
当t中含有质因子 $p_1$ 时，
$\phi (t)=t(1-\frac{1}{p_1})(1-\frac{1}{p_2})(1-\frac{1}{p_3})...$
容易看出$\phi (x)=\phi (t)\ast p_1$
当t中不含有质因子 $p_1$ 时，
$\phi (t)=t(1-\frac{1}{p_2})(1-\frac{1}{p_3})...$
得$\phi (x)=\phi (t)\ast p_1\ast (1-\frac{1}{p_1})=\phi (t)\ast (p_1-1)$
此处的 $t$ 即为以下代码中的 $i$

int p[N], phi[N], cntp;
bool isp[N];
void ol(int mx)
{
    memset(isp, true, sizeof(bool) * (mx + 1));
    isp[0] = isp[1] = false;
    phi[1] = 1;
    for (int i = 2; i <= mx; i++)
    {
        if (isp[i])
        {
            p[++cntp] = i;
            phi[i] = i - 1;
        }
        for (int j = 1; p[j] <= mx / i; j++)
        {
            int x = p[j] * i;
            isp[x] = false;
            if (i % p[j] == 0)
            {
                phi[x] = phi[i] * p[j];
                break;
            }
            else phi[x] = phi[i] * (p[j] - 1);
        }
    }
}

应用：
对于互质的两个正整数 $a,n$ 都有 $a^{\phi (n)}\equiv 1(modn)$
证明：设1到n中与n互质的数为
$b_1,b_2,...,b_{\phi (n)}$
将其中每个数乘以a并取模n后得到新的数组
$a{b}_1\%n,a{b}_2\%n,...,a{b}_{\phi (n)}\%n$

新生成的数组中没有相同的数，这个用反证法证明，
假设有两个数相同即
$a{b}_i\equiv a{b}_j$
得到$b_i\equiv b_j$，与前提矛盾，得证。

对于新的数组，其中没有相同的数，同时所有数还是与n互质的，所以两个数组是相同的集合，就只是排列顺序不同，所以两个数组的乘积同余，即
${\prod }_{i=1}^{\phi (n)}{b}_i\equiv {\prod }_{i=1}^{\phi (n)}a{b}_i(modn)$
${\prod }_{i=1}^{\phi (n)}{b}_i\equiv a^{\phi (n)}{\prod }_{i=1}^{\phi (n)}{b}_i(modn)$
所以： $a^{\phi (n)}\equiv 1(modn)$
证毕。

如有错误欢迎指正。