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;
}