线性筛——莫比乌斯函数
如果不会线性筛素数的话,建议先看这篇博客了解一下线性筛素数。
莫比乌斯函数函数(积性函数都可以线性筛)主要是在线性筛素数的基础上得到的
我们知道:
若
则
所以:
1、当 是质数时,
对于 2和3 设 其中 为 的最小质因子
2、当 是 的某个质因子时, 则
3、当 与 互质时,
good luck and have fun!!!
附上代码:
int mu[MAXN];
void Mobius(int n)
{
memset(vis,0,sizeof(vis));
mu[1]=1;
prime[0]=0;
for(int i=2;i<=n;i++)
{
if(!vis[i])
{
prime[++prime[0]]=i;
mu[i]=-1;
}
for(int j=1;j<=prime[0]&&i<=n/prime[j];j++)
{
vis[i*prime[j]]=1;
if(i%prime[j]==0)
{
mu[i*prime[j]]=0;
break;
}
mu[i*prime[j]]=-mu[i];
}
}
}
good luck and have fun!!!
浙公网安备 33010602011771号