欧拉函数(复习自用笔记)

复习资料来源:刘汝佳.<<算法竞赛入门>>系列
欧拉函数描述为$\varphi (n)=小于n且与n互素的数字个数$.
$如何求得欧拉函数呢,首先运用容斥定理求.$
$n=p_1^{a1}\ast p_2^{a2}..p_3^{a3}$
那么先减去包含$p1,p2,p3$倍数的个数.
$n-(n/p1)-(n/p2)..-(n/p_k)$
再加上同时包含$p1p2,p1p3,p2p3..等含两个因子的个数$
$+(n/(p1\ast p2))+(n/(p1\ast p3)..$
然后再减去含有三个因子的个数…
$\varphi (n)={\sum }_{S\subseteq p1,p2..pk}(-1{)}^{|S|}\ast \frac{n}{{\prod }_{pi\subset s}pi}$.
紫书上提供的式子,看起来十分别扭,也难以计算.
考虑重新组合上述式子.实际上类似于二项式计算.可以认为每一个$p$都是由$(1-\frac{1}{p_i})相乘组合而成$
所以有:
$$\varphi (x)=n\ast \prod (1-\frac{1}{p_i})$$
实际代码的时候需要我们把上述括号通分了再计算.$\frac{(pi-1)}{p}$,那么对于单个数字x的$\varphi (x),找到一个x\%p==0的数字p,直接让n包含p的因子全部去除,这样我们保证了,每次找到一个p，必然是唯一分解中对应的p.因为我们每次找到一个p,是把它对应的贡献全部去掉,那么这个p必然是我们要找到的唯一分解的p$.
单个欧拉函数的计算,类似于试除法.
注意,我们是在枚举$\sqrt{n}$范围内的素数,并没有讨论$n$本身是素数的情况,那种情况我们放在函数末尾处理,即$x>1$的情况.
单个$\varphi (x)的计算情况$

int phi(int x){
	int ans = n;//phi(x) = n / {(1-1/p1)*(1-1/p2)...} 
	for(int i=2;i<=(int)(sqrt(x+0.5));i++){
		if(x%i==0){
			ans = ans / i * (i-1);
		}
		while(x%i==0) x/=i;  
	}
	//如果x还不是1,说明x当前就是一个素数.而我们枚举的素数的范围
	//是[1,根号x].没有枚举自己.
	if(x>1) ans = ans/n*(n-1); 
}

与素数筛法类似,我们也能以几乎线性的时间求出欧拉函数.
先给代码再解释原理:

void phi_table(int n){
	for(int i=2;i<=n;i++) phi[i]=0;
	phi[1]=1;
	for(int i=2;i<=n;i++) {
		if(!phi[i]){
			for(int j=i;j<=n;j+=i){
				if(!phi[j]) phi[j] = j;
				phi[j] = phi[j] / i *(i-1);
			}
		}
	}
}

原理:模仿素数筛,我们想每次枚举出一个素数$i$,然后把包含该素数的所有数字对应的贡献都计算完毕,这样时间复杂度就能优化.
我们枚举该素数$i$的倍数$j$,若当前$phi[j]$没有赋值,就让$phi[j]=j$.这相当于单个$phi$函数中的
然后我们知道,对于每一个枚举出来的,含有$i$因子的数字,都要统计该素数$i$的贡献.
故有$phi[j]=phi[j]/i\ast (i-1)$与单个$phi$函数一致.
最后,代码外层循环
保证每个素数$i$不会被重复计算贡献.
该笔记完,接下来施工练习题.