HDU -- 3844 Mining Your Own Business（点双连通）

题目大意： 一座地下稀有金属矿由n条隧道和一些连接点组成，其中每条隧道连接两个点，任意两个连接点之间只有一条隧道，为降低矿工的危险，在一些连接点出安装逃生装置，使得不管哪个连接点倒塌，不在此连接点的所有矿工都能到达逃生装置，为节约成本，选择经量少的连接点安装逃生装置，计算出逃生装置的数目和最小时的安装方案总数。

思路分析：

①：一个双连通分量有多个割点时，不需要安装逃生装置；

①：双连通分量只有一个割点时，把逃生装置安装在割点是不明智的，一旦割点发生倒塌，图就会变得支离破碎，所以安装一个逃生装置，安装的方案为此双连通分量点的个数减一；

③：不同的双连通之间的方案数是相乘的关系；

④：当图中只有一个双连通分量时，选择两个点来安装逃生装置，方案数为V*（V-1）/ 2，其中V是顶点个数。

代码实现：

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<stack>
#include<cmath>
#include<algorithm>
#include<iostream>
#pragma comment(linker, "/STACK:102400000,102400000")
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define MEM(a) (memset((a),0,sizeof(a)))
#define MEME(a) (memset((a),-1,sizeof(a)))
#define MEMX(a) (memset((a),0x3f,sizeof(a)))
using namespace std;
const int N=100005;
int low[N],dfn[N],tim,bccn,top,n,m,cnt[N];
bool iscut[N];
stack<int> st;
vector<int> bcc[N];
struct Edge{
    int v;
    Edge *next;
}*head[N],e[N];
void Addedge(int from,int to){
    Edge *p=&e[top++];
    p->v=to;
    p->next=head[from];
    head[from]=p;
}
void Tarjan_bcc(int u,int fa){
    dfn[u]=low[u]=++tim;
    st.push(u);
    int v,tmp;
    for(Edge *p=head[u];p;p=p->next){
        v=p->v;
        if(v==fa) continue;
        if(!dfn[v]){
            Tarjan_bcc(v,u);
            low[u]=Min(low[u],low[v]);
            if(low[v]>=dfn[u]){
                ++bccn;
                do{
                    tmp=st.top();
                    st.pop();
                    bcc[bccn].push_back(tmp);
                    cnt[tmp]++;
                }while(tmp!=v);
                bcc[bccn].push_back(u);//一个割点属于多个块
                cnt[u]++;
            }
        }else low[u]=Min(low[u],dfn[v]);
    }
}
void Solve(){
    long long ans1=0,ans2=1;
    for(int i=1;i<=bccn;++i){
        int cnt=0;
        if(bcc[i].size()==1&&iscut[bcc[i][0]]) continue;
        for(int j=0;j<bcc[i].size();++j)
            if(iscut[bcc[i][j]]) cnt++;
        if(cnt==1) ans1++,ans2*=(long long)(bcc[i].size()-1);
    }
    if(bccn==1) ans1=2,ans2=(long long)(bcc[1].size()*(bcc[1].size()-1)/2);
    printf("%I64d %I64d\n",ans1,ans2);
}
int main(){
    int Case=0;
    while(~scanf("%d",&m),m){
        MEM(head),MEM(dfn),MEM(low),MEM(iscut),MEM(cnt);
        bccn=top=tim=n=0;
        int u,v;
        while(m--){
            scanf("%d%d",&u,&v);
            n=Max(n,u),n=Max(n,v);
            Addedge(u,v);
            Addedge(v,u);
        }
        for(int i=1;i<=n;++i) bcc[i].clear();
        for(int i=1;i<=n;++i) if(!dfn[i]) Tarjan_bcc(i,-1);
        for(int i=1;i<=n;++i) if(cnt[i]>1) iscut[i]=1;
        //cout<<"bccn="<<bccn<<endl;
        printf("Case %d: ",++Case);
        Solve();
    }
}