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

Prepare

1、\([P]\)表示当\(P\)为真时\([P]\)为\(1\),否则为\(0\)。

2、\(a|b\)指\(b\)被\(a\)整除。

3、一些奇怪常见的函数:

\(1(n)=1\)

\(id(n)=n\)

\(\sigma(n)=n的约数和\)

\(d(n)=n的约数个数\)

\(\epsilon(n)=[n==1]\)

狄利克雷卷积

数论函数

数论函数指一类定义域是正整数,值域是一个数集的函数。

加法:逐项相加就可以辣\((f+g)(x)=f(x)+g(x)\)

数乘:用一个常数乘\((xf)(n)=x*f(n)\)

狄利克雷卷积

定义两个数论函数的狄利克雷卷积\(*\):

若\(t=f*g\)则

\[t(n)=\sum_{i|n}f(i)g(\frac{n}{i})
\]

等价于

\[t(n)=\sum_{ij=n}f(i)g(j)
\]

狄利克雷卷积有以下性质(两个数论函数相等,是指两个函数的每一项都相等):

1、交换律\(f*g=g*f\)

2、结合律\(f*(g*h)=(f*g)*h\)

3、分配律\(f*h+g*h=(f+g)*h\)

4、没有名字\((xf)*g=x(f*g)\)

5、单位元\(\epsilon*f=f\),其中\(\epsilon(n)=[n==1]\)

6、逆元:对于每一个\(f(1)\neq 0\)的函数\(f\),都有\(f*g=\epsilon\)

讨论一下第六个结论,如何求一个函数的逆呢?

只需要定义

\[g(n)=\frac{1}{f(1)}\left([n==1]-\sum_{i|n,i\neq 1}f(i)g(\frac{n}{i})\right)
\]

这样的话

\[\sum_{i|n}f(i)g(\frac{n}{i})=f(1)g(n)+\sum_{i|n,i\neq1}f(i)g(\frac{n}{i})=[n==1]
\]

积性函数

如果一个数论函数\(f\)有当\(gcd(n,m)==1\)时

\[f(nm)=f(n)f(m)
\]

就称\(f\)为积性函数。

一些常见的积性函数:

\(\epsilon(n)=[n==1]\),\(id(n)=n\),\(id^{k}(n)=n^k\)

事实上他们也满足完全积性(即当\(gcd(n,m)\neq1\)时,也有\(f(nm)=f(n)f(m)\))

特殊的,我们令\(id^0(n)=1(n)=1\)

还有两个普通的积性函数

\(d(n)=n的约数和\)、\(\varphi(n)=[1,n]中与n互质的数的个数\)

还有两个重要结论:

两个积性函数的狄利克雷卷积是积性函数。

积性函数的逆是积性函数。

积性函数有什么用呢?

它可以线性筛

然而还有更有用的\(---\)

莫比乌斯反演

一些理论

我们定义\(1\)的逆是\(\mu\)

这样的话,如果\(g=f*1\),就有\(f=f*1*\mu=g*\mu\)

换句话说,就是

\[g(n)=\sum_{d|n}f(d)\Leftrightarrow f(n)=\sum_{d|n}\mu(\frac{n}{d})g(d)
\]

也可以这样子

\[g(d)=\sum_{d|n}f(n)\Leftrightarrow f(d)=\sum_{d|n}\mu(\frac{n}{d})*g(n)
\]

例子

怎么用呢?举几个例子(以下情况默认\(n\leq m\))


\(Eg1\)

\[\sum_{i=1}^n\sum_{i=1}^m[gcd(i,j)==1]
\]

然后怎么办呢?

\[f(x)=\sum_{i=1}^n\sum_{i=1}^m[gcd(i,j)==x] \\
g(x)=\sum_{x|d}f(d)
\]

\[f(1)=\sum_{1|d}\mu(\frac{d}{1})g(d) \\
f(1)=\sum_{i=1}^n\mu(i)g(i)
\]

考虑\(g(x)\)是什么

\[g(x)=\sum_{x|d}\sum_{i=1}^n\sum_{i=1}^m[gcd(i,j)==d]
\]

\[g(x)=\sum_{i=1}^n\sum_{i=1}^m[x|gcd(i,j)] \\
g(x)=\sum_{i=1}^{\left\lfloor\frac nx\right\rfloor}\sum_{i=1}^{\left\lfloor\frac mx\right\rfloor}[1|gcd(i,j)]\\
g(x)=\left\lfloor\frac nx\right\rfloor\left\lfloor\frac mx\right\rfloor
\]

带回\(f(1)\)

\[Ans=\sum_{x=1}^n\mu(x)\left\lfloor\frac nx\right\rfloor\left\lfloor\frac mx\right\rfloor
\]

