第十一届河南省赛 E 1367-物流配送 //最小费用流(dijkstar实现)
1367-物流配送
- 内存限制:128MB 时间限制:8000ms 特判: No
- 通过数:14 提交数:32 难度:4
题目描述:
物流配送是物流活动中一种非单一的业务形式,它与物品流动、资金流动紧密结合。备货是配送的准备工作或基础工作,备货工作包括筹集货源、订货或购货、集货、进货及有关的质量检查、结算、交接等。配送的优势之一,就是可以集中用户的需求进行一定规模的备货。备货是决定配送成败的初期工作,如果备货成本太高,会大大降低配送的效益。配送中的储存有储备及暂存两种形态。配送储备是按一定时期的配送经营要求,形成的对配送的资源保证。这种类型的储备数量较大,储备结构也较完善,视货源及到货情况,可以有计划地确定周转储备及保险储备结构及数量。
Dr. Kong 所在的研究团队准备为Hai-E集团开发一个物流配送管理系统。已知Hai-E集团已经在全国各地建立了n个货物仓库基地,任意两个基地的货物可以相互调配。现在需要根据用户订货要求,来重新调配每个基地的货物数量。为了节流开源,希望对整个物流配送体系实行统一的货物管理和调度,能够提供一个全面完善的物流仓储配送解决方案,以减少物流配送过程中成本、人力、时间。
输入描述:
第一行: n (1 ≤ n ≤ 1000) 第2行: a1 a2 …… an 表示n个基地当前的物品数量 (0≤ ai ≤ 106 ) 第3行: b1 b2 …… bn 表示调配后,每个基地i应不少于bi个物品 (0≤ bi ≤ 106) 接下来n-1行,每行三个整数: i j k 表示从第i基地调配一个物品到第j基地需要花费为k,或 从第j基地调配一个物品到第i基地需要花费为k。(0≤ k ≤ 10^6)
输出描述:
输出配送后的最小费用。 已知: a1+a2+…+an >=b1+b2+…+bn
样例输入:
复制
6 0 1 2 2 0 0 0 0 1 1 1 1 1 2 2 1 3 5 1 4 1 2 5 5 2 6 1
样例输出:
9
最小费用流板子题,直接跑(套)板子
将所有点连到一个超级源点(下标为0),加边的时候,容量是a[i],费用是0;将所有点再连到一个超级汇点(下标为n),加边的时候,容量是b[i],费用也是0;
双向边,加边的时候直接可以加双向,不会有其他的bug。
这样就是求从源点到汇点的流量为∑(所有的b[i])的最小费用流。就a过啦。
只会套板子
代码已经define int =long long ,中间会炸。
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int max_n=1005;
struct no{int to,cap,cost,rev;};
struct qno{int v,d;}; //队列节点
bool operator<(const qno &x,const qno &y) {return x.d>y.d;} //队列cmp重写
int n; //定点数
vector<no>g[max_n];
int dis[max_n],h[max_n]; //最短距离,顶点的势
int prevv[max_n],preve[max_n]; //在vector下记录路径
int a[max_n],b[max_n];
int cob=0; //∑b[i]
void addedge(int from,int to,int cap,int cost){
g[from].push_back((no){to,cap,cost,g[to].size()});
g[to].push_back((no){from,0,-cost,g[from].size()-1});
}
int min_cost_flow(int s,int t,int f){
int res=0; //answer
for(int i=0;i<=n;i++) h[i]=0; //初始化h
while(f>0){
priority_queue<qno>q;
for(int i=0;i<=n;i++) dis[i]=INT_MAX/2;
dis[s]=0;
q.push((qno){s,0});
while(!q.empty()){
qno now=q.top();q.pop();
int v=now.v;
if(dis[v]<now.d) continue;
for(int i=0;i<g[v].size();i++){
no &e=g[v][i];
if(e.cap>0&&dis[e.to]>dis[v]+e.cost+h[v]-h[e.to]){ //加入啦势
dis[e.to]=dis[v]+e.cost+h[v]-h[e.to];
prevv[e.to]=v;
preve[e.to]=i;
q.push((qno){e.to,dis[e.to]});
}
}
}
if(dis[t]==INT_MAX/2) return -1;
for(int v=0;v<=n;v++) h[v]+=dis[v];
int d=f; //找增广路
for(int v=t;v!=s;v=prevv[v]) d=min(d,g[prevv[v]][preve[v]].cap);
f-=d;res+=d*h[t];
for(int v=t;v!=s;v=prevv[v]){
no &e=g[prevv[v]][preve[v]];
e.cap-=d;
g[v][e.rev].cap+=d;
}
}
return res;
}
void init(){
cin>>n;
for(int i=1;i<=n;i++) cin>>a[i];
for(int i=1;i<=n;i++) cin>>b[i],cob+=b[i];
for(int i=1;i<n;i++) {
int u,v,c;
cin>>u>>v>>c;
addedge(u,v,INT_MAX,c);
addedge(v,u,INT_MAX,c);
}
//连接超级源点
for(int i=1;i<=n;i++) addedge(0,i,a[i],0);
for(int i=1;i<=n;i++) addedge(i,n+1,b[i],0);
n++;
}
signed main(){
ios::sync_with_stdio(false);cin.tie(0);
init();
cout<<min_cost_flow(0,n,cob);
return 0;
}