【题意】

Description

"Help Jimmy" 是在下图所示的场景上完成的游戏。

场景中包括多个长度和高度各不相同的平台。地面是最低的平台,高度为零,长度无限。

Jimmy老鼠在时刻0从高于所有平台的某处开始下落,它的下落速度始终为1米/秒。当Jimmy落到某个平台上时,游戏者选择让它向左还是向右跑,它跑动的速度也是1米/秒。当Jimmy跑到平台的边缘时,开始继续下落。Jimmy每次下落的高度不能超过MAX米,不然就会摔死,游戏也会结束。

设计一个程序,计算Jimmy到底地面时可能的最早时间。

Input

第一行是测试数据的组数t(0 <= t <= 20)。每组测试数据的第一行是四个整数N,X,Y,MAX,用空格分隔。N是平台的数目(不包括地面),X和Y是Jimmy开始下落的位置的横竖坐标,MAX是一次下落的最大高度。接下来的N行每行描述一个平台,包括三个整数,X1[i],X2[i]和H[i]。H[i]表示平台的高度,X1[i]和X2[i]表示平台左右端点的横坐标。1 <= N <= 1000,-20000 <= X, X1[i], X2[i] <= 20000,0 < H[i] < Y <= 20000(i = 1..N)。所有坐标的单位都是米。

Jimmy的大小和平台的厚度均忽略不计。如果Jimmy恰好落在某个平台的边缘,被视为落在平台上。所有的平台均不重叠或相连。测试数据保证问题一定有解。

Output

对输入的每组测试数据,输出一个整数,Jimmy到底地面时可能的最早时间。

Sample Input

1
3 8 17 20
0 10 8
0 10 13
4 14 3

Sample Output

23

【分析】

当Jimmy落在一个平台上后有两种选择(向左走或向右走),而Jimmy走到平台左边和右边的时间很容易计算,如果我们得到了以平台左边为起点及以平台右边为起点到地面的最短时间,那么选择往左走还是往右走就很容易了。这样,原问题就分解为两个子问题这两个子问题和原问题的形式是一致的了,也就找到了“状态”dp[i][j], j = 0, 1 (dp[i][0]表示以i号平台左边为起点到地面的最短时间,dp[i][1]]表示以i号平台右边为起点到地面的最短时间),而“状态转移方程”如下:

dp[i][0] = H[i] - H[m] + Min (dp[m][0] + X1[i] - X1[m], dp[m][1] + X2[m] - X1[i]);  m为i左边下面的平台的编号

dp[i][1] = H[i] - H[m] + Min (dp[m][0] + X2[i] - X1[m], dp[m][1] + X2[m] - X2[i]);  m为i右边下面的平台的编号

【AC代码】
#include <stdio.h>
#include <ctype.h>
#include <string.h>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 1010;
const int inf = 1<<30;
struct node{
   int x1,x2,h;
   node(){}
   node(int x1,int x2,int h):x1(x1),x2(x2),h(h){}
   bool operator<(const node &rhs)const
   {
       return h<rhs.h;
   }
}pla[maxn];
int dp[maxn][2];
int N,X,Y,MAX;
void work_left(int i)
{
    int k=i-1;
    while(k>0&&pla[i].h-pla[k].h<=MAX)
    {
        if(pla[i].x1>=pla[k].x1&&pla[k].x2>=pla[i].x1)
        {
            dp[i][0] = pla[i].h-pla[k].h+min(pla[i].x1-pla[k].x1+dp[k][0],pla[k].x2-pla[i].x1+dp[k][1]);
            return ;
        }
        else
        {
            --k;
        }
    }
    if(pla[i].h-pla[k].h>MAX)dp[i][0]=inf;
    else dp[i][0] = pla[i].h;
}
void work_right(int i)
{
    int k=i-1;
    while(k>0&&pla[i].h-pla[k].h<=MAX)
    {
        if(pla[i].x2<=pla[k].x2&&pla[i].x2>=pla[k].x1)
        {
            dp[i][1] = pla[i].h-pla[k].h+min(pla[i].x2-pla[k].x1+dp[k][0],pla[k].x2-pla[i].x2+dp[k][1]);
            return ;
        }
        else
        {
            --k;
        }
    }
    if(pla[i].h-pla[k].h>MAX)dp[i][1]=inf;
    else dp[i][1] = pla[i].h;
}
int work()
{
    for(int i=1;i<=N+1;i++)
    {
        work_left(i);
        work_right(i);
    }
    return min(dp[N+1][0],dp[N+1][1]);
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d%d",&N,&X,&Y,&MAX);
        for(int i=1;i<=N;i++)
        {
            scanf("%d%d%d",&pla[i].x1,&pla[i].x2,&pla[i].h);
        }
        pla[0]=node(-20000,20000,0);
        pla[N+1]=node(X,X,Y);
        sort(pla,pla+N+2);
        int ans=work();
        printf("%d\n",ans);
    }
    return 0;
}