[BZOJ4596][Shoi2016]黑暗前的幻想乡（容斥+矩阵树定理）

Description

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

Input

第一行包含一个正整数 $N(N\le17)$ ，表示城市个数。
接下来 $N-1$ 行，其中第 $i$ 行表示第 $i$ 个建筑公司可以修建的路的列表：
以一个非负数 $m_i$ 开头，表示其可以修建 $m_i$ 条路，接下来有 $m_i$ 对数，
每对数表示一条边的两个端点。其中不会出现重复的边，也不会出现自环。

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

Solution

显然是一个生成树计数问题。
但生成树中选出的边是有限制的：给定 $n-1$ 组边，规定每条边都必须选每一组边中的一条。
于是就考虑容斥：
$2^{n-1}$ 暴力枚举 $\{1,2,...,n-1\}$ 的一个子集 $S$ ，并强制 $S$ 内的元素对应的边集不能被作为生成树上的边，每次求一遍生成树个数。
这样，答案就为：

$$|S|=0的生成树个数-|S|=1的生成树个数+|S|=2的生成树个数-...$$
求生成树个数，可以用 Matrix-Tree 矩阵树定理求得。
复杂度 $O(2^nn^3\log)$ 。

Code

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define For(i, a, b) for (i = a; i <= b; i++)
using namespace std;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
const int N = 20, M = N * N, PYZ = 1e9 + 7;
int n, p[N][M], q[N][M], ans, mat[N][N];
bool sel[N];
int det(int n) {
    int i, j, k; bool op = 0;
    For (i, 1, n) For (j, 1, n) if (mat[i][j] < 0) mat[i][j] += PYZ;
    For (i, 1, n) For (j, i + 1, n) {
        int a = mat[i][i], b = mat[j][i]; while (b) {
            int tmp = a / b; a %= b; swap(a, b);
            For (k, i, n) mat[i][k] = (mat[i][k] - 1ll * tmp * mat[j][k]
                % PYZ + PYZ) % PYZ;
            For (k, i, n) swap(mat[i][k], mat[j][k]); op ^= 1;
        }
    }
    int ans = 1; For (i, 1, n) ans = 1ll * ans * mat[i][i] % PYZ;
    return op ? (PYZ - ans) % PYZ : ans;
}
void calc(int cnt) {
    int i, j; For (i, 1, n) For (j, 1, n) mat[i][j] = 0;
    For (i, 1, n - 1) if (sel[i]) For (j, 1, p[i][0])
        mat[p[i][j]][q[i][j]]--, mat[q[i][j]][p[i][j]]--;
    For (i, 1, n) For (j, 1, n) if (i != j) mat[i][i] -= mat[i][j];
    int res = det(n - 1); if (cnt & 1) ans = (ans - res + PYZ) % PYZ;
    else ans = (ans + res) % PYZ;
}
void dfs(int dep, int cnt) {
    if (dep == n) return calc(cnt);
    sel[dep] = 1; dfs(dep + 1, cnt);
    sel[dep] = 0; dfs(dep + 1, cnt + 1);
}
int main() {
    int i, j; n = read(); For (i, 1, n - 1) {
        p[i][0] = read(); For (j, 1, p[i][0]) p[i][j] = read(), q[i][j] = read();
    }
    dfs(1, 0); cout << ans << endl; return 0;
}