例题5.6 井下矿工 LA5135

1.题目描述:点击打开链接

2.解题思路:本题利用BCC的性质解决。本题的图论模型是:在一个无向图上选择尽量少的点涂黑,使得任意删除一个点后,每个连通分量至少有个一黑点。通过尝试可以发现如下的规律:

(1)把割顶涂黑是不划算的,而且在一个点-双连通分量中涂黑2个黑点也是不划算的。

(2)当一个连通分量中只有一个割顶时候才需要涂,而且任选一个非割顶涂黑即可。

(3)当整个图没有割顶时候,此时需要涂黑2个点,方案总数是V(V-1)/2。

本题输入时候有一个细节:需要给每个顶点分配一个编号。这样才能保证所有顶点都是连续的。而且是从0开始的。由于整个图是连通的,因此只需要调用dfs(0,-1)即可。

3.代码:

#include<iostream>
#include<algorithm>
#include<cassert>
#include<string>
#include<sstream>
#include<set>
#include<bitset>
#include<vector>
#include<stack>
#include<map>
#include<queue>
#include<deque>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<cctype>
#include<list>
#include<complex>
#include<functional>
using namespace std;

#define me(s)  memset(s,0,sizeof(s))
#define rep(i,n) for(int i=0;i<(n);i++)
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair <int,int> P;


const int N=100000+10;
vector<int>g[N];
int pre[N],iscut[N],bccno[N];
vector<int>bcc[N];
struct Edge{int u,v;};
stack<Edge>s;
int dfs_clock,bcc_cnt;

int dfs(int u,int fa)
{
    int lowu=pre[u]=++dfs_clock;
    int child=0;
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        Edge e=Edge{u,v};
        if(!pre[v])
        {
            s.push(e);
            child++;
            int lowv=dfs(v,u);
            lowu=min(lowu,lowv);
            if(lowv>=pre[u])
            {
                iscut[u]=1;
                bcc[++bcc_cnt].clear();
                for(;;)
                {
                    Edge x=s.top();s.pop();
                    if(bccno[x.u]!=bcc_cnt){bcc[bcc_cnt].push_back(x.u);bccno[x.u]=bcc_cnt;}
                    if(bccno[x.v]!=bcc_cnt){bcc[bcc_cnt].push_back(x.v);bccno[x.v]=bcc_cnt;}
                    if(x.u==u&&x.v==v)break;
                }
            }
        }
        else if(pre[v]<pre[u]&&v!=fa)
        {
            s.push(e);
            lowu=min(lowu,pre[v]);
        }
    }
    if(fa<0&&child==1)iscut[u]=0;
    return lowu;
}

void find_bcc()
{
    dfs_clock=bcc_cnt=0;
    memset(pre,0,sizeof(pre));
    memset(iscut,0,sizeof(iscut));
    memset(bccno,0,sizeof(bccno));
    dfs(0,-1);
}

struct ID  //给每个顶点按顺序分配编号
{
    map<int,int>m;
    int cnt;
    ID():cnt(0){}
    int get(int x)
    {
        if(!m.count(x))m[x]=cnt++;
        return m[x];
    }
};

int main()
{
    int n;
    int kase=0;
    while(~scanf("%d",&n)&&n)
    {
        int u,v;
        me(g);
        ID id;
        for(int i=0;i<n;i++)
        {
            scanf("%d%d",&u,&v);
            u=id.get(u);v=id.get(v);
            g[u].push_back(v);
            g[v].push_back(u);
        }
        find_bcc();
        ll ans1=0,ans2=1;
        for(int i=1;i<=bcc_cnt;i++)
        {
            int cut_cnt=0;
            for(int j=0;j<bcc[i].size();j++)
                if(iscut[bcc[i][j]])cut_cnt++;
            if(cut_cnt==1)
            {
                ans1++;
                ans2*=(ll)(bcc[i].size()-cut_cnt);
            }
        }
        if(bcc_cnt==1) //注意特判!
            ans1=2,ans2=bcc[1].size()*(bcc[1].size()-1)/2;
        printf("Case %d: %lld %lld\n",++kase,ans1,ans2);
    }
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值