[网络流] bzoj2424: [HAOI2010]订货

bzoj2424: [HAOI2010]订货:http://www.lydsy.com/JudgeOnline/problem.php?id=2424

一拿到题认定为是DP的我被网上铺天盖地的费用流题解吓到了
st 连向 每一天 费用为单价 容量为inf
每一天 连向 ed 费用为0 容量为当天需要的数量
相邻两天之间 费用为m 容量为s
因为spfa走出来的是一条路径路径路径!!
所以可以查最小流量(也就是当天库存量)
这样建边那么查最小流量的时候要么是s 要么是当天需要的数量
更新容量的时候顺便搞一搞ans就可以了
表示不会费用流
来自旁边神犇的嘲笑:我用的DP哦 证明单调性

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
struct node
{
    int x,y,c,cost,next,other;
}a[510];
int last[60],len;
int list[1100000],f[60];
int u[60],d[60];
int st,ed,ans=0;
bool v[60];
int pre[60]; 
void build(int x,int y,int c,int cost)
{
    int k1,k2;
    len++;k1=len;
    a[len].x=x;a[len].y=y;a[len].c=c;a[len].cost=cost;a[len].next=last[x];last[x]=len;
    len++;k2=len;
    a[len].x=y;a[len].y=x;a[len].c=0;a[len].cost=-cost;a[len].next=last[y];last[y]=len;
    a[k1].other=k2;a[k2].other=k1;
}
bool spfa()
{
    memset(f,63,sizeof(f));
    memset(v,false,sizeof(v));
    memset(pre,0,sizeof(pre));
    int head=1,tail=1;
    list[1]=st; f[st]=0; v[st]=true;
    while(head<=tail)
    {
        int x=list[head];
        for (int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if (a[k].c>0&&f[y]>f[x]+a[k].cost)
            {
                f[y]=f[x]+a[k].cost;
                pre[y]=k;
                if (v[y]==false)
                {
                    v[y]=true;
                    list[++tail]=y;
                }
            }
        }
        head++;
        v[x]=false; 
    }
    if (f[ed]>=999999999) return false;
    else 
    {
        //spfa出来的是一条路径路径路径 
        int mc=999999999;//路径中的最小流量 即这条路径能走多少 
        int x=ed;
        while (x!=st)
        {
            int k=pre[x];
            mc=min(mc,a[k].c);
            x=a[k].x;
        }
        x=ed;
        while (x!=st)
        {
            int k=pre[x];
            a[k].c-=mc;
            a[a[k].other].c+=mc;
            ans+=a[k].cost*mc;
            x=a[k].x;
        }
        return true;
    }

}
int main()
{
    //st 连向 每一天 费用为单价 容量 inf
    //每一天 连向 ed 费用为0 容量 当天需要的数量
    //相邻两天之间 费用为m 容量 s 
    int n,m,s;
    scanf("%d%d%d",&n,&m,&s);
    for (int i=1;i<=n;i++) scanf("%d",&u[i]);
    for (int i=1;i<=n;i++) scanf("%d",&d[i]);
    st=n+1;ed=st+1;
    for (int i=1;i<=n;i++)
    {
        build(st,i,999999999,d[i]);
        build(i,ed,u[i],0);
        if (i!=1) build(i-1,i,s,m);
    }
    while (spfa()){}
    printf("%d\n",ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值