P3225 [HNOI2012]矿场搭建

题目描述

煤矿工地可以看成是由隧道连接挖煤点组成的无向图。为安全起见,希望在工地发生事故时所有挖煤点的工人都能有一条出路逃到救援出口处。于是矿主决定在某些挖煤点设立救援出口,使得无论哪一个挖煤点坍塌之后,其他挖煤点的工人都有一条道路通向救援出口。

请写一个程序,用来计算至少需要设置几个救援出口,以及不同最少救援出口的设置方案总数。

点双入门。。。 话说题意是真的晦涩。

一眼求割点,求完割点怎么办?

因为我们把割点去掉之后,图中会剩下一个一个的块,这些块就是点双。

那么限制工人逃跑的就是割点,所以我们在乎的是每个点双中有多少个割点。

那么我们就可以通过Tarjan给所有割点打上一个标记。

然后再去DFS出所有的点双,这个过程需要记录点双中的点和点双中的割点数。有些细节。

然后就可以愉快地分情况讨论了。

(一),当前点双内没有割点,那么这个点双是独立的也就是不和其他的块联通,那么对于这个点双内的点我们就需要建立两个救援出口,为什么呢,因为有一个炸了还有另外一个。所以\(ans1+=2\)

然后方案数通过乘法原理和组合数搞一搞就行了,也就是\(ans2*=C(tot,2)\)。

(二),当前点双内只有一个割点,那么可以知道如果这个割点炸了我们就翻车了,所以我们可以继承(一)的思想,只是把这个已有的割点看成一个必然的出口,所以我们还需要建1个出口,方案数为\(ans2*=C(tot,1)\)。

(三),当前点双内有两个割点,还是一样的思想,就会发现我们不用再建出口了,因为可以把两个割点看成两个出口,任意一个炸了我们还有另一个。

ok,那就可以开始愉快地写代码了喵喵喵。。。

code:

#include <iostream>
#include <cstdio>
#include <cstring> using namespace std; #define int long long const int wx=517; inline int read(){
int sum=0,f=1; char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
while(ch>='0'&&ch<='9'){sum=(sum<<1)+(sum<<3)+ch-'0'; ch=getchar();}
return sum*f;
} int n,m,ans1,ans2,num,sjc,tot,zmj,cnt;
int dfn[wx],head[wx],ok[wx],vis[wx],low[wx]; struct e{
int nxt,to;
}edge[wx*2]; void add(int from,int to){
edge[++num].nxt=head[from];
edge[num].to=to;
head[from]=num;
} void Tarjan(int u,int f){
dfn[u]=low[u]=++tot; int child=0;
for(int i=head[u];i;i=edge[i].nxt){
int v=edge[i].to;
if(v==f)continue;
if(!dfn[v]){
Tarjan(v,u);
child++;
low[u]=min(low[u],low[v]);
if(low[v]>=dfn[u])ok[u]=1;
}
else if(dfn[v]<dfn[u]&&v!=f){
low[u]=min(low[u],dfn[v]);
}
}
if(child<=1&&f==-1)ok[u]=0;
} void dfs(int u){
vis[u]=tot; zmj++;
for(int i=head[u];i;i=edge[i].nxt){
int v=edge[i].to;
if(vis[v]!=vis[u]&&ok[v])cnt++,vis[v]=tot;
if(!vis[v]&&!ok[v])dfs(v);
}
} void clear(){
num=0,tot=0;
memset(head,0,sizeof head);
memset(edge,0,sizeof edge);
memset(vis,0,sizeof vis);
memset(dfn,0,sizeof dfn);
memset(ok,0,sizeof ok);
} signed main(){
while(1){
m=read(); if(!m)break;
n=0; sjc++; ans1=0; ans2=1;
clear();
for(int i=1;i<=m;i++){
int x,y;
x=read(); y=read();
add(x,y); add(y,x); n=max(n,max(x,y));
}
for(int i=1;i<=n;i++)if(!dfn[i])Tarjan(i,-1);
for(int i=1;i<=n;i++){
if(ok[i]||vis[i])continue;
tot++; zmj=0,cnt=0; dfs(i);
if(cnt==0)ans1+=2,ans2*=(zmj-1)*zmj/2;
if(cnt==1)ans1++,ans2*=zmj;
}
printf("Case %lld: %lld %lld\n",sjc,ans1,ans2);
}
return 0;
}

最新文章

  1. 一个C#语法高亮插件
  2. lsm-tree
  3. JS比较好用的一些方法搜集
  4. ABAP后台JOB数量控制
  5. windows分屏
  6. C# sogou地图API应用总结
  7. hbase-site.xml 参数设置
  8. (转)C#中的泛型
  9. hdu 4762 Cut the Cake概率公式
  10. win7怎么安装消息队列 MSMQ
  11. boost::asio设置同步连接超时
  12. 在打包程序中自动安装SQL Server数据库 .
  13. Activity进程和线程之间的关系
  14. Ajax提交Form表单的一种方法
  15. GlusterFS最佳实践
  16. [NOIp 2013]货车运输
  17. 第十一章 图像之2D(2)
  18. 痞子衡嵌入式:常用的数据差错控制技术(2)- 奇偶校验(Parity Check)
  19. Debian/Ubuntu安装WPS (转)
  20. CentOS7 使用ntp设置系统时间,开机自动设置时间,

热门文章

  1. SqlServer——for xml path
  2. leetcode872
  3. Spring Cloud Eureka 4 (高可用服务注册中心)
  4. Redis搭建(二):主从复制
  5. 本博文将一步步带领你实现抽屉官网的各种功能:包括登陆、注册、发送邮箱验证码、登陆验证码、页面登陆验证、发布文章、上传图片、form验证、点赞、评论、文章分页处理以及基于tronado的后端和ajax的前端数据处理。
  6. 浅谈块元素绝对定位的margin属性
  7. java中链表的数据(对象)位置交换
  8. ubuntu 下正确安装android手机驱动
  9. python日期,时间函数
  10. flex 布局的深入研究