HDU 2460 Network(Tarjan求解桥+LCA)

题意:求出桥数,然后有Q次填边,每次填边后输出剩余的桥数。

一开始想的是先求解所有的桥,然后双连通缩点建树,再添边求解桥。。感觉很麻烦。。后来发现不用这么麻烦,Tarjan搜索一次后,记录到了所有的dfn值和fa值,即所有点的访问次序以及父亲值,那么借用LCA的思想,图中的边无非是树边与回退边,对每对Q填边操作的u和v,进行一次LCA,搜索树边与回退边。对树边的话,u和v的dfn值一大一小,对大的那个不断递归其父亲值,最终u,v一定会相等,即达到最小公共祖先。对这种情况,添边(u,v)后会形成环,那么环内原先所有是桥的树边都变为非桥,标记一下,从总数中减掉。对回退边,那就需要同时对u,v操作,但是检查的方式还是一样的,仍然是检测是否有桥,检测到后标记一下,从总桥数中减去。如果u,v本来就在一个双连通分量里呢?。。那显然是检测不到桥的。。所以没有什么影响。

代码:

// Header.
#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <string>
#include <queue>
#include <stack>
#include <cmath>
#include <set>
#include <map>
using namespace std;

typedef long long LL;
#define mem(a, n) memset(a, n, sizeof(a))
#define rep(i, n) for(int i = 0; i < (n); i ++)
#define REP(i, t, n) for(int i = (t); i < (n); i ++)
#define FOR(i, t, n) for(int i = (t); i <= (n); i ++)
#define ALL(v) v.begin(), v.end()
#define si(a) scanf("%d", &a)
#define sii(a, b) scanf("%d%d", &a, &b)
#define siii(a, b, c) scanf("%d%d%d", &a, &b, &c)
#define pb push_back
const int inf = 0x3f3f3f3f, N = 1e5 + 5, MOD = 1e9 + 7;

int T, cas = 0;
int n, m, q;
int dfn[N], low[N], fa[N], head[N];
bool brige[N], vis[N];
struct edge {
	int v, next;
}e[4 * N];
int brigeNum = 0, ne = 0, dfsNum = 0;
void addEdge(int u, int v) {
	e[ne].v = v;
	e[ne].next = head[u];
	head[u] = ne ++;
}
// Imp
void init() {
	mem(dfn, 0);
	mem(low, 0);
	mem(fa, -1);
	mem(vis, 0);
	mem(brige, 0);
	mem(head, -1);
    dfsNum = ne = brigeNum = 0;
}

void Tarjan(int u, int f) {
	dfn[u] = low[u] = ++ dfsNum;
	for(int i = head[u]; i != -1; i = e[i].next) {
		int v = e[i].v;
		if(!dfn[v]) {
			vis[v] = 1;
			fa[v] = u;
			Tarjan(v, u);
			low[u] = min(low[u], low[v]);
			if(dfn[u] < low[v]) {
				brigeNum ++;
				brige[v] = 1;
			}
		} else if(v != f) low[u] = min(low[u], dfn[v]);
	}
}

void work(int u) {
	if(brige[u]) brigeNum --, brige[u] = 0;
}

void LCA(int u, int v) {
	if(dfn[u] > dfn[v]) swap(u, v);
	while(dfn[u] < dfn[v]) {
		work(v);
		v = fa[v];
	}
	while(v != u) {
		work(u);
		work(v);
		u = fa[u], v = fa[v];
	}
}

int main(){
#ifdef LOCAL
    freopen("/Users/apple/input.txt", "r", stdin);
//	freopen("/Users/apple/out.txt", "w", stdout);
#endif
	
    while(sii(n, m) != EOF, n + m) {
    	init();
    	int u, v;
    	rep(i, m) {
    		sii(u, v);
    		addEdge(u, v), addEdge(v, u);
    	}
    	FOR(i, 1, n + 1) fa[i] = i;
    	Tarjan(1, -1);
    	printf("Case %d:\n", ++ cas);
    	si(q);
    	rep(i, q) {
    		sii(u, v);
    		LCA(u, v);
    		printf("%d\n", brigeNum);
    	} puts("");

    }
    
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值