Usaco 5.4.5

原创 2012年03月26日 20:33:57

这道题要求删去最少的点,使得网络中给定两点不连通。

这个最少的点在图论中有具体的定义,有名称为点连通度。

那么,如何求这个点连通度呢?

我们可以考虑从前面的知识进行转化,前面有一道题目是让我们求边连通度。

那么从边连通度转化为点连通度是否可行呢?答案是可行的。

可以把在网络中的每个点p(这里p代表其中一个点),拆成两个点,p1,p2,在p1,p2之间连一条边,其容量为1.

那么对于原来就在网络中相连的两点(又因为该题为无向图),假设这两点为p,q,拆点后形成p1,p2,q1,q2,

我们可以在p2和q1之间连一条容量为INF(maxint),在q2和p1之间也连一条容量为INF的边。

做一遍最大流,就可以得出点连通度。


另外,题目还要求我们输出字典序最小的组成大小为点连通度的点割集。

这其实也可以借鉴前面那一道求边连通度的题目。

边连通度的题目,我们尝试按升序去删除每一条边,如果改变删除后,网络中的最大流减小,说名改边存在于边割集中。

同理,因为这道题目我们拆点了,并且拆点后的两个子点之间有边相连,所以我们可以尝试按升序去删除每一条拆点拆出来的边,

如果做最大流后流量减少了1(因为拆点出来的边容量为1),那么该点便属于点割集,一直做到最大流为0为止。


洋洋洒洒写了这么多,只是希望大家能懂。微笑

/*
ID: volz.kz.g
PROB: telecow
LANG: C++
*/
#include <iostream>
#include <fstream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#define INF 0x7FFFFFFF
using namespace std;
ifstream fin("telecow.in");
ofstream fout("telecow.out");

int n,m,s,t,ans;
struct edge_node{
    int v,c,next;
}e[6000];
int edge_num;
int d[6000],c[6000],v[6000];
int g[501][501];
bool cut[501];

void init(){
    memset(e,0,sizeof(e));
    memset(v,0,sizeof(v));
    edge_num=0;
}
void insert(int a,int b,int c){
    //fout << a << " " << b << " " << c << endl;
    ++edge_num;
    e[edge_num].v=b;e[edge_num].c=c;e[edge_num].next=v[a];v[a]=edge_num;
    ++edge_num;
    e[edge_num].v=a;e[edge_num].c=0;e[edge_num].next=v[b];v[b]=edge_num;
}
void makegraph(){
    init();
    for (int i=1;i<=n;++i)
	for (int j=1;j<=n;++j)
	    if (g[i][j]){
		insert(i+n,j,INF);
		insert(j+n,i,INF);
	    }
    for (int i=1;i<=n;++i)
	if (cut[i])
	    insert(i,i+n,1);
	else 
	    insert(i,i+n,0);
}
int opp(int x){
    return (x%2==1)?x+1:x-1;
}
int maxflow(int x,int lim){
    int ret,min_num,max_flow,i,j;
    if (x==t) return lim;
    min_num=2*n;max_flow=0;
    j=v[x];
    while (j!=0){
	if (e[j].c>0){
	    if (d[e[j].v]+1==d[x]){
		ret=min(lim,e[j].c);
		ret=maxflow(e[j].v,ret);
		e[j].c-=ret;e[opp(j)].c+=ret;
		lim-=ret;max_flow+=ret;
		if (lim==0 || d[s]>=2*n) return max_flow;
	    }
	    min_num=min(min_num,d[e[j].v]+1);
	}
	j=e[j].next;
    }
    if (min_num==d[x]) return max_flow;
    --c[d[x]];
    if (c[d[x]]==0) d[s]=2*n;
    d[x]=min_num;c[d[x]]++;
    return max_flow;
}
int main(){
    memset(g,false,sizeof(g));
    memset(cut,true,sizeof(cut));
    fin >> n >> m >> s >>t;
    for (int i=1;i<=m;++i){
	int a,b;fin >> a >> b;
	g[a][b]=true;g[b][a]=true;
    }
    makegraph();
    c[0]=2*n;
    while (d[s]<2*n) ans+=maxflow(s+n,INF);
    fout << ans << endl;//输出点连通度
    int print[201],print_num=0;
    for (int i=1;i<=n;++i){
	if (i!=s && i!=t){
	    cut[i]=false;
	    makegraph();memset(c,0,sizeof(c));memset(d,0,sizeof(d));
	    int tmp=0;c[0]=2*n;
	    while (d[s]<2*n) tmp+=maxflow(s+n,INF);
	    if (tmp==ans-1){
		ans--;
		print[++print_num]=i;
	    }
	    else{
		cut[i]=true;
	    }
	    if (!ans) break;
	}
    }
    for (int i=1;i<print_num;++i) fout << print[i] << " ";fout << print[print_num] << endl;
    return 0;
}


相关文章推荐

[usaco]5.4.5 Big Barn题解

5   Big Barn A Special Treat Farmer John wants to place a big square barn on his square farm. He...

【USACO5.4.5】奶牛的电信(BSOI2155)

【USACO5.4.5】奶牛的电信 Description   农夫约翰的奶牛们喜欢通过电邮保持联系,于是她们建立了一个奶牛电脑网络,以便互相交流。这些机器用如下的方式发送电邮:如果存在一个由c...
  • hwzzyr
  • hwzzyr
  • 2017年02月20日 18:09
  • 181

基于CDH5.4.5(ha)的Hbase 3节点搭建

hbase安装节点: cdhnode3 192.168.80.138 Hmaster cdhnode4 192.168.80.139 HRegionServer cdhno...

在CDH5.4.5版本中使用Scala+sbt+Eclipse+Java遇到的问题总结

在CDH5.4.5中使用Scala编写代码,sbt打包到集群中spark执行任务遇到的一些问题
  • px_dn
  • px_dn
  • 2017年03月02日 15:05
  • 352

hive-jdbc-1.1.0-cdh5.4.5-standalone.jar

  • 2017年12月07日 15:28
  • 18.42MB
  • 下载

CDH5.4.5手动安装hive-1.1.0-cdh5.4.5

本来想在单节点上面玩玩hive和hbase的整合,结果由于单节点上面的hive和hbase版本不兼容,整合的时候报错:org.apache.hadoop.hbase.HTableDescriptor....

PHP-5.4.5.tar.gz

  • 2014年12月02日 21:36
  • 52KB
  • 下载

USACO ORZ(hdu4277 深搜剪枝)

题目链接 题意:  给出n个数, 用这些数组成三个数, 问能够拼成三角形的数目 思路: 深搜, 将这n个数分成三堆再判断是否组成三角形, 但是肯定会爆掉, 所以中间加上set判重剪枝。 虽然在中间...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Usaco 5.4.5
举报原因:
原因补充:

(最多只允许输入30个字)