莫比乌斯与欧拉

#include
using namespace std;
const int N=1e6+10;
int mu[N],prime[N],cnt,ol[N];
bool vis[N];
void ola(int n){
	vis[1]=mu[1]=ol[1]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			prime[++cnt]=i;
			mu[i]=-1;
			ol[i]=i-1;
		}
		for(int j=1;j<=cnt&&i*prime[j]<=n;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0){
				mu[i*prime[j]]=0;
				ol[i*prime[j]]=ol[i]*prime[j];
				break;
			}
			mu[i*prime[j]]=-mu[i];
			ol[i*prime[j]]=ol[i]*ol[prime[j]];
		}
	}
}
int main(){
	int n,T;
	cin>>n>>T;
	ola(n);
	while(T--){
		int opt,x;
		scanf("%d%d",&opt,&x);
		if(opt==1){
			printf("%d\n",mu[x]);
		}
		else{
			printf("%d\n",ol[x]);
		}
	}
	return 0;
}

你可能感兴趣的:(算法,图论,数据结构)