P5933 [清华集训2012]串珠子(连通性 状压 计数)

P5933 [清华集训2012]串珠子

\(\color{yellow}{\bigstar\texttt{Trick}}\):遇到连通性题可以暂时忽略是否联通

\(g_s\) 表示集合为 \(s\) 的点的子图有多少个,可以不连通;\(f_s\) 表示答案。

那么 \(g_s\) 可以直接计算 \(g_s=\prod_{(i<j)\in s}(a_{i,j}+1)\)

\(f_s\) 可以用所有图的数量减去一定不连通的数量,如何计算一定不连通的方案数:我们选出一个一定包含 \(p\) 节点的联通集合,将这些点和其他的分割开来,其他点之间随便连接。对于每个集合都求出这样的答案,和就是一定不连通的方案数。

\(\color{yellow}{\bigstar\texttt{Trick}}\):上面一定包含 \(p\) 的信息非常重要,这样保证了不会计算重复方案。

#define Maxn 18
#define Maxsta 150005
#define mod 1000000007
int n,All;
int c[Maxn][Maxn],f[Maxsta],g[Maxsta];
int main()
{
	n=rd(),All=(1<<n)-1;
	for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) c[i][j]=rd();
	for(int i=0;i<=All;i++) g[i]=1;
	for(int s=0;s<=All;s++) for(int i=1;i<=n;i++) for(int j=i+1;j<=n;j++)
		if((s & (1<<(i-1))) && (s & (1<<(j-1))))
			g[s]=1ll*g[s]*(c[i][j]+1ll)%mod;
	for(int s=1,t;s<=All;s++)
	{
		t=s-(s&(-s)),f[s]=g[s];
		for(int i=t;i;i=(i-1)&t)
			f[s]=(f[s]-1ll*g[i]*f[s^i]%mod+mod)%mod;
	}
	printf("%d\n",f[All]);
	return 0;
}
posted @ 2022-08-14 17:36  EricQian06  阅读(26)  评论(0编辑  收藏  举报