poj 1523 SPF

　　这是割顶题中的极为基础的一道，就是用了《数据结构》中的关节点(P178)，或者说是刘汝佳的《算法艺术与信息学竞赛》中的割顶(P285)。我的理解和解释，就是如果当前结点的某个子结点的祖先的最低遍历编号不小于当前结点的遍历编号（这个解释必须这么长，因为无论是哪本书，对此的解释都是差不多的），那么当前结点就是所要找的割顶，而且割的是之前提及的某个子结点与在当前结点前被遍历的祖先的联系。

　　这题可以直接套《数据结构》的FindArticul、DFSArticul函数，但是还要做的是思考怎样统计割顶的把整个图割成子连通分量的个数，因为题目要求求出割顶的同时，也要输出割顶把图分成多少部分。如果没有，就输出没有。

　　这题要注意的是如果只输入一个0，这个不是一组数据，是不应该被处理的。我一开始wa的几次就只是因为这个，于是后来只好找题解帮忙。可是在我看到解题代码前，我找到一个标准数据，直接发现这个无聊的，但是poj又没有提及的注意事项。

我的代码除了头文件比较长，其他的是模仿《数据结构》的算法的，所以应该还是挺容易理解的：

#include "cstdio"
#include "cstdlib"
#include "cstring"
#include "cmath"
#include "cctype"
#include "vector"
#include "set"
#include "map"
#include "string"
#include "algorithm"
#include "stack"
#include "queue"

#define INF 0x7fffffff
#define reset(a) memset(a, 0, sizeof(a))
#define copy(a, b) memcpy(a, b, sizeof(b))
#define PI acos(-1)
#define FMAX (1E300)
#define MAX 1000000000
#define feq(a, b) (fabs((a)-(b))<1E-6)
#define flq(a, b) ((a)<(b)||feq(a, b))
#define MAXN 10005
#define BASE 137
#define PASS puts("pass")
#define filein freopen("test.in", "r", stdin)
#define fileout freopen("test.out", "w", stdout)

using namespace std;

int e[1001][1001];
int cnt, low[1001], vis[1001];
int spf[1001];
set<int> S;

void dfs(int node){
    int min = ++cnt;

    vis[node] = cnt;
    for (int i = 1; i <= e[node][0]; i++){
        if (!vis[e[node][i]]){
            dfs(e[node][i]);
            if (min > low[e[node][i]])
                min = low[e[node][i]];
            if (low[e[node][i]] >= vis[node])
                spf[node]++;
        }
        else if (min > vis[e[node][i]])
            min = vis[e[node][i]];
    }
    low[node] = min;
}

int main(){
    int a, b, CASE = 1;

    //filein;
    //fileout;
    reset(e);
    S.clear();
    while (1){
        if (!~scanf("%d", &a))
            return 0;
        if (a){
            scanf("%d", &b);
            e[a][++e[a][0]] = b;
            e[b][++e[b][0]] = a;
            S.insert(a);
            S.insert(b);
        }
        else{
            if (!S.size())
                continue;
            if (CASE != 1)
                printf("\n");
            reset(vis);
            reset(spf);
            reset(low);
            bool pnt = false;
            cnt = 1;
            printf("Network #%d\n", CASE++);

            int begin;

            for (begin = 1; begin <= 1000; begin++)
                if (S.count(begin))
                    break;
            low[begin] = vis[begin] = 1;
            for (int i = 1; i <= e[begin][0]; i++){
                if (!vis[e[begin][i]]){
                    dfs(e[begin][i]);
                    if (cnt < S.size())
                        spf[begin]++;
                }
            }
            for (int i = 1; i <= 1000; i++){
                if (spf[i]){
                    printf("  SPF node %d leaves %d subnets\n", i, spf[i] + 1);
                    pnt = true;
                }
            }
            if (!pnt)
                printf("  No SPF nodes\n");
            reset(e);
            S.clear();
        }
    }
}

——written by Lyon