【背诵手记】Matrix Tree定理和一些推广

结论

对于一个无向图\(G=(V,E)\),暂时钦定他是简单图,定义以下矩阵:

  • (入)度数矩阵\(D\),其中\(D_{ii}=deg_i\)。其他=0
  • 邻接矩阵\(A\),其中\(A_{ij}=[\exist e=(i,j)]\)。其他=0 (*******wait!*******)
  • 关联矩阵\(B\),其中\(B_{ij}=[\exist e_i=(a,b)](-1)^{[a>b]}\)。其他=0(后面会用到)
  • 拉普拉斯矩阵\(L=D-A\)

那么,这个图的生成树个数是\(L\)矩阵的任意一个n-1阶主子式的绝对值。

证明

关于拉普拉斯矩阵的性质

  • 拉普拉斯矩阵所有的代数余子式值相同

    证明:

    一个显然的事实是拉普拉斯矩阵任何一行/列的元素和\(=0\)

    我将证明一行内代数余子式\(R_{ij}=R_{ik}\)

    先只看对应位置的余子式,可以发现他们的不同在于一个有第\(j\)列没有第\(k\)列,另一个相反。

    根据之前发现的简单的事实\(L_{ik}+L_{ij}+\text{others}=0\) 。我们对于\(C'_{ik}\),其他列直接乘\(-1\)加到第\(k\)列对应列,此时第\(k\)列对应列就变成了 \(C_{ij}\)的第一列。

    然后考虑把这两个余子式特殊的那一列(j,k),列交换(用冒泡排序的方法,因为我们不能破坏其他列的相对关系)到第一列,每次交换是要乘上一个-1的,所以最终交换完成后,有一个类似\((-1)^{?}\)的系数。此时除了符号,两个行列式相等了。

    两个余子式都要乘上\((-1)^{?}\),而这里的系数正好就可以对应上代数余子式需要乘上的\((-1)^{i+j}\)。

    行内检验无异,列就自然成立,稍微推广QED,由上可知自闭证毕

  • 拉普拉斯矩阵的行列式的值\(=0\)

    根据行列式按_展开法则证明显然

  • 记\(R=R_{11}\),来证明它一些性质

    • \(|R|=0\)时当且仅当这个图不连通

      证明:

      通过交换行列将联通块排在一起,此时的行列式的值是原来的\((-1)^?\)倍,此时整个矩阵变成了这样:
      \[
      L=\begin{pmatrix}
      L_1 &O&...\\
      O&L_2&...\\
      ...&...&...
      \end{pmatrix}
      \]
      根据拉普拉斯矩阵性质\(|L_i|=0\)以及R的定义可知\(|R|=0\)。

    • \(||R||=1\)当且仅当\(G\)是一棵树

      证明:

      我暂时不会

      循环证明:

      显然一棵树的生成树只有一颗


重要补充:
\[
B^TB=L
\]

显然,\(B^TB\)相当于\(B\)中选两列\(e_i,e_j\)列向量然后一个转置后再乘起来,可以验证\(L_{ij}=e_i^Te_j\)。


重要补充:

柯西-比内公式(CB)

对于矩阵\(A[n][m]\)和\(B[m][n]\)有:
\[
\det(AB)=\sum_{|S|=n,S\subseteq\{1\dots m\}}\det\left(A_{\{1...n\},S}\right)\det\left(B_{S,\{1...n\}}\right)
\]
证明

看起来组合意义很强...

这一步好像是理解Matrix Tree本质/性质/拓展性 的关键一步,可惜我看不懂....


回到正题= =

我们把\(B\)的第一列剔除,显然\(R=B'^TB'\) ,套公式有
\[
|R|=\sum_{|S|=n-1,S\subseteq\{1\dots m\}}\det\left(B'^T_{\{1...n-1\},S}\right)\det\left(B'_{S,\{1...n-1\}}\right)
\]
激动人心的时刻到了,根据矩阵转置的性质可以发现
\[
|R|=\sum_{|S|=n-1,S\subseteq\{1\dots m\}}\det\left(B'_{S,\{1...n-1\}}\right)^2
\]
好诱人的组合意义啊。不过我们做一点准备工作,搞个记号\(N=\{1...n-1\}\)。

  • \(\det(B'_{S,N})^2 \in \{0,1\}\)

    可以发现的是,根据\(B'^T_{S,N}B'_{N,S}=R\),当且仅当选出的边不连通时,可以知道\(\det (B'^T_{S,N})=\det(B'_{N,S})=0\)。

    现在要证明选出的边联通时,\(||B'_{N,S}||=1\)。此时由于\(|S|=n-1\)那么一定是一颗树。

    口胡警告

    观察\(B'_{N,S}\)这个矩阵的形式,应该是由初等基本变换变换出来的,然后我们又有\(B'^T_{S,N}B'_{N,S}=R,||R||=1\),基本初等变换得到的矩阵相乘仍然可以看做基本初等变换,由于基本初等变化不会影响行列式的值,就有那么\(||B'_{N,S}||=1\) 。

现在回过头来考虑CB公式的算法,它相当于在\(O({m\choose n-1})\)时间内枚举每一个边集,然后判断是否联通,若联通就\(+1\)的贡献。显然最终累和的值就是生成树的个数。

也就是说|R|就等于这个图的生成树个数。

R是一个行列式,算行列式自然有好算的办法

啦啦啦终于证完了

下面的推广可能完全是假的。

推广

根据行列式的\(O(n!)\)的式子,可以知道\(|R|\)的值是关于\(|R_{ij}|\)的基本初等函数(加减乘除)。根据上面的推导,我们发现\(|R|\)是一堆\(1,-1\)的乘积的和。上面的推导中,不存在除法,虽然存在乘法但是是简单的取相反数运算。(注意,有平方运算!)我们如果将邻接矩阵的01换成边的边权\(w\),那么\(|R|\)就等于所有生成树的边权积(注意,边权\(\in R^+\))的和。

因为行列式本质上是只是一个多元函数\(F(\{\})=\sum \prod (-1)^{?}...\),我们甚至令w=一个生成函数甚至边权!!1...(毒瘤警告)

实际上我们计的数是,在\(m\)个三元组\((i,j,w)\)中选择\(n-1\)个出来,并且保证\((i,j)\)构成一棵树。我们求的是所有满足树条件的每棵树的\(\prod w\) 的和。

口胡警告

更抽象地讲,任何和域\(<R,+,\times>\)有着相同的代数结构的域\(<C,\and,*>\)都满足这个定理,可以得到\(|R|=\and *w\)。

注意到不论我们带什么进去都是好实现的,因为求行列式有\(O(n^3)\)的算法,而且...

这就触及到我的知识盲区了,理性证明等我大学看能不能学会,但是现在可以拍脑袋地讲,假如你有可以对应到加(减)乘(除)的运算法则的元素,那么你带进去没错。这是很对的。

所以:

  • 假如让边权\(w\)等于这条边重边数量(度数也要相应更新),那么\(|R|\)就等于考虑重边生成树个数。
  • 有向图的生成树的话,\(D\)就要钦定是入/出度矩阵,\(A\)就要钦定是\(A_{i,j/j,i}=cnt[e<i,j>]\),得到的\(|R_{ii}|\)是以\(i\)为根的内向树/外向树个数。正确性如果通过上面的证明理解的话貌似是完全错误的,因为就拉普拉斯矩阵的第一个性质就不满足...可能有别的构造出的组合意义?

至此,熟练背诵的Matrix Tree定理就出来了....

【证明与推广与背诵】Matrix Tree定理和一些推广的更多相关文章

  1. 数学-Matrix Tree定理证明

    老久没更了,冬令营也延期了(延期后岂不是志愿者得上学了?) 最近把之前欠了好久的债,诸如FFT和Matrix-Tree等的搞清楚了(啊我承认之前只会用,没有理解证明--),FFT老多人写,而Matri ...

  2. @总结 - 7@ 生成树计数 —— matrix - tree 定理(矩阵树定理)与 prüfer 序列

    目录 @0 - 参考资料@ @0.5 - 你所需要了解的线性代数知识@ @1 - 矩阵树定理主体@ @证明 part - 1@ @证明 part - 2@ @证明 part - 3@ @证明 part ...

  3. [bzoj1016][JSOI2008]最小生成树计数 (Kruskal + Matrix Tree 定理)

    Description 现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树中至少有一条边不同,则这两个最小生成树就是不同的 ...

  4. BZOJ.4031.[HEOI2015]小Z的房间(Matrix Tree定理 辗转相除)

    题目链接 辗转相除解行列式的具体实现? 行列式的基本性质. //864kb 64ms //裸的Matrix Tree定理.练习一下用辗转相除解行列式.(因为模数不是质数,所以不能直接乘逆元来高斯消元. ...

  5. SPOJ.104.Highways([模板]Matrix Tree定理 生成树计数)

    题目链接 \(Description\) 一个国家有1~n座城市,其中一些城市之间可以修建高速公路(无自环和重边). 求有多少种方案,选择修建一些高速公路,组成一个交通网络,使得任意两座城市之间恰好只 ...

  6. HDU 4305 Lightning Matrix Tree定理

    题目链接:https://vjudge.net/problem/HDU-4305 解法:首先是根据两点的距离不大于R,而且中间没有点建立一个图.之后就是求生成树计数了. Matrix-Tree定理(K ...

  7. BZOJ.4894.天赋(Matrix Tree定理 辗转相除)

    题目链接 有向图生成树个数.矩阵树定理,复习下. 和无向图不同的是,度数矩阵改为入度矩阵/出度矩阵,分别对应外向树/内向树. 删掉第i行第i列表示以i为根节点的生成树个数,所以必须删掉第1行第1列. ...

  8. BZOJ.1016.[JSOI2008]最小生成树计数(Matrix Tree定理 Kruskal)

    题目链接 最小生成树有两个性质: 1.在不同的MST中某种权值的边出现的次数是一定的. 2.在不同的MST中,连接完某种权值的边后,形成的连通块的状态是一样的. \(Solution1\) 由这两个性 ...

  9. [模板]Matrix Tree定理

    结论:一个图的生成树个数等于它的度数矩阵减邻接矩阵得到的矩阵(基尔霍夫矩阵)的任意一个n-1阶主子式的行列式的绝对值 证明:不会 求法:高斯消元 例题:[HEOI2013]小Z的房间 #include ...

随机推荐

  1. DataTable CAST 成集合后,进行自定义排序再转换回DataTable

    dt = dt.Rows.Cast<DataRow>().OrderBy(r => Convert.ToInt32(r["数量"])==0?1:0).ThenBy ...

  2. @noi.ac - 507@ 二分图最大权匹配

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 有一天你学了一个能解决二分图最大权匹配的算法,你决定将这个算法应 ...

  3. 弹性FLEX布局

    页面布局一直都是web应用样式设计的重点 我们传统的布局方式都是基于盒模型的 利用display.position.float来布局有一定局限性 比如说实现自适应垂直居中 随着响应式布局的流行,CSS ...

  4. HDU 1026 BSF+优先队列+记录路径、

    #include<iostream> #include<cmath> #include<cstring> #include<cstdio> #inclu ...

  5. 2014年NOIP普及组复赛题解

    题目涉及算法: 珠心算测验:枚举: 比例简化:枚举: 螺旋矩阵:模拟: 子矩阵:状态压缩/枚举/动态规划 珠心算测验 题目链接:https://www.luogu.org/problem/P2141 ...

  6. Codeforces Round #180 (Div. 1 + Div. 2)

    A. Snow Footprints 如果只有L或者只有R,那么起点和终点都在边界上,否则在两者的边界. B. Sail 每次根据移动后的曼哈顿距离来判断是否移动. C. Parity Game 如果 ...

  7. 获取exe和dll里面的资源

    有时候需要仿照另一个程序实现一些对话框,比较笨的办法是打开那个程序,照着样子自己在VC里面画啊画.这样的效率实在有点低. 现在有很多工具可以从exe和dll里面取出图片.图片.字符串.对话框等资源.比 ...

  8. 【React】 百度地图API

    百度地图 开发文档 :http://lbsyun.baidu.com/index.php?title=jspopular 调用接口 需要 内置加载一个 百度api文件    使用自己的ak  申请一个 ...

  9. C# 从零开始写 SharpDx 应用 笔刷

    本文告诉大家如何在 SharpDx 里面使用笔刷,包括纯色笔刷.渐变笔刷和图片笔刷 本文属于 SharpDx 系列 博客,建议从头开始读 初始化 本文将会在 C# 从零开始写 SharpDx 应用 初 ...

  10. 2019.12.15 QLU and SNDU期末联赛

    题目列表: 1582.柳予欣的舔狗行为 1587.柳予欣的女朋友们在分享水果 1585.柳予欣和她女朋友的购物计划 1579.FFFFFunctions 1588.Zeckendorf 1586.柳予 ...