#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<stack>
#include<vector>
#include<cstring>
#define N 100005
using namespace std;
vector<int> prime;
bool pri[N];
int mu[N];
void init()
{
memset(pri,true,sizeof(pri));
mu[1]=1;
for(int i=2;i<N;i++)
{
if(pri[i])
{
prime.push_back(i);
mu[i]=-1;
}
for(int j=0;j<prime.size()&&i*prime[j]<N;j++)
{
pri[i*prime[j]]=false;
if(i%prime[j])
mu[i*prime[j]]=-mu[i];
else
{
mu[i*prime[j]]=0;
break;
}
}
}
}
莫比乌斯函数模板
最新推荐文章于 2022-06-29 18:24:06 发布