欧拉函数
求phi(n)
//phi(1)==1
int phi(int k){
int ans=k,a=k;
for(int i=2;i*i<=a;i++)
if(a%i==0){
ans-=ans/i;
while(a%i==0)a/=i;//算数基本定理
}
if(a>1)ans-=ans/a;//存在大于sqrt(k)的质因子
return ans;
}
欧拉线性筛
筛质数
优秀!O(n)时间复杂度
每个合数仅会被它的最小质因数筛去1次
为什么是最小质因数?
假设A=x*y(x是质数,y是合数)
,y=m*n(m是质数,显然m<y)
,那么A= x*m*n= m*Z
。可得,A可以表示为一个更大的合数和一个更小的质数的乘积,所以A一定会被它最小的质因数筛掉。
比如24,它只会被2*12筛掉。
void getPrime(int n){//筛n以内的质数
int cnt=0;
for(int i=2;i<=n;i++){
if(Mark[i]==false)Prime[++cnt]=i;//没被筛过,是质数
for(int j=1;j<=cnt && i*Prime[j]<=n;j++){
Mark[i*Prime[j]]=true;//标记为已经被筛
if(i%Prime[j]==0)break;//保证后面的数不会重复被筛
}
}
}
求phi(1)+phi(2)+…+phi(n)
欧拉函数的性质(p为质数):
1.phi(p)=p-1
;
2.(k%p=0)phi(k*p)=p*phi(k)
;
3.(k%p!=0,即二者互质)phi(k*p)=(p-1)*phi(k)
(欧拉函数是积性函数).
优秀!同样O(n)的时间复杂度
void phi(int n){
int cnt=0;
for(int i=2;i<=n;i++){
if(Mark[i]==false)Prime[++cnt]=i,Phi[i]=i-1;//性质1
for(int j=1;j<=cnt && i*Prime[j]<=n;j++){
Mark[i*Prime[j]]=true;
if(i%Prime[j]==0){//性质2
Phi[i*Prime[j]]=Phi[i]*Prime[j];break;
}else Phi[i*Prime[j]]=Phi[i]*(Prime[j]-1);//性质3
}
}
}