poj2516----最小费用流(模版+构图)
来源:互联网 发布:数据分析报告 编辑:程序博客网 时间:2024/05/16 06:36
【题意】
有n个商店m个供应商k种商品
给各商店的各种商品的需求量和各供应商的供应量以及供应商和商店之间的运输费用,求满足需求的最小费用
【输入】
多组数据
每组数据第一行n,m,k(均小于50)
接下来n行每行k个数,表示需求量(0到3)
接下来m行每行k个数,表示供应量(0到3)
接下来k个n*m矩阵表示商店和供应商之间的运输费用
输入由第一行三个0结束
【输出】
每组数据输出一个数,表示最小费用
//在多重循环中一定注意下标容易出错,i,j,k!!!
#include<cstdlib>#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include <algorithm>#include<queue>#include<set>#define LL long long#define inf 0x7fffffff#define E 1e-9#define M 100#define N 105using namespace std;int nead[N][N],have[N][N];int n,m,k;int cost[N][N];int flow[N][N],cap[N][N];int d[N],p[N];bool inq[N];int s,t;int num;int getflow(int s,int t){ queue<int> q; memset(flow,0,sizeof(flow)); int c=0,f=0; for(;;) { memset(d,0x7f,sizeof(d));// memset(inq,0,sizeof(inq)); d[s]=0; q.push(s); while(!q.empty()) { int u=q.front(); q.pop(); inq[u]=0; for(int i=1; i<=num; i++) if(cap[u][i]>flow[u][i]&&d[i]>d[u]+cost[u][i]) { d[i]=d[u]+cost[u][i]; p[i]=u; if(!inq[i]) { q.push(i); inq[i]=1; } } } if(d[t]==0x7f7f7f7f)////0x7f7f7f7f break; int a=inf; for(int i=t; i!=s; i=p[i]) { a=min(a,cap[p[i]][i]-flow[p[i]][i]); } for(int i=t; i!=s; i=p[i]) { flow[p[i]][i]+=a; flow[i][p[i]]-=a; } f+=a; c+=d[t]*a; } return c;}int main(){#ifndef ONLINE_JUDGE freopen("ex.in","r",stdin);#endif while(scanf("%d%d%d",&n,&m,&k)) { int ans=0; if(!n&&!m&&!k) break; for(int i=1; i<=n; i++) for(int j=1; j<=k; j++) scanf("%d",&nead[i][j]); for(int i=1; i<=m; i++) for(int j=1; j<=k; j++) scanf("%d",&have[i][j]); int flag=0; for(int i=1; i<=k; i++) { int nead1=0; for(int j=1; j<=n; j++) nead1+=nead[j][i]; int have1=0; for(int j=1; j<=m; j++) have1+=have[j][i]; if(have1<nead1) { flag=1; break; } } for(int p=1; p<=k; p++)//在多重循环中一定注意下标容易出错,i,j,k!!! { memset(cost,0,sizeof(cost)); memset(cap,0,sizeof(cap)); for(int j=1; j<=n; j++) for(int k=1; k<=m; k++) { scanf("%d",&cost[j][k+n]); cost[k+n][j]=-cost[j][k+n]; cap[j][k+n]=nead[j][p];// cap[k+n][j]=0; } if(flag) continue; for(int j=1; j<=n; j++) { cap[0][j]=nead[j][p]; } num=n+m+1; for(int j=1; j<=m; j++) { cap[j+n][num]=have[j][p]; } s=0; t=num; ans+=getflow(s,t); } if(flag) printf("-1\n"); else printf("%d\n",ans); } return 0;}
- poj2516----最小费用流(模版+构图)
- poj2516(最小费用最大流模版)
- poj2516 最小费用流
- 【POJ2516】【最小费用流】
- poj2516 最小费用最大流
- poj2516 最小费用最大流
- poj2516 最小费用最大流
- 最小费用最大流poj2516
- (最小费用最大流)POJ2516 Minimum Cost
- poj2516——Minimum Cost(最小费用最大流)
- POJ2516 Minimum Cost 最小费用最大流
- POJ2516--Minimum Cost最小费用最大流
- POj2516 Minimum Cost 最小费用最大流
- POJ2516-Minimum Cost 最小费用最大流
- poj2516 minimum cost 最小费用流
- 最小费用流模版
- uva 11381(神奇的构图、最小费用最大流)
- 最小费用最大流粗解 poj2516
- linux的sh脚本编程
- Android文件存取
- [HDU-2602] 0-1背包
- 最近笔试面试没答全的题目总结
- 在xp下安装iis5.1注意事项
- poj2516----最小费用流(模版+构图)
- 创建Unicode版本的VC工程
- Fedora 17 安装后记
- hdu - 4324 - Triangle LOVE - 深搜
- 如何判断是iphone还是ipad
- 二叉查找树
- sql 查询本周本月问题
- Android传感器---Motion Sensor(二)
- C#中制作启动窗体的方法和问题