poj 1661 左右分情况 DP

题目地址：http://poj.org/problem?id=1661

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

13 8 17 200 10 80 10 134 14 3

Sample Output

23

Source

POJ Monthly--2004.05.15 CEOI 2000

当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右边下面的平台的编号

//china no.1#pragma comment(linker, "/STACK:1024000000,1024000000")#include <vector>#include <iostream>#include <string>#include <map>#include <stack>#include <cstring>#include <queue>#include <list>#include <stdio.h>#include <set>#include <algorithm>#include <cstdlib>#include <cmath>#include <iomanip>#include <cctype>#include <sstream>#include <functional>#include <stdlib.h>#include <time.h>#include <bitset>using namespace std;#define pi acos(-1)#define PI acos(-1)#define endl '\n'#define srand() srand(time(0));#define me(x,y) memset(x,y,sizeof(x));#define foreach(it,a) for(__typeof((a).begin()) it=(a).begin();it!=(a).end();it++)#define close() ios::sync_with_stdio(0); cin.tie(0);#define FOR(x,n,i) for(int i=x;i<=n;i++)#define FOr(x,n,i) for(int i=x;i<n;i++)#define W while#define sgn(x) ((x) < 0 ? -1 : (x) > 0)#define bug printf("***********\n");#define db double#define ll long longtypedef long long LL;const int INF=0x3f3f3f3f;const LL LINF=0x3f3f3f3f3f3f3f3fLL;const int dx[]={-1,0,1,0,1,-1,-1,1};const int dy[]={0,1,0,-1,-1,1,-1,1};const int maxn=1e3+10;const int maxx=1e6+100;const double EPS=1e-8;const double eps=1e-8;const int mod=1e9+7;template<class T>inline T min(T a,T b,T c) { return min(min(a,b),c);}template<class T>inline T max(T a,T b,T c) { return max(max(a,b),c);}template<class T>inline T min(T a,T b,T c,T d) { return min(min(a,b),min(c,d));}template<class T>inline T max(T a,T b,T c,T d) { return max(max(a,b),max(c,d));}template<class T>inline bool scan_d(T &ret){char c;int sgn;if (c = getchar(), c == EOF){return 0;}while (c != '-' && (c < '0' || c > '9')){c = getchar();}sgn = (c == '-') ? -1 : 1;ret = (c == '-') ? 0 : (c - '0');while (c = getchar(), c >= '0' && c <= '9'){ret = ret * 10 + (c - '0');}ret *= sgn;return 1;}inline bool scan_lf(double &num){char in;double Dec=0.1;bool IsN=false,IsD=false;in=getchar();if(in==EOF) return false;while(in!='-'&&in!='.'&&(in<'0'||in>'9'))in=getchar();if(in=='-'){IsN=true;num=0;}else if(in=='.'){IsD=true;num=0;}else num=in-'0';if(!IsD){while(in=getchar(),in>='0'&&in<='9'){num*=10;num+=in-'0';}}if(in!='.'){if(IsN) num=-num;return true;}else{while(in=getchar(),in>='0'&&in<='9'){num+=Dec*(in-'0');Dec*=0.1;}}if(IsN) num=-num;return true;}void Out(LL a){if(a < 0) { putchar('-'); a = -a; }if(a >= 10) Out(a / 10);putchar(a % 10 + '0');}void print(LL a){ Out(a),puts("");}//freopen( "in.txt" , "r" , stdin );//freopen( "data.txt" , "w" , stdout );//cerr << "run time is " << clock() << endl;struct node{    int h,x1,x2;}Q[maxn];int t,n,x,y,maxc;bool cmp(node a,node b){    return a.h<b.h;}int dp[maxn][3];void Left(int i){    //计算从平台i左边到地面的最短时间    int k = i - 1;    while (k > 0 && Q[i].h - Q[k].h <= maxc)    {        //如果平台i左边下面有平台，且两者相距不超过MAX        if (Q[i].x1 >= Q[k].x1 && Q[i].x1 <= Q[k].x2)        {            dp[i][0] = Q[i].h - Q[k].h +min(Q[i].x1 - Q[k].x1 + dp[k][0],            Q[k].x2 - Q[i].x1 + dp[k][1]);            return;        }        else            --k;    }    //如果平台i左边下面没有平台，或者两者相距超过了MAX    if (Q[i].h - Q[k].h > maxc)        dp[i][0] = INF;    else        dp[i][0] = Q[i].h;}void Right(int i){   //计算从平台i右边到地面的最短时间    int k = i - 1;    while (k > 0 && Q[i].h - Q[k].h <= maxc)    {        //如果平台i右边下面有平台，且两者相距不超过MAX        if (Q[i].x2 >= Q[k].x1 && Q[i].x2 <= Q[k].x2)        {            dp[i][1] = Q[i].h - Q[k].h +min (Q[i].x2 - Q[k].x1 + dp[k][0],            Q[k].x2 - Q[i].x2 + dp[k][1]);            return;        }        else            --k;    }    //如果平台i右边下面没有平台，或者两者相距超过了MAX    if (Q[i].h - Q[k].h > maxc)        dp[i][1] = INF;    else        dp[i][1] = Q[i].h;}void solve(){    FOR(1,n+1,i)    {        Left(i);Right(i);    }    print(min(dp[n+1][1],dp[n+1][0]));}int main(){    scan_d(t);    W(t--)    {        cin>>n>>x>>y>>maxc;        FOR(1,n,i)            cin>>Q[i].x1>>Q[i].x2>>Q[i].h;        Q[0].h=0;Q[0].x1=-20000;Q[0].x2=20000;        Q[n+1].h=y;        Q[n+1].x1=Q[n+1].x2=x;        sort(Q,Q+n+2,cmp);        solve();    }}