洛谷P3225 [HNOI2012]矿场搭建（割点+点双连通分量+组合数学）（点双的两种求法模板）

链接：https://www.luogu.org/problem/P3225

题意：多组样例，有若干个挖煤点，由n条隧道相连。现在要在某些挖煤点修救援出口，使得某个挖煤点倒塌后，其它挖煤点的人可以从救援出口出去。要求的是最少的救援出口数和在保证最少救援出口数的条件下，不同的方案数。（注意，挖煤点的最大编号就是挖煤点的数量，千万不要理解成有n个挖煤点。）

思路：很显然这个题和割点密切相关。先介绍一些前置知识。割点是什么？把该点割掉后，图不连通。点双是什么？点双内任意两点之间都有两条经过不同点的路径，也就是点双内没有割点。而割点肯定是某些点双的公共点，但某些点双中可能没有割点。现在考虑题目要求的最少救援出口，我们肯定要先求出所有点双（因为点双之内的点能经过不同的点互达，这和题目要求的一样。）首先，在一个点双内肯定要建两个救援出口（炸了一个还有另一个嘛）；然后，因为割点是某些点双的公共点，也就是建在割点的救援出口，可以对若个点双都贡献一个救援出口，题目要求救援出口竟可能的少，那肯定有割点的话就要建在割点上，现在按割点的个数分类讨论。

1、如果该点双内没有割点，那么它就是一个独立的点双和非该点双的点没有联系，那么我们要建两个救援出口，救援出口数+2，方案数就是从点双内任选两个点的方案数，也就是。（假设该点双内的点数为num。）

2、如果该点双内只有1个割点，那么一个出口肯定建在割点上，另一个出口从该点双内不是割点的点中任选一个，救援出口数+1，方案数也就是。

3、如果该点双内有超过1个割点，该点双的贡献已被计算，不用再处理。

最后根据乘法原理求出总方案数即可。

PS：贴上两种求点双的方法留模板。

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 510;
struct node
{
    int to,nxt;
    bool use;
}g[N<<1];
int head[N],cnt;
int dfn[N],low[N],id;
bool vis[N],iscut[N];
int bcc[N],bccnum;
int n,m;
int num,ans1,cut;
ll ans2;
void Init()
{
    ans1=0;
    ans2=1;
    bccnum=cnt=id=0;
    memset(iscut,0,sizeof(iscut));
    memset(vis,0,sizeof(vis));
    memset(dfn,0,sizeof(dfn));
    memset(bcc,0,sizeof(bcc));
    memset(head,-1,sizeof(head));
}
void add(int u,int v)
{
    g[cnt].to=v; g[cnt].use=0; g[cnt].nxt=head[u]; head[u]=cnt++;
}
void tarjan(int u,int fa)
{
    int v,child=0;
    dfn[u]=low[u]=++id;
    for(int i=head[u];~i;i=g[i].nxt)
    {
        if(g[i].use) continue;
        g[i].use=g[i^1].use=1;
        v=g[i].to;
        if(!dfn[v])
        {
            child++;
            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(fa==0&&child<2)  iscut[u]=0;
}
void dfs(int u)
{
    int v;
    vis[u]=1;
    bcc[u]=bccnum;
    num++;
    for(int i=head[u];~i;i=g[i].nxt)
    {
        v=g[i].to;
        if(iscut[v]&&bcc[v]!=bccnum)
        {
            cut++,bcc[v]=bccnum;
        }
        if(!iscut[v]&&!vis[v]) dfs(v);
    }
}
int main()
{
    int tt=0,u,v;
    while(~scanf("%d",&m)&&m)
    {
        Init();
        n=0;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&u,&v);
            add(u,v),add(v,u);
            n=max(max(u,v),n);
        }
        for(int i=1;i<=n;i++)
            if(!dfn[i]) tarjan(i,0);
        for(int i=1;i<=n;i++)
        {
            if(vis[i]||iscut[i]) continue;
            num=cut=0;
            bccnum++;
            dfs(i);
            if(!cut) ans1+=2,ans2*=(1LL*num*(num-1)/2);
            if(cut==1) ans1++,ans2*=1LL*num;
        }
        printf("Case %d: %d %lld\n",++tt,ans1,ans2);
    }
    return 0;
}

 

#include <bits/stdc++.h>
#include <vector>
#define ll long long
using namespace std;
const int N = 510;
struct node
{
    int from,to,nxt;
    bool use;
}g[N<<1];
int head[N],cnt;
int dfn[N],low[N],id,sta[N<<1],top;
bool vis[N],iscut[N];
int bccnum;
vector <int> bcc[N];
int n,m;
int num,ans1,cut;
ll ans2;
void Init()
{
    ans1=0;
    ans2=1;
    bccnum=cnt=id=top=0;
    memset(iscut,0,sizeof(iscut));
    memset(vis,0,sizeof(vis));
    memset(dfn,0,sizeof(dfn));
    memset(head,-1,sizeof(head));
}
void add(int u,int v)
{
    g[cnt].from=u; g[cnt].to=v; g[cnt].use=0; g[cnt].nxt=head[u]; head[u]=cnt++;
}
void tarjan(int u,int fa)
{
    int v,child=0;
    dfn[u]=low[u]=++id;
    for(int i=head[u];~i;i=g[i].nxt)
    {
        if(g[i].use) continue;
        g[i].use=g[i^1].use=1;
        v=g[i].to;
        if(!dfn[v])
        {
            sta[++top]=v;
            child++;
            tarjan(v,u);
            low[u]=min(low[v],low[u]);
            if(low[v]>=dfn[u])
            {
                iscut[u]=1;
                bcc[++bccnum].clear();
                while(sta[top]!=v)
                    bcc[bccnum].push_back(sta[top--]);
                bcc[bccnum].push_back(sta[top--]);
                bcc[bccnum].push_back(u);
            }

        }
        else low[u]=min(low[u],dfn[v]);
    }
    if(fa==0&&child<2)  iscut[u]=0;
}

int main()
{
    int tt=0,u,v;
    while(~scanf("%d",&m)&&m)
    {
        Init();
        n=0;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&u,&v);
            add(u,v),add(v,u);
            n=max(max(u,v),n);
        }
        for(int i=1;i<=n;i++)
            if(!dfn[i]) tarjan(i,0);
        for(int i=1;i<=bccnum;i++)
        {
            cut=0;
            num=bcc[i].size();
            for(int j=0;j<num;j++)
                if(iscut[bcc[i][j]]) cut++;
            if(!cut) ans1+=2,ans2*=(1LL*num*(num-1)/2);
            if(cut==1) num--,ans1++,ans2*=1LL*num;
        }
        printf("Case %d: %d %lld\n",++tt,ans1,ans2);
    }
    return 0;
}