P3225 [HNOI2012]矿场搭建

题目描述

煤矿工地可以看成是由隧道连接挖煤点组成的无向图。为安全起见，希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。于是矿主决定在某些挖煤点设立救援出口，使得无论哪一个挖煤点坍塌之后，其他挖煤点的工人都有一条道路通向救援出口。

请写一个程序，用来计算至少需要设置几个救援出口，以及不同最少救援出口的设置方案总数。

输入格式

输入文件有若干组数据，每组数据的第一行是一个正整数 N（N<=500），表示工地的隧道数，接下来的 N 行每行是用空格隔开的两个整数 S 和 T，表示挖 S 与挖煤点 T 由隧道直接连接。输入数据以 0 结尾。

输出格式

输入文件中有多少组数据，输出文件 output.txt 中就有多少行。每行对应一组输入数据的 结果。其中第 i 行以 Case i: 开始（注意大小写，Case 与 i 之间有空格，i 与:之间无空格，: 之后有空格），其后是用空格隔开的两个正整数，第一个正整数表示对于第 i 组输入数据至少需 要设置几个救援出口，第二个正整数表示对于第 i 组输入数据不同最少救援出口的设置方案总 数。输入数据保证答案小于 2^64。输出格式参照以下输入输出样例。

思路：用Tarjan跑出割点，然后DFS搜索所有的联通块，每个连通块是有割点分割，计算连通块的割点数以及总点数，分三种情况：

1、没有割点，说明是一个点双连通图，要两个安全出口，防止只要一个的情况下安全出口坍塌。

2、只有一个割点，那就修一个安全出口，修在非割点位置上，若这个安全出口坏了，还能去割点另一半的安全出口。

3、若有两个以上的割点，则不必修安全出口，因为可以去任何一个割点外的安全出口。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 1e3 + 10;
int n, h[N], cnt, low[N], dfn[N], ind, isc[N], kk, pp, gp;
struct node {
    int v, nt;
} no[N];
void add(int u, int v) {
    no[cnt] = node{v, h[u]};
    h[u] = cnt++;
}
void tarjan(int u, int root) {
    dfn[u] = low[u] = ++ind;
    int son = 0;
    for(int i = h[u]; ~i; i = no[i].nt) {
        int v = no[i].v;
        if(!dfn[v]) {
            tarjan(v, root);
            low[u] = min(low[u], low[v]);
            if(u == root)
                son++;
            if(low[v] >= dfn[u] && u != root)
                isc[u] = 1;
        }
        low[u] = min(low[u], dfn[v]);
    }
    if(u == root && son >= 2)
        isc[root] = 1;
}
void dfs(int u) {
    dfn[u] = gp, pp++;
    for(int i = h[u]; ~i; i = no[i].nt) {
        int v = no[i].v;
        if(dfn[v] != gp && isc[v])
            dfn[v] = gp, kk++;
        if(!dfn[v])
            dfs(v);
    }
}
int main() {
    int t = 0;
    while(~scanf("%d", &n) && n) {
        int ans = 0, r = 0;
        memset(h, -1, sizeof h), cnt = 0, ind = 0;
        memset(isc, 0, sizeof isc);
        memset(dfn, 0, sizeof dfn);
        for(int x, y, i = 1; i <= n; i++) {
            scanf("%d%d", &x, &y);
            add(x, y), add(y, x);
            r = max(r, y), r = max(r, x);
        }
        for(int i = 1; i <= r; i++)
            if(!dfn[i])
                ind = 0, tarjan(i, i);
        ll sum = 1;
        memset(dfn, 0, sizeof dfn);
        for(int i = 1; i <= r; i++)
            if(!dfn[i] && !isc[i]) {
                kk = 0, pp = 0, ++gp;//gp是来判断割点的
                dfs(i);
                if(kk == 0)
                    ans += 2, sum *= (pp * (pp - 1) / 2ll);
                if(kk == 1)
                    ans += 1, sum *= (pp * 1ll);
            }
        printf("Case %d: %d %lld\n", ++t, ans, sum);
    }
    return 0;
}