首页 > 其他分享 >莫比乌斯反演常见的三个模型

莫比乌斯反演常见的三个模型

时间:2023-05-11 16:44:10浏览次数:36  
标签:lfloor frac gcd 乌斯 sum rfloor 反演 莫比 td

莫比乌斯反演常见模型

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

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

模型2:\(\sum_{i=1}^n\sum_{j=1}^mij[gcd(i,j)=t]\)

\[\begin{aligned} \sum_{i=1}^n\sum_{j=1}^mij[gcd(i,j)=t]&= \sum_{d=1}^{min(n,m)}\mu(d)\sum_{d|i}^{\lfloor\frac{n}{t}\rfloor}\sum_{d|j}^{\lfloor\frac{m}{t}\rfloor}ij\\ &=\frac{\sum_{d=1}^{min(n,m)}\lfloor\frac{n}{td}\rfloor\lfloor\frac{m}{td}\rfloor(\lfloor\frac{n}{td}\rfloor+1)(\lfloor\frac{m}{td}\rfloor+1)}{4} \end{aligned} \]

至于推导的话就是等差数列求和而已

模型3:\(\sum_{i=1}^n\sum_{j=1}^mf[gcd(i,j)=t]\)

这个的推导就总结向上面一样照葫芦画瓢推就行了

答案如下:

\[\sum_{u=1}^{min(n,m)}\lfloor\frac{n}{u}\rfloor\lfloor\frac{m}{u}\rfloor\sum_{d|u}f(d)\mu(\frac{u}{d}) \]

注意所有的求和都要用到整除分块

标签:lfloor,frac,gcd,乌斯,sum,rfloor,反演,莫比,td
From: https://www.cnblogs.com/Zhangrx-/p/17391543.html

相关文章

  • 莫比乌斯反演学习笔记(内涵反演详细推导和证明!)
    莫比乌斯反演前言记得上次学这个东西已经是一年前了,题到没有做过几道……一些有趣的东西莫比乌斯这个人还是很nb的,相信大家都听说过莫比乌斯带,就是他发现的,跟所有数学大佬一样的,他还喜欢天文学和物理废话不扯多了前置知识整除分块一个在数论计算中异常常见的东西,一般和......
  • 拉格朗日反演公式(lagrange inversion)组合证明
    Thereisasimplecombinatorialproof.Theoriginalformis\[[t^n]w^k=\frac{k}{n}[t^{n-k}]\phi^k\]where\(w=t\phi(w)\)consider\(w\)asegf.ofthewaysofsometrees.\(\phi\)asageneratingruleconcerningdegree.\[n![x^n]\frac{w^k}{k......
  • 单窗算法的地表温度反演:谷歌地球引擎GEE代码
      本文介绍在GEE中基于Landsat遥感影像实现地表温度(LST)单窗算法反演的代码。1背景知识  基于遥感数据的地表温度(LST)反演目前得到了广泛的应用,尤其是面向大尺度、长时间范围的温度数据需求,遥感方法更是可以凸显其优势。目前,基于各类遥感数据源的地表温度反演方法不断得以改......
  • 二项式反演
    反演就是对于两个整数函数\(f\)和\(g\),从用\(g\)表示\(f\)转化为用\(f\)表示\(g\)。简而言之,\(f(n)\)是\(g(0),g(1),\cdots,g(n)\)的一个线性组合,那么很明显,有\(f(n)=\sum_{i=0}^na_{n,i}g(i)\)。如果把\(g(i),f(i)\)用向量\(G,F\)表示,那么\(F=\{a_{i,j}\}*G......
  • 莫比乌斯反演
    反演我们再来看看容斥原理实质上发生了什么——根据定义我们有\[N_\geq(S)=\sum\limits_{S\subseteqJ}N_=(J)\]而容斥原理(一般形式)表明\[N_=(S)=\sum\limits_{J:S\subseteqJ\subseteq\mathscr{A}}(-1)^{|J|-|S|}N_\geq(J)\]也就是说,容斥原理其实是由\((1)\)式解出了\((......
  • 二项式反演
    若\[g_n=\sum_{i=0}^n\dbinom{n}{i}f_n\]有\[f_n=\sum_{i=1}^{n}\dbinom{n}{i}g_n\]若\[g_i=\sum_{j=i}^{n}\dbinom{j}{i}f_j\]则\[f_i=\sum_{j=i}^{n}\dbinom{j}{i}(-1)^{j-i}g_j\]P4859已经没有什么好害怕的了给两个数列\(a\),\(b\),要求\(a_i,b_i\)两两匹配,......
  • 线性筛,整除分块,欧拉函数与莫比乌斯反演
    埃氏筛法说到筛质数,就不得不提到大名鼎鼎的埃氏筛法了思想非常简单,就是对于每一个素数pri[i],我们都把它的倍数筛去,譬如对于素数2来说,我们就把\(2*2,2*3,2*4,2*5....2*n\)的数全部筛掉代码voidzhishu(intn){ for(inti=2;i<=n;i++){ if(p[i]==0) for(intj=i*2;......
  • 莫比乌斯反演
    说到筛质数,就不得不提到大名鼎鼎的埃氏筛法了埃氏筛法思想非常简单,就是对于每一个素数pri[i],我们都把它的倍数筛去,譬如对于素数2来说,我们就把\(2*2,2*3,2*4,2*5....2*n\)的数全部筛掉代码voidzhishu(intn){ for(inti=2;i<=n;i++){ if(p[i]==0) for(intj=i*2;......
  • 有关莫比乌斯函数
    前置知识:整除分块结论:集合\(\{x|\existsi\in\mathbbN^*,\lfloor\fracni\rfloor=x\}\)中的元素数量为\(\sqrtn\)级别个。具体来说,不超过\(2\sqrtn\)个。证......
  • 容斥与反演技巧(二)
    年更大作FJOI2017矩阵填数\(\max=v\)拆成\(\lev\)和存在一个\(=v\),对第二个条件容斥发现变成\(<v\),离散化后对每个点求出上界直接算。复杂度\(O(n^32^n)\)。......