关闭

Usaco 5.4.5

标签: insert网络c
336人阅读 评论(0) 收藏 举报
分类:

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

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

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

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

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

可以把在网络中的每个点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;
}


0
0

猜你在找
【直播】计算机视觉原理及实战—屈教授
【套餐】深度学习入门视频课程—唐宇迪
【套餐】Hadoop生态系统零基础入门--侯勇蛟
【套餐】嵌入式Linux C编程基础--朱有鹏
【套餐】2017软考系统集成项目——任铄
【套餐】Android 5.x顶级视频课程——李宁
【直播】广义线性模型及其应用——李科
【直播】从0到1 区块链的概念到实践
【直播】机器学习之凸优化——马博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:32100次
    • 积分:669
    • 等级:
    • 排名:千里之外
    • 原创:34篇
    • 转载:0篇
    • 译文:0篇
    • 评论:6条
    最新评论