【bzoj2424】[HAOI2010]订货

Description

某公司估计市场在第i个月对某产品的需求量为Ui，已知在第i月该产品的订货单价为di，上个月月底未销完的单位产品要付存贮费用m，假定第一月月初的库存量为零，第n月月底的库存量也为零，问如何安排这n个月订购计划，才能使成本最低？每月月初订购，订购后产品立即到货，进库并供应市场，于当月被售掉则不必付存贮费。假设仓库容量为S。

Input

第1行：n, m, S (0<=n<=50, 0<=m<=10, 0<=S<=10000)

第2行：U1 , U2 , … , Ui , … , Un (0<=Ui<=10000)

第3行：d1 , d2 , …, di , … , dn (0<=di<=100)

Output

只有1行，一个整数，代表最低成本

Sample Input

3 1 1000
2 4 8
1 2 4
Sample Output

34
题解
建图方法：把每个月抽象成一个点
S向每个点连边，费用为当月购进价，容量无限大；每个点向T连边，费用为0，容量为当月卖出量
每个点向后一个月的点连边，费用为仓库储存费用m，容量为仓库的容量s

代码（第一道费用流）

#include<bits/stdc++.h>
#define N 500005
#define ll long long
#define inf 0x7fffffff
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,s,T;
int tot=1,ans;
int from[2005],q[2005],dis[2005],head[2005];
bool inq[2005];
struct data{int from,to,next,v,c;}e[1000001];
void ins(int u,int v,int w,int c)
{
    tot++;
    e[tot].to=v;e[tot].from=u;
    e[tot].v=w;e[tot].c=c;
    e[tot].next=head[u];head[u]=tot;
}
void insert(int u,int v,int w,int c){ins(u,v,w,c);ins(v,u,0,-c);}
bool spfa()
{
    for (int i=1;i<=n+1;i++) dis[i]=inf;
    int t=0,w=1,i,now;
    dis[0]=q[0]=0;inq[0]=1;
    while (t!=w)
    {
        now=q[t];t++;if(t==2001)t=0;
        for(int i=head[now];i;i=e[i].next)
        {
            if(e[i].v&&dis[e[i].to]>dis[now]+e[i].c)
            {
                from[e[i].to]=i;
                dis[e[i].to]=dis[now]+e[i].c;
                if(!inq[e[i].to])
                {
                    inq[e[i].to]=1;
                    q[w++]=e[i].to;
                    if(w==2001)w=0;
                }
            }
        }
        inq[now]=0;
    }
    if (dis[T]==inf) return 0;return 1;
}
void mcf()
{
    int i,x=inf;
    i=from[T];
    while(i)
    {
        x=min(e[i].v,x);
        i=from[e[i].from];
    }
    i=from[T];
    while(i)
    {
        e[i].v-=x;
        e[i^1].v+=x;
        ans+=x*e[i].c;
        i=from[e[i].from];
    }
}
int main()
{
    n=read();m=read();s=read();T=n+1;
    for (int i=1;i<=n;i++)
    {
        int x=read();
        insert(i,T,x,0);
    }
    for (int i=1;i<=n;i++)
    {
        int x=read();
        insert(0,i,inf,x);
    }
    for (int i=1;i<n;i++) insert(i,i+1,s,m);
    while (spfa())mcf();
    printf("%d",ans);
    return 0;
}