[BZOJ2424/Luogu2517][HAOI2010]订货

题目链接：

BZOJ2424

Luogu2517

费用流大水题。

源点与每月连边，容量\(+\infty\)（可以生产无数件商品），费用\(d_i\)（商品单位成本）。

每月与下一月连边，容量\(S\)（仓库容量），费用\(m\)（单位物品存储费用）。

每月与汇点连边，容量\(U_i\)（需求量），费用\(0\)。

#include <queue>
#include <cstdio>
#include <cstring>

inline int Min(int a,int b){return a<b?a:b;}

int n,m,s,St,Ed,MaxFlow,MinCost;
int Head[60],Next[405],To[405],Val[405],Cos[405],En=1;
int Dis[60],Pre[60],Ref[60],Inq[60];

inline void Add(int x,int y,int z,int w)
{
	Next[++En]=Head[x],To[Head[x]=En]=y,Val[En]=z,Cos[En]=+w;
	Next[++En]=Head[y],To[Head[y]=En]=x,Val[En]=0,Cos[En]=-w;
}

bool SPFA()
{
	std::queue<int> q;
	memset(Dis,0x3f,sizeof Dis);
	q.push(St),Dis[St]=0,Ref[St]=1<<30;
	for(int x,y;!q.empty();q.pop(),Inq[x]=0)
		for(int i=Head[x=q.front()];i;i=Next[i])
			if(Val[i]&&Dis[y=To[i]]>Dis[x]+Cos[i])
			{
				Dis[y]=Dis[x]+Cos[Pre[y]=i];
				Ref[y]=Min(Ref[x],Val[i]);
				if(!Inq[y])Inq[y]=1,q.push(y);
			}
	if(Dis[Ed]==0x3f3f3f3f)return false;
	MaxFlow+=Ref[Ed],MinCost+=Ref[Ed]*Dis[Ed];
	for(int x=Ed,i;x!=St;x=To[i^1])
		Val[i=Pre[x]]-=Ref[Ed],Val[i^1]+=Ref[Ed];
	return true;
}

int main()
{
	scanf("%d%d%d",&n,&m,&s),St=n+1,Ed=St+1;
	for(int i=1,u;i<=n;++i)
	{
		scanf("%d",&u);
		Add(i,Ed,u,0);
		if(i!=n)Add(i,i+1,s,m);
	}
	for(int i=1,d;i<=n;++i)
	{
		scanf("%d",&d);
		Add(St,i,0x3f3f3f3f,d);
	}
	while(SPFA());
	printf("%d\n",MinCost);
	return 0;
}