HDU 6736【2019秦皇岛ccpc】 F. Forest Program(点的双连通分量)

21 篇文章 0 订阅
7 篇文章 0 订阅

HDU 6736 Forest Program

题目

给你一个无向图,保证每条边最多在一个简单环里。问你有多少删边方式使图不包括环。mod 998244353

分析

首先对于一个环(至少三个点)里的边,一定要删去一个,那么如果环里有 n 条边,就有 2 n − 1 2^n - 1 2n1
再考虑不在环里的边,那么删不删不所谓,如果不再环里的边一共有 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);
	}
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值