bzoj2730/UVALIVE5135 矿场搭建 求割点

前（sai）记（hua）

第一次见到这道题是在某大佬的博客里（那时它还叫UVALIVE5135），以为它是点双连通的裸题（然而我最后没有用点双），没调出来，未果。
在班上大佬讲课时又遇见它【然后它变成了bzoj的题？】了，调出来了

解法

在一个连通图有多个连通分量时

该连通图矿井个数为割点为1的连通图的个数。
方案数为割点为1的连通图的点数-1（割点）相乘

在一个连通图只有一个连通分量时

该连通图的矿井个数为2
方案数为c(2,总点数）

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=5e4+5,M=5e4+5;
bool iscut[N],vis[N];
int low[N],dfn[N],idc;
int head[N],nxt[M],to[M],etot;
int u[N],v[N],c[N],cut[N],tot; 
void adde(int u,int v)
{
    to[++etot]=v;
    nxt[etot]=head[u];
    head[u]=etot;
}
void tarjan(int u,int fa)
{
    int son=0;
    low[u]=dfn[u]=++idc;
    for(int i=head[u];i;i=nxt[i]){
        int v=to[i];
        if(v==fa) continue;
        if(!dfn[v]){
            son++;
            tarjan(v,u);
            low[u]=min(low[v],low[u]);
            if(low[v]>=dfn[u]) iscut[u]=1;
        }   
        else low[u]=min(low[u],dfn[v]); 
    }
    if(u==fa&&son<=1) iscut[u]=0;  
}
int cnt;
int dfs(int u,int sta)
{
    int ans=0;
    vis[u]=1;cnt++;
    for(int i=head[u];i;i=nxt[i]){
        int v=to[i];
        if(vis[v]) continue;
        if(iscut[v]){
            if(cut[v]!=sta) {cut[v]=sta;ans++;}
            continue;
        }
        ans+=dfs(v,sta);
    }
    return ans;
}
void init()
{
    etot=0;idc=0;tot=0;
    memset(vis,0,sizeof(vis));
    memset(head,0,sizeof(head));
    memset(low,0,sizeof(low));
    memset(dfn,0,sizeof(dfn));
    memset(iscut,0,sizeof(iscut));
    memset(cut,0,sizeof(cut));
}
int main()
{
    int n,k=0;
    while(scanf("%d",&n)&&n){
        k++;
        init();
        for(int i=1;i<=n;i++){
            scanf("%d%d",&u[i],&v[i]);
            c[++tot]=u[i],c[++tot]=v[i];
        }
        sort(c+1,c+1+n*2);
        int sum=unique(c+1,c+1+n*2)-c-1;
        for(int i=1;i<=n;i++){
            u[i]=lower_bound(c+1,c+1+sum,u[i])-c;
            v[i]=lower_bound(c+1,c+1+sum,v[i])-c;
            adde(u[i],v[i]);
            adde(v[i],u[i]);
        }
        tarjan(1,1);
        int exit=0;long long ans=1;
        for(int i=1;i<=sum;i++)
        if(!vis[i]&&!iscut[i]){
            cnt=0;
            int x=dfs(i,i);
            if(x==1) {
                exit++;
                ans=(long long)ans*cnt;
            }
            if(!x){
                exit+=2;
                ans=(long long)ans*cnt*(cnt-1)/2; 
            }  
        }
        printf("Case %d: %d %lld\n",k,exit,ans);    
        //bzoj要改成lld？ 
    }
    return 0;
}