poj 1459 Power Network （测模板）

最新推荐文章于 2019-10-01 10:52:14 发布
最新推荐文章于 2019-10-01 10:52:14 发布
阅读量532 收藏
点赞数
分类专栏： 网络流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/gotoac/article/details/8138859
版权

题目链接:http://poj.org/problem?id=1459


题目大意:

多源多汇网络流.


题目思路:

学习网络流的第一题,也是拿来测模板.

定一个虚源点指向所有源点,容量为每个原点的Pmax.

定一个虚汇点使所有汇点指向它,容量为每个汇点的Cmax.

流一下得到最大流即可.

该题用邻接表写速度可能会变慢,因为邻接表里,会出现重边,但是不影响最终结果.


代码:

Edmonds-Karp


#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;

#define ll long long
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle (l+r)>>1
#define eps (1e-9)
#define type int
#define clr_all(x,c) memset(x,c,sizeof(x))
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n+1))
#define MOD 1000000007
#define inf 0x3f3f3f3f
#define pi acos(-1.0)

template <class T> void _swap(T &x,T &y){T t=x;x=y;y=t;}
template <class T> T _max(T x,T y){return x>y? x:y;}
template <class T> T _min(T x,T y){return x<y? x:y;}
int test,cas;
const int M=100 +5;

int n,m,np,nc;
int cap[M][M];
int flow[M][M],fa[M];
int q[M],head,tail;
char str[22];

int maxFlow(int s,int t){
	clr_all(flow,0);
	int f=0;
	for(;;){
		clr_all(fa,-1);
		for(q[head=tail=0]=s;head<=tail;){
			int u=q[head++];
			for(int v=0;v<=n+2;v++){
				if(fa[v]==-1 && flow[u][v]<cap[u][v]){
					fa[v]=u,q[++tail]=v;
				}
			}
		}
		if(fa[t]==-1) break;
		int a=inf;
		for(int u=t;u!=s;u=fa[u])
			a=_min(a,cap[fa[u]][u]-flow[fa[u]][u]);
		for(int u=t;u!=s;u=fa[u]){
			flow[fa[u]][u]+=a;
			flow[u][fa[u]]-=a;
		}
		f+=a;
	}
	return f;
}

void run(){
	int i,j;
	clr_all(cap,0);
	int u,v,cp;
	while(m--){
		scanf("%s",str);
		sscanf(str,"(%d,%d)%d",&u,&v,&cp);
		++u,++v;
		cap[u][v]=cp;
	}
	while(np--){
		scanf("%s",str);
		sscanf(str,"(%d)%d",&u,&cp);
		++u;
		cap[0][u]=cp;
	}
	while(nc--){
		scanf("%s",str);
		sscanf(str,"(%d)%d",&u,&cp);
		++u;
		cap[u][n+2]=cp;
	}
	printf("%d\n",maxFlow(0,n+2));
}

void preSof(){
}

int main(){
	//freopen("input.txt","r",stdin);
	//freopen("output.txt","w",stdout);
	preSof();
	//run();
	while(~scanf("%d%d%d%d",&n,&np,&nc,&m)) run();
	//for(scanf("%d",&test),cas=1;cas<=test;cas++) run();
	return 0;
}


Dinic(邻接表)


#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;

#define ll long long
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle (l+r)>>1
#define eps (1e-9)
#define type int
#define clr_all(x,c) memset(x,c,sizeof(x))
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n+1))
#define MOD 1000000007
#define inf 0x3f3f3f3f
#define pi acos(-1.0)

template <class T> void _swap(T &x,T &y){T t=x;x=y;y=t;}
template <class T> T _max(T x,T y){return x>y? x:y;}
template <class T> T _min(T x,T y){return x<y? x:y;}
int test,cas;
const int M=100 +5;

int n,m,np,nc;
char str[22];

struct edge{
	int from,to,cap,flow;
	edge(){}
	edge(int u,int v,int c,int f){from=u,to=v,cap=c,flow=f;}
};

