狄利克雷卷积简介

卷积这名字听起来挺学究的,今天学了之后发现其实挺朴实hhh。

卷积:

“(n)”表示到n的一个范围。
设\(f,g\)是两个数论函数(也就是说,以自然数集为定义域的复数值函数),则卷积运算\(f\ast g\)定义为
\[(f\ast g)(n) = \sum_{ij=n}{f(i)g(j)}\]
另一种写法就是:
\[(f\ast g)(n) = \sum_{d\mid n}{f(d)g(\frac{n}{d})}\]
这里给一段数论函数的定义:

数论函数亦称算术函数,一类重要的函数,指定义在正整数集上的实值或复值函数,更一般地,也可把数论函数看做是某一整数集上定义的函数。

一些数论函数

首先最简单的就是常数函数,直接映射到一个正整数,比如\(f(x)=1,f(x)=2\)这样的。
再有就是一些整数域数列的通项公式函数,例如\(f(x)=x\)这样的。
还有就是\(\phi(x)\)欧拉函数,表示因数个数。
另外就是元函数e,写成表达式就是\(e(x)=[x=1]\).
还有特殊的常数函数,把所有的数字映射成1的\(u(x)=1\)
莫比乌斯函数:通常,莫比乌斯函数\(\mu\)定义为
\(\mu(1)=1;\)
\(\mu(n)=(-1)^k\)如果n能写成k个不同素数之积;
\(\mu(n)=0\),其他情况。

以下参考如何证明莫比乌斯反演?by Syu Gau

一些简单性质

交换律

根据\[(f\ast g)(n) = \sum_{ij=n}{f(i)g(j)}\]
这个定义,结论是显然的了。

结合律

只要证明\((f*g)*h=f*(g*h)\)就可以了。
于是左边就是
\[
\begin{align} ((f\ast g)\ast h)(n) &= \sum_{lk=n}(f\ast g)(l)h(k) \\ &= \sum_{lk=n}\left(\sum_{ij=l}f(i)g(j)\right)h(k)\\ &= \sum_{ijk=n} f(i)g(j)h(k) \end{align}
\]
右边是
\[
\begin{align} (f\ast (g\ast h))(n) &= \sum_{il=n}f(i)(g\ast h)(l) \\ &= \sum_{il=n}f(i)\left(\sum_{jk=l}g(j)h(k)\right)\\ &= \sum_{ijk=n} f(i)g(j)h(k) \end{align}
\]
得证。

加法的结合律

看不懂网上的证明,简单贴一下。

存在单位元\(\iota\) 使得\(\iota\ast f=f\)
我们需要
\[(\iota\ast f)(n)=\sum_{ij=n}\iota(i)f(j)=f(n)\]
故不难猜到\(\iota\) 应该定义为\(\iota(n)=\) \begin{cases} 1&n=1\ 0&n\neq1 \end{cases}
事实上,直接验证可得
\[(\iota\ast f)(n)=\sum_{ij=n}\delta_{i,1}f(j)=f(n)\]
以上说明数论函数在卷积意义下构成一个交换群。

卷积差不多就这些。。。

莫比乌斯反演证明

\(\mu\)函数的性质

为什么要发明这个函数呢,肯定是有道理的。
我们一般把\(\mu\)看做是\(u(x)=1\)在卷积意义下的逆元。就是说它满足:
\[\mu\ast u=e\]
1就是函数f(n)=1。展开来写就是
\[\sum_{d\mid n}\mu(d)*1\]
当 \(n=1\)时,显然成立。
当 \(n>1\)时,根据唯一分解定理我们可以把n拆成\(n=p^{k_1}_1*p^{k_2}_2*\cdots*p^{k_n}_n\)
当 \(\exists k_x=1\)时,\(\mu\)值肯定为0,所以我们把\(k_x\)都看作1。
而d枚举的就是n的因子,其实就是在n的质因子集合里取走任意个。所以这个式子可以写成这个样子:
\[
\begin{align} \sum_{d\mid n}\mu(d) =& \mu(1)+\mu(p_1)+\mu(p_2)+\cdots+\mu(p_k)+\mu(p_1p_2)+\cdots+\mu(p_1p_2\cdots p_k) \\ =& \binom{k}{0}+\binom{k}{1}(-1)+\binom{k}{2}(-1)^2+\cdots+\binom{k}{k}(-1)^k \\ =&(1-1)^k=0 \end{align}
\]
那么\[\sum_{d\mid n}\mu(d)=1\]就得证了。

