首页 > 其他分享 >莫比乌斯反演学习笔记

莫比乌斯反演学习笔记

时间:2024-12-21 14:08:54浏览次数:3  
标签:lfloor geq frac 函数 乌斯 sum rfloor 反演 莫比

前置知识


一、整除分块

即按照 \(\lfloor \frac{n}{i} \rfloor\) 的值域进行分块,块数 $ \leq \lfloor 2\sqrt n \rfloor $。分 $i \leq \lfloor \sqrt n \rfloor $ ,$ i > \lfloor \sqrt n \rfloor $ 讨论。

\(i\) 所在块的右端点为 $ \lfloor \frac{n}{ \lfloor \frac{n}{i} \rfloor } \rfloor $ ,代码即
$ r=n/(n/l) $

这样,我们可以把 \(O(n)\) 的求 $ \sum _{i=1}^{n} \lfloor \frac{n}{i} \rfloor * f(i)$ 值过程变为 \(O(\sqrt n)\) ( \(f(i)\)用前缀和处理 )

  • 无关的小知识

\(1^2+2^2+3^2+...+n^2=\frac{n(n+1)(2n+1)}{6}\)

\(1^3+2^3+3^3+...+n^3=\frac{n^2(n+1)^2}{4}\)

另外,在做题的时候时刻关注的是式子的形式和范围,除数不为0,r不大于n

二、普通生成函数

形如 \(f(x)=\sum _n a_nx^n\),此处a可以有限也可以无限

如序列 \(a=1,3,5,7...\) (编号从0开始) 的生成函数是 $ \sum _{n \geq 0} (2n+1)x^n$ ,如果a有通项公式,则它的普通生成函数的系数就是通项公式

基本运算

考虑 \(a,b\) 的普通生成函数 \(f(x),g(x)\)

加法,$ f(x) \pm g(x) =\sum _n (a_n \pm b_n)x^n$ ,即 $ f(x) \pm g(x) $ 是序列 $ q_n = (a_n \pm b_n) $ 的普通生成函数

乘法,或卷积,$ f(x)g(x)= \sum_n x^n \sum_{i=0}^n a_i b_{n-i} $,即 \(f(x)g(x)\) 是序列 \(q_n=\sum_{i=0}^n a_i b_{n-i}\) 的普通生成函数

关键在构造,一般指数是题中给的限制,某项系数是最终答案

可以解决多重集组合数

三、指数生成函数

形如 $f(x)=\sum_{n \geq 0} a_n \frac{x^n}{n!} $

这里,如果原数列是有限的,就把大于数列上界的 \(a_n\) 设为 0 即可

little question: 封闭形式如何得到?(坑)

加法同上

乘法(卷积),

\(f(x)g(x)=\sum_{i \geq 0} a_i \frac{x^i}{i!} \sum_{j \geq 0} b_j \frac{x^j}{j!}\)

\(=\sum_{n \geq 0} x^n \sum_{i=0}^n a_i b_{n-i} \frac{1}{i!(n-i)!}\)

$=\sum_{n \geq 0} \frac{x^n}{n!} \sum_{i=0}^{n} \frac{n!}{i!(n-i)!} a_i b_{n-i} $

\(=\sum_{n \geq 0} \frac{x^n}{n!} \sum_{i=0}^{n} C_n^i a_i b_{n-i}\)

即 \(f(x)g(x)\) 是序列 \(q_n=\sum_{i=0}^{n} C_n^i a_i b_{n-i}\)

可以解决多重集排列数

四、狄利克雷生成函数

形如 \(f(x)=\sum_{n>0} \frac{a_n}{n^x}\)

乘法(卷积),

\(\sum_{i>0} \frac{a_i}{i^x} \sum_{j>0} \frac{b_j}{j^x}\)

\(=(\frac{a_1}{1^x}+\frac{a_2}{2^x}+\frac{a_3}{3^x}+\frac{a_4}{4^x}+...)(\frac{b_1}{1^x}+\frac{b_2}{2^x}+\frac{b_3}{3^x}+\frac{b_4}{4^x}+...)\)

\(=\frac{a_1b_1}{1^x}+\frac{a_1b_2+a_2b_1}{2^x}+\frac{a_1b_3+a_3b_1}{3^x}+\frac{a_1b_4+a_2b_2+a_4b_1}{4^x}+...\)

\(=\sum_{n>0} \frac{1}{n^x} \sum_{d|n} a_db_{\frac{n}{d}}\)

五、积性函数

\(f(1)=1\),当\(gcd(a,b)=1\)时,有\(f(ab)=f(a)f(b)\),则\(f(n)\)为积性函数

欧拉函数和莫比乌斯函数都是积性函数

欧拉函数

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

又 \(\varphi(n)=n(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_s})\),其中p为n的质因数

即n以内,和n互质的数的个数

性质:\(\sum_{d|n} \varphi(d)=n\)

证明:将以n为分母,且比1小的非负分数写出来,化简后按分母归类,发现每一类 \(d\) 的个数都是 \(\varphi(d)\)

正主出场:莫比乌斯函数

定义:

