zoukankan      html  css  js  c++  java
  • 莫比乌斯反演理解

    学了下莫比乌斯反演,实质上mu函数是从容斥原理推出来的

    详细的原理解读+例子讲解可以看博客 https://www.cnblogs.com/chenyang920/p/4811995.html

    好博客 https://blog.csdn.net/tomandjake_/article/details/81082837

    首先要了解狄利克雷卷积公式 

      (fg)(n)=d|nf(d)g(n/d)

    一些常用的数论函数(积性函数):

      
    积性函数

      第一个是狄利克雷单位元,相当于常数1(可以参考线性代数里的单位矩阵)

      1(n) = 1 是 单位函数

      Id(n) = n 是 不变函数

    关于一些性质和公式 https://www.cnblogs.com/Mychael/p/8744633.html

    莫比乌斯反演

     (1)    F(n)=sum_{d|n}f(d)    ====》      f(n)=sum_{d|n}mu(d)F(frac{n}{d})  

     (2) 

    莫比乌斯的逆元:单位函数:u(n) * 1(n) = € (n)

    证明 1:F(n)=(f*1)(n)  ==> (F*u)(n) = (f*1*u)(n) = f(n) ==> f(n)=sum_{d|n}mu(d)F(frac{n}{d})

    另外还有一种证明方法,个人认为更加直接硬核

      

    其中第三个等号那里用了一下置换法

    引理:∑d|n u(d)=  € (n)

    证明2:略 ,基本上是魔方上面第二种证明方法

    关于u(n)函数的一些直观理解

      g(1)=f(1)

      g(p1)=f(1)+f(p2)

      g(p2)=f(1)+f(p2)

      g(p1p2)=f(1)+f(p1)+f(p2)+f(p1p2)

      那么可以往回推出 f(p1p2)=g(p1p2)-g(p1)-g(p2)+g(1)

      来看每项系数,g(p1p2)=1,因为 f(p1p2)仅仅存在g(p1p2)中,所以第一项系数为1,对应莫比乌斯函数u(1)=1

      然后看第二项和第三项 ,因为f(p1)和f(p2)在g(p1p2)中,所以需要减掉,所以第二项和第三项的系数为-1 对应莫比乌斯函数 u(p1)=u(p2)=-1

      同理,第四项系数为1 对应莫比乌斯函数u(p1p2)=1

    然后推广到n=p1p2p3p4...pk,即n没有平方因子的情况 

      第一项g(n)的系数必定是1

      g(n/pi)的系数必定是 -1=u(pi)

      g(n/(pi*pj))的系数必定是 1=u(pi*pj)  

      利用容斥原理推导公式,依次类推,当u(n),n为没有平方因子的数,k为n的质因子个数,u(n)=(-1)^n

    再推广到有平方因子的情况

      我们只考虑平方因子在的情况

      g(n/(pi*pi))的系数是 0=u(pi*pi)

      g(n/(pi*pi*pj))的系数是 0=u(pi*pi*pj)

      根据容斥原理推导公式,当u(n),n有平方因子,则u(n)=0

     

    关于莫比乌斯函数的两条性质

                  用二项式定理可证,实质还是容斥原理(莫比乌斯函数的逆元是单位函数也是这么来的)

                  不知道什么东西。。

  • 相关阅读:
    hdu5358 推公式+在一个区间内的尺取+枚举法
    poj3349 散列表(hash)
    hdu3282 链表或者对顶堆
    hdu5178 尺取
    hdu5672 尺取
    hdu3244完全背包+二分答案 留坑
    hdu5256 二分求LIS+思维
    hdu5646数学构造+二分
    hdu4190 二分答案
    Python中Scapy网络嗅探模块的使用
  • 原文地址:https://www.cnblogs.com/zsben991126/p/10732968.html
Copyright © 2011-2022 走看看