[HAOI2010] 订货

题目描述：

大大大雾。

题目分析：

开始把仓库容量想成了购买+储存的容量上限，还傻敷敷的拆点（哇，什么脑回路…
建模很好搞，就不分析了

题目链接：

Luogu 2517
BZOJ 2424

Ac 代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
const int inf=0x7ffffff;
const int maxm=2e6+100;
int head[maxm],to[maxm<<1],net[maxm<<1],cost[maxm<<1],cap[maxm<<1];
int cnt=1;
int n,m,s;
int u[maxm],costx[maxm];
inline void addedge(int u,int v,int c1,int c2)
{
    cnt++;
    to[cnt]=v,cap[cnt]=c1,cost[cnt]=c2,net[cnt]=head[u],head[u]=cnt;
    cnt++;
    to[cnt]=u,cap[cnt]=0,cost[cnt]=-c2,net[cnt]=head[v],head[v]=cnt;
}
namespace MFMC{
    int pre[maxm],id[maxm],dis[maxm],flow[maxm];
    int mincost,maxflow;
    bool vis[maxm];
    std::queue <int> dl;
    inline bool SPFA(int s,int t)
    {
        memset(pre,-1,sizeof(pre));
        memset(dis,127/3,sizeof(dis));
        dis[s]=0,pre[s]=0,flow[s]=inf,vis[s]=1;
        dl.push(s);
        while(!dl.empty())
        {
            int now=dl.front();
            dl.pop();
            vis[now]=0;
            for(int i=head[now];i;i=net[i])
            if(cap[i]&&dis[to[i]]>dis[now]+cost[i])
            {
                dis[to[i]]=dis[now]+cost[i];
                pre[to[i]]=now;
                id[to[i]]=i;
                flow[to[i]]=std::min(flow[now],cap[i]);
                if(!vis[to[i]]) vis[to[i]]=1,dl.push(to[i]);
            }
        }
        return pre[t]!=-1;
    }
    inline void change_cap(int s,int t,int x)
    {
        int now=t;
        while(now!=s)
        {
            cap[id[now]]-=x,cap[id[now]^1]+=x;
            now=pre[now];
        }
    }
    inline int mfmc(int s,int t)
    {
        maxflow=mincost=0;
        while(SPFA(s,t))
        {
            maxflow+=flow[t];
            mincost+=flow[t]*dis[t];
            change_cap(s,t,flow[t]);
        }
        return mincost;
    }
} 
int main()
{
    scanf("%d%d%d",&n,&m,&s);
    int S=0,T=2e6+1;
    for(int i=1;i<=n;i++)
     scanf("%d",&u[i]);
    for(int i=1;i<=n;i++)
     scanf("%d",&costx[i]);
    for(int i=1;i<=n;i++)
    {
        addedge(S,i,inf,costx[i]),addedge(i,T,u[i],0);
        if(i!=n) addedge(i,i+1,s,m);
    } 
    printf("%d\n",MFMC::mfmc(S,T));
    return 0;
}