P5933 [清华集训2012]串珠子

2022-07-29 12:02:34 阅读：148 来源： 互联网

标签：联通 P5933 rep 个数子集 2012 集训 define

题意

给定一张 \(n\) 个点的图，其中 \(i\) 和 \(j\) 两点间有 \(c_{i,j}\) 种边可以连。求把这 \(n\) 个点连成连通块的方案数是多少。

Solution

还是考虑拍在脸上的状压。

令 \(f_S\) 表示点集 \(S\) 中的点联通图的个数。如果我们考虑 \(c_{i,j}=1\)，那么容易想到这就是考虑有多少联通图，那么以前做过，可以用生成函数。当时是令 \(g_S\) 表示点集 \(S\) 构成不一定联通的图的个数。这题也差不多。首先 \(g\) 随便求：

\[g_S=\prod_{i,j\in S}(c_{i,j}+1) \]

然后在求 \(f_S\) 的时候，你先求 \(S\) 中点不联通的图的个数，然后用 \(g_S\) 减去它。后者可以先枚举出 \(S\) 的一个子集 \(T\)，然后假设这个子集是联通的，那么方案数是 \(f_T\)。但是接下来你会发现这里算重了！

\(\color{red}{\bigstar}\) 为了保证不算重，可以固定一个点 \(p\) 在集合 \(T\) 中，也就是说我们只枚举包含 \(p\) 的子集。这样子就不会有重复的情况了。

那对于 \(S-T\) 的部分就可以用 \(g_{S-T}\) 了。所以有：

\[f_S=g_{S}-\sum_{p\in T\subset S}f_Tg_{S-T} \]

Code

// Problem: P5933 [清华集训2012]串珠子
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P5933
// Memory Limit: 500 MB
// Time Limit: 1000 ms

#include<bits/stdc++.h>
// #define ll long long
#define inf (1<<30)
#define INF (1ll<<60)
#define pb emplace_back
#define pii pair<int,int>
#define mkp make_pair
#define fi first
#define se second
#define all(a) a.begin(),a.end()
#define rep(i,j,k) for(int i=(j);i<=(k);i++)
#define per(i,j,k) for(int i=(j);i>=(k);i--)
#define pt(a) cerr<<#a<<'='<<a<<' '
#define pts(a) cerr<<#a<<'='<<a<<'\n'
#define int long long
using namespace std;
const int MOD=1e9+7;
const int MAXN=1e6+10;
int c[20][20],f[MAXN],g[MAXN];
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	int n;cin>>n;
	rep(i,1,n) rep(j,1,n) cin>>c[i][j];
	rep(s,0,(1<<n)-1){
		g[s]=1;
		rep(i,1,n){
			if(!(s&(1<<(i-1)))) continue;
			rep(j,i+1,n)
				if(s&(1<<(j-1)))
					g[s]=g[s]*(c[i][j]+1)%MOD;
		}
	}
	rep(s,1,(1<<n)-1){
		int p=s&(-s),sub=0;
		for(int t=(s-1)&s;t;t=(t-1)&s){
			if(!(t&p)) continue;
			sub=(sub+f[t]*g[s^t]%MOD)%MOD;
		}f[s]=(g[s]-sub+MOD)%MOD;
	}
	cout<<f[(1<<n)-1]<<'\n';
	return 0;
}

标签：联通,P5933,rep,个数,子集,2012,集训,define
来源： https://www.cnblogs.com/ZCETHAN/p/16531791.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

P5933 [清华集训2012]串珠子

题意

Solution

Code