首页 > 其他分享 >二项式反演

二项式反演

时间:2023-04-25 20:45:31浏览次数:32  
标签:BF sum 矩阵 反演 binom 二项式

反演就是对于两个整数函数 \(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\)

而反演的实质就是想要找到一个下三角矩阵 \(B=\{b_{i,j}\}\) 满足 \(G=\{b_{i,j}\}*F\),也就是 \(g(n)=\sum_{i=0}^nb_{n,i}f(i)\)

记 \(A=\{a_{i,j}\}\),那么有 \(F=AG,G=BF\),提一下,\(F=A(BF)\),\(F=(AB)F\),所以 \(AB=E\),所以 \(B=A^{-1}\)。

所以实际上我们就是要找到 \(A\) 的逆矩阵。

求取这个逆矩阵其实很简单,我们只需要在想要的 \(a_{i,j}\) 下,把特殊的 \(f(n)\) 和 \(g(n)\) 带入进去,得到 \(b_{i,j}\) 就可以得到对普适的 \(f(n),g(n)\) 的结论。

例如,二项式反演,\(f(n)=\sum_{i=0}^n\binom{n}{i}g(i)\),那么我们令 \(g(i)=(x-1)^i\),那么根据二项式定理,\(f(n)=x^n\),接下来,\((x-1)^n=\sum_{i=0}^n b_{n,i}x^i\),很明显,因为 \((x-1)^n=\sum_{i=0}^n (-1)^{n-i}\binom{n}{i}x^i\),所以 \(b_{n,i}=(-1)^{n-i}\binom{n}{i}\)。

我们就得到了 \(g(n)=\sum_{i=0}^n(-1)^{n-i}\binom{n}{i}f(i)\)。

标签:BF,sum,矩阵,反演,binom,二项式
From: https://www.cnblogs.com/jucason-xu/p/17353783.html

相关文章

  • 莫比乌斯反演
    反演我们再来看看容斥原理实质上发生了什么——根据定义我们有\[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;......
  • 容斥与反演技巧(二)
    年更大作FJOI2017矩阵填数\(\max=v\)拆成\(\lev\)和存在一个\(=v\),对第二个条件容斥发现变成\(<v\),离散化后对每个点求出上界直接算。复杂度\(O(n^32^n)\)。......
  • 莫比乌斯反演
    莫比乌斯反演引入莫比乌斯反演用处:对于一些函数\(f(n)\),如果比较难以求出它的值,但容易求出其倍数和或约束和\(g(n)\),则可以通过莫比乌斯反演简化运算。莫尼乌斯函数......
  • 二项式反演
    学习参照cmd的博客,知乎,oi-wiki,某神仙的博客组合恒等式\[\binom{n}{k}=\binom{n-1}{k}+\binom{n-1}{k-1}\\\binom{n}{n_1,n_2,\ldots,n_k}=\frac{(n_1+n_2+\ldots+n_k......
  • 【应用】Lagrange 反演应用
    证明鸽了,所以先开始应用篇。对于一元多项式\(F,G\)我们有Lagrange反演公式:\[n[x^n]F^k=k[w^{-k}]G^{-n}\]绝大多数情况我们都取\(k=1\)。其中多项式\(G\)为\(F......
  • 【洛谷】P2257 YY的GCD(莫比乌斯反演)
    原题链接题意\(T\)组询问,每次询问求:\[\sum_{i=1}^{n}\sum_{i=1}^{m}[\gcd(i,j)\inprime]\]\(T=10^4,n,m\leq10^7\)。思路不难想到枚举质数,将原式化简为:\[\sum......
  • 莫比乌斯反演 & 狄利克雷卷积
    大家好,我不会数学实锤了。文章内容较杂,分章节叙述了的大部分有关内容。为什么把这俩放一起?我不知道。积性函数积性函数:\(\foralla,b\),\(a\perpb\),如果一个函数\(f\)......