Description

Input

Output

Sample Input

1
3 3
1 2
2 3
1 3

Sample Output

6

题目大意：求所有生成树的prufer序列和（prufer中有重复序列的只算一次！！！）

题解：

这样的话，对于每一颗生成树，我们可以把所有的点全都加进去，然后再减去叶子结点的和。

我们不可能找到所有的生成树然后一个一个的计算，因此我们用矩阵树定理来做。

我们先计算图所有的点的和，并且乘以生成树的数量，把他们放在sum里。然后再把所有的叶子结点减去，就好了

如果一个叶子节点出现在一颗子树里，那么把这个点去掉，仍然可一得到图的该生成树，而如果这个点是内部节点就不行了。

注意：如果这个叶子节点的度不为1，那么要用这个叶子节点的度数乘以生成树的数量，才是这个叶子节点对应生成树的个数。

sum -= 去掉该节点生成树的数量*该节点的度*该节点的值。

最后得到的sum就是答案

代码：

#include<iostream>。
#include<cmath>
#include <cstring>
using namespace std;
#define zero(x)((x>0? x:-x)<1e-15)
#define int long long
int const MAXN = 105;
const int mod = 1e9 + 7;
int a[MAXN][MAXN];
int b[MAXN][MAXN];
int g[103][103];
int d[105];
int n, m;
int det(int a[MAXN][MAXN], int n){ int s=0;for(int i=0; i<n; i++){ int r=i; for(; r<n; r++)    // error-prone if(a[r][i]) break; if(r==n+1) return 0; if(r!=i){ s^=1; for(int j=i; j<n; j++) swap(a[i][j], a[r][j]); } for(int j=i+1; j<n; j++){ int x=i, y=j; for(; a[y][i]; ){ // print(a, n); int t=a[y][i]/a[x][i]; if(t){ for(int k=i; k<n; k++){ a[y][k] -= t*a[x][k]%mod; a[y][k] %= mod; } if(a[y][i]==0) break; } swap(x, y); } if(x!=i){ for(int k=i; k<n; ++k) swap(a[x][k], a[y][k]); s ^= 1; } } } int res=1; for(int i=0; i<n; i++) res*=a[i][i], res%=mod; if(s){ res=-res; }if(res < 0) res+=mod; return res; 
} void prep(int n,int x)
{for(int i = 0;i < n;i++){for(int j = 0;j < n;j++){a[i][j] = (i == j)?d[i]-g[i][x]:-g[i][j];}}if(x + 1){for(int i = 0;i < n;i++) a[x][i] = a[i][x] = 0;a[x][x] = 1;}
}main() 
{	int cas;scanf("%lld", &cas);while (cas--) {memset(g,0,sizeof(g));memset(d,0,sizeof(d));memset(a,0,sizeof(a));memset(b,0,sizeof(b));scanf("%lld%lld", &n,&m);for(int i = 0;i < m;i++){int u,v;scanf("%lld%lld",&u,&v);u--,v--;d[u]++;d[v]++;g[u][v] = g[v][u] = 1;}prep(n,-1);int sum = det(a,n-1)*((1+n)*n/2) % mod;for(int i = 0;i < n-1;i++){prep(n,i);sum = (sum - (i+1)*d[i] % mod * det(a,n-1) % mod + mod)%mod;//cout<<':'<<sum<<endl;}prep(n,n-1);sum = (sum - n*d[n-1] % mod * det(a,n-2) % mod + mod)%mod;cout<<sum<<endl;}return 0;
}