[HNOI2012]矿场搭建

[HNOI2012]矿场搭建

题目描述
煤矿工地可以看成是由隧道连接挖煤点组成的无向图。为安全起见，希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。于是矿主决定在某些挖煤点设立救援出口，使得无论哪一个挖煤点坍塌之后，其他挖煤点的工人都有一条道路通向救援出口。

请写一个程序，用来计算至少需要设置几个救援出口，以及不同最少救援出口的设置方案总数。

输入格式
输入文件有若干组数据，每组数据的第一行是一个正整数 N（N<=500），表示工地的隧道数，接下来的 N 行每行是用空格隔开的两个整数 S 和 T，表示挖 S 与挖煤点 T 由隧道直接连接。输入数据以 0 结尾。

输出格式
输入文件中有多少组数据，输出文件 output.txt 中就有多少行。每行对应一组输入数据的 结果。其中第 i 行以 Case i: 开始（注意大小写，Case 与 i 之间有空格，i 与:之间无空格，: 之后有空格），其后是用空格隔开的两个正整数，第一个正整数表示对于第 i 组输入数据至少需 要设置几个救援出口，第二个正整数表示对于第 i 组输入数据不同最少救援出口的设置方案总 数。输入数据保证答案小于 2^64。输出格式参照以下输入输出样例。

输入输出样例
输入 #1复制
9
1 3
4 1
3 5
1 2
2 6
1 5
6 3
1 6
3 2
6
1 2
1 3
2 4
2 5
3 6
3 7
0
输出 #1复制
Case 1: 2 4
Case 2: 4 1
说明/提示
Case 1 的四组解分别是(2,4),(3,4),(4,5),(4,6)；

Case 2 的一组解为(4,5,6,7)。

【题目分析】
割点+乱搞。一道提高难度的题目。
首先不难想到，关键信息是只会坍塌一个点。所以我们必须找到割点，因为割点坍塌的话，会被分为多块，每一块都要有逃生出口。

如图所示，红色的是割点，蓝色区域就是需要设置逃生出口的区域，根据乘法原理，假设有N块逃生区域，每块数量为Fi，则总方案数为${\prod }_1^N$Fi。

接下来就是逃生区域的划分，不难发现，不是所有割点划出的区域都需要设置逃生出口，比如中间的两个区域就不需要（因为可以向左跑到逃生出口，如果左边的点坍塌不能向左，根据只有1个点坍塌的条件，说明其它点没有坍塌，就可以跑到其它的区域）。
这种逃生区域看起来很像一个连通的叶子，就管它叫做叶子连通块吧。
所以分两步:
①求出割点
②求出叶子连通块
至于求法，直接dfs，判断是否只是碰到1个割点，若是，说明是叶子连通块，若不是，说明是中间的某个连通块。

还有一种特殊情况。当整图是一个环（或者没有割点），就考虑设置两个逃生出口（因为要考虑逃生出口点坍塌的情况）。

#include<cstdio>
#include<cstring>
#define min(a,b) ((a)<(b))?a:b
#define max(a,b) ((a)>(b))?a:b
#define maxn 539
using namespace std;
typedef long long ll;
struct EDGE{
	int to, nex;
}edges[maxn];
int head[maxn], k;
void add(int from, int to){edges[k].to = to;edges[k].nex = head[from];head[from] = k++;}
int dfn[maxn], low[maxn], num, gd[maxn], zs[maxn], f[maxn], g[maxn], flag, cnt, t;
ll ans, sum;
void tarjan(int u, int pre){
	dfn[u] = low[u] = ++num;
	for(int i = head[u]; i!=-1; i=edges[i].nex){
		int v = edges[i].to;
		if(v!=pre){
			if(dfn[v]){
				low[u] = min(dfn[v], low[u]);
			}
			else{
				zs[u]++;
				tarjan(v, u);
				low[u] = min(low[u], low[v]);
				if(low[v]>=dfn[u]){
					gd[u] = 1;
				}
			}
		}
	}
	return;
}
void dfs(int u, int pre){
	if(gd[u])if(f[u]!=t)flag++, f[u] = t;
	if(f[u]==t)return;
	cnt++;
	f[u] = t;
	for(int i = head[u]; i!=-1; i=edges[i].nex){
		int v = edges[i].to;
		if(v!=pre&&f[v]!=t){
			dfs(v, u);
		}
	}
	return;
}
ll cal(int x){
	if(x==1)return 1;
	return (ll)x*(x-1)/2;
}
int N, M;
int main(){
	//freopen("1.in", "r", stdin);
	int T = 0;
	while(1){
		memset(head, -1, sizeof(head)); k = 0;
		memset(dfn, 0, sizeof(dfn));
		memset(low, 0, sizeof(low));
		memset(gd, 0, sizeof(gd));
		memset(zs, 0, sizeof(zs));
		memset(f, 0, sizeof(f));
		scanf("%d", &M);
		if(!M)break;
		N = 0;
		int a, b;
		while(M--){scanf("%d%d", &a, &b);add(a, b);add(b, a);N = max(max(a, N), b);}
		for(int i = 1; i < N+1; i++){
			if(!dfn[i]){
				tarjan(i, -1);
				if(zs[i]>1)gd[i] = 1;else gd[i] = 0;
			}
		}
		ans = 1, sum = 0;
		for(int i = 1; i < N+1; i++)
			if(!f[i]&&!gd[i]){
				cnt = 0;
				t++;
				flag = 0;
				dfs(i, -1);
				if(flag==1)ans *= cnt, sum++;
				if(!flag){
					ans = cal(cnt);
					sum = 2;
					break;
				}
			}
		printf("Case %d: %d ", ++T, sum);
		printf("%lld\n", ans);
	}
 	return 0;
}