题目大意:将长城看成一直线段,有n个需要修缮的点。对于第i个点用三个参数x,c,d表述。x表示位置,c表示立即修缮的费用,d表示单位时间增加的修缮费用。修缮时间忽略不计。机器人初始位置在x处,移动速度恒定为v,找到修缮所有点的最小费用
将所有点按照位置排序。由于修复的点一定是一段连续的区间,因此我们可以考虑用区间dp解决
f[i][j][k]表示修复i~j区间所有点最终停在k点时(k=0表示在i,k=1表示在j)所需最小费用。但是这样做的问题是,不知道当前经过的时间,就无法计算后续产生的费用。
我们可以事先把还没有发生但是肯定会产生的费用计算下来,然后把t归零后进行状态转移。
于是f[i][j][k]表示修复完i~j已经发生的总费用和未来一定会产生的费用总和。这里只有两个决策:往左走和往右走。状态转移方程在代码中有体现,这里就不再写出。
注意:在输入每组数据时while(scanf("%d%lf%lf",&n,&v,&x0))莫名其妙会超时,全改成cin就能过了。特别神奇……
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1005;
const double inf=1e9+5;
struct Point{
double x,c,d;
bool operator <(const Point &h)const{
return x<h.x;
}
}s[maxn];
int n,bp,tur,cnt[maxn][maxn][2];
double v,x0,sum0;
double f[maxn][maxn][2],sum[maxn];
void clear(){
memset(s,0,sizeof(s));
memset(sum,0,sizeof(sum));
sum0=0;
}
void spework1(){
double ans=sum0,cur=x0;
for(int i=1;i<=n;i++){
ans+=(s[i].x-cur)/v*(sum[n]-sum[i-1]);
cur=s[i].x;
}
printf("%d\n",(int)ans);
}
void spework2(){
double ans=sum0,cur=x0;
for(int i=n;i;i--){
ans+=(cur-s[i].x)/v*(sum[i]-sum[0]);
cur=s[i].x;
}
printf("%d\n",(int)ans);
}
double getf(int i,int j,int k){
if(cnt[i][j][k]==tur)return f[i][j][k];
if(i==j&&i!=bp&&i!=bp+1)return inf;
cnt[i][j][k]=tur;
f[i][j][k]=inf;
if(k==1){
f[i][j][k]=min(f[i][j][k],getf(i,j-1,1)+(s[j].x-s[j-1].x)/v*(sum[i-1]+sum[n]-sum[j-1]));
f[i][j][k]=min(f[i][j][k],getf(i,j-1,0)+(s[j].x-s[i].x)/v*(sum[i-1]+sum[n]-sum[j-1]));
}
else {
f[i][j][k]=min(f[i][j][k],getf(i+1,j,0)+(s[i+1].x-s[i].x)/v*(sum[i]+sum[n]-sum[j]));
f[i][j][k]=min(f[i][j][k],getf(i+1,j,1)+(s[j].x-s[i].x)/v*(sum[i]+sum[n]-sum[j]));
}
return f[i][j][k];
}
int main(){
while(cin>>n>>v>>x0){
if(n==0&&v==0&&x0==0)break;
tur++;
clear();
for(int i=1;i<=n;i++)scanf("%lf%lf%lf",&s[i].x,&s[i].c,&s[i].d);
sort(s+1,s+1+n);
for(int i=1;i<=n;i++){
sum[i]=sum[i-1]+s[i].d;
sum0+=s[i].c;
}
if(x0<=s[1].x){
spework1();
continue;
}
if(x0>=s[n].x){
spework2();
continue;
}
for(int i=1;i<=n;i++)
if(s[i].x<=x0&&s[i+1].x>=x0){
f[i][i][1]=f[i][i][0]=(x0-s[i].x)/v*sum[n]+sum0;
f[i+1][i+1][1]=f[i+1][i+1][0]=(s[i+1].x-x0)/v*sum[n]+sum0;
cnt[i][i][0]=cnt[i][i][1]=cnt[i+1][i+1][0]=cnt[i+1][i+1][1]=tur;
bp=i;
}
printf("%d\n",min((int)getf(1,n,0),(int)getf(1,n,1)));
}
}