[poj] 1523 SPF

原题

这是一道点双联通分量。
对于一个割点来说，在搜索他所分割的块与他的连边一定会使这个割点满足割点条件，也就是这时候讲这个割点所分割出来的块+1即可。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define M 500010
#define N 1010
using namespace std;
int n,m,a,b,cnt,sub,head[N],dfn[N],low[N],t,is_AP[N],num;
bool vis[N];
struct hhh
{
    int to,next;
}edge[2*M];

int read()
{
    int ans=0,fu=1;
    char j=getchar();
    for (;(j<'0' || j>'9') && j!='-';j=getchar()) ;
    if (j=='-') j=getchar(),fu=-1;
    for (;j>='0' && j<='9';j=getchar()) ans*=10,ans+=j-'0';
    return ans*fu;
}

void add(int u,int v)
{
    edge[cnt].to=v;
    edge[cnt].next=head[u];
    head[u]=cnt++;
}

void Tarjan(int u,int f)
{
    int child=0;
    dfn[u]=low[u]=++t;
    for (int i=head[u],v;i;i=edge[i].next)
    {
	v=edge[i].to;
	if (v==f) continue;
	if (!dfn[v])
	{
	    child++;
	    Tarjan(v,u);
	    low[u]=min(low[u],low[v]);
	    if (low[v]>=dfn[u] && f!=0)
	    {
		if (!is_AP[u]) num++;
		is_AP[u]++;
	    }
	}
	else low[u]=min(low[u],dfn[v]);
    }
    if (f==0 && child>1) is_AP[u]=child,num++;
    else if (is_AP[u]) is_AP[u]++;
}

int main()
{
    while (~scanf("%d",&a) && a!=0)
    {
	memset(dfn,0,sizeof(dfn));
	memset(vis,0,sizeof(vis));
	memset(head,0,sizeof(head));
	memset(is_AP,0,sizeof(is_AP));
	cnt=1;
	t=0;
	sub++;
	b=read();
	add(a,b);
	add(b,a);
	num=0;
	vis[a]=vis[b]=1;
	n=2;
	while (~scanf("%d",&a) && a!=0)
	{
	    b=read();
	    add(a,b);
	    add(b,a);
	    if (!vis[a]) n++,vis[a]=1;
	    if (!vis[b]) n++,vis[b]=1;
	}
	for (int i=1;i<=n;i++)
	    if (!dfn[i]) Tarjan(i,0);
	printf("Network #%d
",sub);
	if (!num) printf("  No SPF nodes
");
	for (int i=1;i<=n;i++)
	    if (is_AP[i]) printf("  SPF node %d leaves %d subnets
",i,is_AP[i]);
	printf("
");
    }
    return 0;
}