La5135 Mining your own Business

题目大意:

给你n条边，让你添加逃生通道，使得不论删除某一个点之后，任何一个点都能到达至少一个逃生通道。

分析：

若任意一个双连通分量中有两个割顶，因为割顶是与其他双连通分量的交点，所以我们只需要考虑只有一个割顶的双连通分量。对于最坏的情况显然是将割顶删去，所以逃生通道不能建在割顶上，可以证明：若逃生通道不在割顶上，即使它被删去，也仍然可以通过割顶到达另一个双连通分量的逃生通道，所以只需特判整个图都是一个简单环的特殊情况即可。

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<stack>
using namespace std;
const int maxm = 50000 + 10;
const int maxn = 100000 + 10;
typedef long long LL;

struct Edge { int u, v; };

int Begin[maxn], Next[maxm*2], To[maxm*2], E;
void Add(int x, int y) {
To[++E] = y;
Next[E] = Begin[x];
Begin[x] = E;
}

stack<Edge> S;
vector<int> bcc[maxn];
bool iscut[maxn];
int pre[maxn], low[maxn], bccno[maxn], bcc_cnt, dfs_clock;

int dfs(int u, int f) {
int child = 0;
int lowu = pre[u] = ++dfs_clock;
for(int i=Begin[u]; i; i=Next[i]) {
int v = To[i];
Edge e = (Edge){u, v};
if(!pre[v]) {
S.push(e);
child++;
int lowv = dfs(v, u);
lowu = min(lowu, lowv);
if(lowv >= pre[u]) {
iscut[u] = true;
bcc[++bcc_cnt].clear();
while(1) {
Edge x = S.top(); S.pop();
if(bccno[x.u] != bcc_cnt) bcc[bcc_cnt].push_back(x.u), bccno[x.u] = bcc_cnt;
if(bccno[x.v] != bcc_cnt) bcc[bcc_cnt].push_back(x.v), bccno[x.v] = bcc_cnt;
if(x.u == u && x.v == v) break;
}
}
}else if(pre[v] < pre[u] && v != f) {
S.push(e);
lowu = min(lowu, pre[v]);
}
}
if(f < 0 && child == 1) iscut[u] = 0;
return lowu;
}

void find_bcc(int n) {
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);
}

int n, m, u, v, Case;
void init() {
n = 0; E = 0;
memset(Begin, 0, sizeof(Begin));
}
int main() {
#ifndef ONLINE_JUDGE
freopen("data.txt", "r", stdin);
freopen("ans.txt", "w", stdout);
#endif
while(scanf("%d", &m) == 1 && m) {
init();
for(int i=1; i<=m; i++) {
scanf("%d%d", &u, &v);
Add(u, v), Add(v, u); n = max(u, max(v, n));
}
find_bcc(n);
LL ans1 = 0, ans2 = 1;
for(int i=1; i<=bcc_cnt; i++) {
int cut_cnt = 0;
for(int j=0; j<bcc[i].size(); j++)
if(iscut[bcc[i][j]]) cut_cnt++;
if(cut_cnt == 1) {
ans1++; ans2 *= (LL)(bcc[i].size() - 1);
}
}
if(bcc_cnt == 1) {
ans1 = 2; ans2 = (LL) bcc[1].size() * (bcc[1].size() - 1) / 2;
}
printf("Case %d: %lld %lld\n", ++Case, ans1, ans2);
}
return 0;
}