G - Connectivity 2
题意:给你一张无向图,有 n n n 个顶点和 m m m 条边,存在 2 m 2^m 2m 种不同的子图。求子图的数量满足 1 1 1 和 i i i 联通。 n < = 17 n<=17 n<=17。
Solution:
状压 dp 神题。
- 1 1 1 和 i i i 联通等价于 1 1 1 和 i i i 在同一连通块中。
- 看到总状态数
2
m
2^m
2m ,于是想到容斥,不过本题的条件是 至少存在一条
1->i
的路径 ,所以很遗憾暂时没有看到容斥的做法。 - 回到连通块的思考中。启发我们 枚举终止状态中 1 1 1 和 i i i 所处连通块集合,此时大概是一个连通块,内部点没有向外部点连的边。
- 那么我们只要计算出 f ( S ) f(S) f(S) 表示使集合 S S S 联通的图的方案数即可。令 c n t ( S ) cnt(S) cnt(S) 表示集合 S S S 内边 ( x , y ) (x,y) (x,y) 的数量,满足 x , y ∈ S x,y\in S x,y∈S 。那么答案为: a n s ( x ) = ∑ 1 ∈ S , x ∈ S f ( S ) 2 c n t ( U − S ) ans(x)=\sum_{1\in S,x\in S}f(S)2^{cnt(U-S)} ans(x)=1∈S,x∈S∑f(S)2cnt(U−S)
- f ( S ) f(S) f(S) 的状态转移为:(这里的方案数只包括 c n t ( S ) cnt(S) cnt(S) 中边的选择方案) f ( S ) = 2 m − ∑ x ∈ T , T ⊂ S f ( T ) 2 c n t ( S − T ) f(S)=2^{m}-\sum_{x\in T,T\subset S}f(T)2^{cnt(S-T)} f(S)=2m−x∈T,T⊂S∑f(T)2cnt(S−T)
注意到 c n t ( S ) cnt(S) cnt(S) 可以在 O ( 2 n m ) O(2^nm) O(2nm) 内求出, f ( S ) f(S) f(S) 本质上是一个子集枚举。
时间复杂度 O ( 3 n + 2 n m ) O(3^n+2^nm) O(3n+2nm) 。总结:本题贵在状态的设计做到了不重不漏,在转移时巧妙运用了正难则反的思想,因为直接枚举子集会造成图的割边不唯一的情况。
一遍过 qwq… 其实当时不是没想到连通块的思路,就是不知道怎么算联通块的方案数。其实只要把总方案数减去小的连通块方案数就行了。
#include<bits/stdc++.h>
#define INF 1e9
#define ll long long
#define PII pair<ll,int>
#define All(a) a.begin(),a.end()
using namespace std;
const int mx=17;
const int mod=998244353;
inline int read() {
int x=0,f=1; char c=getchar();
while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
while(c>='0'&&c<='9') {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
return x*f;
}
int n,m,u[1<<mx],v[1<<mx],fac[1<<mx];
ll res,cnt[1<<mx],f[1<<mx],ans[mx];
int main() {
scanf("%d%d",&n,&m);
fac[0]=1; for(int i=1;i<=m;i++) fac[i]=fac[i-1]*2%mod;
for(int i=1;i<=m;i++) {
scanf("%d%d",&u[i],&v[i]);
u[i]--,v[i]--;
}
for(int i=0;i<1<<n;i++) {
for(int j=1;j<=m;j++) {
if(i>>u[j]&1 && i>>v[j]&1) {
cnt[i]++;
}
}
}
for(int i=1;i<1<<n;i++) {
f[i]=fac[cnt[i]];
int x;
for(int j=0;j<n;j++) {
if(i>>j&1) {
x=j;
break;
}
}
for(int j=(i-1)&i;j;j=(j-1)&i) {
if((j>>x&1)==0) continue;
f[i]=(f[i]-f[j]*fac[cnt[i-j]]%mod)%mod;
}
}
for(int i=0;i<1<<n;i++) {
if(i%2==0) continue;
ll tmp=f[i]*fac[cnt[(1<<n)-1-i]]%mod;
for(int j=0;j<n;j++) {
if((i>>j&1)==0) continue;
ans[j]=(ans[j]+tmp)%mod;
}
}
for(int i=1;i<n;i++) {
if(ans[i]<0) ans[i]+=mod;
printf("%lld\n",ans[i]);
}
}