struct dinic{
	vector<edge>edg;      //边表
	vector<int>g[M];      //邻接表:g[i][j]表示结点i连接的第j条边在edg中的序号
	int n,m,s,t;          //顶点数,边数(包括反向弧),源点,汇点
	int dis[M],cur[M];    //从起点到i的距离,当前弧下标
	int q[M],head,tail;
	void clear(){
		edg.clear();
		for(int i=0;i<=n;i++) g[i].clear();
	}
	void insert(int u,int v,int c){
		edg.push_back(edge(u,v,c,0));
		edg.push_back(edge(v,u,0,0));
		m=edg.size();
		g[u].push_back(m-2);
		g[u].push_back(m-1);
	}
	int bfs(){
		clr(dis,0,n);dis[s]=1;
		for(q[head=tail=0]=s;head<=tail;){
			int u=q[head++];
			for(int i=0;i<g[u].size();i++){
				edge& e=edg[g[u][i]];
				if(!dis[e.to] && e.cap>e.flow){
					dis[e.to]=dis[e.from]+1;
					q[++tail]=e.to;
				}
			}
		}
		return dis[t];
	}
	int dfs(int u,int a){
		if(u==t || a==0) return a;
		int flow=0,f;
		for(int& i=cur[u];i<g[u].size();i++){
			edge& e=edg[g[u][i]];
			if(dis[e.to]==dis[u]+1 && (f=dfs(e.to,_min(a,e.cap-e.flow)))>0){
				e.flow+=f,edg[g[u][i]^1].flow-=f;
				flow+=f,a-=f;
				if(a==0) break;
			}
		}
		return flow;
	}
	int maxFlow(int _s,int _t){
		s=_s,t=_t;
		int flow=0;
		while(bfs()){
			clr(cur,0,n);
			flow+=dfs(s,inf);
		}
		return flow;
	}
}p;

void run(){
	int i,j;
	int u,v,cp;
	p.clear();
	p.n=n+2;
	while(m--){
		scanf("%s",str);
		sscanf(str,"(%d,%d)%d",&u,&v,&cp);
		p.insert(u+1,v+1,cp);
	}
	while(np--){
		scanf("%s",str);
		sscanf(str,"(%d)%d",&u,&cp);
		p.insert(0,u+1,cp);
	}
	while(nc--){
		scanf("%s",str);
		sscanf(str,"(%d)%d",&u,&cp);
		p.insert(u+1,n+2,cp);
	}
	printf("%d\n",p.maxFlow(0,n+2));
}

void preSof(){
}

int main(){
	//freopen("input.txt","r",stdin);
	//freopen("output.txt","w",stdout);
	preSof();
	//run();
	while(~scanf("%d%d%d%d",&n,&np,&nc,&m)) run();
	//for(scanf("%d",&test),cas=1;cas<=test;cas++) run();
	return 0;
}

SAP(邻接表)


#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <ctype.h>
#include <math.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;

#define ll long long
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle (l+r)>>1
#define eps (1e-9)
#define type int
#define clr_all(x,c) memset(x,c,sizeof(x))
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n+1))
#define MOD 1000000007
#define inf 0x3f3f3f3f
#define pi acos(-1.0)
#define _max(x,y) (((x)>(y))? (x):(y))
#define _min(x,y) (((x)<(y))? (x):(y))
#define _abs(x) ((x)<0? (-(x)):(x))
#define getmin(x,y) (x= (x==-1 || (y)<x)? (y):x)
#define getmax(x,y) (x= ((y)>x)? (y):x)
template <class T> void _swap(T &x,T &y){T t=x;x=y;y=t;}
int test,cas;
const int M=100 +5;

int n,m,np,nc;
char str[22];

