【spoj104】highways

来源:互联网 发布:mac 查询端口号 编辑:程序博客网 时间:2024/06/05 17:09

生成树计数


Description
一个有n座城市的组成国家,城市1至n编号,其中一些城市之间可以修建高速公路。现在,需要有选择的修建一些高速公路,从而组成一个交通网络。你的任务是计算有多少种方案,使得任意两座城市之间恰好只有一条路径?
Input
第一行两个整数n和m,分别表示点数和边数;
接下来m行,每行包含两个整数x,y,表示城市x和城市y之间的有公路相连
Output
输出一行,表示方案的个数。
Sample Input
5 6
1 2
1 3
2 3
2 4
4 5
3 5
Sample Output
11


这有专门的算法
什么什么构造一个拉普拉斯矩阵呀
再算一下行列式啊
以下是行列式计算的模板


#include<cstdio>#include<iostream>#include<cstring>#include<cstdlib>#include<cmath>#include<algorithm>using namespace std;int n,m,g[101][101],d[101][101];int a[101][101];#define zero(x) ((x>0? x:-x)<1e-15)double det(){    int i,j,k,sign=0;    double b[101][101],ret=1,t;    for(i=0;i<n;i++)      for(j=0;j<n;j++)         b[i][j]=a[i][j];    for(i=0;i<n;i++)    {       if(zero(b[i][i]))       {        for(j=i+1;j<n;j++)         if(!zero(b[j][i])) break;       if(j==n) return 0;       for(k=i;k<n;k++) swap(b[i][k],b[j][k]);       sign++;       }       ret*=b[i][i];       for(k=i+1;k<n;k++) b[i][k]/=b[i][i];       for(j=i+1;j<n;j++)         for(k=i+1;k<n;k++)           b[j][k]-=b[j][i]*b[i][k];     }     return fabs(ret);}int main(){    int i,j,x,y;   cin>>n>>m;   for(i=1;i<=m;i++)    {     cin>>x>>y;     g[x-1][y-1]=g[y-1][x-1]=1;   }   for(i=0;i<n;i++)   {     int dd=0;     for(j=0;j<n;j++)        if(g[i][j]) dd++;     a[i][i]=dd;   }   for(i=0;i<n;i++) for(j=0;j<n;j++) if(g[i][j]) a[i][j]=-1;   n--;   printf("%0.0lf",det());   return 0;}
0 0