BZOJ 2730 Tarjan求割点

        首先我们考虑只有一个割点的情况，那么一定是在割点分开的这两个连通块中各建一个，如果有多个割点，那么我们一定是在只被一个割点控制的那个连通块内建一个，因为如果一个连通块被多个割点控制，任何一个割点坍塌，这个连通块的人都能通过其他割点到别的连通块去，特殊的，如果一张图里没有只被一个割点控制的连通块，那么我们只需随便建两个就行，建两个是因为只建一个的话如果建的那个坍塌别人就没地方去了，至于方案数可以用乘法原理来算。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
using namespace std;
#define maxn 50000
int last[maxn],pre[maxn],other[maxn];
int low[maxn],dfn[maxn],cut[maxn],scc[maxn],num[maxn],siz[maxn];
int n,m,dfs_clock,l,tot;
long long ans1,ans2;
bool flag[maxn],f[maxn];

void init(void)
{
	memset(last,0,sizeof last);
	memset(dfn,0,sizeof dfn);
	memset(cut,0,sizeof cut);
	memset(scc,0,sizeof scc);
	memset(num,0,sizeof num);
	memset(siz,0,sizeof siz);
	memset(flag,0,sizeof flag);	
	memset(f,0,sizeof f);
	memset(siz,0,sizeof siz);
	dfs_clock=l=n=ans1=0;
	ans2=1;
}

void connect(int x,int y)
{
	l++;
	pre[l]=last[x];
	last[x]=l;
	other[l]=y;	
}

void dfs(int u,int fa)
{
	int cnt=0;
	dfn[u]=low[u]=++dfs_clock;
	for (int p=last[u];p;p=pre[p]) 
	{
		int v=other[p];
		if (v==fa) continue;
		if (!dfn[v]) 
		{
			dfs(v,u);
			cnt++;
			low[u]=min(low[u],low[v]);
			if (dfn[u]<=low[v]&&fa!=-1) cut[u]=1;
		}
		else low[u]=min(low[u],dfn[v]);
	}	
	if (fa==-1&&cnt>1) cut[u]=1;
}	

void NB_dfs(int u,int fa) 
{
	for (int p=last[u];p;p=pre[p]) 
	{
		int v=other[p];
		if (v==fa||cut[v]) continue;
		if (!scc[v]) 
		{
			scc[v]=scc[u];
			NB_dfs(v,u);
		}
	}
}

int main()
{
	while (1) 
	{
		++tot;
		init();	
		scanf("%d",&m);
		if (m==0) break;
		for (int i=1;i<=m;i++) 
		{
			int x,y;
			scanf("%d%d",&x,&y);
			connect(x,y);
			connect(y,x);
			f[x]=f[y]=1;
			n=max(n,x);
			n=max(n,y);
		}	
		for (int i=1;i<=n;i++) 
			if (!dfn[i]&&f[i]) dfs(i,-1);
		for (int i=1;i<=n;i++) 
			if (!scc[i]&&!cut[i]) 
			{
				scc[i]=i;
				NB_dfs(i,-1);
			}
		for (int i=1;i<=n;i++) 
			if (scc[i]) siz[scc[i]]++;
		for (int i=1;i<=n;i++) 
			if (cut[i]) 
			{
				memset(flag,0,sizeof flag);
				for (int p=last[i];p;p=pre[p]) 
				{
					int v=other[p];
					if (!flag[scc[v]]) 
					{
						num[scc[v]]++;
						flag[scc[v]]=1;
					}
				}
			}
		for (int i=1;i<=n;i++) 
			if (num[i]==1) 
			{
				ans1++;
				ans2=ans2*(long long)siz[i];
			}
		if (ans1==0) 
		{
			ans2=0;
			for (int i=1;i<=n;i++) 
				if (f[i]) ans2++;
		}
		if (ans1) printf("Case %d: %lld %lld\n",tot,ans1,ans2);
		else printf("Case %d: 2 %lld\n",tot,ans2*(ans2-1)/2);
	}
	return 0;	
}