生成树的计数与基尔霍夫矩阵

来源:互联网 发布:mac桌面卸载图标 编辑:程序博客网 时间:2024/05/29 15:41

     jsoi08巨额奖金不会,看了07年集训队论文《生成树的计数与运用》,好像明白了点东西,不过清楚地发现我真是数学着急呀,行列式的计算与高斯消元完全被虐-_-! 

       在此仅写写他的操作吧。

    何谓基尔霍夫矩阵?我不明了,只晓得他可以计算任意一张图的生成树个数

    对于一张图首先构造他的基尔霍夫矩阵:

    Cij:     1.i==j,c[i][j]等于i的度数;

              2.i!=j:

                      1.i到j有k条边相连时,c[i][j]=-k;

                      2.i到j没有边相连时,c[i][j]=0;

     基尔霍夫矩阵的一个简单构图方式:用他的度数矩阵减去他的邻接矩阵(邻接矩阵要记录度数);

     构造完了就是计算了:

     将矩阵都第R行与R列同时去掉,得到一个n-1阶的主子式,他的值的绝对值就是原图生成树的个数。

     行列式的计算可以使用高斯消元,将其转化为三角形行列式,他主对角线的乘积即为行列式的值(被蹂躏呀-_-!!!)

        贴一段程序吧:

inline double hls(){ double a[200][200];  for (int i=1;i<=n-1;i++){  for (int j=i;j<=n;j++)   if (a[j][i]!=0){    for (int k=1;k<=n;k++)     swap(a[i][k],a[j][k]);    break;    }  for (int j=i+1;j<=n;j++){   double t=a[j][i]/a[i][i];   for (int k=i;k<=n;k++)    a[j][k]=a[j][k]-a[i][k]*t;   }  } double ans=1; for (int i=1;i<=n;i++)  ans*=a[i][i]; return ans;}

原创粉丝点击