HNOI 2012 矿场搭建 题解

题目传送门

题目大意： 给出一张图，现在让你在某些点上建救援出口，要求在随便一个点坍塌（即被删除）后剩余的点都可以到达一个救援出口，问至少建几个救援出口以及方案数。

题解

很显然能够联想到点双，我们考虑一个点双里面的割点数量：

0个时 说明这个点双与其他的点双之间没有联系，那么这个点双内部至少要建两个救援出口，因为是点双，所以建在哪里都可以，不能只建一个，因为如果恰好是那个点塌了就完了。

1个时 说明这个点双与另一个点双之间有联系，那么这个点双里面就只需要建一个救援出口，如果这个出口塌了，剩下的点也可以去另一个点双里面找救援出口，如果那个割点塌了，那么这个点双内的点去这个点双里面的救援出口就好了。

大于等于2个时 无论哪一个点塌了，这个点双里面的点都可以去其他点双里面找救援出口，所以这个点双里面不需要建救援出口。

然后还有一个要注意的，不能在割点上建救援出口，想想就明白了。然后剩下的方案数随便搞搞就出来了，并不是重点。

代码：

#include <cstdio>
#include <cstring>
#include <map>
#include <algorithm>
using namespace std;
#define maxn 510

int n,m,s[maxn<<1],t;
map<int,int> mp;
struct rd{int x,y;};
rd road[maxn];
struct edge{int x,y,next;};
edge e[maxn<<1];
int first[maxn],len;
void buildroad(int x,int y)
{
	e[++len]=(edge){x,y,first[x]};
	first[x]=len;
}
int dfn[maxn],low[maxn],id;
bool cut[maxn];
void dfs1(int x,int fa)
{
	dfn[x]=low[x]=++id; int son=0;
	for(int i=first[x];i;i=e[i].next)
	{
		int y=e[i].y;
		if(y==fa)continue;
		if(!dfn[y])
		{
			dfs1(y,x); son++;
			if(low[y]<low[x])low[x]=low[y];
			if(low[y]>=dfn[x])cut[x]=true;
		}
		else if(dfn[y]<low[x])low[x]=dfn[y];
	}
	if(fa==-1&&son==1)cut[x]=false;
}
long long ans1,ans2;
struct par{
	int size,point;//记录非割点数和割点数
	par(int x,int y):size(x),point(y){}
	par operator +(const par &b){return par(size+b.size,point+b.point);}
};
par dfs2(int x,int fa)
{
	if(cut[x])return par(0,1);
	par re(1,0);
	for(int i=first[x];i;i=e[i].next)
	{
		if(e[i].y==fa)continue;
		if(dfn[e[i].y]!=dfn[x])dfn[e[i].y]=dfn[x],re=re+dfs2(e[i].y,x);
	}
	return re;
}
int tot;

int main()
{
	while(scanf("%d",&m),m!=0)
	{
		tot++;n=t=0;
		memset(first,0,sizeof(first));len=1;
		for(int i=1;i<=m;i++)
		scanf("%d %d",&road[i].x,&road[i].y),s[++t]=road[i].x,s[++t]=road[i].y;
		sort(s+1,s+t+1);
		for(int i=1;i<=t;i++)
		if(s[i]!=s[i-1])mp[s[i]]=++n;
		for(int i=1;i<=m;i++)
		buildroad(mp[road[i].x],mp[road[i].y]),buildroad(mp[road[i].y],mp[road[i].x]);
		
		memset(dfn,0,sizeof(dfn));
		memset(low,0,sizeof(low));id=0;
		memset(cut,false,sizeof(cut));
		for(int i=1;i<=n;i++)
		if(dfn[i]==0)dfs1(i,-1);
		ans1=0;ans2=1;
		memset(dfn,0,sizeof(dfn));
		for(int i=1;i<=n;i++)
		if(dfn[i]==0&&!cut[i])
		{
			dfn[i]=i; par p=dfs2(i,-1);
			if(p.point==0)ans1+=min(p.size,2),ans2*=max(p.size*(p.size-1ll)/2ll,1ll);
			if(p.point==1)ans1+=1,ans2*=p.size;
		}
		printf("Case %d: %lld %lld\n",tot,ans1,ans2);
	}
}