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

17

分析

首先我们要知道一个叫做矩阵树定理的东西，其实矩阵树定理说的是什么呢？矩阵树定理就是说：对于一个无向图 G ，它的生成树个数等于其基尔霍夫Kirchhoff矩阵任何一个N-1阶主子式的行列式的绝对值。所谓的N-1阶主子式就是对于一个任意的一个 r ，将矩阵的第 r 行和第 r 列同时删去得到的新矩阵。
基尔霍夫Kirchhoff矩阵的一种求法：
基尔霍夫Kirchhoff矩阵 K =度数矩阵 D - 邻接矩阵 A
那等我们知道这个东西了之后令 f[i] 表示选出 (N-1)-i 个公司来修路(即有i个公司一定不修)，且不管每个公司只能修一条路这一限制的方案数。
那么 答案 ANS=0个公司不修的方案数 - 1个公司不修的方案数 +2个公司不修的翻案数 …
即 ANS= f[0] - f[1] +f[2] - … + (-1)i*f[i]
f[i]的求法呢，就是先O(2N)枚举公司集合情况，
然后用矩阵树定理 O(N3) 求出当前情况下的生成树方案数。就可以了

CODE

#include <bits/stdc++.h>

#define add(x,y) (((1ll*(x)+(y))%mod+mod)%mod)
#define mul(x,y) (((1ll*(x)*(y))%mod+mod)%mod)
using namespace std;

const int mod = 1000000007;

struct Matrix
{
    int Val[20][20],*X[20],R,C;
    void Init(int r,int c)
    {
        R = 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;
    }
}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;
}