acwing-396. 矿场搭建(Tarjan点双连通分量)

煤矿工地可以看成是由隧道连接挖煤点组成的无向图。

为安全起见，希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。

于是矿主决定在某些挖煤点设立救援出口，使得无论哪一个挖煤点坍塌之后，其他挖煤点的工人都有一条道路通向救援出口。

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

输入格式
输入文件有若干组数据，每组数据的第一行是一个正整数 N，表示工地的隧道数。

接下来的 N 行每行是用空格隔开的两个整数 S 和 T，表示挖煤点 S 与挖煤点 T 由隧道直接连接。

注意，每组数据的挖煤点的编号为 1∼Max，其中 Max 表示由隧道连接的挖煤点中，编号最大的挖煤点的编号，可能存在没有被隧道连接的挖煤点。

输入数据以 0 结尾。

输出格式
每组数据输出结果占一行。

其中第 i 行以 Case i: 开始（注意大小写，Case 与 i 之间有空格，i 与 : 之间无空格，: 之后有空格）。

其后是用空格隔开的两个正整数，第一个正整数表示对于第 i 组输入数据至少需要设置几个救援出口，第二个正整数表示对于第 i 组输入数据不同最少救援出口的设置方案总数。

输入数据保证答案小于 264，输出格式参照以下输入输出样例。

数据范围
1≤N≤500，
1≤Max≤1000

输入样例：
9
1  3
4  1
3  5
1  2
2  6
1  5
6  3
1  6
3  2
6
1  2
1  3
2  4
2  5
3  6
3  7
0
输出样例：
Case 1: 2 4
Case 2: 4 1

#include<bits/stdc++.h>
using namespace std;
const int N = 1e3 + 10;
const int M = 1e3 + 10;
typedef unsigned long long ULL;
struct Edge{
    int v,next;
}edge[M];
int head[N],cnt;
void add(int u,int v){
    edge[cnt].v = v;
    edge[cnt].next = head[u];
    head[u] = cnt ++;
}
int in[N];
int sccno;
vector<int>g[N];
int low[N],num[N],dfn;
int sta[N],top;
int is_solate[N],cut[N];
//点双连通分量
void Tarjan(int u,int fa){
    num[u] = low[u] =  ++dfn;
    if(fa == -1 && head[u] == -1){      //孤立节点特判
        sccno ++;
        g[sccno].push_back(u);
        is_solate[sccno] = true;
        return;
    }
    sta[top ++] = u;
    int child = 0;
    for(int i = head[u];~i;i = edge[i].next){
        int v = edge[i].v;
        if(v == fa)continue;
        if(!num[v]){
            child ++;
            Tarjan(v,u);
            low[u] = min(low[u],low[v]);
            if(low[v] >= num[u]){
                if(fa != -1)cut[u] = true;
                int t;
                sccno ++;
                while(true){
                    t = sta[-- top];
                    g[sccno].push_back(t);
                    if(t == v)break;
                }
                g[sccno].push_back(u);
            }
        }
        else low[u] = min(low[u],num[v]);
    }
    if(fa == -1 && child >= 2)cut[u] = true;
}
int main(){
    int m = 0;
    int x,y;
    int n;
    int T = 0;
    while(cin>>m,m){
        n = 0;
        memset(head,-1,sizeof head);
        memset(num,0,sizeof num);
        memset(low,0,sizeof low);
        memset(cut,0,sizeof cut);
        dfn = 0,cnt = 0;
        for(int i = 1;i <= sccno;i ++)g[i].clear();
        sccno = 0;
        for(int i = 0;i < m;i ++){
            cin>>x>>y;
            add(x,y),add(y,x);
            n = max(n,max(x,y));
        }
        for(int i = 1;i <= n;i ++){
            if(!num[i])Tarjan(i,-1);
        }
        ULL res = 0;
        ULL c = 1;
        for(int i = 1;i <= sccno;i ++){
            int d = 0;
            for(auto &a : g[i]){
                if(cut[a])d ++;
            }
            int a = g[i].size() - d;
            if(d == 0){
                if(a == 1)res += 1;
                else res += 2,c *= (a * (a - 1) / 2);
            }
            else if(d == 1){
                if(a == 1)res += 1;
                else res += 1,c *= a;
            }
        }
        printf("Case %lld: %lld %lld\n",++T,res,c);
    }
    return 0;
}