矩阵树 Matrix-Tree 定理实现模板(高斯消元求解行列式)
生活随笔
收集整理的這篇文章主要介紹了
矩阵树 Matrix-Tree 定理实现模板(高斯消元求解行列式)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
大佬1博客:https://www.cnblogs.com/zj75211/p/8039443.html
大佬2博客:https://www.cnblogs.com/yangsongyi/p/10697176.html
三個定理:
對于上述三個定理,首先都是需要構造基爾霍夫Kirchhoff矩陣,構造方法如下:
基爾霍夫Kirchhoff矩陣 K?= 度數矩陣 D - 鄰接矩陣 A
對于上述三個定理,分別求解:
最后掛個高斯消元求解行列式的模板,時間復雜度為 n^3 ,空間復雜度為 n^2,允許取模
代碼:
const int N=110;const int mod=998244353;LL a[N][N];LL q_pow(LL a,LL b) {LL ans=1;while(b){if(b&1)ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans; }LL Gauss(int n) {LL ans=1;for(int i=2;i<=n;i++){for(int j=i+1;j<=n;j++)if(!a[i][i]&&a[j][i]) {ans=-ans;swap(a[i],a[j]);break;}LL inv=q_pow(a[i][i],mod-2);for(int j=i+1;j<=n;j++){LL temp=a[j][i]*inv%mod;for(int k=i;k<=n;k++)a[j][k]=(a[j][k]-a[i][k]*temp%mod+mod)%mod;}ans=ans*a[i][i]%mod;}return ans; }void addedge(int x,int y) {a[x][x]=(a[x][x]+1)%mod;a[y][y]=(a[y][y]+1)%mod;a[x][y]=(a[x][y]-1+mod)%mod;a[y][x]=(a[y][x]-1+mod)%mod; }?
總結
以上是生活随笔為你收集整理的矩阵树 Matrix-Tree 定理实现模板(高斯消元求解行列式)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: HDU多校6 - 6831 Fragra
- 下一篇: HDU多校6 - 6836 Expect