首页 > 其他分享 >莫比乌斯反演

莫比乌斯反演

时间:2024-02-27 18:00:28浏览次数:19  
标签:frac 函数 积性 乌斯 sum mu 反演 莫比 aligned

积性函数:

若函数 \(f(x)\) 满足:\(f(1)=1\) 且 \(∀x,y∈N_+,\gcd(x,y)=1\) 都有 \(f(xy)=f(x)f(y)\),则称它为积性函数。
若函数 \(f(x)\) 满足:\(f(1)=1\) 且 \(∀x,y∈N_+\) 都有 \(f(xy)=f(x)f(y)\),则称它为完全积性函数。

性质:若 \(f(x),g(x)\) 均为积性函数,那么则以下函数也为积性函数:

\[\begin{aligned} h(x)&=f(x^p)\\ h(x)&=f^p(x)\\ h(x)&=f(x)g(x)\\ h(x)&=∑_{d|x}f(d)g(\frac{x}{d}) \end{aligned} \]

积性函数的例子:

  • 单位函数:\(ε(n)=[n=i]\) (完全积性)
  • 恒等函数:\(id_k(n)=n^k\),通常地,我们将 \(id_1(n)\) 记作 \(id(n)\) 。(完全积性)
  • 常数函数:\(1(n)=1\)(完全积性)
  • 除数函数:\(σ_k(n)=∑_{d|n}d^k\) ,通常地,我们将 \(σ0(n)\) 记作 \(d(n)\) 或者 \(τ(n)\),\(σ_1(n)\) 记作\(σ(n)\) 。
  • 欧拉函数:\(φ(n)=∑^n_{i=1}[gcd(i,n)=1]\)
  • 莫比乌斯函数:\(\mu(n)=\begin{cases}\begin{aligned}&1& &n = 1\\&0& &\exists \ d>1,d^2|n \\&(-1)^{\omega(n)} && \text{otherwise}\end{aligned}\end{cases}\)

迪利克雷卷积(Dirichlet卷积)

定义

卷积专指的是函数见乘法运算,迪利克雷卷积是其中一种。对于两个数论函数 \(f(x)\)
和 \(g(x)\),将他们进行迪利克雷卷积,结果 \(h(x)\) 定义为:

\[h(x)=∑_{d|x}f(d)g(\frac{x}{d})=∑_{ab=x}f(a)g(b) \]

可以简写成:

\[h=f*g \]

性质

  • 交换律:\(f*g=g*f\)
  • 结合律:\((f*g)*h=f*(g*h)\)
  • 分配律:\((f+g)*h=f*h+g*h\)
  • 等式的性质:\(f=g\) 的充要条件:\(f*h=g*h\),其中 \(h(x)\) 满足 \(h(1)≠0\)。

单位元:单位函数 \(ε\) 是迪利克雷卷积运算中的单位元,即:对于任意函数 \(f\),存在 \(ε*f=f\)。
逆元:对于任意一个满足 \(f(x)≠0\) 的数论函数,假设存在一个 \(g(x)\),满足 \(g*f=ε\),则称\(g(x)\) 是 \(f(x)\) 的逆元。并且有且仅有 \(1\) 个。

莫比乌斯反演

定义 \(μ\) 为莫比乌斯函数,定义:

\[μ= \begin{cases} \begin{aligned} &1 & &n = 1\\ &0 & &\text{n 含有平方因子}\\ &(-1)^k & & \text{k 为 n 本质不同质因子个数} \end{aligned} \end{cases} \]

性质

\[∑_{d|n}μ(d)=\begin{cases}1 &n=1\\ 0 &n\ne 1 \end{cases} \]

即:\(∑_{d|n}μ(d)=ε(n),μ*1=ε\)

\[[\gcd(i,j)=1]=∑_{d|\gcd(i,j)}μ(d) \]

\(\text{Description}\):

已知 \(F_{m} = \sum_{a_1=1}^{m} \sum_{a_2=1}^{m} \cdots \sum_{a_n=1}^{m} [\gcd(a_1,a_2,\cdots a_n)=1]\) ,求 \(\sum_{i=1}^{k} F_i\),\(1\le k,n\le 2\times 10^6\)

\(\text{Solution}\):

考虑:

\[\begin{aligned} S &= \sum_{i=1}^n \sum_{j=1}^m [\gcd(i,j)=1]\\ &= \sum_{i=1}^n \sum_{j=1}^m \sum_{d|\gcd(i,j)} \mu(d) \\ &= \sum_{i=1}^{n} \sum_{d|i} \mu(d) \times \lfloor \frac{m}{d} \rfloor\\ &= \sum_{d=1}^{n} \mu(d) \times \lfloor \frac{m}{d} \rfloor \times \lfloor \frac{n}{d} \rfloor \end{aligned} \]

所以有

\[F_k = \sum_{d=1}^k \mu(d) \times {\lfloor \frac{k}{d} \rfloor} ^ n \]

直接整除分块是 \(k \times \sqrt k \log n\) 的,注意到:

\[\lfloor \frac{k}{d} \rfloor \ne \lfloor \frac{k-1}{d} \rfloor \Leftrightarrow d | k \]

我们发现 \(F_k\) 在每 \(d\) 个段都会发生变化,而在段内不变。考虑维护差分数组:

\[\begin{aligned} \Delta f &= \sum_{d|k}^k \mu(d) \times ({\lfloor \frac{k}{d} \rfloor}^n - {\lfloor \frac{k-1}{d} \rfloor}^n) \end{aligned} \]

