nyoj 1251 山区修路（DP）

山区修路

时间限制： 1000 ms  |  内存限制： 65535 KB

难度： 2

描述

某山区的孩子们上学必须经过一条凹凸不平的土路，每当下雨天，孩子们非常艰难。现在村里走出来的Dr. Kong决定募捐资金重新修建着条路。由于资金有限，为了降低成本，对修好后的路面高度只能做到单调上升或单调下降。

为了便于修路，我们将整个土路分成了N段，每段路面的高度分别A1，A2,….,An。由于将每一段路垫高或挖低一个单位的花费成本相同，修路的总费用与路面的高低成正比。

现在Dr. Kong希望找到一个恰好含N个元素的不上升或不下降序列B1，B2,….,Bn，作为修过的路路段的高度。要求：

 

         | A1-B1| + | A2–B2| + ... + | An-Bn|------>最小

输入

第一行： K 表示有多少组测试数据。
接下来对每组测试数据：
第1行: N 表示整个土路分成了N段
第2~N+1行： A1 A2 ……AN 表示每段路面的高度
2≤k≤10 0≤Ai≤10000000 0≤N≤500 (i=1,…, N)
所有数据都是整数。 数据之间有一个空格。
数据保证| A1-B1|+| A2-B2|+ ... +| An-Bn|的最小值不会超过109

输出

对于每组测试数据，输出占一行：| A1-B1|+| A2-B2|+ ... +| An-Bn|的最小值。

样例输入

271 3 2 4 5 3 958 6 5 6 2

样例输出

31

来源

第七届河南省程序设计大赛

因为最终结果每一位一定是原500个数中的一位，所以我们枚举每一位上取500位中的哪一位更优。

代码：

 #include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <queue>
#define mem(p,k) memset(p,k,sizeof(p));
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define inf 0x6fffffff
#define LL long long
using namespace std;
struct node{
    int l,r,le;

}kk[100100],c;
//vector<int> b[100100];
char s[11000];
LL a[1100],dp[510][510],b[1100];
LL m(LL a,LL b){
    return a>b?b:a;
}
int main()
{
    int t,n;
    cin>>t;
    while(t--){
        cin>>n;
        LL minn=2147483647;
        for(int i=0;i<n;i++){
                cin>>a[i];b[i]=a[i];
        }
        sort(b,b+n);
        mem(dp,0);
        for(int i=0;i<n;i++)
            dp[0][i]=abs(b[i]-a[0]);
        for(int i=1;i<n;i++){
            LL minn=dp[i-1][0];
            for(int j=0;j<n;j++){
                minn=m(minn,dp[i-1][j]);
                dp[i][j]=minn+abs(a[i]-b[j]);
                //cout<<i<<j<<'='<<dp[i][j]<<endl;
            }
        }

        for(int i=0;i<n;i++){
                minn=m(minn,dp[n-1][i]);
        }
        reverse(b,b+n);
        mem(dp,0);
        for(int i=0;i<n;i++){
            dp[0][i]=abs(b[i]-a[0]);//cout<<dp[0][i]<<endl;
        }

        for(int i=1;i<n;i++){
            LL minn=dp[i-1][0];
            for(int j=0;j<n;j++){
                minn=m(minn,dp[i-1][j]);
                dp[i][j]=minn+abs(a[i]-b[j]);
                //cout<<i<<j<<'='<<dp[i][j]<<endl;
            }
        }

        for(int i=0;i<n;i++)minn=m(minn,dp[n-1][i]);
        cout<<minn<<endl;
    }

    return 0;
}