[Luogu] P3225 [HNOI2012]矿场搭建

题目描述

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

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

题目解析

tarjan找双联通分量，找割点。

在一个分量里有 ≥2 个割点，这个分量就怎么都能跑出去，不需要建出口。

在一个分量里有 1 个割点，这个分量就要防止割点塌了，需要建1个出口。

在一个分量里没有割点，说明它不和别的分量连通，为了防止出口塌掉，要建两个出口。

看起来有些坑，值得注意的是，一个点只会在一个强连通分量，但可能同时处于多个双联通分量

Code

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

const int MAXN =  + ;

struct Edge {
    int nxt;
    int to;
} l[MAXN<<];

int n,m,res,T;
int tot,stamp;
int head[MAXN],cnt;
int low[MAXN],dfn[MAXN];
int sum,num;
int root,deg;
int vis[MAXN],col[MAXN];
bool cut[MAXN];
long long ans1,ans2 = ;

void tarjan(int x,int from) {
    dfn[x] = low[x] = ++stamp;
    for(int i = head[x];i;i = l[i].nxt) {
        if(!dfn[l[i].to]) {
            tarjan(l[i].to,x);
            low[x] = min(low[x],low[l[i].to]);
            if(low[l[i].to] >= dfn[x]) {
                if(x == root) deg++;
                else cut[x] = true;
            }
        } else if(l[i].to != from) low[x] = min(low[x],dfn[l[i].to]);
    }
    return;
}

void dfs(int x) {
    vis[x] = tot;
    if(cut[x]) return;
    sum++;
    for(int i = head[x];i;i = l[i].nxt) {
        if(cut[l[i].to] && vis[l[i].to] != tot) num++,vis[l[i].to] = tot;
        if(!vis[l[i].to]) dfs(l[i].to);
    }
}

inline void add(int x,int y) {
    cnt++;
    l[cnt].nxt = head[x];
    l[cnt].to = y;
    head[x] = cnt;
    return;
}

inline void clean() {
    memset(low,,sizeof(low));
    memset(dfn,,sizeof(dfn));
    memset(head,,sizeof(head));
    memset(col,,sizeof(col));
    memset(cut,,sizeof(cut));
    memset(vis,,sizeof(vis));
    num = sum = ;
    ans1 = n = tot = stamp = cnt = ;
    ans2 = ;
}

int main() {
    while(~scanf("%d",&m) && m) {
        clean();
        int x,y;
        for(int i = ;i <= m;i++) {
            scanf("%d%d",&x,&y);
            n = max(n,max(x,y));
            add(x,y),add(y,x);
        }
        for(int i = ;i <= n;i++) {
            if(!dfn[i]) {
                deg = ;
                tarjan(root = i,);
                if(deg >= ) cut[root] = true;
            }
        }
        for(int i = ;i <= n;i++) {
            if(!vis[i] && !cut[i]) {
                tot++;
                sum = num = ;
                dfs(i);
                if(!num && sum > ) ans1 += , ans2 *= sum*(sum-)/;
                else if(num == ) ans1++, ans2 *= sum;
            }
        }
        printf("Case %d: %lld %lld\n",++T,ans1,ans2);
    }
    return ;
}