[BZOJ] 2730: [HNOI2012]矿场搭建

2730: [HNOI2012]矿场搭建

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2317  Solved: 1076
[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)。

 

Source

day1

Analysis

qwq调个水题调错了近十次，，，

那么我们根据数据可以建图，这个图我们可以得出这些性质：

1. 割点不是出口，如果割点是出口的话，塌掉好像，，，就塌掉了，不是最优解。正解应该是在点双连通分量里建一个出口。

2. 有两个割点的连通块就没必要建出口了：因为可以跑到两个割点相连接的其他连通块

3. 没有割点的连通块至少有两个出口：万一塌掉一个出口，，，，

第三点尤为重要。

那么解法就是：

首先求双连通分量顺便把割点数目计算出来

然后把每个连通块内建出口的方案乘起来

注意只有一个割点和无割点的情况才需要计算

还有一件事，n不是点的数量

还有一件事，初始化

还有一件事，需要输出 Case %d: 

Code

#include<cstdio>
#include<iostream>
#include<vector>
#include<stack>
#include<cstring>
#define maxn 600
using namespace std;

int dfn[maxn],TIM,size[maxn],bcc_cnt,iscut[maxn],bccno[maxn],n,buck[maxn],m;

struct edge{ int from; int u; int v; }e[maxn]; int tot,first[maxn];
void insert(int u,int v){ tot++; e[tot].from = first[u]; e[tot].u = u; e[tot].v = v; first[u] = tot; }

vector<int> bcc[maxn];
stack<edge> S;

int dfs(int u,int fa){
    int lowu = dfn[u] = ++TIM;
    int child = 0;
    for(int i = first[u];i;i = e[i].from){
        int v = e[i].v;
        edge e = (edge){0,u,v};
        if(!dfn[v]){
            S.push(e);
            child++;
            int lowv = dfs(v,u);
            lowu = min(lowu,lowv);
            if(lowv >= dfn[u]){
                iscut[u] = true;
                bcc_cnt++; bcc[bcc_cnt].clear();
                for(;;){
                    edge x = S.top();S.pop();
                    if(bccno[x.u] != bcc_cnt){
                        if(bccno[x.u]) buck[bccno[x.u]]++;
                        bcc[bcc_cnt].push_back(x.u); bccno[x.u] = bcc_cnt;
                    }
                    if(bccno[x.v] != bcc_cnt){
                        if(bccno[x.v]) buck[bccno[x.v]]++;
                        bcc[bcc_cnt].push_back(x.v); bccno[x.v] = bcc_cnt;
                    }
                    if(x.u == u && x.v == v) break;
                }
            }
        }else if(dfn[v] < dfn[u] && v != fa) {
            S.push(e); lowu = min(lowu,dfn[v]);
        }
    }if(fa < 0 && child == 1) iscut[u] = false;
    return lowu;
}

void Build(){
    memset(dfn,0,sizeof(dfn));
    memset(iscut,0,sizeof(iscut));
    memset(bccno,0,sizeof(bccno));
    memset(buck,0,sizeof(buck));
    TIM = bcc_cnt = 0;
    
    for(int i = 1;i <= n;i++)
        if(!dfn[i]) dfs(i,-1);
}

int main(){
    int T = 1;
    scanf("%d",&m);
    while(m){
        n = 0;
        for(int i = 1;i <= m;i++){
            int x,y; scanf("%d%d",&x,&y);
            n = max(n,max(x,y));
            insert(x,y); insert(y,x);
        }Build();
        for(int i = 1;i <= n;i++)
            if(iscut[i]) buck[bccno[i]]++;
            
        long long ans = 1;
        int tot = 0;
        for(int i = 1;i <= bcc_cnt;i++){
            if(buck[i] == 1) ans *= max(((int)bcc[i].size()-1),1),tot++;
            if(!buck[i]) ans *= bcc[i].size()*(bcc[i].size()-1)/2,tot += 2;
//            cout << "ans" << ans << endl;
        }
        printf("Case %d: %d %lld\n",T,tot,ans); T++;
//        cout << ans << endl;
//        cout << "buck: ";
//        for(int i = 1;i <= n;i++) printf("%d ",buck[i]); cout << endl;
//        cout << "bcc: ";
//        for(int i = 1;i <= bcc_cnt;i++) printf("%d ",bcc[i].size()); cout << endl;
        memset(first,0,sizeof(first));
        tot = 0;
        scanf("%d",&m);
    }return 0;
}