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;
}