欧拉函数

在说欧拉函数前，先提一下一个数学关系：互质。

    互质，是一种数学关系，它表示两个数除了1以外没有其他公因子。即对于两
个数a, b而言, 若gcd(a, b) = 1,则称a、b互质。

欧拉函数($\phi$)

定义

  $1$~$N$中与$N$互质的数的个数被称为欧拉函数记为$\phi (N)$。

1.计算式

$$\phi (N)=N\times \frac{p_1-1}{p_1}\times \frac{p_2-1}{p_2}\times \cdots \frac{p_m-1}{p_m}$$
该式可化简为
$$\phi (N)=N\times \prod _{质数p|N}(1-\frac{1}{p})$$
其中，$p$为$N$通过算术基本定理得出的质因子。

  事实上，该式用计算机写的时候并不是特别优，所以我们求$\phi (N)$时并不是用该式来求的。

2.性质

1.$\forall n$ > 1,1~$n$中与$n$互质的数的和为 $n\times \phi (n)÷$ 2

  因为$gcd(n,x)=gcd(n,n-x)$，所以与$n$不互质的数是成对出现的且若$x$与$n$不互质，则$n-x$与$n$也不互质。所以与$n$不互质的数的平均数为$n/2$，所以与$n$互质的数的平均数为$n/2$，进而得到上述性质。

2.若$a，b$互质，则$\phi (ab)=\phi (a)\phi (b)$

  该性质还可以表示为欧拉函数为积性函数。

3.若$f$是积性函数，且在算术基本定理中$n={\prod }_{i=1}^m{p}_i^{c_i}$则$f(n)={\prod }_{i=1}^{m}f(p_i^{c_i})$。

  这其实是所有积性函数的共有性质，实际上就是积性函数的定义的另一种表现。

4.设$p$ 为质数，若$p\mid n$且$p^2\mid n$，则$\phi (n)=\phi (n/p)\times p$

  显然，如果我们将$n$分解质因数。那么，对于满足上述条件的$p$而言，$n$与$n/p$一定有相同的质因子，且仅$p$的指数不同。因此，$\phi (n)$与$\phi (n/p)$的商为$p$。所以该性质成立。

5.设$p$为质数，若$p\mid n$ 且$p^2\nmid n$，则$\phi (n)=\phi (n/p)\times (p-1)$

  因为$p\mid n$且$p^2\nmid n$所以$n$与$n/p$互质。因为欧拉函数为积性函数，所以$\phi (n)=\phi (n/p)\times \phi (p)$，又因为$p$是质数，易知$\phi (p)=p-1$。所以该性质成立。

3.求解

  当我们求$\phi (N)$时，如果用计算式求，实际上限制了求的次数。事实上，我们可以利用欧拉函数的几个性质来得出求解方法。

  由欧拉函数性质4和性质5可以想到，如果我们就出小于$N$的所有质数，我们就可以求出$\phi (x)(x\in [1,N))$的值。因此，我们可以利用素数筛，每当找到一个素数，就可以利用欧拉函数的性质4和性质5来求出一系列的$\phi (N)(N=p\times k)$，而某一个$N=p\times k$又正好是欧拉筛计算时的过程量。

代码如下

int f[Size], q[Size] , phi[Size];

void make_q(int n)
{
    int top = 0;
    for(int i = 2 ; i <= n ; i++)
    {
    	if(!f[i])
    	{
    	    f[i] = i;
    	    q[top++] = i;
    	}
    	for(int j = 0 ; q[j] <= f[i] && j < top && q[j] < n / i ; j++)
    	{
    	    f[i * q[j]] = q[j];
    	    if(i % q[j] == 0)//性质4
    	    {
    	    	phi[i * q[j]] = phi[i] * q[j];
    	    }
    	    else//性质5
    	    {
    	    	phi[i *q[j]] = phi[i] * phi[q[j]];
    	    }
    	}
    }
}