矩阵树定理

【行列式】

【性质】

1.单位矩阵的行列式为1

2.交换两行要变号

3.有两行一样的话行列式为0

用上面的第二条,将这两行交换,加了符号实际上却未变,就是0

4.线性性

5.用矩阵的一行加上另一行的倍数,行列式不变

【求法】

利用高斯消元的方式,再加上上面的性质加以辅助即可

我们对于 A 进行消元,然后记录操作对于行列式的影响,最后得到上三角矩阵,行列式就是对角线乘积,如果消不出来那么返回 0

【矩阵树定理】

矩阵树定理是用来求一个无向图有多少生成树的方法

设这个图的邻接矩阵为G,度数矩阵为D(只有对角线上有值,为d[i])

定义一个图的Kirchhoff矩阵K=D−A

可以证明这个矩阵去掉任意一行和一列,行列式的值不变

且这个值就是生成树的个数

【例题1】P4111 [HEOI2015]小 Z 的房间

sol.

【习题1】P4336 [SHOI2016]黑暗前的幻想乡

这个题是无向图的板子+简单的一点容斥 sol.

【带权的矩阵树定理】

正常的带权的矩阵树求的是所有生成树的边权乘积之和

构建方法就是把邻接矩阵的1改成对应边权即可

【例题2】P3317 [SDOI2014]重建

这道题目就是普通的利用带权的性质,构造了一下邻接矩阵的值

sol.

【有向图矩阵树定理】

这个就是简单改造一下,边的邻接矩阵不变

度数矩阵分类讨论,外向树度数只统计入度,内向图度数只统计出度

然后注意这里这就不能任意去掉一行一列了,一定要去掉起始点(root)的那一行和列

【例题3】P4455 [CQOI2018]社交网络

sol.

【变形——带权矩阵树定理】

这是一类特别的问题,求所有生成树的边权和的和

不同于一般的生成树带权的积的和,我们无法直接利用普通的行列式计算了

这时,我们就可以考虑利用计算贡献的方式

先算出整张图的生成树个数tot

然后枚举m条边,删除它求一下生成树个数cnt,我们知道当前边的贡献为$(tot-sum)*w[i]$

这样的做时间复杂度较大,为$O(n^3m)$

我们继续考虑优化

我们可以自行构造行列式内记录的信息来维护,我们给行列式每个位置赋值(x,y)

表示方案数为x,所有方案的边权的和的和为y

然后我们依然是计算乘法

(a,b)*(c,d)=(a*c,a*d+b*c)

这样我们第二维度的数就是答案

考虑我们在高斯消元的时候需要除法,根据乘法定义(a,b)*(c,d)=(a*c,a*d+b*c),推出(a*c,a*d+b*c)/(a,b)=(c,d)然后计算就行了,这个的没什么实际意义,就相当于一个逆运算

【例题4】P6624 [省选联考 2020 A 卷] 作业题

这道题需要简单莫比乌斯反演,枚举一下因数,然后就是裸的了

提交记录

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
矩阵定理是一种基于图的行列式计算算法,可以用来判断两个图是否同构。下面是使用矩阵定理判断两个图是否同构的具体流程: 1. 对于两个图G1和G2,先计算它们的拉普拉斯矩阵L1和L2。 2. 对于拉普拉斯矩阵L1和L2,分别选择一个相同的大小的行列子矩阵,称为L1'和L2'。 3. 计算L1'和L2'的行列式,如果它们相等,则认为G1和G2同构。否则,认为它们不同构。 下面是使用C语言代码实现矩阵定理判断两个图是否同构的示例: ``` #include <stdio.h> #include <stdlib.h> #include <string.h> // 定义矩阵大小 #define MAX_N 100 #define MAX_M 5000 // 定义拉普拉斯矩阵 int L1[MAX_N][MAX_N], L2[MAX_N][MAX_N]; // 定义行列式计算函数 int determinant(int n, int mat[MAX_N][MAX_N]) { int det = 1; for (int i = 0; i < n; i++) { for (int j = i + 1; j < n; j++) { while (mat[j][i]) { int t = mat[i][i] / mat[j][i]; for (int k = i; k < n; k++) { mat[i][k] = (mat[i][k] - t * mat[j][k]); // 交换i,j两行 int temp; temp = mat[i][k]; mat[i][k] = mat[j][k]; mat[j][k] = temp; } det = -det; } } det *= mat[i][i]; } return det; } int main() { // 读入两个图 int n, m1, m2; scanf("%d%d%d", &n, &m1, &m2); memset(L1, 0, sizeof(L1)); memset(L2, 0, sizeof(L2)); for (int i = 0; i < m1; i++) { int u, v; scanf("%d%d", &u, &v); L1[u - 1][u - 1]++; L1[v - 1][v - 1]++; L1[u - 1][v - 1]--; L1[v - 1][u - 1]--; } for (int i = 0; i < m2; i++) { int u, v; scanf("%d%d", &u, &v); L2[u - 1][u - 1]++; L2[v - 1][v - 1]++; L2[u - 1][v - 1]--; L2[v - 1][u - 1]--; } // 计算行列式 int k = n / 2; int det1 = determinant(k, L1); int det2 = determinant(k, L2); // 判断是否同构 if (det1 == det2) { printf("Same\n"); } else { printf("Different\n"); } return 0; } ``` 在上面的代码中,我们先输入两个图的节点数、边数和边的信息,然后构造它们的拉普拉斯矩阵。接着,我们选择相同大小的行列子矩阵,用行列式计算函数计算它们的行列式,最后比较两个行列式是否相等,从而判断它们是否同构。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值