Buying Feed, 2010 Nov (单调队列优化DP)

时间:2022-08-31 17:04:27

约翰开车回家,又准备顺路买点饲料了(咦?为啥要说“又”字?)回家的路程一共有 E 公里,
这一路上会经过 K 家商店,第 i 家店里有 Fi 吨饲料,售价为每吨 Ci 元。约翰打算买 N 吨饲料,他
知道商家的库存是足够的,至少所有店的库存总和不会少于 N。除了购买饲料要钱,运送饲料也是
要花油钱的,约翰的卡车上如果装着 X 吨饲料,那么他行驶一公里会花掉 X 2 元,行驶 D 公里需要
D X 2 元。已知第 i 家店距约翰所在的起点有 Xi 公里,那么约翰在哪些商店买饲料运回家,才能做到
最省钱呢?

输入格式
• 第一行:三个整数 KE N, 1 K 10000 , 1 E 500 , 1 N 500
• 第二行到第 N + 1 行:第 i + 1 行有三个整数 XiFi Ci, 0 < Xi < E, 1 Fi 10000, 1
Ci 107

输出格式
• 单个整数:表示购买及运送饲料的最小费用

样例输入
2 5 3
3 1 2
4 1 2
1 1 1

样例输出
9

解释
在离家较近的两家商店里各购买一吨饲料,
则花在路上的钱是 1 + 4 = 5,花在店里的钱是
2 + 2 = 4

【分析】

  嗯,啊,还是好笨,想了挺久。

  先列DP,f[i][x]=min(f[j][k]+(x-k)^2*(d[i]-d[j])+(x-k)*c[i]) d[i][x]表示走到i,一共买了x个东西的最小费用。

  但是这样列的话很难降维,因为答案跟d[j]有关,所以可以用 计算未来费用的思想,就是买的时候直接算他运到终点了。

  f[i][x]=min(f[j][k]+(x-k)*c[i]+(x^2-k^2)*(s-d[i])) 这样就可以降维了。

  f[x]=min(f[k]+(x-k)*c[i]+(x^2-k^2)*(s-d[i])) i直接for,不过要注意一点是要用的是i之前算出的f而不能是i时计算出的f

  如果没有限制的话,这样的方程当然存一个最优解就好了,但是有限制,就要看限制的单调性,我们要x-k<=sm[i] 即 k>=x-sm[i]

  x按顺序枚举的话就有单调性了。

  啊,又是一道限制为主的单调队列ORZ、、、

  

代码如下:

 #include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
#define Maxn 510
#define Maxm 200010
#define LL long long struct node
{
LL d,sm,w;
}t[Maxn]; LL mymin(LL x,LL y) {return x<y?x:y;}
LL mymax(LL x,LL y) {return x>y?x:y;} bool cmp(node x,node y) {return x.d<y.d;} LL q[Maxm],st[Maxm],f[Maxm]; int main()
{
LL v,s,n;
scanf("%lld%lld%lld",&v,&s,&n);
for(LL i=;i<=n;i++) scanf("%lld%lld%lld",&t[i].d,&t[i].sm,&t[i].w);
sort(t+,t++n,cmp);
for(LL i=;i<=n;i++) t[i].d=s-t[i].d;
memset(f,,sizeof(f));
f[]=;
int ql,qr;
for(LL i=;i<=n;i++)
{
ql=qr=;q[qr]=;st[qr]=;
for(LL j=;j<=v;j++)
{
while(ql<qr&&(j-st[ql])>t[i].sm) ql++;
LL now=f[j];
f[j]=mymin(f[j],q[ql]+t[i].d*j*j+t[i].w*j);
while(now-j*j*t[i].d-t[i].w*j<=q[qr]&&qr>=ql) qr--;
q[++qr]=now-j*j*t[i].d-t[i].w*j;st[qr]=j;
}
}
printf("%lld\n",f[v]);
return ;
}

2016-10-20 09:14:21