bzoj2730 [HNOI2012]矿场搭建 （UVAlive5135 Mining Your Own Business）

2730: [HNOI2012]矿场搭建

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1147  Solved: 528
[Submit][Status][Discuss]

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)。

【思路】

点-双连通分量。

求出bcc，对于每个bcc而言，最优的方案就是在每一个只有一个割点的bcc中安置一个逃生装置。

【代码】

 #include<cstdio>
 #include<cstring>
 #include<stack>
 #include<vector>
 #include<iostream>
 #include<cstdlib>
 #define FOR(a,b,c) for(int a=(b);a<(c);a++)
 using namespace std;

 typedef long long LL;
 const int maxn = +;

 struct Edge{ int u,v;
 };

 int pre[maxn],iscut[maxn],bccno[maxn],dfs_clock,bcc_cnt;
 vector<int> G[maxn],bcc[maxn];

 stack<Edge> S;

 int dfs(int u,int fa) {
     int lowu=pre[u]=++dfs_clock;
     int ch=;
     for(int i=;i<G[u].size();i++) {
         int v=G[u][i];
         Edge e=(Edge) {u,v};
         if(!pre[v]) {
             S.push(e);
             ch++;
             int lowv=dfs(v,u);
             lowu=min(lowu,lowv);
             if(lowv>=pre[u]) {
                 iscut[u]=;
                 bcc_cnt++; bcc[bcc_cnt].clear();
                 for(;;) {
                     Edge x=S.top(); S.pop();
                     if(bccno[x.u]!=bcc_cnt) bcc[bcc_cnt].push_back(x.u),bccno[x.u]=bcc_cnt;
                     if(bccno[x.v]!=bcc_cnt) bcc[bcc_cnt].push_back(x.v),bccno[x.v]=bcc_cnt;
                     if(x.u==u && x.v==v) break;
                 }
             }
         }
         else if(pre[v]<pre[u] && v!=fa) {
             S.push(e);  lowu=min(lowu,pre[v]);
         }
     }
     if(fa< && ch==) iscut[u]=;
     return lowu;
 }
 void find_bcc(int n) {
     memset(pre,,sizeof(pre));
     memset(iscut,,sizeof(iscut));
     memset(bccno,,sizeof(bccno));
     dfs_clock=bcc_cnt=;
     for(int i=;i<n;i++)
         if(!pre[i]) dfs(i,-);
 }

 int n,m,kase;

 int main() {
     while(scanf("%d",&m) && m) {
         int u,v;
         n=;
         for(int i=;i<*m;i++) G[i].clear();
         for(int i=;i<m;i++) {
             scanf("%d%d",&u,&v);
             u--,v--;
             n=max(n,max(u,v));
             G[u].push_back(v);
             G[v].push_back(u);
         }
         while(!S.empty()) S.pop();

         find_bcc(n);

         LL ans1=,ans2=;
         for(int i=;i<=bcc_cnt;i++) {
             int cnt=;
             for(int j=;j<bcc[i].size();j++)
                 if(iscut[bcc[i][j]]) cnt++;
             if(cnt==) {
                 ans1++; ans2*=(long long)(bcc[i].size()-);
             }
         }
         if(bcc_cnt==) {
             ans1=; ans2=bcc[].size()*(bcc[].size()-)/;
         }
         printf("Case %d: %lld %lld\n",++kase,ans1,ans2);
     }
     return ;
 }