HDU 6736 Forest Program
题目
给你一个无向图,保证每条边最多在一个简单环里。问你有多少删边方式使图不包括环。mod 998244353
分析
首先对于一个环(至少三个点)里的边,一定要删去一个,那么如果环里有 n 条边,就有
2
n
−
1
2^n - 1
2n−1 。
再考虑不在环里的边,那么删不删不所谓,如果不再环里的边一共有 m 条边,就有
2
m
2^m
2m。
根据乘法法则,将上面的乘起来就是最终答案。
由于点的双连通分量内部没有割点,因此所有的环就是点双之后点大于三的分量。
套个点双模板即可
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define fuck(x) cout<<x<<endl
const int N = 1e6 + 5;
const ll mod = 998244353;
int n, m;
int ne[N], h[N], to[N], fr[N], idx; // 存图
int pre[N], iscut[N], bccno[N], low[N], dfs_clock, bcc_cnt;
// 开始时间戳、割点、bccno[i] 表示点 i 所在的分量(vector)
vector<int> bcc[N]; // 每个 vector 保存一个点双分量
stack<int> s;
void add(int u, int v){
fr[idx] = u, to[idx] = v, ne[idx] = h[u], h[u] = idx++;
}
void dfs(int u, int fa){
low[u] = pre[u] = ++dfs_clock;
int son = 0;
for(int i = h[u]; ~i; i = ne[i]){
int v = to[i];
if(!pre[v]){
s.push(i), son++;
dfs(v, u);
low[u] = min(low[u], low[v]);
if(low[v] >= pre[u]){ // u 是割点,栈里面的边都是一个点双分量
iscut[u] = 1;
bcc[++bcc_cnt].clear();
while(1){
int x = s.top();
s.pop();
if(bccno[fr[x]] != bcc_cnt){
bcc[bcc_cnt].push_back(fr[x]);
bccno[fr[x]] = bcc_cnt;
}
if(bccno[to[x]] != bcc_cnt){
bcc[bcc_cnt].push_back(to[x]);
bccno[to[x]] = bcc_cnt;
}
if(x == i) break;
}
}
}else if(v != fa && pre[v] < low[u]){
s.push(i);
low[u] = pre[v];
}
}
if(fa < 0 && son == 1)
iscut[u] = 0;
}
void find_bcc(int n){ // 点双模板
memset(low, 0, sizeof(low));
memset(pre, 0, sizeof(pre));
memset(iscut, 0, sizeof(iscut));
memset(bccno, 0, sizeof(bccno));
dfs_clock = bcc_cnt = 0;
for(int i = 1; i <= n; i++){
if(!pre[i]){
dfs(i, -1);
}
}
}
ll pp(ll x, ll y){
ll res = 1;
while(y){
if(y&1)
res = res * x % mod;
x = x * x % mod;
y >>= 1;
}
return res;
}
int main(){
while(~scanf("%d%d", &n, &m)){
memset(h, -1, sizeof(h)), idx = 0;
for(int i = 0, u, v; i < m; i++){
scanf("%d%d", &u, &v);
add(u, v);
add(v, u);
}
find_bcc(n);
ll ans = 1, num = 0;
for(int i = 1; i <= bcc_cnt; i++){
if(bcc[i].size() >= 3){ // 当这个分量有超过三个点,就是一个环
num += bcc[i].size();
ans = ans * (pp(2, bcc[i].size()) - 1) % mod;
}
}
printf("%lld\n", ans * pp(2, m - num) % mod);
}
}