UVA1336 Fixing the great wall

题目大意:将长城看成一直线段，有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)));
	}
}