[bzoj4596][SHOI2016]黑暗前的幻想乡

题目描述

四年一度的幻想乡大选开始了，最近幻想乡最大的问题是很多来历不明的妖
怪涌入了幻想乡，扰乱了幻想乡昔日的秩序。但是幻想乡的建制派妖怪（人类）
博丽灵梦和八云紫等人整日高谈所有妖怪平等，幻想乡多元化等等，对于幻想乡
目前面临的种种大问题却给不出合适的解决方案。
风间幽香是幻想乡里少有的意识到了问题的严重性的大妖怪。她这次勇敢的
站了出来参加幻想乡大选。提出包括在幻想乡边境建墙（并让人类出钱），大力
开展基础设施建设挽回失业率等一系列方案，成为了大选年出人意料的黑马并顺
利的当上了幻想乡的大统领。
幽香上台以后，第一项措施就是要修建幻想乡的公路。幻想乡有 N 个城市，
之间原来没有任何路。幽香向选民承诺要减税，所以她打算只修 N- 1 条路将
这些城市连接起来。但是幻想乡有正好 N- 1 个建筑公司，每个建筑公司都想
在修路的过程中获得一些好处。
虽然这些建筑公司在选举前没有给幽香钱，幽香还是打算和他们搞好关系，
因为她还指望他们帮她建墙。所以她打算让每个建筑公司都负责一条路来修。
每个建筑公司都告诉了幽香自己有能力负责修建的路是哪些城市之间的。所
以幽香打算选择 N-1 条能够连接幻想乡所有城市的边，然后每条边都交给一
个能够负责该边的建筑公司修建，并且每个建筑公司都恰好修一条边。
幽香现在想要知道一共有多少种可能的方案呢？两个方案不同当且仅当它
们要么修的边的集合不同，要么边的分配方式不同。

矩阵树定理

我们可以发现有n-1个条件，第i个条件可以看做第i个公司要修边。
同时满足这n-1个条件并造出生成树就是合法的。
那么可以指数级枚举去枚举哪些条件不满足，剩余的做矩阵树定理求生成树个数即可。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int maxn=17+10,mo=1000000007;
int cr[maxn][maxn],c[maxn][maxn],d[maxn][maxn],a[maxn][maxn];
int edge[maxn][maxn*maxn][2],cnt[maxn];
bool bz[maxn];
int i,j,k,l,t,n,m,ans;
int det(){
    int i,j,k,l=1,t;
    fo(i,1,n-1){
        fo(j,i+1,n-1)
            while (cr[j][i]){
                t=cr[i][i]/cr[j][i];
                fo(k,i,n) cr[i][k]=(cr[i][k]-(ll)t*cr[j][k]%mo)%mo;
                fo(k,i,n) swap(cr[i][k],cr[j][k]);
                l=-l;
            }
        l=(ll)l*cr[i][i]%mo;
    }
    return l;
}
int work(){
    int i,j,k,l;
    fo(i,1,n)
        fo(j,1,n)
            d[i][j]=a[i][j]=0;
    fo(i,1,n-1)
        if (!bz[i])
            fo(j,1,cnt[i]){
                k=edge[i][j][0];l=edge[i][j][1];
                d[k][k]++;d[l][l]++;
                a[k][l]++;a[l][k]++;
            }
    fo(i,1,n)
        fo(j,1,n)
            c[i][j]=d[i][j]-a[i][j];
    fo(i,1,n-1)
        fo(j,1,n-1)
            cr[i][j]=c[i+1][j+1];
    return det();
}
void dfs(int x,int y){
    if (x==n){
        ans=(ans+y*work())%mo;
        return;
    }
    bz[x]=0;
    dfs(x+1,y);
    bz[x]=1;
    dfs(x+1,-y);
}
int main(){
    scanf("%d",&n);
    fo(i,1,n-1){
        scanf("%d",&cnt[i]);
        fo(j,1,cnt[i]) scanf("%d%d",&edge[i][j][0],&edge[i][j][1]);
    }
    dfs(1,1);
    (ans+=mo)%=mo;
    printf("%d\n",ans);
}