HDU3844.Mining Your Own Business井下矿工——双连通分量

http://acm.hdu.edu.cn/showproblem.php?pid=3844

题目描述：
有一座地下的稀有金属矿由n条隧道和一些连接点组成，其中每条隧道连接两个连接点，任意两个连接点之间最多只有一条隧道。为了降低矿工的危险，你的任务是在一些连接点处安装太平井和相应的逃生装置，使得不管哪个连接点倒塌，不在此连接点的所有矿工都能到达太平井逃生（假定除了倒塌的连接点不能通行外，其他隧道和连接点完好无损）。为了节约成本，比应当在尽量少的连接点安装太平井。还需要计算出当太平井的数目最小时的总方案数。

分析：
本题的模型：在一个无向图上选择尽量少的点涂黑，使得任意删除一个点后，每个连通分量至少有一个黑点。不难发现，把割点涂黑是不划算的，而且在同一个点双连通分量中涂两个黑点也是不划算的。进一步分析发现，当一个点双连通分量只有一个割点时才需要涂，而且是任选一个非割点涂黑

一个特殊情况是整个图没有割点，那么此时只需涂两个点，方案数为n*(n-1)/2 n个数中取两个

这个无限WA了两天，终于发现n*(n-1)溢出了
原来这样写的ans=(ll)(n*(n-1)/2);
应该这样ans=(ll)n*(n-1)/2;

要手动扩栈

//530MS 10132K  2468 B
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<algorithm>
typedef long long ll;
const int MAXN=100100;
const int MAXM=100010;
using namespace std;
int low[MAXN],dfn[MAXN],Stack[MAXN];
int dfs_clock,top,block;
bool iscut[MAXN];
vector<int> bcc[MAXN];
struct Edge {
    int to,next;
} edge[MAXM<<1];int head[MAXN],tot;
void addedge(int u,int v){
    edge[tot].to=v;edge[tot].next=head[u];head[u]=tot++;
}
void init(){
    tot=0;
    memset(head,0xff,sizeof(head));
}
void Tarjan(int u,int pre){
    int v;
    low[u]=dfn[u]=++dfs_clock;
    Stack[top++]=u;
    int son=0;
    for(int i=head[u]; i!=-1; i=edge[i].next) {
        v=edge[i].to;
        if(v==pre) continue;
        if(!dfn[v]) {
            son++;
            Tarjan(v,u);
            low[u]=min(low[u],low[v]);
            if(low[v]>=dfn[u]) {
                iscut[u]=1;
                block++;
                bcc[block].clear();
                int vn;
                for(;;) {
                    vn=Stack[--top];
                    bcc[block].push_back(vn);
                    if(vn==v) break;
                }
                bcc[block].push_back(u);
            }
        } else if(low[u]>dfn[v]) {
            low[u]=dfn[v];
        }
    }
    if(pre<0&&son==1) iscut[u]=false;
}
void solve(int n)
{
    memset(dfn,0,sizeof(dfn));
    memset(iscut,false,sizeof(iscut));
    dfs_clock=top=block=0;
    for(int i=1;i<=n;++i)
        if(!dfn[i]) Tarjan(i,-1);
    ll ans1=0,ans2=1;
    if(block==1) {
        ans1=2ll;
        ans2=(ll)n*(n-1)/2;
    } else {
        for(int i=1; i<=block; ++i) {
            int z=bcc[i].size();
            int cut_cnt=0;
            for(int j=0; j<z; ++j) {
                if(iscut[bcc[i][j]]) cut_cnt++;
            }
            if(cut_cnt==1) {
                ans1++;
                ans2*=(ll)(z-1);
            }
        }
    }
    printf("%lld %lld\n",ans1,ans2);
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.cpp","r",stdin);
#endif // ONLINE_JUDGE
    int m,u,v,cas=1;
    while(scanf("%d",&m)&&m) {
        int n=0;
        init();
        while(m--) {
            scanf("%d%d",&u,&v);
            addedge(u,v);
            addedge(v,u);
            n=max(n,max(u,v));
        }
        printf("Case %d: ",cas++);
        solve(n);
    }
    return 0;
}