bzoj2730 [HNOI2012]矿场搭建(点双)

先tarjan求一下所有的点双。
然后对于每一个点双,分类讨论:
1、只有一个割点,必须选一个非割点。
2、有>=2个割点,不用选
3、有0个割点,必须选俩。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 510
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,h[N],num=0,dfn[N],low[N],dfnum=0,scc=0,tst=0;
bool cut[N];
struct edge{
    int to,next;
}data[N];
vector<int>a[N];
stack<int>qq;
inline void tarjan(int x,int Fa){
    dfn[x]=low[x]=++dfnum;qq.push(x);int sz=0;
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(y==Fa) continue;
        if(!dfn[y]){
            tarjan(y,x);low[x]=min(low[x],low[y]);++sz;
            if(!Fa){
                ++scc;while(1){
                    int z=qq.top();qq.pop();
                    a[scc].push_back(z);if(z==y) break;
                }a[scc].push_back(x);continue;
            }if(low[y]>=dfn[x]){
                ++scc;while(1){
                    int z=qq.top();qq.pop();
                    a[scc].push_back(z);if(z==y) break;
                }a[scc].push_back(x);cut[x]=1;
            }
        }else low[x]=min(low[x],dfn[y]);
    }if(!Fa&&sz>1) cut[x]=1;if(!Fa) qq.pop();
}
int main(){
//  freopen("a.in","r",stdin);
    while(1){
        m=read();if(!m) break;n=0;++tst;
        memset(h,0,sizeof(h));memset(cut,0,sizeof(cut));
        num=0;dfnum=0;scc=0;memset(dfn,0,sizeof(dfn));
        while(m--){
            int x=read(),y=read();n=max(n,max(x,y));
            data[++num].to=y;data[num].next=h[x];h[x]=num;
            data[++num].to=x;data[num].next=h[y];h[y]=num;
        }for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i,0);
        ll ans=1;int cnt=0;
        for(int i=1;i<=scc;++i){
            int tmp=0;
            for(int j=0;j<a[i].size();++j) tmp+=cut[a[i][j]];
            if(!tmp){ans*=(a[i].size()-1)*a[i].size()/2;cnt+=2;continue;}
            if(tmp>1) continue;++cnt;ans*=(a[i].size()-1);
        }
        printf("Case %d: %d %lld\n",tst,cnt,ans);
        for(int i=1;i<=scc;++i) a[i].clear();
    }return 0;
}

猜你喜欢

转载自blog.csdn.net/Icefox_zhx/article/details/80526361
今日推荐