部分经典反演的推导(附矩阵形式)

反演:已知有对向量 u 的线性变换:Au=v。现在已知 v 要反解 u 的过程称为反演。

反演相当于做线性变换:A1v=u,于是一般我们需要知道 A1 的形式。

部分转置问题:由于 (AT)1=(A1)T,于是 ATu=v(A1)Tv=u

:下文 Ai,j1 表示 (A1)i,j(Ai,j)1 表示求逆。部分除法用 /

写成数列形式:fi=jAi,jgjgi=jAi,j1fj

转置形式:fi=jAj,igjgi=jAj,i1fj

或者你可以理解为下面要对一些特殊矩阵写出其求逆后的形式。

子集反演#

用这个简单点的引入。下面集合符号均在二进制下考虑,线性变换矩阵:Ai,j=[ji],其中 i,j[0,2n)

注意到这其实就是 or 卷积的 DFT 矩阵。反演矩阵 A1 就自然是 IDFT 矩阵了。

Ai,j1=(1)|i||j|[ji],其中 |i|=popcount(i)

写成数列形式应该更好理解:fS=TSgTgS=TS(1)|S||T|fT,就是一个简答的容斥。

转置形式:Ai,j=[ij]Ai,j1=(1)|j||i[ij]fS=STgTgS=ST(1)|T||S|fT

莫比乌斯反演#

这个数列形式极其常见,先写:fn=dng(d)gn=dnμ(n/d)f(d)

太典啦,去 OI-wiki 自己学吧,懒得写啦。

矩阵形式:Ai,j=[ji]Ai,j1=μ(i/j)[ji]。转置形式略。

二项式反演#

线性变换矩阵为 Ai,j=(ij),即变换 gffn=i=0n(ni)gi

直接亮结论:逆矩阵 Ai,j1=(1)ij(ij)。数列形式:fn=i=0n(ni)gign=i=0n(1)ni(ni)fi。转置形式略。

代数证法#

(A×A1)i,j=k(ik)(kj)(1)kj=k=0ij(ik+j)(k+jj)(1)k=(ij)k=0ij(ijk)(1)k={1(i=j)0(ij)

于是 A×A1=I

多项式证法#

已知:fn=i=0n(ni)gi,写成 EGF 形式:fnn!=i=0ngii!×1(ni)!

设出 EGF:令 f(x)=i=0nfixii!,g(x)=i=0ngixii!,则有 f(x)=g(x)×(i0xii!)=exg(x)

于是 g(x)=exf(x)gn=i=0n(1)ni(ni)fi

斯特林反演#

线性变换矩阵为 Ai,j={ij},即变换 gffn=i=0n{ni}gi

直接亮结论:逆矩阵 Ai,j1=(1)ij[ij]。数列形式:fn=i=0n{ni}gign=i=0n(1)ni[ni]fi

类似的,有:Ai,j=[ij]Ai,j1=(1)ij{ij},数列形式:fn=i=0n[ni]gign=i=0n(1)ni{ni}fi

由于两者形式相似下文会证明第一种(有的是证明第一种转置形式)。

下面证明参考 dwt 博客mol dwt

多项式证法 1#

首先要介绍几个经典的斯特林数恒等式。首先考虑组合意义,有:

xn=k=0n{nk}xk_ xn¯=k=0n[nk]xk


同时做代换:xn¯=(1)nxn_,有:

xn=k{nk}(1)nkxk¯xn_=k[nk](1)nkxk


考虑转置命题:已知 fn=in{in}gi,考虑下降幂多项式:f(x)=i=0nfixi_,有:

f(x)=i=0nji{ji}gjxi_=j=0ngjij{ji}xi_=j=0ngjxj

提取两边 xn 系数,有:gn=in(1)in[in]fi。就证明了转置形式。

多项式证法 2#

考虑斯特林数 EGF,有:{nm}=[xnn!](ex1)mm![nm]=[xnn!](ln(1x))mm!

已知 fn=i=0n{ni}gi,设 F,G 分别为 f,gEGF,有:

F(x)=kfkxkk!=kxkk!i=0k{ki}gi=igiki{ki}xkk!=i=0gi(ex1)ii!=G(ex1)

G(x)=F(ln(1x))=i0filn(1x)ii!=i0(1)ifi(ln(1x))ii!=i0(1)ifiji1j![ji]xj=jxjj!ij(1)ifi[ji]

提取两边 [xnn!] 系数,有:(1)ngn=i=0n(1)i[ni]fi,把 (1)n 移项即可。

代数证法#

请拜读 dwt 博客,有空补,多项式证法好。

练习#

由于本人没有系统总结过,于是请参考 cmd 博客,也附一下 dwt 反演博客

posted @   HaHeHyt  阅读(289)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· AI与.NET技术实操系列(五):向量存储与相似性搜索在 .NET 中的实现
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 上周热点回顾(3.3-3.9)
点击右上角即可分享
微信分享提示
主题色彩
主题色彩