【BZOJ2730】[HNOI2012]矿场搭建 Tarjan

【BZOJ2730】[HNOI2012]矿场搭建

Description

煤矿工地可以看成是由隧道连接挖煤点组成的无向图。为安全起见，希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。于是矿主决定在某些挖煤点设立救援出口，使得无论哪一个挖煤点坍塌之后，其他挖煤点的工人都有一条道路通向救援出口。请写一个程序，用来计算至少需要设置几个救援出口，以及不同最少救援出口的设置方案总数。

Input

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

Output

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

Sample Input

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

Sample Output

Case 1: 2 4
Case 2: 4 1

HINT

Case 1 的四组解分别是(2,4),(3,4),(4,5),(4,6)；

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

题解：先Tarjan求点双，然后得到一个森林，对于一个森林，我们其实只需要在它的每个叶子中都放一个（不能是割点）。但是注意：如果一个连通块中只有一个点双，那么这个点双里要放两个。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int n,top,tot,cnt,Cnt,sum,a1,cas;
int To[1010],Next[1010],Head[1010],dep[1010],low[1010],to[2010],next[2010],head[2010],siz[2010],sta[1010],size[2010];
unsigned long long a2;
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+gc-'0',gc=getchar();
	return ret*f;
}
void Add(int a,int b)
{
	To[Cnt]=b,Next[Cnt]=Head[a],Head[a]=Cnt++;
}
void add(int a,int b)
{
	to[cnt]=b,next[cnt]=head[a],head[a]=cnt++;
}
void tarjan(int x)
{
	dep[x]=low[x]=++tot,sta[++top]=x;
	for(int i=Head[x],y,t;i!=-1;i=Next[i])
	{
		y=To[i];
		if(!dep[y])
		{
			tarjan(y),low[x]=min(low[x],low[y]);
			if(low[y]>=dep[x])
			{
				sum++;
				do
				{
					t=sta[top--],add(sum,t),add(t,sum),size[sum]++;
				}while(t!=y);
				add(sum,x),add(x,sum),size[sum]++;
			}
		}
		else	low[x]=min(low[x],dep[y]);
	}
}
void dfs(int x,int fa)
{
	siz[x]=(x>1000);
	int i,tmp=0;
	for(i=head[x];i!=-1;i=next[i])	if(to[i]!=fa)
	{
		dfs(to[i],x),siz[x]+=siz[to[i]];
		if(siz[to[i]])	siz[x]+=siz[to[i]],tmp++;
	}
	if(x>1000)
	{
		if(!fa)
		{
			if(siz[x]==1)
			{
				if(size[x]==1)	a1++;
				else	a1+=2,a2*=size[x]*(size[x]-1)>>1;
			}
			else	if(tmp==1)	a1++,a2*=size[x]-1;
		}
		else	if(!tmp)	a1++,a2*=size[x]-1;
	}
}
int main()
{
	while(1)
	{
		n=rd();
		if(!n)	return 0;
		int i,a,b;
		memset(head,-1,sizeof(head)),memset(Head,-1,sizeof(Head)),cnt=Cnt=top=tot=0;
		memset(dep,0,sizeof(dep)),memset(low,0,sizeof(low)),memset(size,0,sizeof(size)),memset(siz,0,sizeof(siz));
		for(i=1;i<=n;i++)	a=rd(),b=rd(),Add(a,b),Add(b,a);
		sum=1000,a1=0,a2=1;
		for(i=1;i<=1000;i++)	if(!dep[i]&&Head[i]!=-1)	tarjan(i);
		for(i=1001;i<=sum;i++)	if(!siz[i])	dfs(i,0);
		printf("Case %d: %d %llu\n",++cas,a1,a2);
	}
}