[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

17

这题直接Matrix-Tree肯定是不行的，因为我们要把n-1条边分给n-1个不同公司

直接计算出来的方案包括有些公司未参与的修建方案

令f[i]为禁止i个公司参与的方案数，利用容斥原理

$ans=f[0]-f[1]+f[2]-f[3]....(-1)^{n}*f[n]$

可以枚举二进制状态，参加为1，禁止参加为0

然后建出相应的基尔霍夫矩阵，根据Matrix-Tree定理算出行列式的值

根据0的数量决定加减

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
struct ZYYS
{
  int u,v;
}edge[18][501];
int n,a[18][18],Mod=1e9+7,len[18],ans;
int guass()
{int i,j,k,s;
  n--;
  for (i=1;i<=n;i++)
    {
      for (j=1;j<=n;j++)
    {
      a[i][j]=(a[i][j]+Mod);
      if (a[i][j]>=Mod) a[i][j]-=Mod;
    }
    }
  s=1;
  for (i=1;i<=n;i++)
    {
      for (j=i+1;j<=n;j++)
    while (a[j][i])
      {
        int t=a[i][i]/a[j][i];
        for (k=i;k<=n;k++)
          {
        a[i][k]=(a[i][k]-1ll*t*a[j][k]%Mod+Mod);
        if (a[i][k]>=Mod) a[i][k]-=Mod;
        swap(a[i][k],a[j][k]);
          }
        s*=-1;
      }
      s=1ll*s*a[i][i]%Mod;
      if (!s)
    {
      n++;
      return 0;
    }
    }
  n++;
  return (s+Mod)%Mod;
}
int main()
{int i,j,k,as,cnt;
  cin>>n;
  for (i=1;i<=n-1;i++)
    {
      scanf("%d",&len[i]);
      for (j=1;j<=len[i];j++)
    {
      scanf("%d%d",&edge[i][j].u,&edge[i][j].v);
    }
    }
  for (i=0;i<=(1<<n-1)-1;i++)
    {
      cnt=0;
      memset(a,0,sizeof(a));
      for (j=1;j<=n-1;j++)
    if ((i&(1<<j-1)))
      {
        for (k=1;k<=len[j];k++)
          {
        int x=edge[j][k].u,y=edge[j][k].v;
        a[x][y]--;a[y][x]--;
        a[x][x]++;a[y][y]++;
          }
      }
    else cnt++;
      as=guass();
      if (cnt%2)
      ans=(ans-as+Mod);
      else ans=(ans+as);
      if (ans>=Mod) ans-=Mod; 
    }
  cout<<(ans+Mod)%Mod;
}