POj2516 Minimum Cost 最小费用最大流
来源:互联网 发布:中国印度边境冲突 知乎 编辑:程序博客网 时间:2024/05/16 08:55
题目链接:
poj2516
题意:
Dearboy 是一个优秀的食品供应商,他现在面临一个大问题,需要你的帮忙。在他的销售地区,有N 个店主(编号从1~N)帮他销售食品。Dearboy 有M 仓库(编号从1~M),每个仓库第6 章 网络流问题可以提供K 种不同的食品(编号从1~K)。一旦有店主向他订食品,Dearboy 应该安排哪个仓库、向该店主提供多少食品,以减少总的运输费用?
现已知道,从不同的仓库向不同的店主运输不同种类的单位重量食品所需的费用是不同的。给定每个仓库K 种食品格子的储藏量,N 个店主对K 种食品的订量,以及从不同的仓库运输不同种食品到不同的店主的所需费用,你的任务是安排每个仓库的各种食品供应量,以减少总的运输费用。
解题思路:
题目中有 需求数量(容量) , 运输费用(费用) 很容易想到最小费用最大流
但是 题目还给出了一个K 表示品种数量,同时每个品种的需求数量,运输费用各不相同;也就是说如果要构建网络,每条边要抽象为k条边,每条边表示一种费用.orz
重新浏览题目,注意到 题目给出的k很大(50) ,但给出的点很少(50+50+2),
这样我们就可以针对每一个品种,建一次图,跑一次最小最小费用最大流即可.
最后将所有品种的最小费用加起来即可.
考虑答案为-1的情况:仅当某类产品供不应求时才会出现 , 因为给出的k个费用矩阵是全部有值的.
代码:
#include<iostream>#include<cstdio>#include<cstring>#define maxn 205#define INF 0x3f3f3f3fusing namespace std;struct node{ int to,c,f,w,next;} edge[maxn*maxn];int head[maxn],ss;int supply[55],need[55];int a[55][55],b[55][55];int cost[55][55][55];int dis[maxn],vis[maxn],pre[maxn];int n,m,k,s,t;void init(){ memset(head,-1,sizeof(head)); ss=0;}void addedge(int a,int b,int c,int d){ edge[ss]=(node){b,c,0,d,head[a]}; head[a]=ss++; edge[ss]=(node){a,c,c,-d,head[b]}; head[b]=ss++;}int spfa(){ int q[maxn*10]; memset(vis,0,sizeof(vis)); memset(dis,INF,sizeof(dis)); int u,v,head1=0,tail=0; q[tail++]=s; dis[s]=0; while(head1<tail) { u=q[head1++]; vis[u]=0; for(int i=head[u]; i!=-1; i=edge[i].next) { v=edge[i].to; if(edge[i].c-edge[i].f>0 && dis[v]>dis[u]+edge[i].w) { dis[v]=dis[u]+edge[i].w; pre[v]=i; if(!vis[v]) { vis[v]=1; q[tail++]=v; } } } } if(dis[t]<INF) return 1; return 0;}int MCMF(){ int minn; int sum=0; while(spfa()) { minn=INF; //最多能买minn个 for(int i=t; i!=s; i=edge[pre[i]^1].to) { if(edge[pre[i]].c-edge[pre[i]].f<minn) minn=edge[pre[i]].c-edge[pre[i]].f; } for(int i=t; i!=s; i=edge[pre[i]^1].to) { edge[pre[i]].f+=minn; edge[pre[i]^1].f-=minn; } sum+=minn*dis[t]; } return sum;}void build(int kind){ init(); for(int i=1; i<=m; i++) addedge(s,i,b[i][kind],0); for(int j=1; j<=n; j++) addedge(j+m,t,a[j][kind],0); for(int i=1; i<=n; i++) for(int j=1; j<=m; j++) addedge(j,m+i,INF,cost[kind][i][j]);}int main(){// freopen("in.txt","r",stdin); while(scanf("%d%d%d",&n,&m,&k)&&(n+m+k)) { s=0,t=n+m+1; int flag=0,ans=0; memset(supply,0,sizeof(supply)); memset(need,0,sizeof(need)); for(int i=1; i<=n; i++) for(int j=1; j<=k; j++) { scanf("%d",&a[i][j]); need[j]+=a[i][j]; } for(int i=1; i<=m; i++) for(int j=1; j<=k; j++) { scanf("%d",&b[i][j]); supply[j]+=b[i][j]; } for(int i=1; i<=k; i++) if(supply[i]<need[i]) //供不应求 flag=1; for(int kk=1; kk<=k; kk++) for(int i=1; i<=n; i++) for(int j=1; j<=m; j++) scanf("%d",&cost[kk][i][j]); if(flag) { printf("-1\n"); continue; } for(int i=1; i<=k; i++) //建k次图 跑k次最小费用最大流 orz { build(i); ans+=MCMF(); } printf("%d\n",ans); } return 0;}
1 0
- POJ2516 Minimum Cost 最小费用最大流
- POJ2516--Minimum Cost最小费用最大流
- POj2516 Minimum Cost 最小费用最大流
- POJ2516-Minimum Cost 最小费用最大流
- POJ2516 Minimum Cost 最大流最小费用典型题目
- (最小费用最大流)POJ2516 Minimum Cost
- poj2516——Minimum Cost(最小费用最大流)
- poj2516 minimum cost 最小费用流
- POJ2516 Minimum Cost(K次费用流)
- poj2516 Minimum Cost(费用流)
- PKU2516 Minimum Cost - 最小费用最大流
- Minimum Cost-最小费用最大流
- poj2516 最小费用最大流
- poj2516 最小费用最大流
- poj2516 最小费用最大流
- 最小费用最大流poj2516
- Minimum Cost----最小费用流
- poj2516--Minimum Cost(费用流,分别建图)
- C#高级编程----反射的小结
- 线程的三种创建方法
- 基本命令(按学习进度随机添加)
- ORM进阶之 ORM简介
- poj 2752 Seek the Name, Seek the Fame 【在原串中 找所有既是原串前缀又是原串后缀的子串长度】
- POj2516 Minimum Cost 最小费用最大流
- HDOJ Zipper (DFS)
- style-visibility
- java 深入解析FileInputStream和FileOutputStream
- java 序列化与反序列化
- 杭电 2087 剪花布条 kmp水题
- c++/c中流的概念
- UVa 10806 Dijkstra, Dijkstra (无向图+最小费用最大流+最短路:)
- HDOJ 5348 MZL's endless loop