hnoi矿场搭建——Tarjan割点

Description

煤矿工地可以看成是由隧道连接挖煤点组成的无向图。为安全起见，希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。于是矿主决定在某些挖煤点设立救援出口，使得无论哪一个挖煤点坍塌之后，其他挖煤点的工人都有一条道路通向救援出口。请写一个程序，用来计算至少需要设置几个救援出口，以及不同最少救援出口的设置方案总数。
Input

输入文件有若干组数据，每组数据的第一行是一个正整数 N（N≤500），表示工地的隧道数，接下来的 N 行每行是用空格隔开的两个整数 S 和 T，表示挖 S 与挖煤点 T 由隧道直接连接。输入数据以 0 结尾。

Output

输入文件中有多少组数据，输出文件 output.txt 中就有多少行。每行对应一组输入数据的 结果。其中第 i 行以 Case i: 开始（注意大小写，Case 与 i 之间有空格，i 与:之间无空格，: 之后有空格），其后是用空格隔开的两个正整数，第一个正整数表示对于第 i 组输入数据至少需 要设置几个救援出口，第二个正整数表示对于第 i 组输入数据不同最少救援出口的设置方案总 数。输入数据保证答案小于 2^64。输出格式参照以下输入输出样例。

Sample Input

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

Sample Output

Case 1: 2 4

Case 2: 4 1

HINT

Case 1 的四组解分别是(2,4),(3,4),(4,5),(4,6)；

Case 2 的一组解为(4,5,6,7)。

---

我们先建出第一个样例的图，如下图：

我们发现这个图中有且只有一个割点，即点1。如果我们删去这个割点以及和其有关的边，那么该图剩下两个联通块，即（4）和（2,3,5,6），我们在这两个连通块里各建造一个营救点即可。

我们再建造样例二的图例，发现和图一的原理也一样（即删去和割点有关的边，剩下每个连通块建一个营救点）。

所以我们猜想：求出割点并删去，然后求连通块的个数即可。

但这个猜想是不是正确的呢？其实这个猜想存在一个漏洞：如果一个连通块与两个或以上个割点相连通时，就不需要建造营救点了。因为即使与它相连的某个割点坍塌，它还可以从别的与它相连的割点跑走。

所以我们求出割点并删去，然后找出每个连通块，统计与它相连的割点数，按照不同的点数进行讨论即可。

那么割点怎么求呢？我找了个blog，里面写的很详细，我在这里就概述一下：我们用类似Tarjan的方法求取，如果对于某条边上的两点u,v,存在low[u]>=dfn[v],则v为割点

Code：

#include<bits/stdc++.h>
#include<vector>
#define ll long long
using namespace std;
ll read(){
    char c;ll x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
vector<ll> f[1005];
struct node{
    ll sum,tot;
    ll re[1005];
}col[1005];
ll n,p,co,ta,low[1005],dfn[1005],ans,k,result;
ll vis[1005],fa[1005],cut[1005];
ll most(ll x,ll y,ll z){ll ans=0;ans=max(x,y);ans=max(ans,z);return ans;}
void tarjan(ll x){
    dfn[x]=low[x]=ta++;vis[x]=1;
    ll res=0;
    for(ll i=0;i<f[x].size();i++){
        ll to=f[x][i];
        if(!vis[to]){
            res++;fa[to]=x;tarjan(to);
            low[x]=min(low[x],low[to]);
            if(!fa[x]&&res>=2) cut[x]=1;
            else if(fa[x]&&low[to]>=dfn[x]) cut[x]=1;
        }
        else if(to!=fa[x]){
            low[x]=min(low[x],dfn[to]);
        }
    }
}
void dfs(ll x){
    col[co].sum++;vis[x]=1;
    for(ll i=0;i<f[x].size();i++){
        ll to=f[x][i];
        if(cut[to]){
            if(!col[co].re[to]){
                col[co].re[to]=1;
                col[co].tot++;
            }
        }
        else if(!vis[to]){
            dfs(to);
        }
    }
}
int main()
{
    while(1){
        k++;
        n=read();ans=1;result=0;
        if(!n) break;
        for(ll i=1;i<=co;i++){
            col[i].sum=col[i].tot=0;
            memset(col[i].re,0,sizeof(col[i].re));
        }
        memset(vis,0,sizeof(vis));memset(fa,0,sizeof(fa));memset(cut,0,sizeof(cut));
        for(ll i=1;i<=p;i++) f[i].clear();
        p=co=ta=0;
        for(ll i=1;i<=n;i++){
            ll x=read(),y=read();
            f[x].push_back(y);
            f[y].push_back(x);
            p=most(p,x,y);
        }
        for(ll i=1;i<=p;i++)
            if(!dfn[i]) tarjan(i);
        memset(vis,0,sizeof(vis));
        for(ll i=1;i<=p;i++)
            if(!vis[i]&&!cut[i]) co++,dfs(i);
        for(ll i=1;i<=co;i++){
            if(col[i].tot>=2) continue;
            if(col[i].tot==1) ans*=col[i].sum,result+=1;;
            if(!col[i].tot) ans*=(col[i].sum*(col[i].sum-1))/2,result+=2;
        }
        printf("Case %d: ",k);
        printf("%lld %lld\n",result,ans);
    }
    return 0;
}