struct sap{
	struct edge{
		int from,to,flow;
		edge(){}
		edge(int u,int v,int f){from=u,to=v,flow=f;}
	};
	int n,m,s,t;          //顶点数,边数(包括反向弧),源点,汇点
	int dis[M],cur[M];    //从起点到i的距离,当前所在的弧
	int gap[M],pre[M];    //距离标号计数,可增广路上的上一条弧的序号
	vector<edge>edg;      //边表
	vector<int>g[M];      //g[i][j]表示点i连接的第j条边的序号
	void clear(){edg.clear();for(int i=0;i<n;i++) g[i].clear();}
	void insert(int u,int v,int c){
		edg.push_back(edge(u,v,c)),edg.push_back(edge(v,u,0));
		m=edg.size();
		g[u].push_back(m-2),g[v].push_back(m-1);
	}
	int maxFlow(int _s,int _t){
		s=_s,t=_t;
		int res=0;
		clr(dis,0,n),clr(gap,0,n),clr(cur,0,n);
		gap[0]=n;
		for(int u=s,a;dis[s]<n;){
			if(u==t){//增广
				for(a=-1;u!=s;u=edg[pre[u]].from)
					getmin(a,edg[pre[u]].flow);
				for(u=t;u!=s;u=edg[pre[u]].from)
					edg[pre[u]].flow-=a,edg[pre[u]^1].flow+=a;
				res+=a;
			}
			bool ok=0;
			for(int i=cur[u];i<g[u].size();i++){//前推
				edge& e=edg[g[u][i]];
				if(dis[u]==dis[e.to]+1 && e.flow){
					pre[e.to]=g[u][i],cur[u]=i,u=e.to,ok=1;break;
				}
			}
			if(!ok){//回退
				int mindis=n-1;
				for(int i=0;i<g[u].size();i++){
					edge& e=edg[g[u][i]];
					if(e.flow && mindis>dis[e.to])
						mindis=dis[e.to],cur[u]=i;
				}
				if(--gap[dis[u]]==0) break;
				gap[dis[u]=mindis+1]++;
				if(u!=s) u=edg[pre[u]].from;
			}
		}
		return res;
	}
}p;

void run(){
	int i,j;
	int u,v,cp;
	p.clear();
	p.n=n+3;
	while(m--){
		scanf("%s",str);
		sscanf(str,"(%d,%d)%d",&u,&v,&cp);
		p.insert(u+1,v+1,cp);
	}
	while(np--){
		scanf("%s",str);
		sscanf(str,"(%d)%d",&u,&cp);
		p.insert(0,u+1,cp);
	}
	while(nc--){
		scanf("%s",str);
		sscanf(str,"(%d)%d",&u,&cp);
		p.insert(u+1,n+2,cp);
	}
	printf("%d\n",p.maxFlow(0,n+2));
}

void preSof(){
}

int main(){
	//freopen("input.txt","r",stdin);
	//freopen("output.txt","w",stdout);
	preSof();
	//run();
	while(~scanf("%d%d%d%d",&n,&np,&nc,&m)) run();
	//for(scanf("%d",&test),cas=1;cas<=test;cas++) run();
	return 0;
}


