该算法在可在线性时间内筛素数的同时求出所有数的欧拉函数。
需要用到如下性质(p为质数):
1. phi(p)=p-1 因为质数p除了1以外的因数只有p,故1至p的整数只有p与p不互质
2. 如果i mod p = 0, 那么phi(i * p)=p * phi(i) 证明如下
(上述证明存在bug。。感谢@PrimaryOIer指教)
上面的过程证明了从区间[1,i]->[i+1,i+i],若整数n不与i互质,n+i依然与i不互质。下面给出另一个证明:若整数n与i互质,n+i与i依然互质
3.若i mod p ≠0, 那么phi(i * p)=phi(i) * (p-1)
i mod p 不为0且p为质数, 所以i与p互质, 那么根据欧拉函数的积性phi(i * p)=phi(i) * phi(p) 其中phi(p)=p-1即第一条性质
PS:我的代码如下:
void get_phi_and_primes(int n){
for(int i = 1;i <= n;i++) phi[i] = 0;
phi[1] = 1;
for(int i = 2;i <= n;i++){
if(!phi[i]) {
primes[++pcnt] = i;phi[i] = i - 1;
}
for(int j = 1;j <= pcnt && primes[j] * i <= n;j++){
int t = primes[j];
if(i % primes[j] == 0){
phi[i * t] = phi[i] * t;
break;
}
else phi[i * t] = phi[i] * (t - 1);
}
}
}