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

     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;
}