第十一届河南省赛 E 1367-物流配送 //最小费用流(dijkstar实现）

http://nyoj.top/problem/1367

第十一届河南省赛 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;
}