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

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

Time Limit: 20 Sec   Memory Limit: 256 MB
Submit: 212   Solved: 135
[ Submit][ Status][ Discuss]

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

HINT

Source

By 佚名上传

[ Submit][ Status][ Discuss]
看上去问题很复杂。。先简单化 如果给一张无向图G，如何快速求其生成树个数？？？ 引入矩阵树定理，这里有vfleaking巨神的证明传送门 有兴趣的朋友可以看看啦~但是苟蒻博主只是学了结论= = 嗯。。。大概是这样的 定义无限图G的基尔霍夫矩阵为C[G]，度数矩阵为A[G]，邻接矩阵为B[G] 在A[G]中， 在B[G]中， 然后C[G] = A[G] - B[G] 此时图G的生成树为C[G]的行列式的任意一个n - 1阶主子式的值 而计算行列式的话，，暴力法是O(n!)的 由行列式的计算法则，任意一行的任意倍加到另一行后，行列式的值不变 如果交换任意不同的两行，行列式的值取相反数 这样可以将行列式进行高斯消元，变换成上三角行列式 行列式的值就是对角线元素乘积了 消元的话利用扩展欧几里得算法类似的步骤，即(a,b) -> (b,a%b) 复杂度多个log
如果题中没有建筑公司的限制，到这里就完了。。。还剩最后一点点小步骤 因为题中的限制，n - 1条边由n - 1个不同的建筑公司分别承包 如果直接统计方案太难。。但是，如果有Matrix-tree定理，已经可以解决总方案数的问题了 那么，ans = 总方案数 - 不合法方案数，又一个公司要且必须要承包一条边 那么ans = 总方案数 - 一个公司不修路的方案 + 两个公司不修路的方案...... 这里n很小，完全可以O(2^n)枚举所有公司的情况，就用容斥原理统计方案，，问题解决

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<bitset>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

typedef long long LL;
const LL mo = 1000000007;

struct E{
	int x,y; E(){}
	E(int x,int y): x(x),y(y){}
};

int Ans,n,ti,A[20][20];

vector <E> v[20];

void Swap(int x,int y,int k) {for (int i = k; i < n; i++) swap(A[x][i],A[y][i]);}
int Inc(LL x,LL y) {return (x + y) % mo;}
int Dec(LL x,LL y) {return (x - y + mo) % mo;}
int Mul(LL x,LL y) {return x*y%mo;}

void GCD(int x,int y)
{
	if (!A[y][x]) return;
	int t = A[x][x] / A[y][x];
	for (int i = x; i < n; i++) 
		A[x][i] = Dec(A[x][i],Mul(t,A[y][i]));
	Swap(x,y,x); ++ti; 
	GCD(x,y);
}

void Gauss(int x)
{
	if (x == n - 1) return;
	if (!A[x][x])
		for (int j = x + 1; j < n; j++)
			if (A[j][x]) {Swap(x,j,x); ++ti; break;}
	if (!A[x][x]) return;
	for (int j = x + 1; j < n; j++)
		if (A[j][x]) GCD(x,j);
	Gauss(x + 1);
}

int Work()
{
	int ret = 1; ti = 0; Gauss(1);
	for (int i = 1; i < n; i++) ret = Mul(ret,A[i][i]);
	return (ti&1)?mo - ret:ret;
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	cin >> n;
	for (int i = 0; i < n - 1; i++)
	{
		int k; scanf("%d",&k);
		while (k--)
		{
			int x,y; scanf("%d%d",&x,&y);
			v[i].push_back(E(--x,--y));
		}
	}
	
	for (int i = 0; i < (1<<(n-1)); i++)
	{
		int tot = 0; memset(A,0,sizeof(A));
		for (int j = 0; j < n - 1; j++)
			if (!(i&(1<<j)))
			{
				for (int k = 0; k < v[j].size(); k++)
				{
					E e = v[j][k];
					++A[e.x][e.x];
					++A[e.y][e.y];
					A[e.x][e.y] = Dec(A[e.x][e.y],1);
					A[e.y][e.x] = Dec(A[e.y][e.x],1);
				}
			}
			else ++tot;
		if (tot&1) Ans = Dec(Ans,Work());
		else Ans = Inc(Ans,Work());
	}
	cout << Ans;
	return 0;
}