poj 3436 最大流的增广路算法

题目大意:

有N台机器,每台机器有P部分,每部分有自己的输入、输出,因此每台机器有2*P+1种参数:第一个参数Q:该机器的容量;接下来P个参数S:该机器各部分的输入条件;接下来P个参数D:该机器各部分的输出情况。

其中输入有三种情况:0,1,2

0:该部分必须不要;1:该部分必须要;2:该部分可有可无

输出有2种情况:0,1  

0:该部分不存在;1:该部分存在。

第一行输出最大流,所经过的路径数

第二行开始每行输出所经过的每条路径的三个数据,分别是开始点、结束点、流量。

解题思路:

每个点都有输入和输出的条件,因此可以把每个点拆成一个入点和一个出点。

每次拆点时都和其它点进行一次判断,如果出点能满足入点的条件就在两点之间加一条路径,路径的流量上限为两点上限的最小值。

最后设置两个源点,一个连接条件中没有1的入点,另一个连接输出全是1的出点,路径的流量上限为连接点的上限。


#include<iostream>
#include<stdio.h>
#include<math.h>
#include<string.h>
#include<vector>
#include<queue>
#include<stack>
#include<queue>
#include<algorithm>
using namespace std;
struct line //把边定义为结构体,每条边中间储存出发点、终点、容量、流量
{
	int from;
	int to;
	int cap;
	int flow;
	line(int f,int t,int c,int ff):from(f),to(t),cap(c),flow(ff){}
};
int min(int a,int b)
{
	if(a<b)return a;
	else return b;
}
int P,N,Q[205];
int maps[205][205];//maps中存点i的所有可到达边的序号maps[i][j]
vector<line>lines;//lines中存所有边
int a[205];//a中存出发点到i的可改进量
int p[205];//最短路上p[i]的入弧编号
void addline(int from,int to,int cap)//增加一条边
{
	lines.push_back(line(from,to,cap,0));
	lines.push_back(line(to,from,0,0));
	int k=lines.size();
	maps[from][maps[from][0]]=k-2;
	maps[from][0]++;
	maps[to][maps[to][0]]=k-1;
	maps[to][0]++;
}
int maxflow(int s,int e)//用bfs求最大流
{
	int flow=0;
	while(true)
	{
		memset(a,0,sizeof(a));
		queue<int>q;
		q.push(s);
		a[s]=10005;
		while(!q.empty())
		{
			int t=q.front();
			q.pop();
			for(int i=1;i<maps[t][0];i++)
			{
				line l=lines[maps[t][i]];
				if(!a[l.to]&&l.flow<l.cap)
				{
					a[l.to]=min(a[t],l.cap-l.flow);
					p[l.to]=maps[t][i];
					q.push(l.to);
				}
			}
			if(a[e])
				break;
		}
		if(!a[e])
			break;
		for(int i=e;i!=s;i=lines[p[i]].from)
		{
			lines[p[i]].flow+=a[e];
			lines[p[i]^1].flow-=a[e];
		}
		flow+=a[e];
	}
	return flow;
}
int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	while(scanf("%d%d",&P,&N)!=EOF)
	{
		int point[200][15];
		int mark[105][105];//标记两点间是否存在路径
		memset(a,0,sizeof(a));
		memset(p,0,sizeof(p));
		memset(point,0,sizeof(point));
		lines.clear();
		memset(mark,0,sizeof(mark));
		memset(maps,0,sizeof(maps));
		memset(Q,0,sizeof(Q));
		for(int i=0;i<205;i++)
			maps[i][0]=1;
		int num=0;
		for(int i=0;i<P;i++)//设置源点
		{
			point[2*N][i]=0;
			point[2*N+1][i]=1;
		}
		for(int i=0;i<2*N;i+=2)
		{
			scanf("%d",&Q[i/2]);
			int in=i,out=i+1;
			for(int j=0;j<P;j++)
				scanf("%d",&point[i][j]);
			for(int j=0;j<P;j++)
				scanf("%d",&point[i+1][j]);
			int flag1=0;
			for(int j=0;j<P;j++)
				if(point[in][j]!=point[out][j])//排除入点出点相同的点
				{
					flag1=1;
					break;
				}
				if(flag1==1)
				{
			addline(in,out,Q[i/2]);//在入点出点间建立路径
			for(int k=1;k<i;k+=2)
			{
				int flag=0;
				for(int j=0;j<P;j++)//出点满足入点条件
					if((point[in][j]!=point[k][j]&&point[in][j]!=2)||mark[k/2][in])
					{
						flag=1;
						break;
					}
				if(flag==0)
				{
					addline(k,in,min(Q[k/2],Q[in/2]));
					mark[k/2][in/2]=1;
				}
			}
			for(int k=0;k<i;k+=2)
			{
				int flag=0;
				for(int j=0;j<P;j++)
					if((point[out][j]!=point[k][j]&&point[k][j]!=2)||mark[out/2][k/2])
					{
						flag=1;
						break;
					}
				if(flag==0)
				{
					addline(out,k,min(Q[k/2],Q[out/2]));
					mark[out/2][k/2]=1;
				}
			}
			}
		}
		for(int i=0;i<2*N;i+=2)//建立源点与入点的路径
		{
			int flag=0;
			for(int j=0;j<P;j++)
				if(point[i][j]==1)
				{
					flag=1;
					break;
				}
				if(flag==0)
					addline(2*N,i,Q[i/2]);
		}
		for(int i=1;i<2*N;i+=2)//建立出点与源点的路径
		{
			int flag=0;
			for(int j=0;j<P;j++)
				if(point[i][j]!=1)
				{
					flag=1;
					break;
				}
				if(flag==0)
					addline(i,2*N+1,Q[i/2]);
		}
		int flow=maxflow(2*N,2*N+1);
		int sum=0;
		int mem[10005][3];
			for(int i=0;i<lines.size();i++)
				if(lines[i].flow>0)//查找路径
					{
						if(lines[i].from/2==lines[i].to/2||lines[i].from>=2*N||lines[i].to>=2*N)continue;
						else
						{
						mem[sum][0]=(lines[i].from/2)+1;
						mem[sum][1]=(lines[i].to/2)+1;
						mem[sum][2]=lines[i].flow;
						sum++;
						}
				}
		printf("%d %d\n",flow,sum);
		if(sum!=0)
			for(int i=0;i<sum;i++)
				printf("%d %d %d\n",mem[i][0],mem[i][1],mem[i][2]);
		}
	return 0;
}