反演形式1证明

法1

莫比乌斯反演形式1就是,如果\(f(n)=\sum_{d\mid n}g(d)\),则\(g(n)=\sum_{d\mid n}\mu\left(\frac{n}{d}\right)f(d)\)
写成卷积的形式就是,如果\(f=g*e\),则\(g=f*\mu\)。
这样写就比原来哪样要好记而且简介多了。
有了之前的铺垫,接下来就很容易了。
把原方程两边乘一个\(\mu\)
\[f*\mu=g*e*\mu\]
\[f*\mu=g*(e*\mu)\]
由于之前有证明\(\mu*e=1\)所以就有\(f*\mu=g\)于是得证。
感觉这种方法非常巧妙啊。。

法2

听知乎上大佬讲的,莫比乌斯反演其实就是偏序集上的容斥,简单理解了一下大概是这样的。
我们知道容斥定理的公式是
\[g(S)=\sum_{V\subset S}f(V)\implies f(S)=\sum_{V\subset S}(-1)^{\mid S\mid-\mid V\mid}*g(V)\]
用叉姐的话将就是:n 个坏事都不发生的概率,可以通过 2 n 个同时发生的概率计算,定义一个由数映射到它质因子集合的映射,映射关系显然是整除,V看做是S的质因子但不是V的质因子的乘积,那么莫比乌斯反演定理就和容斥的式子长得一模一样了。\(\mu\)就是\((-1)^{\mid S\mid-\mid V\mid}\)

反演形式2证明

以后再填坑。。感觉效率好低QAQ

\(\phi\)和μ的关系

有一个经典公式就是:
\[\sum_{d\mid n}\phi(d)=n\]
这个公式怎么证明呢?
我们可以把它简记为
\[\phi*e=id\]
然后两边乘一个\(\mu\)
\[\phi*(\mu*e)=id*\mu\]
\[\phi=id*\mu\]
再化回来
\[\phi(n)=\sum_{d\mid n}\mu(d)*\lfloor\frac{n}{d}\rfloor\]
μ只有在d质因数分解之后各个质因子个数为1的时候才有贡献,为奇数个因子的时候-,偶数为+,这不就是一个容斥求因子个数么。。于是左边等于右边,得证

