预处理欧拉函数值

#include<bits/stdc++.h>
using namespace std;
const int maxn=2333333;
bitset<maxn>notprime;
int phi[maxn],prime[maxn],cnt=0;
void pre(){
	phi[1]=1; 
	for(int i=2;i<=maxn;i++){
        if(!notprime[i]){
            prime[++cnt]=i;
            phi[i]=i-1;//i为素数时,phi[i]=i-1
        }
        for(int j=1;j<=cnt&&prime[j]*i<=maxn;j++){
            notprime[i*prime[j]]=1;
            if(i%prime[j]==0){//每个数只被它的最小质因数给筛掉
                phi[i*prime[j]]=phi[i]*prime[j];
                //当a与b互质时,满足phi(a∗b)=phi(a)∗phi(b),积性函数
                break;
            }
            else phi[i*prime[j]]=phi[i]*(prime[j]-1);
            //phi[i∗prime[j]]=phi[i]∗phi[prime[j]]=phi[i]∗(prime[j]−1);
        }
    }
}
int main(){
	pre();
}

if(i%prime[j]==0) break;

这句话让复杂度变为O(n)     orz%%

猜你喜欢

转载自blog.csdn.net/g21wcr/article/details/82290872
今日推荐