这个用整除分块可以做到\(O(\sqrt n)\)


\(Eg2\)

\[\sum_{i=1}^n\sum_{j=1}^mgcd(i,j)
\]

可化为

\[\sum_{d=1}^n\sum_{i=1}^n\sum_{j=1}^md[gcd(i,j)==d] \\
\sum_{d=1}^nd\sum_{i=1}^{\left\lfloor\frac nd\right\rfloor}\sum_{j=1}^{\left\lfloor\frac md\right\rfloor}[gcd(i,j)==1]
\]

\[x=\left\lfloor\frac nd\right\rfloor, y=\left\lfloor\frac md\right\rfloor
\]

\[\sum_{d=1}^nd\sum_{i=1}^x\sum_{j=1}^y[gcd(i,j)==1]
\]

套入我们刚才在\(Eg1\)求得的

\[\sum_{d=1}^nd\sum_{i=1}^x\mu(i)\left\lfloor\frac xi\right\rfloor\left\lfloor\frac yi\right\rfloor\\
\sum_{d=1}^nd\sum_{i=1}^{n/d}\mu(i)\left\lfloor\frac {x}{id}\right\rfloor\left\lfloor\frac {y}{id}\right\rfloor
\]

化到现在是\(O(n)\)的,因为前后都可以数论分块

但是我们能做得更好

令\(T=id\)

原式化为

\[\sum_{d=1}^nd\sum_{i=1}^n\mu(i)\left\lfloor\frac {x}{T}\right\rfloor\left\lfloor\frac {y}{T}\right\rfloor\\
\sum_{T=1}^n\left\lfloor\frac {x}{T}\right\rfloor\left\lfloor\frac {y}{T}\right\rfloor\sum_{d|T}d\mu(\frac Td)
\]

乍一看还是\(O(n)\)的呀,但是对于后面那一坨

\[f(T)=\sum_{d|T}d\mu(\frac Td)\\
\Leftrightarrow\\
f(T)=\sum_{ij=T}id(i)\mu(j)
\]

两个积性函数相乘,可以线性筛呀!!

所以复杂度被我们压到了\(O(\sqrt n)\)

\(upd\;on\;2019.3.9:\)

发现以前没有讲线性筛,导致现在自己都不知道是怎么搞得了。。。

线性筛

我们筛\(\mu\)的函数是长这样的(自动认为有模数):

int prime[MAX_N], mu[MAX_N]
bool nprime[MAX_N];
void sieve() {
mu[1] = 1;
for (int i = 2; i <= N; i++) {
if (!nprime[i]) prime[++tot] = i, mu[i] = Mod - 1;
for (int j = 1; i * prime[j] <= N; j++) {
nprime[i * prime[j]] = 1;
if (i % prime[j] == 0) break;
mu[i * prime[j]] = Mod - mu[i];
}
}
}

其实就是和筛素数的是一样的,

if (i % prime[j] == 0) break;

这句话保证了复杂度,因为你存的素数是递增的,

如果\(i\)被\(prime[j]\)整除后,\(i*prime[j+k](k>0)\)一定可以被\(prime[j]*x\)的形式表示出来。

那么就有我们下面的一个问题:


\(Eg3\)

给定你一个数组\(f\),求

\[\sum_{i=1}^n\sum_{j=1}^mf[gcd(i,j)]
\]

其中\(n,m\leq 10^7\),数据组数\(T\leq 10^4\)。

由我们上面推的东西,将\(f\)看作一个数论函数,可以知道只要求出一个函数\(g=\mu * f\)的前缀和,

这个问题就解决了。

一下是解决这个问题的几种方法(蒯的):

void get_g_1(int N, const int *f, int *g) {
for (int i = 1; i <= N; i++) g[i] = 0;
for (int i = 1; i <= N; i++)
for (int j = 1; i * j <= N; j++)
g[i * j] = (g[i * j] + mu[i] * f[j]) % mod;
} // 依照定义,O(nlogn) void get_g_2(int N, const int *f, int *g) {
for (int i = 1; i <= N; i++) g[i] = f[i];
for (int i = 1; i <= N; i++)
for (int j = 2; i * j <= N; j++)
g[i * j] = (g[i * j] - g[i]) % mod;
} // 类似求狄利克雷卷积逆的方式,不需要线性筛 mu ,O(nlogn) void get_g_3(int N, const int *f, int *g) {
for (int i = 1; i <= N; i++) g[i] = f[i];
for (int i = 0; i < prime_count; i++)
for (int j = N / prime[i]; j >= 1; j--)
g[j * prime[i]] = (g[j * prime[i]] - g[j]) % mod;
} // Magic! O(nloglogn)

对于最后一种方法,理解成dp:

