#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%%