河南省ACM nyoj 1367 物流配送

题目描述:

物流配送是物流活动中一种非单一的业务形式，它与物品流动、资金流动紧密结合。备货是配送的准备工作或基础工作，备货工作包括筹集货源、订货或购货、集货、进货及有关的质量检查、结算、交接等。配送的优势之一，就是可以集中用户的需求进行一定规模的备货。备货是决定配送成败的初期工作，如果备货成本太高，会大大降低配送的效益。配送中的储存有储备及暂存两种形态。配送储备是按一定时期的配送经营要求，形成的对配送的资源保证。这种类型的储备数量较大，储备结构也较完善，视货源及到货情况，可以有计划地确定周转储备及保险储备结构及数量。

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

基本就是费用流裸题

建立一个虚拟源点，和所有点连边，容量是a[i]，费用是0

建立一个虚拟汇点，和所有点连边，容量为b[i]，费用是0

注意longlong

 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=5e3+5,M=1e5+5;
const int INF=0x3f3f3f3f;
int n,tot=1,lnk[N],cur[N],ter[M],nxt[M],cap[M],cost[M],dis[N];
ll ret;
bool vis[N];

void add(int u,int v,int w,int c) {
    ter[++tot]=v,nxt[tot]=lnk[u],lnk[u]=tot,cap[tot]=w,cost[tot]=c;
}
void addedge(int u,int v,int w,int c) {
    add(u,v,w,c),add(v,u,0,-c);
}
bool spfa(int s,int t) {
    memset(dis,0x3f,sizeof(dis));
    memcpy(cur,lnk,sizeof(lnk));
    std::queue<int> q;
    q.push(s),dis[s]=0,vis[s]=1;
    while(!q.empty()) {
        int u=q.front(); q.pop(),vis[u]=0;
        for(int i=lnk[u];i;i=nxt[i]) {
            int v=ter[i];
            if(cap[i]&&dis[v]>dis[u]+cost[i]) {
                dis[v]=dis[u]+cost[i];
                if(!vis[v]) q.push(v),vis[v]=1;
            }
        }
    }
    return dis[t]!=INF;
}
int dfs(int u,int t,int flow) {
    if(u==t) return flow;
    vis[u]=1;
    int ans=0;
    for(int &i=cur[u];i&&ans<flow;i=nxt[i]) {
        int v=ter[i];
        if(!vis[v]&&cap[i]&&dis[v]==dis[u]+cost[i]) {
            int x=dfs(v,t,std::min(cap[i],flow-ans));
            if(x) ret+=x*cost[i],cap[i]-=x,cap[i^1]+=x,ans+=x;
        }
    }
    vis[u]=0;
    return ans;
}
int mcmf(int s,int t) {
    int ans=0;
    while(spfa(s,t)) {
        int x;
        while((x=dfs(s,t,INF))) ans+=x;
    }
    return ans;
}
int main()
{
    int s,t;
    scanf("%d",&n);
    s=0,t=n+1;     //虚拟源点，虚拟汇点
    for(int i=1;i<=n;i++)
    {
    	int x;scanf("%d",&x);
    	addedge(s,i,x,0);
    }
    for(int i=1;i<=n;i++)
    {
    	int x;scanf("%d",&x);
    	addedge(i,t,x,0);
    }
    for(int i=1;i<n;i++)
    {
    	int u,v,c;scanf("%d%d%d",&u,&v,&c);
    	addedge(u,v,INF,c);
    	addedge(v,u,INF,c);
    }
    ll ans=mcmf(s,t);
    printf("%lld\n",ret);
    return 0;
}