学习笔记：欧拉函数

欧拉函数定义：

简单推理：

当$p$是质数时，显然有$\phi (p)=p-1$。

当$p$，$q$为质数时，有$$\phi (pq)=(p-1)(q-1)=pq(1-\frac{1}{p})(1-\frac{1}{q})$$
首先，对于$1$ ~ $pq$ 这$pq$个数，我们考虑有多少个数与$pq$是不互质的。
有， $p,2p,3p,...,(q-1)p$共$q-1$个 和 $q,2q,3q,...,(p-1)q$共p-1个，最后在加上$pq$这一个。
所以，与$pq$互质的数有$pq-(p-1)-(q-1)-1$个。
$pq-(p-1)-(q-1)-1=pq-p-q+1=(p-1)(q-1)$。

对于$p^k$，其中$p$是质数，有
$$\phi (p^k)=p^k-p^{k-1}=p^k\ast (1-\frac{1}{p})$$
同样的，我们只考虑在这$p^k$个数中有多少个是与$p^k$不互质的。
有$p,2p,3p,...,p^{k-1}p$共$p^{k-1}$个。
所以，与$p^k$互质的数有$p^k-p^{k-1}$个。

同时，欧拉函数还有一个性质

欧拉函数是积性函数,但不是完全积性函数。
当且只当$n$可以分解成两个互质的整数之积，$n=p_1\times p_2$，则$\phi (n)=\phi (p_1{p}_2)=\phi (p_1)\phi (p_2)$

---

有了以上的知识，我们就可以快速的利用公式计算任意一个正整数$n$的欧拉函数。

由算术基本定理可以知道，任意一个正整数$x$都可以唯一的分解成若干个质数的乘积。
即$$x=p_1^{k_1}\ast p_2^{k_2}\ast p_3^{k_3}\ast ...\ast p_n^{k_n}$$

所以，
$$\begin{gathered} \phi (x)=\phi (p_1^{k_1})\ast \phi (p_2^{k_2})\ast \phi (p_3^{k_3})\ast ...\ast \phi (p_n^{k_n}) \\ =p_1^{k_1}\ast (1-\frac{1}{p_1})\ast ...\ast p_n^{k_n}\ast (1-\frac{1}{p_n}) \\ =x\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})...\ast (1-\frac{1}{p_n}) \end{gathered}$$

---

单个欧拉函数求法：

int euler(int x)
{
    int res = x;
    for(int i = 2; i <= x / i; i ++)
        if(x % i == 0)
        {
            res = res / i * (i - 1);
            while(x % i == 0)   x /= i;
        }
    if(x > 1)   res = res / x * (x - 1);
    return res;
}

区间欧拉函数筛法：
（求1~n的欧拉函数值）
1、埃氏筛思想：

const int maxn=1e6+5;
int phi[maxn]; //phi[i]是i的欧拉函数值
void Euler(int n)
{
    //首先将每个数的欧拉函数值都初始化为它自己 
	for(int i=1;i<=n;i++)	phi[i]=i;
	
	for(int i=2;i<=n;i++)
	{
		//对于每一个质数p，对2p,3p,...的欧拉函数值/p*(p-1) 
		if(phi[i]==i)
		{
			for(int j=i;j<=n;j+=i)
				phi[j]=phi[j]/i*(i-1);
		}
	}
}
int main()
{
	int n; cin>>n;
	Euler(n);
	return 0;
}

2、线性筛过程中求欧拉函数：

const int N = 1e6+5;

int prime[N],cnt,phi[N];
bool vis[N];

void get_euler(int n)
{
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            prime[cnt++]=i;
            phi[i]=i-1;
        }
        for(int j=0;prime[j]<=n/i;j++)
        {
            vis[prime[j]*i]=1;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i]=phi[i]*prime[j];
                break;
            }
            phi[prime[j]*i]=phi[i]*(prime[j]-1);
        }
    }
}

int main()
{
    int n; cin>>n;
    get_euler(n);
    return 0;
}