【连通分量】BZOJ 2730: [HNOI2012]矿场搭建

2730: [HNOI2012]矿场搭建

---

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

---

Solution

先求割点
然后推一下可得如果一个双连通分量连1个割点那这里要设一个出口
2个以上不用连
如果只有一个双联通那么之间输出2 n*(n-1)/2
求割点由于数据太小可以n^2
还有一种更神的O（n）待学习..
然后对于一个点与几个割点相连可以枚举割点染色,染2个以上为-1
注意输出..冒号后面有空格

Code

#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>using namespace std;inline int in(){    int x=0;char ch=getchar();    while(ch<'0'||ch>'9')ch=getchar();    while(ch<='9'&&ch>='0')x=x*10+ch-'0',ch=getchar();    return x;}#define maxn 501int n,cnt,col[maxn],tot,last[maxn],num[maxn];struct edge{    int to,lst;}e[maxn*2];void add(int u,int v){    e[++tot]=(edge){v,last[u]};last[u]=tot;    e[++tot]=(edge){u,last[v]};last[v]=tot;}bool vis[maxn],cut[maxn];void dfs(int poi,int del){    vis[poi]=1;    for(int i=last[poi];i;i=e[i].lst)    {        int u=e[i].to;        if(vis[u]||u==del)continue;        dfs(e[i].to,del);    }}void dfs2(int poi){    vis[poi]=1;    if(!col[poi])col[poi]=cnt;    else col[poi]=-1;    for(int i=last[poi];i;i=e[i].lst)        if(!cut[e[i].to]&&!vis[e[i].to])            dfs2(e[i].to);}int main(){    int m,u,v,T=0;    while(1)    {        n=0;        memset(cut,0,sizeof(cut));        memset(last,0,sizeof(last));        memset(col,0,sizeof(col));        memset(num,0,sizeof(num));        tot=0;cnt=0;        m=in();        if(!m)return 0;        for(int i=1;i<=m;i++)        {            u=in(),v=in();            add(u,v);            n=max(u,n);            n=max(v,n);        }        for(int i=1;i<=n;i++)        {            memset(vis,0,sizeof(vis));            dfs(i%n+1,i);            for(int j=1;j<=n;j++)                if(!vis[j]&&j!=i)                {                    cut[i]=1;                    break;                }        }        for(int i=1;i<=n;i++)            if(cut[i])            {                memset(vis,0,sizeof(vis));                for(int j=last[i];j;j=e[j].lst)                    if(!vis[e[j].to]&&!cut[e[j].to])                        cnt++,dfs2(e[j].to);            }        int sum=0;        T++;printf("Case %d: ",T);        for(int i=1;i<=n;i++)            if(col[i])            {                if(col[i]==-1)continue;                num[col[i]]++;            }        long long ans=1;        for(int i=1;i<=cnt;i++)            if(num[i])sum++,ans*=num[i];        if(!sum)printf("2 %d\n",n*(n-1)/2);        else printf("%d %lld\n",sum,ans);    }}