题链：

http://www.lydsy.com/JudgeOnline/problem.php?id=4596

题解：

容斥，矩阵树定理，矩阵行列式

先说说容斥：(一共有 N-1个公司)
令 f[i] 表示选出 (N-1)-i 个公司来修路(即有i个公司一定不修)，且不管每个公司只能修一条路这一限制的方案数。
那么 答案 ANS=0个公司不修的方案数 - 1个公司不修的方案数 +2个公司不修的翻案数 ...
即 ANS= f[0] - f[1] +f[2] - ... + (-1)ⁱ*f[i]
f[i]的求法呢，就是先O(2^N)枚举公司集合情况，
然后用矩阵树定理 O(N³) 求出当前情况下的生成树方案数。
另外再本题中，在构造上三角矩阵用以求行列式时，
既要避免小数，还要不影响原矩阵的行列式的值，
所以采用辗转相除的高斯消元法去构造上三角矩阵。复杂度多一个(logN)
由矩阵行列式的性质可知，这样辗转相除的高斯消元法不会影响行列式的绝对值，
只会影响符号位的正负，所以统计一下正负号的变化就好了。
(还不会矩阵树定理，看看这里，入门一波。)
所以总的时间复杂度为 O(2^N*N³*logN)。
(都是这个复杂度，不晓得为什么我的代码跑到这么慢，都垫底了......)

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#define add(x,y) (((1ll*(x)+(y))%mod+mod)%mod)
#define mul(x,y) (((1ll*(x)*(y))%mod+mod)%mod)
#define filein(x) freopen(#x".in","r",stdin)
#define fileout(x) freopen(#x".out","w",stdout)
using namespace std;
const int mod=1000000007;
struct Matrix{int Val[20][20],*X[20],R,C;void Init(int r,int c){//r==cR=r; C=c;memset(Val,0,sizeof(Val));for(int i=1;i<=R;i++) X[i]=Val[i];}void Modify(int r,int c,int v){X[r][c]=add(X[r][c],v);}void operator =(const Matrix &rtm){Init(rtm.R,rtm.C);for(int i=1;i<=R;i++)for(int j=1;j<=C;j++)Val[i][j]=rtm.X[i][j];}Matrix operator -(const Matrix & rtm) const{Matrix now; now=*this;for(int i=1;i<=R;i++)for(int j=1;j<=C;j++)now.X[i][j]=add(now.X[i][j],-rtm.X[i][j]);return now;}void Gauss_Euclidean(int p,int &ti){//形成上三角矩阵 if(p==R-1) return;if(!X[p][p]) for(int i=p+1;i<R;i++) if(X[i][p]){swap(X[i],X[p]); ti++; break;}if(!X[p][p]) return;for(int i=p+1;i<R;i++){while(X[i][p]){int t=X[p][p]/X[i][p];for(int j=p;j<R;j++)X[p][j]=add(X[p][j],-mul(X[i][j],t));swap(X[p],X[i]); ti++;}}Gauss_Euclidean(p+1,ti);}int Determinant(){int ti=0,ans=1;Gauss_Euclidean(1,ti);for(int i=1;i<R;i++) ans=mul(ans,X[i][i]);if(ti&1) ans=mul(ans,-1);return ans;}void print(){for(int i=1;(i!=1?printf("\n"):0),i<=R;i++)for(int j=1;j<=R;j++)printf("%d ",X[i][j]);}
}A,B,K;
int cpy[20][500];
int ANS,N,tmp;
void dfs(int p,int num){if(p>=N) return;//选for(int i=1,a,b;i<=2*cpy[p][0];i+=2){a=cpy[p][i]; b=cpy[p][i+1];A.Modify(a,a,1); A.Modify(b,b,1);B.Modify(a,b,1); B.Modify(b,a,1);}K=A-B; tmp=K.Determinant();if(((N-1)-(num+1))&1) tmp=mul(tmp,-1);ANS=add(ANS,tmp);dfs(p+1,num+1);//不选for(int i=1,a,b;i<=2*cpy[p][0];i+=2){a=cpy[p][i]; b=cpy[p][i+1];A.Modify(a,a,-1); A.Modify(b,b,-1);B.Modify(a,b,-1); B.Modify(b,a,-1);}dfs(p+1,num);
}
int main()
{scanf("%d",&N);A.Init(N,N); B.Init(N,N);for(int i=1;i<=N-1;i++){scanf("%d",&cpy[i][0]);for(int j=1;j<=2*cpy[i][0];j++)scanf("%d",&cpy[i][j]);}dfs(1,0);printf("%d",ANS);return 0;
}