【HNOI2012】矿场搭建

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


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


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


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


样例输入 
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
样例输出 
Case 1: 2 4
Case 2: 4 1
提示
Case 1 的四组解分别是(2,4),(3,4),(4,5),(4,6)；

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

分析：先开始看到这个题想到的是缩点，因为第一个样例很具有这方面的暗示性，后来发现这是与割点有关而不是割桥，所以并不适用。

将此题抽象为数学模型就是，要让每个点满足对于一些特定的点联通。

将割点求出后，删除割点，得到了若干个连通块。，dfs求出连通块所连的割点数量。对于每个连通块，只要被2个割点（或2个以上）连接过，就不用管了。

因为一个割点炸了，将它与图的其他部分分开后，还可以从另一个割点通过。所以其实对于每个连通块需要两个割点。

对于没有割点的连通块，选任意两个点建立，组合求出种数，对于只有一个割点的，只需随意选一个建即可。

当然，对于不是割点的点炸了也无所谓（只要不是一个单的点。但是貌似测试数据里面也没有这种情况），因为总有道路可以逃离这一块

#include<bits/stdc++.h>
using namespace std;

#define N 10001

int n,m,cnt,cot,tot,group_num,cut,root,num=0;
long long ans1,ans2;
int first[N],dfn[N],low[N],flag[N],vis[N];

struct email
{
	int u,v;
	int nxt;
}e[N*10];

void add(int u,int v)
{
	e[++cnt].nxt=first[u];first[u]=cnt;
	e[cnt].u=u;e[cnt].v=v;
}

void init()
{
	cnt=0;cot=0;ans1=0;ans2=1;n=0;
	memset(first,0,sizeof(first));
	memset(low,0,sizeof(low));
	memset(dfn,0,sizeof(dfn));
	memset(flag,0,sizeof(flag));
	memset(vis,0,sizeof(vis));
}

void tarjan(int u,int fa)
{
	int son=0;
	dfn[u]=low[u]=++cot;
	for(int i=first[u];i;i=e[i].nxt)
	{
		int v=e[i].v;
		if(!dfn[v])
		{
			son++;
			tarjan(v,u);
			low[u]=min(low[u],low[v]);
			if(low[v]>=dfn[u])
				flag[u]=1;
		}	
		else 
			if(v!=fa)
				low[u]=min(low[u],dfn[v]);
		if(u==root&&son==1)
			flag[u]=0;
	}
}

void dfs(int u)
{
	vis[u]=group_num;
	tot++;
	for(int i=first[u];i;i=e[i].nxt)
	{
		int v=e[i].v;
		if(vis[v]!=group_num&&flag[v])
		{
			cut++;
			vis[v]=group_num;
		}
		if(!vis[v])
			dfs(v);
	}
} 

void calculate(int cut)
{
	if(!cut) 
	{
		ans1+=2;
		ans2*=(tot-1)*tot/2;
	}
	if(cut==1)
	{
		ans1+=1;
		ans2*=tot;
	}	
	return ;	
}

int main()
{
	while(scanf("%d",&m)&&m)
	{
		num++;
		init();
		for(int i=1;i<=m;i++)
		{
			int x,y;
			scanf("%d%d",&x,&y);
			add(x,y);add(y,x);
			n=max(n,max(x,y));
		}
		for(int i=1;i<=n;i++)
		{
			if(!dfn[i]) 
			{
				root=i;
				tarjan(i,i);
			}
		}
		for(int i=1;i<=n;i++)
		{
			if(!vis[i]&&!flag[i])
			{
				group_num++;
				tot=0;cut=0;
				dfs(i);
				calculate(cut);
			} 	
		}
		printf("Case %d: %lld %lld\n",num,ans1,ans2);
	}
	return 0;
}