gotoac
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
POJ 1459 Power Network
u010084308的专栏
04-16 727
链接：http://poj.org/problem?id=1459 题目： Power Network Time Limit: 2000MS   Memory Limit: 32768K Total Submissions: 21928   Accepted: 11453 Description A power netwo
POJ1459 Power Network（网络流）
KeyboardPianist的博客
10-15 415
Power Network Time Limit: 2000MS   Memory Limit: 32768K Total Submissions: 27409   Accepted: 14256 Description A power network consists of nodes (power stations, consumers
poj 1459 Power Network
aizhengcuo9261的博客
10-01 74
A的第一道网络流，嘿嘿，相比较来说蛮合适的。题目不难，套模板，注意加超级源、超级汇。难就难在能把题看完，题目很纠结啊。。 题目大意就是说，输电网络中:p是发电站，c是用户。由Sample Input：输入边、p、c ，求最大流。 多源多汇！？加个源点连接所有p，所有c 连接汇点，又是单源单汇了。 ek 1 #include<stdio....
POJ1459Power Network（dinic模板）
apk6909的博客
12-04 95
Power Network Time Limit:2000MS Memory Limit:32768K Total Submissions:25832 Accepted:13481 Description A power network consists of nodes (power stations, consu...
POJ 1459 Power Network 网络流
Adolphrocs的博客
12-03 182
题目源：https://vjudge.net/problem/POJ-1459 http://poj.org/problem?id=1459 模板源：https://blog.csdn.net/Adolphrocs/article/details/84368661 https://blog.csdn.net/Adolphrocs/article/details/84779691 https...
POJ-1459 Power Network
哇-WA 的博客
12-22 147
Power Network Time Limit: 2000MS   Memory Limit: 32768K Total Submissions: 29271   Accepted: 15192 Description A power network consists of nodes (power stations, consumer
POJ 1459 power network(网络流 dinic模版)
徐伯莱
08-26 134
题意：给几个发电站，给几个消耗站，再给几个转发点。 发电站只发电，消耗站只消耗电，转发点只是转发电，再给各个传送线的传电能力。 问你消耗站能获得的最多电是多少. #include &lt;cstdio&gt; #include &lt;cstring&gt; #include &lt;iostream&gt; #include &lt;string&gt; #include &lt;algor...
poj1459 Power Network题解
xxtz_xt的博客
07-22 214
传送门 A power network consists of nodes (power stations, consumers and dispatchers) connected by power transport lines. A node u may be supplied with an amount s(u) >= 0 of power, may produce an am
poj1459 Power Network
AaronPolaris
12-08 1648
一道网络流题，还是比较简单的
poj 1459 Power Network.md
11-17
poj 1459 Power Network.md
POJ1459-Power Network
07-30
【标题】"POJ1459-Power Network"是一个编程竞赛题目，源自北京大学的在线判题系统POJ（Problem Set of Peking University）。这个题目主要涉及到图论和算法的应用，特别是广度优先搜索（BFS）以及一种优化技巧——...
poj2488.rar_poj24_poj2488_方向模板法
09-14
标题中的“poj2488.rar_poj24_poj2488_方向模板法”指的是一个解决编程竞赛问题POJ2488的压缩文件，其中可能包含了解决该问题的代码示例。POJ是Programming Online Judge的缩写，是一个在线的编程竞赛平台，参与者...
POJ2531-Network Saboteur
07-31
【标题】"POJ2531-Network Saboteur" 是一道来自北京大学在线判题系统POJ（Problem Set）的编程题目。该题目属于算法竞赛中的经典问题，旨在考察参赛者的图论知识和编程能力。 【描述】"北大POJ2531-Network ...
POJ2109-Power of Cryptography
07-29
《POJ2109-Power of Cryptography：解题报告与AC代码解析》 在计算机科学领域，算法和编程竞赛是检验和提升编程技能的重要途径。POJ（Problem Set of Peking University）是由北京大学主办的一个在线编程竞赛平台，...
nbu 2431 星际转移问题
Go.To.AC的专栏
12-13 763
题目链接：http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2431   题目大意： 中文题意，不解释。   题目思路： 源自王晓东《线性规划和网络流24题》。 原题的解法，一种枚举答案+动态建图，一种二分答案+多次建图，因为原题最多30天。 但是这题最多可以是1000多天，所以二分答案+多次建图会TLE。 不过两种解法的思路
nbu 2428 Risk
Go.To.AC的专栏
12-14 575
题目链接:http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2428 题目大意: 有n个区域(1 现在要通过一次调兵(设从a调兵,那么只能调到和a相邻的自己的区域)使得和敌人接壤的区域的最小士兵数最大,且自己的所有区域的兵的数量为正. 题目思路: 假设答案为x. (1)建立一个超级源点S,并把所有的敌人的区域看成一个
nbu 2429 Transfer stations
Go.To.AC的专栏
12-14 568
题目链接:http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2429 题目大意: 现在有一项工程,有n(n 而假设a站点和b站点均建立了,那么可以获利C,这样的关系有m个(m 问建立那些站点可获得最大纯利润. 题目思路: 裸的最大权闭合图. 这种题型的建图有两种方法. 一:(建完后的图的原图完全不同) (1)将
nbu 2427 Pigs
Go.To.AC的专栏
12-14 532
题目链接:http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2427 题目大意: 有m(1)个猪圈,n(1 第i人第i天去买猪,有ki把钥匙,最多买bi头猪. 第i天的时候,钥匙对应的猪圈会打开,这个时候可以对猪圈任意分配猪的数量(和不会超过打开的猪圈的猪的数量之和.) 问最多能卖多少猪. 题目思路: (1)建立一
西工大C语言POJ习题答案.docx
最新发布
11-15
西工大C语言POJ习题答案.docx文件包含了一季度的10道题目的答案，其中第五题的问题已经得到解决并且确认为正确答案。这份答案文件是由西安工业大学的学生们共同完成，并且经过严格的检验与确认。在这份答案文件中，...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值