POJ - 2516 Minimum Cost 每次要跑K次费用流
题意:
有m个仓库,n个买家,k个商品,每个仓库运送不同商品到不同买家的路费是不同的。问为了满足不同买家的订单的最小的花费。
思路:
设立一个源点S和汇点T,从源点S到每个仓库(1~m)连上容量为商品A的库存、费用为0的边,每个仓库再向每个不同的买家连上容量inf,费用为路费的边、每个顾客向汇点连一条容量为自己对商品A的需求个数、费用为0的边。跑一边费用流即可。这只有运送一个商品的费用,对,那我们就对不同商品建不同的图,一共跑K边,累计答案即可。
自己一直在想建图,怎么一次就跑出费用。感觉上面的思路和多开一维的观点差不多。
//#pragma GCC optimize(3)
//#pragma comment(linker, "/STACK:102400000,102400000") //c++
// #pragma GCC diagnostic error "-std=c++11"
// #pragma comment(linker, "/stack:200000000")
// #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
// #pragma GCC optimize("-fdelete-null-pointer-checks,inline-functions-called-once,-funsafe-loop-optimizations,-fexpensive-optimizations,-foptimize-sibling-calls,-ftree-switch-conversion,-finline-small-functions,inline-small-functions,-frerun-cse-after-loop,-fhoist-adjacent-loads,-findirect-inlining,-freorder-functions,no-stack-protector,-fpartial-inlining,-fsched-interblock,-fcse-follow-jumps,-fcse-skip-blocks,-falign-functions,-fstrict-overflow,-fstrict-aliasing,-fschedule-insns2,-ftree-tail-merge,inline-functions,-fschedule-insns,-freorder-blocks,-fwhole-program,-funroll-loops,-fthread-jumps,-fcrossjumping,-fcaller-saves,-fdevirtualize,-falign-labels,-falign-loops,-falign-jumps,unroll-loops,-fsched-spec,-ffast-math,Ofast,inline,-fgcse,-fgcse-lm,-fipa-sra,-ftree-pre,-ftree-vrp,-fpeephole2",3) #include <algorithm>
#include <iterator>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <iomanip>
#include <bitset>
#include <cctype>
#include <cstdio>
#include <string>
#include <vector>
#include <stack>
#include <cmath>
#include <queue>
#include <list>
#include <map>
#include <set>
#include <cassert> using namespace std;
#define lson (l , mid , rt << 1)
#define rson (mid + 1 , r , rt << 1 | 1)
#define debug(x) cerr << #x << " = " << x << "\n";
#define pb push_back
#define pq priority_queue typedef long long ll;
typedef unsigned long long ull;
//typedef __int128 bll;
typedef pair<ll ,ll > pll;
typedef pair<int ,int > pii;
typedef pair<int,pii> p3; //priority_queue<int> q;//这是一个大根堆q
//priority_queue<int,vector<int>,greater<int> >q;//这是一个小根堆q
#define fi first
#define se second
//#define endl '\n' #define OKC ios::sync_with_stdio(false);cin.tie(0)
#define FT(A,B,C) for(int A=B;A <= C;++A) //用来压行
#define REP(i , j , k) for(int i = j ; i < k ; ++i)
#define max3(a,b,c) max(max(a,b), c);
//priority_queue<int ,vector<int>, greater<int> >que; const ll mos = 0x7FFFFFFF; //
const ll nmos = 0x80000000; //-2147483648
const int inf = 0x3f3f3f3f;
const ll inff = 0x3f3f3f3f3f3f3f3f; //
const int mod = 1e9+;
const double esp = 1e-;
const double PI=acos(-1.0);
const double PHI=0.61803399; //黄金分割点
const double tPHI=0.38196601; template<typename T>
inline T read(T&x){
x=;int f=;char ch=getchar();
while (ch<''||ch>'') f|=(ch=='-'),ch=getchar();
while (ch>=''&&ch<='') x=x*+ch-'',ch=getchar();
return x=f?-x:x;
} /*-----------------------showtime----------------------*/
const int maxn = ;
int n,m,k,x;
int req[maxn][maxn],hav[maxn][maxn],fy[maxn][maxn][maxn];
int res[maxn]; const int maxl = 1e6+;
struct Edge{
int to,val,cost,nxt;
}gEdge[maxl]; int h[maxn],gPre[maxn];
int gPath[maxl],gDist[maxn];
bool in[maxn];
int gcount = ; bool spfa(int s,int t){
memset(gPre,-,sizeof(gPre));
memset(gDist, inf,sizeof(gDist));
memset(in, false, sizeof(in)); gDist[s] = ; in[s] = true;
queue<int>q;
q.push(s);
while(!q.empty()){
int u = q.front();
q.pop(); in[u] = false;
for(int e = h[u]; e!=-;e = gEdge[e].nxt){
int v = gEdge[e].to, w = gEdge[e].cost;
if(gEdge[e].val > && gDist[v] > gDist[u] + w){
gDist[v] = gDist[u] + gEdge[e].cost;
gPre[v] = u;
gPath[v] = e;
if(!in[v]){
q.push(v); in[v] = true;
}
}
}
}
return gPre[t] != -;
}
int MinCostFlow(int s,int t,int nd){
int cost = ,flow = ;
while(spfa(s,t)){
int f = inf;
for(int u=t; u!=s; u = gPre[u]){
if(gEdge[gPath[u]].val < f){
f = gEdge[gPath[u]].val;
}
}
flow +=f;
cost += gDist[t]*f;
for(int u=t; u!=s; u=gPre[u]){
gEdge[gPath[u]].val -= f;
gEdge[gPath[u]^].val += f;
}
}
if(flow != nd){
return -;
}
return cost;
} void addedge(int u,int v,int val,int cost){
gEdge[gcount].to = v;
gEdge[gcount].val = val;
gEdge[gcount].cost = cost;
gEdge[gcount].nxt = h[u];
h[u] = gcount++; gEdge[gcount].to = u;
gEdge[gcount].val = ;
gEdge[gcount].cost = -cost;
gEdge[gcount].nxt = h[v];
h[v] = gcount++;
}
int main(){
while(~scanf("%d%d%d", &n, &m, &k) && n+m+k){ memset(req,,sizeof(req));
memset(hav,,sizeof(hav));
memset(res, ,sizeof(res)); for(int i=; i<=n; i++)
for(int j=; j<=k; j++)
scanf("%d", &req[i][j]),res[j] += req[i][j]; for(int i=; i<=m; i++)
for(int j=; j<=k; j++)
{
scanf("%d", &hav[i][j]);
}
for(int i=; i<=k; i++)
for(int t = ; t<=n; t++)
for(int j=; j<=m; j++)
{
scanf("%d",&fy[i][j][t]);
}
int ans = ,kk = k,flag = ;
int s = ,t = n+m+;
for(int kind = ; kind <= k; kind ++){
memset(h,-,sizeof(h));
gcount = ;
for(int i=; i<=m; i++) addedge(s,i,hav[i][kind],);
for(int i=; i<=m; i++)
for(int j=; j<=n; j++){
addedge(i,m+j,inf,fy[kind][i][j]);
}
for(int i=; i<=n; i++){
addedge(m+i,t, req[i][kind],);
}
int tmp = MinCostFlow(s,t,res[kind]);
if(tmp == -) flag = ;
else ans += tmp;
}
if(flag == )puts("-1");
else
printf("%d\n", ans);
}
return ;
}
POJ 2516
POJ - 2516 Minimum Cost 每次要跑K次费用流的更多相关文章
- POJ 2516 Minimum Cost (网络流,最小费用流)
POJ 2516 Minimum Cost (网络流,最小费用流) Description Dearboy, a goods victualer, now comes to a big problem ...
- Poj 2516 Minimum Cost (最小花费最大流)
题目链接: Poj 2516 Minimum Cost 题目描述: 有n个商店,m个仓储,每个商店和仓库都有k种货物.嘛!现在n个商店要开始向m个仓库发出订单了,订单信息为当前商店对每种货物的需求 ...
- POJ 2516 Minimum Cost (最小费用最大流)
POJ 2516 Minimum Cost 链接:http://poj.org/problem?id=2516 题意:有M个仓库.N个商人.K种物品.先输入N,M.K.然后输入N行K个数,每一行代表一 ...
- POJ 2516 Minimum Cost (费用流)
题面 Dearboy, a goods victualer, now comes to a big problem, and he needs your help. In his sale area ...
- POJ 2516 Minimum Cost 最小费用流
题目: 给出n*kk的矩阵,格子a[i][k]表示第i个客户需要第k种货物a[i][k]单位. 给出m*kk的矩阵,格子b[j][k]表示第j个供应商可以提供第k种货物b[j][k]单位. 再给出k个 ...
- POJ 2516 Minimum Cost 最小费用流 难度:1
Minimum Cost Time Limit: 4000MS Memory Limit: 65536K Total Submissions: 13511 Accepted: 4628 Des ...
- POJ 2516 Minimum Cost(最小费用流)
Description Dearboy, a goods victualer, now comes to a big problem, and he needs your help. In his s ...
- POJ - 2516 Minimum Cost(最小费用最大流)
1.K种物品,M个供应商,N个收购商.每种物品从一个供应商运送到一个收购商有一个单位运费.每个收购商都需要K种物品中的若干.求满足所有收购商需求的前提下的最小运费. 2.K种物品拆开来,分别对每种物品 ...
- POJ 2516 Minimum Cost(拆点+KM完备匹配)
题目链接:http://poj.org/problem?id=2516 题目大意: 第一行是N,M,K 接下来N行:第i行有K个数字表示第i个卖场对K种商品的需求情况 接下来M行:第j行有K个数字表示 ...
随机推荐
- O2优化的实质
重点:如果使用多个-O选项(包含或不包含级别编号),则最后一个选项是有效的选项.------------ ------------ ------------例如:#pragma GCC optimiz ...
- 【iOS】手动抛出异常
之前没遇到过需要手动抛出异常的时候,这次见到了,记录一下.示例代码如下: /** 如果调用 [[BNRItemStore alloc] init],就提示应该使用 [BNRItemStore shar ...
- oracle的开窗函数
原创 select * from (select province, commodity, sum(price), ROW_NUMBER() OVER(PARTITION BY province o ...
- Extjs的使用总结笔记
一:Extjs自带验证 1.alpha //只能输入字母,无法输入其他(如数字,特殊符号等) 2.alphanum//只能输入字母和数字,无法输入其他 3.email//email验证,要求的格式是& ...
- 1. 源码分析---SOFARPC可扩展的机制SPI
这几天离职在家,正好没事可以疯狂的输出一下,本来想写DUBBO的源码解析的,但是发现写DUBBO源码的太多了,所以找一个写的不那么多的框架,所以就选中SOFARPC这个框架了. SOFARPC是蚂蚁金 ...
- C程序设计(第四版)课后习题完整版 谭浩强编著
//复习过程中,纯手打,持续更新,觉得好就点个赞吧. 第一章:程序设计和C语言 习题 1.什么是程序?什么是程序设计? 答:程序就是一组计算机能识别和执行的指令.程序设计是指从确定任务到得到结果,写出 ...
- 【JDK】JDK源码分析-CyclicBarrier
概述 CyclicBarrier 是并发包中的一个工具类,它的典型应用场景为:几个线程执行完任务后,执行另一个线程(回调函数,可选),然后继续下一轮,如此往复. 打个通俗的比方,可以把 CyclicB ...
- Qt基于sqlite数据库的管理小软件
闲来无事,写了一个基于sqlite的数据库管理小软件. 先上图 中心思想就是: 创建一个数据库 然后每一个分组对应一个数据表 然后遍历该数据表.将名字以treewidgetItem显示出来.添加删除实 ...
- SpringBoot 集成Jedis操作set
题外话: Redis是个有趣的东西,相信搞java的或多或少都会用到,面试时也总离不开问Redis,之前觉得redis只是用做缓存,飞快!也因为最初在封装底层的时候,使用Redisson,所以大部分都 ...
- openldap介绍和使用
openldap介绍和使用 为什么会有本文? 早期,公司是没有统一认证这个东西的,所以各自玩各自的.于是, confluence一个用户体系,gitlab一个用户体系,Jenkins一个用户体系等等, ...