\(\mu(n)=\left\{ \begin{array}{rcl} 1 & & n=1 \\ (-1)^s & & n=p_1p_2...p_s \\ 0 & & \text{n包含相同质因子}\\ \end{array} \right.\)

性质: $\sum_{d|n} \mu(d)=[n=1] $

简单证明:

\(n>1\) 时,\(n=p_1^{a_1}p_2^{a_2}..p_s^{a_s}\)

令 \(n'=p_1p_2...p_s\)

又 \(\sum_{d|n} \mu(d)=\sum_{d|n'}\mu(d)\)

约数由质因子的成绩构成,又有容斥原理

\(\sum_{d|n'}\mu(d)=C_s^0+(-1)C_s^1+(-1)^2C_s^2+...+(-1)^sC_s^s\)

\(=(1+(-1))^s\)

\(=0\)

与欧拉函数的联系:\(\sum_{d|n} \mu(d) \frac{n}{d} = \varphi(n)\)

证明是把n提出来,因式分解,化成欧拉函数的形式。

性质应用:\([gcd(i,j)=1]=\sum_{d|gcd(i,j)} \mu(d)\)

标签:lfloor,geq,frac,函数,乌斯,sum,rfloor,反演,莫比
From: https://www.cnblogs.com/YYYmoon/p/18620364

相关文章

  • 反演笔记
    反演反演若已知\(f(n)=\sumg(k)\),用\(f\)表示\(g\)的过程就叫“反演”。二项式反演参考一下邓老师的PPT。经典题:\(n\)个元素错排的方案数。要求线性。考虑枚举有\(k\)个人非错排,可以得到这\(k\)个人一共有\(\dbinom{n}{k}\)种选法,剩下的人有\((n-k)\)......
  • 人工智能原理期末速成——消解反演与反演求解
    消解反演:证明真或假反演求解:求解变量值消解反演解法:1.否定命题,将否定后的命题加入前提2.通过前提之间互相组合,得到新的前提3.最终前提与前提互相矛盾,得到NIL,此时证明完成注:在第1步之后第2步之前,还需要将所有的命题转换成子句形式例1:设已知:(1)能阅读的人是识字的(2)海豚不......
  • 遥感技术在生态系统碳储量、碳收支、碳循环以及人为源排放反演等领域的技术发展
    以全球变暖为主要特征的气候变化已成为全球性环境问题,对全球可持续发展带来严峻挑战。2015年多国在《巴黎协定》上明确提出缔约方应尽快实现碳达峰和碳中和目标。2019年第49届 IPCC全会明确增加了基于卫星遥感的排放清单校验方法。随着碳中和目标以及全球碳盘点的现实压力,基于......
  • 子集反演 & sos dp 学习笔记
    子集反演&sosdp学习笔记子集反演设\(g(S)\)表示集合\(S\)的答案,\(f(S)\)为\(S\)的子集的答案和。根据定义:\[f(S)=\sum_{T\inS}g(T)\]子集反演就是:\[g(S)=\sum_{T\inS}(-1)^{|S|-|T|}f(T)\]本质上就是容斥原理,可感性理解,证明略(给你你也记不住)。于是便可以通......
  • 数论 莫比乌斯反演
    前置需求数论分块概念对于一个形如\(\sum_{x=1}^n\lfloor{\frac{n}{x}}\rfloor\)的式子,我们发现对于一部分的\(x\),它们的\(\lfloor{\frac{n}{x}}\rfloor\)值相同,因此我们没必要\(\mathcal{O(n)}\)计算,可以采用数论分块的办法将这一步的复杂度降低至\(\mathcal{O(\sqrt......
  • 二项式反演学习笔记
    前言万字长文!这里有我的一些思考和领会,网络上的教程都太潦草了。并且我发现了新的反演公式!概述二项式反演用于转化两个具有特殊关系的函数\(f\)和\(g\),从而方便求解问题。一般来说,直接计算恰好满足\(n\)个限制的答案不好求,但是可以计算出“至少”/“至多”满足\(n\)......
  • 莫比乌斯反演入门
    来自这位大佬的视频的整理先整理几个重要的数论函数。1.莫比乌斯函数$\mu(n)$当\(n=1\)时取1,当\(n\)存在平方因子的时候取0,否则取\((-1)^k\),其中\(k\)是\(n\)所含的质因子数量。2.欧拉函数\(\phi(n)=\displaystyle\sum_{d=1}^n[gcd(d,n)=1]\),就是小于等于n且与\(n\)互质......
  • 【2】容斥与二项式反演
    【2】容斥与二项式反演1.1容斥原理容斥原理基于的是下面的恒等式:\[\sum\limits_{i=0}^n\dbinom{n}{i}(-1)^i=0^n=[n=0]\]这个式子有什么意义呢?我们考虑一个长度为\(N\)的序列,并且要求其中每个元素都满足某个限制,计算满足这个条件的序列数量。每个元素都满足限制\(\Leftri......
  • 9. 容斥与反演
    9.容斥与反演容斥原理:\[|\bigcup_{i=1}^nP_i|=\sum_{S\subseteqU}(-1)^{|S|-1}|\bigcap_{s\inS}P_s|\]感性理解:\(P_i\):”满足某种性质的元素的集合“;左边:具有任意一种性质的元素的并,右边:至少具备多个性质的元素。证明:考虑一个元素\(x\),设其包含在\(k\)个集合内,那么当......
  • 2024牛客多校第九场 C.Change Matrix 欧拉反演
    这题是欧拉反演的应用,之前没学过欧拉函数和欧拉反演,傻傻对着\(gcd(i,j)\)不知道怎么化简。首先对原来的矩阵进行转化,拆成\(n\)个小矩阵因为\(gcd(i,j)=\sum_{x|i,x|j}\phi(x)\)这是因为对于任意的正整数\(n\)都有\(n=\sum_{d|n}\phi(d)\),证明见oiwiki:https://oi-wi......