BZOJ 4596: [Shoi2016]黑暗前的幻想乡

Description

四年一度的幻想乡大选开始了，最近幻想乡最大的问题是很多来历不明的妖

怪涌入了幻想乡，扰乱了幻想乡昔日的秩序。但是幻想乡的建制派妖怪（人类）

博丽灵梦和八云紫等人整日高谈所有妖怪平等，幻想乡多元化等等，对于幻想乡

目前面临的种种大问题却给不出合适的解决方案。

风间幽香是幻想乡里少有的意识到了问题的严重性的大妖怪。她这次勇敢的

站了出来参加幻想乡大选。提出包括在幻想乡边境建墙（并让人类出钱），大力

开展基础设施建设挽回失业率等一系列方案，成为了大选年出人意料的黑马并顺

利的当上了幻想乡的大统领。

幽香上台以后，第一项措施就是要修建幻想乡的公路。幻想乡有 N 个城市，

之间原来没有任何路。幽香向选民承诺要减税，所以她打算只修 N- 1 条路将

这些城市连接起来。但是幻想乡有正好 N- 1 个建筑公司，每个建筑公司都想

在修路的过程中获得一些好处。

虽然这些建筑公司在选举前没有给幽香钱，幽香还是打算和他们搞好关系，

因为她还指望他们帮她建墙。所以她打算让每个建筑公司都负责一条路来修。

每个建筑公司都告诉了幽香自己有能力负责修建的路是哪些城市之间的。所

以幽香打算选择 N-1 条能够连接幻想乡所有城市的边，然后每条边都交给一

个能够负责该边的建筑公司修建，并且每个建筑公司都恰好修一条边。

幽香现在想要知道一共有多少种可能的方案呢？两个方案不同当且仅当它

们要么修的边的集合不同，要么边的分配方式不同。

Input

第一行包含一个正整数 N(N<=17)，表示城市个数。

接下来 N-1 行，其中第 i行表示第 i个建筑公司可以修建的路的列表：

以一个非负数mi 开头，表示其可以修建 mi 条路，接下来有mi 对数，

每对数表示一条边的两个端点。其中不会出现重复的边，也不会出现自环。

Output

仅一行一个整数，表示所有可能的方案数对 10^9 + 7 取模的结果。

Sample Input

4
2 3 2 4 2
5 2 1 3 1 3 2 4 1 4 3
4 2 1 3 2 4 1 4 2

Sample Output

HINT

Source

By 佚名上传

想法：容斥+矩阵树定理。复杂度$O(2^{n-1}*n^3)$。如果分不清奇加偶减，可以考虑全选的一定是加，那么奇偶性和$n-1$一样就好了。因为模数为$10^9+7$是个质数，可以用辗转相除法。

#include< cstdio > #include< vector > #define gec getchar #define FILE(F) freopen(F".in","r",stdin),freopen(F".out","w",stdout) #define DEBUG fprintf(stderr,"Passing [%s] in Line (%d)\n",__FUNCTION__,__LINE__); typedef long long ll; template inline void read(T&x) { x=0;bool f=0;char c=gec(); for(;c<‘0‘||c>‘9‘;c=gec())f=(c==‘-‘); for(;c>=‘0‘&&c<=‘9‘;c=gec())x=x*10+c-‘0‘; x=f?-x:x; } const int N(20),MP(1e9+7); ll A[N][N],G[N][N];int Ans,n; struct Data{int a,b;}; std::vectorEdge[20];int siz[20]; void swap(ll &a,ll &b){ll t(a);a=b;b=t;} int Det()//删第n行n列 { for(int i=1;i<=n;i++) for(int j=1;j<=n;j++)A[i][j]=G[i][j]; ll ans=1,tmp; for(int i=1;i<n;i++) { for(int j=i+1;j<n;j++) while(A[j][i]) { tmp=A[i][i]/A[j][i]; for(int k=i;k<n;k++) (A[i][k]-=A[j][k]*tmp%MP)%=MP; for(int k=i;k<n;k++) swap(A[i][k],A[j][k]); ans=-ans; } (ans*=A[i][i])%=MP; } ans+=ans<0?MP:0; return ans; } void update(int x,int y) { for(int v=0;v<siz[x];v++) {G[Edge[x][v].a][Edge[x][v].a]+=y;G[Edge[x][v].b][Edge[x][v].b]+=y; G[Edge[x][v].b][Edge[x][v].a]-=y;G[Edge[x][v].a][Edge[x][v].b]-=y;} } void DFS(int x,int sum) { for(int j=x+1;j<n;j++) update(j,1),DFS(j,sum+1),update(j,-1); if(sum%2==(n-1)%2) Ans+=Det();else Ans-=Det(); Ans%=MP; } int main() { #ifndef ONLINE_JUDGE FILE("C"); #endif read(n); for(int i=1;i<n;i++) { read(siz[i]); for(int j=1,a,b;j<=siz[i];j++) read(a),read(b),Edge[i].push_back((Data){a,b}); } DFS(0,0); Ans+=Ans<0?MP:0; printf("%d\n",Ans); return 0; }