【JZOJ 省选模拟】赢家(winner)

题目

Description
PinkRabbit 是一位人赢。
福州市可以抽象成一个n个点m条边的，不包含重边与自环的无向图，PinkRabbit 住在1号
点，而他的妹子住在2号点。
某一天，PinkKitten 施放了一个大魔法，让这个无向图上所有的边都变成了单向边。现在
PinkRabbit 关心的是他是否能够和他的妹子见面。
具体地，PinkRabbit 能和他的妹子见面，当且仅当存在一个点 u，满足新图上从1号点出发能够
到达u，从2号点出发也能到达 。
现在你需要计算出，在把所有 m条边进行定向的所有2^m 种方案中，有多少种方案能让
PinkRabbit 和他的妹子见面。你只需输出其对10^9+7 取模后的结果。

Input
第一行三个正整数，分别为 n、m 、id 。n 和 m为图的点数和边数，id 为子任务编号。
接下来有m 行，每行两个正整数x 和y ，描述一条边。

Output
输出一个整数表示答案。

Sample Input
样例输入 1：
3 2 1
1 3
2 3
样例输入 2：
4 5 1
1 3
2 3
1 4
2 4
3 4
样例输入 3：
5 6 1
1 3
2 3
3 4
3 5
1 4
2 5

Sample Output
样例输出 1：
3
样例输出 2：
30
样例输出 3：
55

Data Constraint

Hint
【样例解释 1】
这个样例的id=1 ，满足子任务1 的m<=20 限制。
合法的方案有3种：
（1）1->3 ，2->3
（2）1->3 ，3->2
（3）3->1 ，2->3

思路

考虑用总方案-不合法方案。

对于不合法方案，枚举1能走到的点集是S,2能走到的点集是T，S∩T=∅。

设Z=总集−S−T

那么Z和S、T之间的边方向确定，S和T之间不能有边。

现在就是求S的方案数（T同理）。

设f[S]表示S的方案数，同样用总-不合法。
不合法就是是S的一个子集S′,f[S]−=f[S′]∗(S′和S−S′中间的边定向)

预处理cnt[S]表示S内的边数就都可以快速计算了。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int mod=1e9+7,N=16;
int n,m,id,x,y,b[N][N],cnt[1<<15];
ll a2[N*N],f[1<<15],g[1<<15];
int main()
{
//	freopen("winner.in","r",stdin); freopen("winner.out","w",stdout);
	scanf("%d%d%d",&n,&m,&id);
	for(int i=1; i<=m; i++) scanf("%d%d",&x,&y),b[x][y]=b[y][x]=1;
	a2[0]=1;
	for(int i=1; i<=n*n; i++)a2[i]=a2[i-1]*2%mod;
	for(int s=0; s<1<<n; s++)
	{
		for(int i=0; i<n; i++) if(s>>i&1)
		{
			cnt[s]=cnt[s^(1<<i)];
			for(int j=0; j<n; j++) if(s>>j&1) cnt[s]+=b[i+1][j+1];
			break;
		}
	}
	for(int s=0; s<1<<n; s++)
	{
		if(s&1)
		{
			f[s]=a2[cnt[s]];
			for(int t=(s-1)&s; t>0; t=(t-1)&s)
				f[s]=(f[s]-f[t]*a2[cnt[s^t]])%mod;
		}
		if(s&2)
		{
			g[s]=a2[cnt[s]];
			for(int t=(s-1)&s; t>0; t=(t-1)&s)
				g[s]=(g[s]-g[t]*a2[cnt[s^t]])%mod;
		}
	}
	ll ans=a2[m];
	for(int s=0; s<1<<n; s++) if(f[s])
	{
		for(int t=s+1; t<(1<<n); t=(t+1)|s) if(g[s^t])
		{
			int z=(1<<n)-1-t;
			if(cnt[s^z]+cnt[s^t^z]-cnt[z]==m)
			{
				ans=(ans-f[s]*g[s^t]%mod*a2[cnt[z]])%mod;
			}
		}
	}
	ans=(ans%mod+mod)%mod;
	printf("%lld\n",ans);
}