狄利克雷卷积&莫比乌斯反演证明的更多相关文章

  1. 狄利克雷卷积&莫比乌斯反演总结

    狄利克雷卷积&莫比乌斯反演总结 Prepare 1.\([P]\)表示当\(P\)为真时\([P]\)为\(1\),否则为\(0\). 2.\(a|b\)指\(b\)被\(a\)整除. 3.一 ...

  2. 狄利克雷卷积&莫比乌斯反演

    昨天刚说完不搞数论了,刚看到一个\(gcd\)的题目dalao用这个做了,虽然比正解麻烦,还是打算学一学了 数论函数: 数论函数的定义: 数论函数亦称算术函数,一类重要的函数,指定义在正整数集上的实值 ...

  3. 狄利克雷卷积 & 莫比乌斯反演

    积性函数与完全积性函数 积性函数 若一个数论函数\(f\)满足当\(gcd(n,m)=1\)时,\(f(nm)=f(n)f(m)\) 则称\(f\)为积性函数 一些常见的积性函数 完全积性函数 若一个 ...

  4. [复习]莫比乌斯反演,杜教筛,min_25筛

    [复习]莫比乌斯反演,杜教筛,min_25筛 莫比乌斯反演 做题的时候的常用形式: \[\begin{aligned}g(n)&=\sum_{n|d}f(d)\\f(n)&=\sum_ ...

  5. 中国剩余定理 & 欧拉函数 & 莫比乌斯反演 & 狄利克雷卷积 & 杜教筛

    ssplaysecond的博客(请使用VPN访问): 中国剩余定理: https://ssplaysecond.blogspot.jp/2017/04/blog-post_6.html 欧拉函数: h ...

  6. [BZOJ 2301] [HAOI 2011] Problem b (莫比乌斯反演)(有证明)

    [BZOJ 2301] [HAOI 2011] Problem b (莫比乌斯反演)(有证明) 题面 T组询问,每次给出a,b,c,d,k,求\(\sum _{i=a}^b\sum _{j=c}^d[ ...

  7. BZOJ3601. 一个人的数论(狄利克雷卷积+高斯消元)及关于「前 $n$ 个正整数的 $k$ 次幂之和是关于 $n$ 的 $k+1$ 次多项式」的证明

    题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=3601 题解 首先还是基本的推式子: \[\begin{aligned}f_d(n) &a ...

  8. 我也不知道什么是"莫比乌斯反演"和"杜教筛"

    我也不知道什么是"莫比乌斯反演"和"杜教筛" Part0 最近一直在搞这些东西 做了将近超过20道题目吧 也算是有感而发 写点东西记录一下自己的感受 如果您真的 ...

  9. 莫比乌斯反演学习笔记+[POI2007]Zap(洛谷P3455,BZOJ1101)

    先看一道例题:[POI2007]Zap BZOJ 洛谷 题目大意:$T$ 组数据,求 $\sum^n_{i=1}\sum^m_{j=1}[gcd(i,j)=k]$ $1\leq T\leq 50000 ...

随机推荐

  1. Ui自动化测试上传文件方法都在这里了

    前言 实施UI自动化测试的时候,经常会遇见上传文件的操作,那么对于上传文件你知道几种方法呢?今天我们就总结一下几种常用的上传文件的方法,并分析一下每个方法的优点和缺点以及哪种方法效率,稳定性更高 被测 ...

  2. excel报表开发-- 根据datatable个数自动生成新sheet

    总结一下很久之前做的报表小程序,今日有问题又调试了一下. DB中存在一个表,记录了ID<自增长>,SqlStatement<sql查询语句>及其他必要字段,比如SheetNam ...

  3. getlasterror() 输出错误信息,

    得自http://bbs.csdn.net/topics/390416234 LPVOID lpMsgBuf;  FormatMessage(      FORMAT_MESSAGE_ALLOCATE ...

  4. Spring_AOP基于AspectJ的注解开发&JDBC的模板使用&事务管理(学习笔记3)

    一:AOP基于AspectJ的注解开发 1,简单的实例: 1)引入相应的jar包 ​ 2)在配置文件里引入相关约束 <beans xmlns="http://www.springfra ...

  5. APP爬虫(2)把小姐姐的图片down下来

    APP爬虫(1)想学新语言,又没有动力,怎么办? 书接上文.使用appium在帖子列表界面模拟上划动作,捕捉不到列表的规律.上划结束后,列表只能获取到屏幕呈现的部分,而且下标还是从0开始的. 根据酸奶 ...

  6. tar 命令详解(持续更新)

    可以用man tar查看tar命令使用的权威解释 Main operation mode: -c: 建立压缩档案 -r:向压缩归档文件末尾追加文件 -t:查看内容 -u:更新原压缩包中的文件 -x:解 ...

  7. jsp数据交互(一).2

    01.什么是JSP内置对象(jsp核心)? Java 内置对象 Java  作用域 解析:jsp内置对象是web容器创建的一组对象.我们都知道tomcat这款软件可以看成是一种web容器,所以我们可以 ...

  8. Guitar Pro如何更改五线谱的符杆方向

    可能有的小伙伴不知道Guitar Pro是什么软件,我先稍微给大家介绍一下~ Guitar Pro是专为帮助所有吉他爱好者学习.绘谱.创作的多功能软件.它包含所有吉他的现有指法和音色,可以帮助我们了解 ...

  9. MOCTF-Crypt-writeup

    MOctf Crypt Writeup记录 都不难,就随便记录记录下. MOCTF平台地址:http://www.moctf.com 0x01 数据库密码 hint:20岁的小刚,自幼热爱信息安全,一 ...

  10. StarUML 3.0 破解方法

    首先在我这里下载  StarUML3.0 破解替换文件app.asar 链接:https://pan.baidu.com/s/1wDMKDQkKrE9D1c0YeXz0xg 密码:y65m 然后参照下 ...