然后做完了,复杂度 \(O(n \times (\log n + \log k))\)


线性筛求 \(\mu(x)\)


int pi[N], Ip[N], mu[N], tot;
void get() {
	mu[1] = 1;
	for(int i = 2; i <= lim; ++i) {
		if(!Ip[i]) {
			pi[++tot] = i;
			mu[i] = -1;
		}
		for(int j = 1; j <= tot; ++j) {
			if(pi[j] * i > lim) break;
			Ip[i * pi[j]] = 1;
			if(i % pi[j] == 0) {
				mu[i * pi[j]] = 0;
				break;
			}
			mu[i * pi[j]] = -mu[i];
		}
	}
}

标签:frac,函数,积性,乌斯,sum,mu,反演,莫比,aligned
From: https://www.cnblogs.com/Saka-Noa/p/18037459

相关文章

  • 莫比乌斯反演学习笔记
    莫比乌斯反演目录莫比乌斯反演反演公式&性质例题[HAOI2011]ProblembYY的GCD于神之怒加强版Crash的数字表格/JZPTAB[SDOI2014]数表[SDOI2015]约数个数和反演公式&性质\[f(n)=\sum_{d|n}g(d)\\g(n)=\sum_{d|n}\mu(d)f(\fracnd)\]感觉我不太会用上面那个我只会用莫比乌斯函......
  • 【数论】卷积反演大集合
    不知道为啥脑抽要学数论,骂声一片中发现数论还没入门(悲)。1.狄利克雷卷积与数论函数1.1数论函数定义:数论函数为值域为整数的函数。简单数论函数:\(I(n)\),恒等函数,恒等为\(1\)。\(e(n)\),元函数,卷积中的单位元,若\(n=1\),\(e(n)=1\)。否则为\(e(n)=0\)。\(id(n)\),单位函数,\(......
  • 单位根反演
    单位根反演通常用于求\(\sum\limits_{i=0}^n[x\midi]f_i\)。形式\[[n\midk]=\frac{1}{n}\sum\limits_{i=0}^{n-1}\omega_n^{ik}\]其中\(\omega_n\)是\(n\)次单位根,模意义下可以被原根替换。证明当\(n\midk\)时,\(\omega_n^{ki}=1\)。所以右边等于\(\frac{1}{n}......
  • 莫比乌斯反演
    数论分块引理\[\lfloor\frac{a}{bc}\rfloor=\lfloor\frac{\lfloor\frac{a}{b}\rfloor}{c}\rfloor\]数论分块对于\(\displaystyleh(i)=\lfloor\frac{n}{i}\rfloor\),设\(f(l)=x\)。则\(\displaystyle\foralli\in[l,\lfloor\frac{n}{\lfloor\frac{n}{l......
  • 单位根反演学习笔记
    单位根反演\[[n|a]=\dfrac{1}{n}\sum\limits_{k=0}^{n-1}w^{ak}_{n}\]证明:当\(i\not\equiv0\pmodn\)时,由等比数列求和公式可得:原式\(=\dfrac{1}{n}\times\dfrac{w^{an}_n-1}{w^a_{n}-1}\),而右边分子为0,分母不为0,因此和为0。当\(i\equiv0\pmodn\)时,有原式\(=\dfrac{1}{n}......
  • 容斥与反演
    目录容斥与反演容斥原理容斥原理广义容斥原理P3813[FJOI2017]矩阵填数ProblemSolution反射容斥2023.11.16T1ProblemSolutionP3266[JLOI2015]骗我呢ProblemSolution集合反演P3349[ZJOI2016]小星星ProblemSolution[AGC005D]~KPermCountingProblemSolutionP5644[PKUWC2018]......
  • 莫比乌斯反演——优美地解决容斥问题
    莫比乌斯反演假设现在有一个函数\(f\),令\(F(n)=\sum\limits_{d|n}f(d)\),如\(F(1)=f(1),F(4)=f(1)+f(2)+f(4)\),现在我们要通过\(F\)反推\(f\),如\(f(1)=F(1),f(4)=F(4)-F(2)\),这就是莫比乌斯反演。可以推出这样的公式:\(F(n)=\sum\limits_{d|n}f(d......
  • <学习笔记> 二项式反演
    二项式反演证明我们设\(g(x)\)为任意\(x\)个集合的交集的大小,\(f(x)\)表示任意\(x\)个集合补集的交集大小。首先有(组合数学6.2)\[|\overline{S_1}\cap\overline{S_2}\cap...\cap\overline{S_{n-1}}\cap\overline{S_n}|=|U|-|{S_1}|-|{S_2}|-...+(-1)^n\times|{S......
  • 莫比乌斯反演
    莫比乌斯反演补了补暑假欠下的账(你怎么寒假才学)推狮子>>写代码。数论函数:定义域为正整数的函数。积性函数,对于一个数论函数,任意两个互质的正整数\(x,y\),都有\(f(xy)=f(x)f(y)\)完全积性函数就是不要求\(x,y\)互质的积性函数。常见的积性函数:单位函数\(\epsilon(n)......
  • 莫比乌斯反演
    前置:积性函数与狄利克雷卷积和整除分块两个基础积性函数:\(\varepsilon(n)=[n=1]\),\(1(n)=1\)。性质:\(\varepsilon*f=f\),\(f\)是任意函数。结论:\(f(n)\)是积性函数\(\iffg(n)=\displaystyle\sum_{d|n}f(d)\)是积性。证明:$\Rightarrow$方向:\(g=f*1\),狄利克雷卷积的性......