莫比乌斯函数(Mobius)的求法每日一遍，算法再见！

最新推荐文章于 2025-02-05 21:53:11 发布

鲜果维他命

最新推荐文章于 2025-02-05 21:53:11 发布

阅读量1.9k

点赞数

分类专栏：数论文章标签：算法

本文链接：https://blog.csdn.net/TheWayForDream/article/details/113871548

版权

数论专栏收录该内容

20 篇文章

订阅专栏

博客介绍了莫比乌斯函数的概念及其性质，重点讲解了如何使用欧拉筛法来高效计算莫比乌斯函数。代码示例展示了如何初始化一个数组mu，通过遍历并标记素数及其倍数，根据莫比乌斯函数的规则更新数组值，从而得到所有整数的莫比乌斯函数值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

莫比乌斯函数

我们用符号u表示莫比乌斯函数，u(x)表示x的莫比乌斯函数
①x=1,u(1)=1.
②x=p1p2p3p4…pk,其中p1~pk是x的素数因子，u(x)= $1)^k$
③其他情况u(x)=0

②和③的情况通俗的讲，就是把x用素数乘积的形式表示，如果存在一个素数的指数大于1那么u(x)=0,否则u(x)= $1)^k$ .

欧拉筛求Mobius函数
代码：

void get_mu()
{
	memset(vis,false,sizeof(vis));
	mu[1]=1;
	vis[1]=true;
	for(int i=2;i<=Maxn;i++)
	{
	    if(!vis[i])
	    {
	    	prime[++prime[0]]=i;
	    	vis[i]=true;
	    	mu[i]=-1;
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<=Maxn;j++)
		{
			vis[i*prime[j]]=true;
			if(i%prime[j]==0)
			{
				mu[i*prime[j]]=0;
				break;
			}
			else
			{
				mu[i*prime[j]]=-mu[i];
			}
		}
	}
	return ;
}