洛谷 P3225 [HNOI2012]矿场搭建

题目描述
煤矿工地可以看成是由隧道连接挖煤点组成的无向图。为安全起见，希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。于是矿主决定在某些挖煤点设立救援出口，使得无论哪一个挖煤点坍塌之后，其他挖煤点的工人都有一条道路通向救援出口。
请写一个程序，用来计算至少需要设置几个救援出口，以及不同最少救援出口的设置方案总数。
输入输出格式
输入格式：
输入文件有若干组数据，每组数据的第一行是一个正整数 N（N<=500），表示工地的隧道数，接下来的 N 行每行是用空格隔开的两个整数 S 和 T，表示挖 S 与挖煤点 T 由隧道直接连接。输入数据以 0 结尾。
输出格式：
输入文件中有多少组数据，输出文件 output.txt 中就有多少行。每行对应一组输入数据的 结果。其中第 i 行以 Case i: 开始（注意大小写，Case 与 i 之间有空格，i 与:之间无空格，: 之后有空格），其后是用空格隔开的两个正整数，第一个正整数表示对于第 i 组输入数据至少需 要设置几个救援出口，第二个正整数表示对于第 i 组输入数据不同最少救援出口的设置方案总 数。输入数据保证答案小于 2^64。输出格式参照以下输入输出样例。
输入输出样例
输入样例#1：
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
输出样例#1：
Case 1: 2 4
Case 2: 4 1
说明
Case 1 的四组解分别是(2,4),(3,4),(4,5),(4,6)；
Case 2 的一组解为(4,5,6,7)。

//首先看到割点就是Tarjan搞 但是怎么搞
//首先假设我们把所有的点双都缩点 那么我们一定可以得到一棵树 然后我们就会发现
//叶子节点(只含有一个割点的点双)必须建 因为叶子节点如果不建 一旦割点被爆就死翘了
//非叶节点(含有两个或两个以上的割点的点双)不用建 因为即使一个割点被爆了也可以沿着另一个割点走到一个叶节点
//还有一种情况就是整个联通块都是点双(即不含割点的点双) 这样我们讨论点双的大小
//如果只有一个点 那么这个点必须建 数据没有卡这个的点所以我没写(其实是我忘写了 然后还过了)
//如果有两个或两个以上的点 那么要建两个 一个被爆了还可以走另一个
//方案数就是乘法原理的问题了 注意叶节点那里出口不能建在割点上
//注：分析来自PoPoQQQ世界大佬
#include<cstdio>
#include<cstring>
#include<iostream>
#include<stack>
#include<algorithm>
using namespace std;
#define MAXN 1010
#define LL long long
struct Edge{ int to,next; }e[MAXN<<4];
int head[MAXN],tot,n,m,cases,ans,visx,root; LL sum;
int dfn[MAXN],low[MAXN],belong[MAXN],cur;

stack<int> st;
bool exist[MAXN],cut[MAXN];
inline void prepare(){
    memset(head,0,sizeof head );
    memset(dfn,0,sizeof dfn );
    memset(low,0,sizeof low );
    tot=0,n=0,ans=0,sum=1;//sum 求 方案数  乘法原理
    memset(belong,0,sizeof belong );
}

inline void Add_Edge(int u,int v){ e[++tot].to=v,e[tot].next=head[u],head[u]=tot; }

void Tarjan(int u){
    dfn[u]=low[u]= ++visx;
    for(int i=head[u];i;i=e[i].next){
        int v=e[i].to;
        if(!dfn[v]){
            Tarjan(v);
            low[u]=min(low[u],low[v]);
            if(low[v]>=dfn[u]) belong[u]++;
        }
        else low[u]=min(low[u],dfn[v]);
    }
}

void Tarjan_(int u){
    dfn[u]=low[u]= ++visx; st.push(u);
    for(int i=head[u];i;i=e[i].next){
        int v=e[i].to;
        if(dfn[v]) low[u]=min(low[u],dfn[v]);
        else {
            Tarjan_(v);
            low[u]=min(low[u],low[v]);
            if(low[v]>=dfn[u]){
                int t,temp=0,size=0;
                do{
                    t=st.top();st.pop();
                    if(belong[t]>=2) ++temp;
                    ++size;
                }while(t!=v);
                t=u;
                if(belong[t]>=2) ++temp;
                ++size;
                if(!temp) ans+=2,sum*=size*(size-1)/2;
                else if(temp==1) ans++,sum*=size-1;
            }
        }
    }
}

int main(){
    while(scanf("%d",&m)==1){
        if(m==0) return 0;
        prepare();
        for(int u,v,i=1;i<=m;++i){
            scanf("%d%d",&u,&v);
            n=max(n,max(u,v));
            Add_Edge(u,v),Add_Edge(v,u);
        }
        for(int i=1;i<=n;++i)
            if(!dfn[i]) Tarjan(i);
            else ++belong[i];
        memset(dfn,0,sizeof dfn );
        for(int i=1;i<=n;++i)
            if(!dfn[i]) Tarjan_(i);
        printf("Case %d: ",++cases),printf("%d %lld\n",ans,sum);
    }
    return 0;
}