日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

矩阵树 Matrix-Tree 定理实现模板(高斯消元求解行列式)

發布時間:2024/4/11 编程问答 43 豆豆
生活随笔 收集整理的這篇文章主要介紹了 矩阵树 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

    對于上述三個定理,分別求解:

  • D矩陣:度數,A矩陣:無向圖,刪除任意第 r 行和第 r 列后的行列式的絕對值即為答案
  • D矩陣:入度,A矩陣:有向圖,假設點為 x ,則刪除第 x 行和第 x 列后的行列式的絕對值即為答案
  • D矩陣:出度,A矩陣:有向圖,假設點為 x ,則刪除第 x 行和第 x 列后的行列式的絕對值即為答案
  • 最后掛個高斯消元求解行列式的模板,時間復雜度為 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 定理实现模板(高斯消元求解行列式)的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。