UVA 1108 - Mining Your Own Business

刘汝佳书上都给出了完整的代码

在这里理一下思路：

由题意知肯定存在一个或者多个双连通分量；

假设某一个双连通分量有割顶。那太平井一定不能打在割顶上。

而是选择割顶之外的随意一个点；

假设没有割顶，则要在该双连通分量上打两个井

至于打井方案。见代码

#include <cstdio>  
#include <cstring>  
#include <vector>  
#include <stack>  
#include <map>  
using namespace std;  
  
const int N = 50005;  
  
struct Edge {  
    int u, v;  
    Edge() {}  
    Edge(int u, int v) {  
        this->u = u;  
        this->v = v;  
    }  
};  
  
int pre[N], bccno[N], dfs_clock, bcc_cnt;  
bool iscut[N];  
  
vector<int> g[N], bcc[N];  
stack<Edge> S;  
  
int dfs_bcc(int u, int fa) {  
    int lowu = pre[u] = ++dfs_clock;  
    int child = 0;  
    for (int i = 0; i < g[u].size(); i++) {  
        int v = g[u][i];  
        Edge e = Edge(u, v);  
        if (!pre[v]) {  
            S.push(e);  
            child++;  
            int lowv = dfs_bcc(v, u);  
            lowu = min(lowu, lowv);  
            if (lowv >= pre[u]) {  
                iscut[u] = true;  
                bcc_cnt++; bcc[bcc_cnt].clear(); //start from 1  
                while(1) {  
                    Edge x = S.top(); S.pop();  
                    if (bccno[x.u] != bcc_cnt) {bcc[bcc_cnt].push_back(x.u); bccno[x.u] = bcc_cnt;}  
                    if (bccno[x.v] != bcc_cnt) {bcc[bcc_cnt].push_back(x.v); bccno[x.v] = bcc_cnt;}  
                    if (x.u == u && x.v == v) break;  
                }  
            }  
        } else if (pre[v] < pre[u] && v != fa) {  
            S.push(e);  
            lowu = min(lowu, pre[v]);  
        }  
    }  
    if (fa < 0 && child == 1) iscut[u] = false;  
    return lowu;  
}  
  
int st;  
  
void find_bcc() {  
    memset(pre, 0, sizeof(pre));  
    memset(iscut, 0, sizeof(iscut));  
    memset(bccno, 0, sizeof(bccno));  
    dfs_clock = bcc_cnt = 0;  
    dfs_bcc(0, -1);  
}  
  
int n, m;  
  
typedef long long ll;  
  
void solve() {  
    ll ans1 = 0, ans2 = 1;  
    for (int i = 1; i <= bcc_cnt; i++) {  
        int cut_cnt = 0;  
        for (int j = 0; j < bcc[i].size(); j++)  
            if (iscut[bcc[i][j]]) cut_cnt++;  
        if (cut_cnt == 1) {  
            ans1++;   
            ans2 *= (ll)(bcc[i].size() - cut_cnt);  
        }  
    }  
    if (bcc_cnt == 1) {  
        ans1 = 2;   
        ans2 = (ll)bcc[1].size() * (bcc[1].size() - 1) / 2;  
    }  
    printf(" %lld %lld\n", ans1, ans2);  
}  
  
int main() {  
    int cas = 0;  
    while (~scanf("%d", &m) && m) {  
        int u, v, Max = 0;  
        while (m--) {  
            scanf("%d%d", &u, &v);  
            u--; v--;  
            g[u].push_back(v);  
            g[v].push_back(u);  
            Max = max(Max, u);  
            Max = max(Max, v);  
        }  
        find_bcc();  
        printf("Case %d:", ++cas);  
        solve();  
        for (int i = 0; i <= Max; i++)  
            g[i].clear();  
    }  
    return 0;  
}