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

莫比乌斯反演学习笔记

时间:2024-04-02 21:23:28浏览次数:15  
标签:乌斯 sum mu 反演 莫比 cases

莫比乌斯反演学习笔记

前言

之前学了一遍,只学了朴素的莫比乌斯反演,现在第二次面对不知道能否有所长进。

性质

莫比乌斯反演是数论中的重要内容。对于一些函数 \(f(n)\),如果难以直接求出它的值,但容易求得其倍数和或约数和 \(g(n)\),那么可以通过莫比乌斯函数反演简化运算,从而求得 \(f(n)\) 的值。

前置知识

数论分块+狄利克雷卷积。

莫比乌斯函数

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

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

即对于 \(n=\prod p_i^{\alpha_i}\):

\[\mu(n)= \begin{cases} 1&n=1\\ 0&\exists \alpha_i\ge 2\\ (-1)^i&\forall\alpha_i=1\\ \end{cases} \]

容易看出,莫比乌斯函数是积性函数,同时满足以下性质:

\[\sum_{d|n}\mu(d)=\begin{cases} 1&n=1\\ 0&n\ne 1 \end{cases} \]

求解

由于莫比乌斯函数是积性函数,所以可以采用线性筛求解。

代码

void get_miu(int N){
	miu[1]=1;
	for(int i=2;i<=N;i++){
		if(!vis[i]){
			p[++cnt]=i;
			miu[i]=-1;
		}
		for(int j=1;p[j]<=N/i;j++){
			vis[i*p[j]]=1; 
			if(i%p[j]==0)break;//说明其含有平方因子,不做处理
			else miu[i*p[j]]=-miu[i];//否则根据积性函数的性质求解
		}
	}
    return ;
}

莫比乌斯变换

有以下两种形式,即文章开头提及的约数和和倍数和的形式。

约数和:如果 \(F(n)=\sum_{d|n}f(d)\),那么 \(f(n)=\sum_{d|n}\mu(\frac{n}{d})F(d)\)。

倍数和:如果 \(F(n)=\sum_{n|d}f(d)\),那么 \(f(n)=\sum_{n|d}\mu(\frac{d}{n})F(d)\)。

标签:乌斯,sum,mu,反演,莫比,cases
From: https://www.cnblogs.com/DycBlog/p/18111525

相关文章

  • 反演问题求解:基于MATLAB的反演问题求解算法实现和应用,包括反演问题数值模拟、反演问题
    鱼弦:公众号【红尘灯塔】,CSDN内容合伙人、CSDN新星导师、全栈领域优质创作者 、51CTO(Top红人+专家博主) 、github开源爱好者(go-zero源码二次开发、游戏后端架构https://github.com/Peakchen)基于MATLAB的反演问题求解:原理、应用、实现与分析反演问题是指由间接观测数......
  • *【莫比乌斯反演】数表[SDOI2014]
    问题有一张\(N\timesN\)的数表(\(N=10^5\)),其第\(i\)行第\(j\)列(\(1\lei\len\),\(1\lej\lem\))的数值为能同时整除\(i\)和\(j\)的所有自然数之和。有T次询问,每次询问给定\(n,m,A\),计算数表(1,1)至(n,m)中不大于\(A\)的数之和(\(|A|\le10^9\))。每组数据输出一行一个整数......
  • 3101: *【莫比乌斯反演:练习】gcd(i,j)=k的对数[POI2007]ZAP-Queries
    问题给出\(n,m,k\)(\(1\leqn,m,k\leq10^5\)),求\(\sum\limits_{i=1}^n\sum\limits_{j=1}^m\lbrack\gcd(i,j)=k\rbrack\),即:满足\(1\leqi\leqn\),\(1\leqj\leqm\),且\(\gcd(i,j)=k\)的二元组\((i,j)\)的数量。题解\(\sum\limits_{i=1}......
  • 2024 年春节集训 _ 第三课 - 莫比乌斯反演
    练习\(5\)\(\color{orange}{\texttt{E->link}}\)求\[\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\]\(n,m\leq10^7,\T\leq10^4\)贴个照片。及其丑陋的照片(我的草稿)如上化简最后可以得到\[\sum_{d=1}^nd\sum_{k=1}^{\left[\dfrac{n}{d}\right]}\mu(k)k^2......
  • 2024 年春节集训 _ 第二课 - 莫比乌斯反演
    练习\(5\)\(\color{orange}{\texttt{E->link}}\)求\[\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\]\(n,m\leq10^7,\T\leq10^4\)贴个照片。及其丑陋的照片(我的草稿)如上化简最后可以得到\[\sum_{d=1}^nd\sum_{k=1}^{\left[\dfrac{n}{d}\right]}\mu(k)k^2\......
  • 单位根反演小记
    核心式子\[[k|n]=\dfrac1k\sum_{i=0}^{k-1}\omega_k^{i\cdotn}\]证明:当\(n\)是\(k\)的因数时,\(\dfrac1k\sum\limits_{i=0}^{k-1}\omega_k^{i\cdotn}=\dfrac1k\sum\limits_{i=0}^{k-1}\omega_k^0=\dfrac1k\cdotk=1\)当\(n\)不是\(k\......
  • 莫比乌斯反演
    积性函数:若函数\(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)\)均为积性函数,那么则以下函数也为......
  • 莫比乌斯反演学习笔记
    莫比乌斯反演目录莫比乌斯反演反演公式&性质例题[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}......