4455: [Zjoi2016]小星星

4455: [Zjoi2016]小星星

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 246   Solved: 161
[ Submit][ Status][ Discuss]

Description

小Y是一个心灵手巧的女孩子，她喜欢手工制作一些小饰品。她有n颗小星星，用m条彩色的细线串了起来，每条细

线连着两颗小星星。有一天她发现，她的饰品被破坏了，很多细线都被拆掉了。这个饰品只剩下了n?1条细线，但

通过这些细线，这颗小星星还是被串在一起，也就是这些小星星通过这些细线形成了树。小Y找到了这个饰品的设

计图纸，她想知道现在饰品中的小星星对应着原来图纸上的哪些小星星。如果现在饰品中两颗小星星有细线相连，

那么要求对应的小星星原来的图纸上也有细线相连。小Y想知道有多少种可能的对应方式。只有你告诉了她正确的

答案，她才会把小饰品做为礼物送给你呢。

Input

第一行包含个2正整数n,m，表示原来的饰品中小星星的个数和细线的条数。

接下来m行，每行包含2个正整数u,v，表示原来的饰品中小星星u和v通过细线连了起来。

这里的小星星从1开始标号。保证u≠v，且每对小星星之间最多只有一条细线相连。

接下来n-1行，每行包含个2正整数u,v，表示现在的饰品中小星星u和v通过细线连了起来。

保证这些小星星通过细线可以串在一起。

n<=17,m<=n*(n-1)/2

Output

输出共1行，包含一个整数表示可能的对应方式的数量。

如果不存在可行的对应方式则输出0。

Sample Input

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

Sample Output

6

HINT

 题解:JudgeOnline/upload/201603/4455.txt

Source

[ Submit][ Status][ Discuss]



显然是要dp统计方案，但怎么统计是个问题。。。可以利用容斥来做

每一种合法方案，一定是树上的点与原图的点建立一一对应的关系

那么，不合法的方案，就是树上多个点对应到原图的一个点了

令f[i][j]:点i对应原图中点j，i为根的子树与原图中的点的对应方案数

转移的话，逐个合并子树信息即可O(n^2)，每个点只在其父亲处被合并，共O(n^3)

每次枚举一个点集，规定不许用这个点集的点转移，这样就能得出一定不合法的某种方案

如此进行容斥即可，，O(2^n*n^3)，常数要写小一点，这个数值有点极限- -

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<cmath>
using namespace std;

const int N = (1<<17);
typedef long long LL;

LL Ans,f[17][17],g[17];
int n,m,tp,h[17];
bool E[17][17],bo[17];

vector <int> v[17];

void Dfs(int x,int fa)
{
	for (int i = 0; i < tp; i++) f[x][h[i]] = 1;
	for (int i = 0; i < v[x].size(); i++)
	{
		int y = v[x][i];
		if (y == fa) continue;
		Dfs(y,x); for (int j = 0; j < tp; j++) g[h[j]] = 0;
		for (int A = 0; A < tp; A++)
			for (int B = 0; B < tp; B++)
				if (E[h[A]][h[B]]) g[h[A]] += f[x][h[A]]*f[y][h[B]];
		for (int j = 0; j < tp; j++) f[x][h[j]] = g[h[j]];
	}
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	cin >> n >> m;
	while (m--)
	{
		int x,y; 
		scanf("%d%d",&x,&y);
		--x; --y;
		E[x][y] = E[y][x] = 1;
	}
	for (int i = 1; i < n; i++)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		--x; --y;
		v[x].push_back(y);
		v[y].push_back(x);
	}
	
	for (int i = 0; i < (1<<n); i++)
	{
		int tot = 0; LL ret = 0; tp = 0;
		for (int j = 0; j < n; j++)
			if (i&(1<<j)) bo[j] = 1,++tot;
			else bo[j] = 0,h[tp++] = j;
		memset(f,0,sizeof(f)); Dfs(0,-1);
		for (int j = 0; j < tp; j++) ret += f[0][h[j]];
		if (tot&1) Ans -= ret; else Ans += ret;
	}
	cout << Ans;
	return 0;
}