欧拉函数

欧拉函数

定义

欧拉函数$\phi (x)$，表示1~x所有数中与x互质的数的个数。

计算公式

欧拉函数有两个神奇的规律：

1.如果x为素数，则$\phi (x)=x-1$

2.如果$gcd(x,y)=1$，则$\phi (x\ast y)=\phi (x)\ast \phi (y)$

证明见Orz zzkksunboy，我的智商已经无法承受这个证明了，自己看。

那么现在要得到欧拉函数的通项公式

对n分解质因数，有$n=p_1^{x_1}\ast p_2^{x_2}\ast \dots \ast p_k^{x_k}$

注意，这k个数都互质，所以$\phi (n)=\phi (p_1^{x_1})\ast \phi (p_2^{x_2})\ast \dots \ast \phi (p_k^{x_k})$

所以问题又变成了求$\phi (p^k)$

而显而易见的是，因为p是素数，所以在1~$p^x$中除了p的倍数外所有数都与$p^x$互质。而有多少个数呢？$\frac{p^k}{p}$个。

所以$\phi (p^k)=p^k-\frac{p^k}{k}=p^k\ast (1-\frac{1}{p})$

所以综上，欧拉函数的通项公式为$\phi (n)=n\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast \dots \dots \ast (1-\frac{1}{p_k})$

求法

对于要求的欧拉函数的次数大致有两种方法。

1.唯一分解法

专门求$\phi (x)$，直接找出所有x的素数。然后通项公式解决。

int getp(int x){
	int sum=x; for (int i=2,S=sqrt(x);i<=S;i++) if (x%i==0) {sum=sum/i*(i-1); while (x%i==0) x/=i;}
	if (x>1) sum=sum/x*(x-1); return sum;
}

2.筛法，如果要求多次欧拉函数，那么可以筛法求出$\phi$(1~n)然后O（1）得到。

这里采用本身效率高的欧拉筛用来求欧拉函数。

void makep(){
	memset(vs,1,sizeof(vs)); vs[0]=vs[1]=0; p[0]=0; phi[1]=1;
	for (int i=2;i<=maxn;i++){
		if (vs[i]) {p[++p[0]]=i; phi[i]=i-1;}
		for (int j=1;j<=p[0]&&i*p[j]<=maxn;j++){
			vs[i*p[j]]=0;
			if (i%p[j]) phi[p[j]*i]=phi[p[j]]*phi[i];
			else {phi[p[j]*i]=p[j]*phi[i]; break;}
		}
	}
}

扩展

同样copy于zzk

1~n与n互质的所有数的和为：
n=1:nφ(n)/2+1/2=1
n>1:nφ(n)/2
其实很好求：如果(x,n)=1，那么(n-x,n)=1。也就是说与n互质的数是成对出现的，那么显然这些数的和为n*φ(n)/2（1比较特殊，需要特判）。

欧拉定理

如果$gcd(a,p)=1$，有$a^{\phi (p)}\equiv 1(modp)$

注意费马小定理：当p为素数时，存在$a^{p-1}\equiv 1(modp)$ 而此时$\phi (p)=p-1$，所以一个有趣的发现：费马小定理其实是欧拉定理的特殊情况