\[g_{i,n}=\sum_{d|n,d只包含前i种质因子}\mu(d)f(n/d)
\]

那么转移:

\[g_{i,n}=\begin{cases}
g_{i-1,n}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;p_i\nmid n\\
g_{i-1,n}-g_{i-1,n/p_i}\;\;\;\;\;p_i\mid n
\end{cases}
\]

复杂度\(O(n\log\log n)\)。


参考文章1

参考文章2

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

  1. 狄利克雷卷积&莫比乌斯反演证明

    狄利克雷卷积简介 卷积这名字听起来挺学究的,今天学了之后发现其实挺朴实hhh. 卷积: "(n)"表示到n的一个范围. 设\(f,g\)是两个数论函数(也就是说,以自然数集为定义域 ...

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

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

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

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

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

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

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

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

  6. 【BZOJ3529】数表(莫比乌斯反演,树状数组)

    [BZOJ3529]数表(莫比乌斯反演,树状数组) 题解 首先不管\(A\)的范围的限制 要求的东西是 \[\sum_{i=1}^n\sum_{j=1}^m\sigma(gcd(i,j))\] 其中\ ...

  7. 【Luogu3768】简单的数学题(莫比乌斯反演,杜教筛)

    [Luogu3768]简单的数学题(莫比乌斯反演,杜教筛) 题面 洛谷 \[求\sum_{i=1}^n\sum_{j=1}^nijgcd(i,j)\] $ n<=10^9$ 题解 很明显的把\( ...

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

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

  9. 【51NOD 1847】奇怪的数学题(莫比乌斯反演,杜教筛,min_25筛,第二类斯特林数)

    [51NOD 1847]奇怪的数学题(莫比乌斯反演,杜教筛,min_25筛,第二类斯特林数) 题面 51NOD \[\sum_{i=1}^n\sum_{j=1}^nsgcd(i,j)^k\] 其中\( ...

随机推荐

  1. loadrunner 脚本优化-检查点设置

    脚本优化-检查点设置 by:授客 QQ:1033553122 VuGen判断脚本是否执行成功是根据服务器返回的状态来确定的,如果服务器返回的是HTTP状态为200 OK,那么VuGen就认为脚本正确地 ...

  2. java数据结构 • 面向对象 • 异常 • 随机数·时间

    • 语法基础 • 控制流 • 数据结构 • 面向对象 • 异常 • 随机数 //String常用的方法: indexOf   charAt   charAt   codePointAt   compa ...

  3. Linux 中提高的 SSH 的安全性

    SSH 是远程登录 Linux 服务器的最常见的方式.且 SSH 登录的时候要验证的,相对来讲会比较安全.那只是相对,下面会介绍一些方式提高 SSH 的安全性 SSH 的验证 而SSH 登录时有两种验 ...

  4. Jmeter和Loadrunner的异同

    1.jmeter的架构跟loadrunner原理一样,都是通过中间代理,监控&收集并发客户端发现的指令,把他们生成脚本,再发送到应用服务器,再监控服务器反馈的结果的一个过程. 2.分布式中间代 ...

  5. Lua不显示小数点0的部分

    我的环境:Unity3D 5.3.7p4 XLua版本v2.1.6 基于Lua5.3 (https://github.com/Tencent/xLua) 在Lua中数字不区分整型或浮点型,所有都是nu ...

  6. 注入攻击(SQL注入)

    注入攻击是web安全领域中一种最为常见的攻击方式.注入攻击的本质,就是把用户输入的数据当做代码执行.这里有两个关键条件,第一是用户能够控制输入,第二个就是原本程序要执行的代码,将用户输入的数据进行了拼 ...

  7. mybatis : ERROR. token : COMMA, pos : 373

    查询的字段符号问题,可能是多了" , "  也可能是少了 " , "  仔细检查, 都能解决

  8. SpringCloud之初识Hystrix熔断器 ----- 程序的保护机制

    在上一篇的-负载均衡Robbin中,我们简单讲解到负债均衡的算法和策略.负载均衡就是分发请求流量到不同的服务器,以减小服务器的压力和访问效率,但是当负载均衡的某个服务器或是服务挂掉之后,那么程序会出现 ...

  9. 【Teradata】grouping和rollup窗口函数

    1.group by后带rollup子句 先按一定的规则产生多种分组,然后返回各个分组所产生的结果集的并集,且没有去掉重复数据(统计出的数据是求和还是最大值还是平均值等这就取决于SELECT后的聚合函 ...

  10. 【Linux基础】Linux常用命令汇总

    3-1文件目录操作命令(cd pwd mkdir rmdir rm) 绝对路径:由根目录(/)开始写起的文件名或目录名称, 例如 /home/dmtsai/.bashrc: 相对路径:相对于目前路径的 ...