zoj 1119 / poj 1523 SPF (典型例题 求割点 Tarjan 算法)

// File Name    :poj1523.cpp
// Author        :Freetion
// Created Time    :2013年09月11日 星期三 22时40分42秒

#define LOCAL  //Please annotate this line when you submit
/********************************************************/
#include <iostream>
#include <stdio.h>
#include <math.h>
#include <algorithm>
#include <string.h>
#include <string>
#include <map>
#define CLE(name) memset(name, 0, sizeof(name))
using namespace std;

const int max_n = 1001;
int edge[max_n][max_n]; //边的邻接矩阵
int vis[max_n];//顶点访问状态
int node; //记录最大的节点号
int tmpdfn;  // 在dfs中记录当前深度优先数
int dfn[max_n]; // 记录每个顶点的深度优先数
int low[max_n]; // 记录每个顶点通过回边或邻边所能达到的最低深度优先数
int son; // 根的子女数
int subnet[max_n]; // 记录每个节点（去掉该节点后）的连通分量个数

//深度优先搜索，求每个节点的low值（根据low值判断是否为割点）
void dfs(int u)
{
for (int v = 1; v <= node; v ++)
{
if (edge[u][v]) // 边存在
{
if (!vis[v]) // 没有被访问过，则v是u的儿子节点
{
vis[v] = 1; //标记
dfn[v] = low[v] = ++ tmpdfn; //dfs值 和 初始的low值
dfs(v); // 递归
low[u] = min(low[u], low[v]); //求low的第二种情况
if (low[v] >= dfn[u])
{
if (u != 1)
subnet[u] ++;
else son ++;
//因为跟节点没有进入的边，只有出的边所以需要单独处理一下
}
}
else low[u] = min(low[u], dfn[v]); // 求low第三种情况的
}
}
}

int main()
{
int i, u, v, find, num = 1;
while (~scanf ("%d", &u) && u)
{
CLE(edge);
node = 0;
scanf ("%d", &v);
if (u > node) node = u;
if (v > node) node = v;
edge[u][v] = edge[v][u] = 1;
while (scanf ("%d", &u) && u)
{
scanf ("%d", &v);
if (u > node)
node = u;
if (v > node)
node = v;
edge[u][v] = edge[v][u] = 1;//标记边存在
}
//以下是初始化
dfn[1] = low[1] = tmpdfn = 1;
son = 0;
CLE(vis);    CLE(subnet);
vis[1] = 1;
//以上是初始化
dfs(1);  // 默认根节点为 1；
if (son > 1) //单独处理根
subnet[1] = son -1;
find = 0;
if (num > 1)
puts("");
printf ("Network #%d\n", num ++);
for (int i = 1; i <= node; i ++)
{
if (subnet[i])
{
find = 1;
printf ("  SPF node %d leaves %d subnets\n", i, subnet[i] +1);
//因为有一条入边所以还要加上一个 1；
}
}
if (find == 0)
printf ("  No SPF nodes\n");
}
return 0;
}