【暑假】[深入动态规划]UVa 12170 Easy Climb

UVa 12170 Easy Climb

题目：

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=24844

思路：

 引别人一个题解琢磨一下：

 

  from：http://blog.csdn.net/glqac/article/details/45257659

代码：

#include<iostream>
#include<algorithm>
#define FOR(a,b,c) for(int a=(b);a<(c);a++)
using namespace std;

typedef long long LL;

const int maxn = 100 + 5;
const int maxx = maxn*maxn*2;
const LL INF = (1LL << 60);

LL h[maxn],x[maxx],dp[2][maxx];

int main() {
    int T; cin>>T;
    while(T--) {
        int n; LL d;
        cin>>n>>d;
        FOR(i,0,n) cin>>h[i];
        if(abs(h[0]-h[n-1])> (n-1)*d) { //input无解 
            cout<<"impossible
";
            continue;
        }
        
        int nx=0;
        FOR(i,0,n) 
          FOR(j,-(n-1),n) //-(n-1)..(n-1)
              x[nx++]=h[i]+j*d;
        sort(x,x+nx);  //从小到大单调增长 
        nx=unique(x,x+nx)-x; //去重且返回最终尾标  //x构造为修改后的所有可能取值 
        
        FOR(i,0,nx) {  //dp_init
            dp[0][i]=INF;                                                             
            if(x[i]==h[0]) dp[0][i]=0; //原本第i个就是h[0] //第0个数不能修改 
        } 
        
        //d[i][j]意味着已经修改i个数其中第i个数修改成为x[j]需要的最小费用 
        
        int t=0;  //滚动数组的空间优化 
        FOR(i,1,n) {
            int k=0;
            FOR(j,0,nx) {
                while(k<nx && x[k]<x[j]-d) k++;  //两者的取值必须不超过d 
                while(k+1<nx && x[k+1]<=x[j]+d && dp[t][k+1]<=dp[t][k])   //在 滑动窗口 中最小的  
                  k++;
                //刷表法 更新 
                if(dp[t][k]==INF) dp[t^1][j]=INF;
                else dp[t^1][j]=dp[t][k]+abs(x[j]-h[i]);
            }
            t^=1;
        }
        
        FOR(i,0,nx) if(x[i]==h[n-1])  //第N-1个数不能修改 
          cout<<dp[t][i]<<"
";
        
    }
    return 0;
}