莫比乌斯反演学习笔记

最新推荐文章于 2019-04-03 21:04:29 发布

dengpailuo0807

最新推荐文章于 2019-04-03 21:04:29 发布

阅读量76

点赞数

原文链接：http://www.cnblogs.com/HenryHuang-Never-Settle/p/10472218.html

版权

前置芝士

~~极高的数学造诣与不怕劳累的精神~~

还有可能会用到的这么个东西 ~~虽然好像真的rbl~~

正文

我们可以非常轻易地手动模拟出\({f(n)}\)与\({g(d)}\)的关系

\[{f(1)=g(1)}\]

\[{f(2)=g(1)+g(2)}\]

\[{...}\]

\[{f(6)=g(1)+g(2)+g(3)+g(6)}\]

与此同时，我们也一样能够从\({f(n)}\)逆推到\({g(d)}\)

\[{g(1)=f(1)}\]

\[{g(2)=f(2)-f(1)}\]

\[{...}\]

\[{g(6)=f(6)-f(3)-f(2)+f(1)}\]

我们会发现，这样的一个式子，并没有从\({f(n)}\)到\({g(d)}\)那么和谐，而是有加有减，似乎不太好处理。

\[{g(n)=\sum_{d|n}^{}}f(d)*? \]

我们会发现，每一个\({g(n)}\)都是由上面这样一个式子推出来的，只不过还多了一个系数

然后我们就要引入一个新的东西

莫比乌斯函数

莫比乌斯函数是一个很有趣的函数，其中有很多很优美，我们也特别喜欢的性质

性质1：莫比乌斯函数是一个积性函数

这个东西应该没什么好说的，自己带一带，YY一下就能出来~~（其实是我自己懒得证明啦QAQ）~~

这对我们有什么好处呢？

还记得欧拉函数\({\varphi (n)}\)吗

所以同样，我们也可以用线性筛来求解莫比乌斯函数

第一种求法：利用定义来求莫比乌斯函数

贴代码吧，自行理解

void mobius()
{
    for(int i=1;i<=maxn;++i)
        miu[i]=1;
    for(int i=2;i<=maxn;++i)
    {
        if(!vis[i])
        {
            miu[i]=-1;
            for(int j=i+i;j<=maxn;j+=i)
            {
                vis[j]=1;
                if((j/i)%i==0) miu[j]=0;
                else miu[j]*=-1;
            }
        }
    }
}