莫比乌斯函数理论基础+实战例题

最新推荐文章于 2023-01-26 16:02:51 发布

H_Z_01

最新推荐文章于 2023-01-26 16:02:51 发布

阅读量239

点赞数

文章标签：算法莫比乌斯反演

本文链接：https://blog.csdn.net/weixin_46783197/article/details/118423778

版权

在这里插入图片描述

求莫比乌斯函数μ(n)模板：

const int Max=1e7+10;
bool ver[Max];
int mo[Max],prim[Max],ans=0;
void set_mobius(int n)
{
    memset(ver,false,sizeof(ver));
    mo[1]=1;
    for(long long i=2;i<n;i++)
    {
        if(!ver[i])
        {
            ver[i]=true;
            prim[ans]=i;
            ans++;
            mo[i]=-1;
        }
        for(long long j=0;j<ans&&prim[j]*i<n;j++)
        {
            if(prim[j]*i<=n)
            {
                ver[prim[j]*i]=true;
                if(i%prim[j]==0)
                {
                    break;
                }
                else
                {
                    mo[i*prim[j]]=-mo[i];
                    //mo[i*prim[j]]=mo[i]*mo[prim[j]];
                }
            }
        }
    }
}

题目地址：https://vjudge.net/contest/436024#problem/F
在这里插入图片描述

ac代码：

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
const int Max=1e7+10;
bool ver[Max];
int mo[Max],prim[Max],ans=0;
void set_mobius(int n)
{
    memset(ver,false,sizeof(ver));
    mo[1]=1;
    for(long long i=2;i<n;i++)
    {
        if(!ver[i])
        {
            ver[i]=true;
            prim[ans]=i;
            ans++;
            mo[i]=-1;
        }
        for(long long j=0;j<ans&&prim[j]*i<n;j++)
        {
            if(prim[j]*i<=n)
            {
                ver[prim[j]*i]=true;
                if(i%prim[j]==0)
                {
                    break;
                }
                else
                {
                    mo[i*prim[j]]=-mo[i];
                    //mo[i*prim[j]]=mo[i]*mo[prim[j]];
                }
            }
        }
    }
}
int main()
{
    set_mobius(Max);
    long long t,n;
    cin>>t;
    long long sum;
    while(t--)
    {
        cin>>n;
        sum=0;
        for(long long i=1;i<=n/i;i++)
            if(mo[i])
                sum+=n/i/i*mo[i];
        cout<<sum<<endl;
    }
    return 0;
}