#include<iostream>#include<stdio.h>#include<math.h>#include<string.h>#include<vector>#include<queue>#include<stack>#include<queue>#include<algorithm>using namespace std;struct line //把边定义为结构体,每条边中间储存出发点、终点、容量、流量{int from;int to;int cap;int flow;line(int f,int t,int c,int ff):from(f),to(t),cap(c),flow(ff){}};int min(int a,int b){if(a<b)return a;else return b;}int P,N,Q[205];int maps[205][205];//maps中存点i的所有可到达边的序号maps[i][j]vector<line>lines;//lines中存所有边int a[205];//a中存出发点到i的可改进量int p[205];//最短路上p[i]的入弧编号void addline(int from,int to,int cap)//增加一条边{lines.push_back(line(from,to,cap,0));lines.push_back(line(to,from,0,0));int k=lines.size();maps[from][maps[from][0]]=k-2;maps[from][0]++;maps[to][maps[to][0]]=k-1;maps[to][0]++;}int maxflow(int s,int e)//用bfs求最大流{int flow=0;while(true){memset(a,0,sizeof(a));queue<int>q;q.push(s);a[s]=10005;while(!q.empty()){int t=q.front();q.pop();for(int i=1;i<maps[t][0];i++){line l=lines[maps[t][i]];if(!a[l.to]&&l.flow<l.cap){a[l.to]=min(a[t],l.cap-l.flow);p[l.to]=maps[t][i];q.push(l.to);}}if(a[e])break;}if(!a[e])break;for(int i=e;i!=s;i=lines[p[i]].from){lines[p[i]].flow+=a[e];lines[p[i]^1].flow-=a[e];}flow+=a[e];}return flow;}int main(){//freopen("in.txt","r",stdin);//freopen("out.txt","w",stdout);while(scanf("%d%d",&P,&N)!=EOF){int point[200][15];int mark[105][105];//标记两点间是否存在路径memset(a,0,sizeof(a));memset(p,0,sizeof(p));memset(point,0,sizeof(point));lines.clear();memset(mark,0,sizeof(mark));memset(maps,0,sizeof(maps));memset(Q,0,sizeof(Q));for(int i=0;i<205;i++)maps[i][0]=1;int num=0;for(int i=0;i<P;i++)//设置源点{point[2*N][i]=0;point[2*N+1][i]=1;}for(int i=0;i<2*N;i+=2){scanf("%d",&Q[i/2]);int in=i,out=i+1;for(int j=0;j<P;j++)scanf("%d",&point[i][j]);for(int j=0;j<P;j++)scanf("%d",&point[i+1][j]);int flag1=0;for(int j=0;j<P;j++)if(point[in][j]!=point[out][j])//排除入点出点相同的点{flag1=1;break;}if(flag1==1){addline(in,out,Q[i/2]);//在入点出点间建立路径for(int k=1;k<i;k+=2){int flag=0;for(int j=0;j<P;j++)//出点满足入点条件if((point[in][j]!=point[k][j]&&point[in][j]!=2)||mark[k/2][in]){flag=1;break;}if(flag==0){addline(k,in,min(Q[k/2],Q[in/2]));mark[k/2][in/2]=1;}}for(int k=0;k<i;k+=2){int flag=0;for(int j=0;j<P;j++)if((point[out][j]!=point[k][j]&&point[k][j]!=2)||mark[out/2][k/2]){flag=1;break;}if(flag==0){addline(out,k,min(Q[k/2],Q[out/2]));mark[out/2][k/2]=1;}}}}for(int i=0;i<2*N;i+=2)//建立源点与入点的路径{int flag=0;for(int j=0;j<P;j++)if(point[i][j]==1){flag=1;break;}if(flag==0)addline(2*N,i,Q[i/2]);}for(int i=1;i<2*N;i+=2)//建立出点与源点的路径{int flag=0;for(int j=0;j<P;j++)if(point[i][j]!=1){flag=1;break;}if(flag==0)addline(i,2*N+1,Q[i/2]);}int flow=maxflow(2*N,2*N+1);int sum=0;int mem[10005][3];for(int i=0;i<lines.size();i++)if(lines[i].flow>0)//查找路径{if(lines[i].from/2==lines[i].to/2||lines[i].from>=2*N||lines[i].to>=2*N)continue;else{mem[sum][0]=(lines[i].from/2)+1;mem[sum][1]=(lines[i].to/2)+1;mem[sum][2]=lines[i].flow;sum++;}}printf("%d %d\n",flow,sum);if(sum!=0)for(int i=0;i<sum;i++)printf("%d %d %d\n",mem[i][0],mem[i][1],mem[i][2]);}return 0;}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值