基尔霍夫矩阵

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

    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阶的主子式，他的值的绝对值就是原图生成树的个数。
 

ll det(int n)    //生成树计数:Matrix-Tree定理   //这个是求行列式的精华啊(实质上这个子函数就是用来求行列式的值的!)(所以也适用于求行列式值的题)
{
    ll ret=1;
    for(int i=2;i<=n;i++)    //计算任意一个n-1阶行列式的值就是答案.
    {
        for(int j=i+1;j<=n;j++)
        {     
            while(a[j][i])
            {
                ll t=a[i][i]/a[j][i];    //如果要模,则这里就一下
                for(int k=i; k<=n; k++)
                    a[i][k]=(a[i][k]-a[j][k]*t);    //加一下
                for(int k=i; k<=n; k++)
                    swap(a[i][k],a[j][k]);      
            }
        }
        if(a[i][i]==0) return 0;     //图无法联通.
        ret*=abs(a[i][i]);     //加一下
    }
    return ret;
}