【DP】[Usaco2008 Feb]Making the Grade 路面修整

【问题描述】

大老板Mr. Yao打算好好修一下公司门口的那条凹凸不平的路。按照Mr. Yao的设想，修好后的路面高度应当单调上升或单调下降，也就是说，高度上升与高度下降的路段不能同时出现在修好的路中。

    整条路被分成了N段，N个整数A_1, ... , A_N (1 <= N <= 2,000)依次描述了每一段路的高度(0 <= A_i <= 1,000,000,000)。Mr. Yao希望找到一个恰好含N个元素的不上升或不下降序列B_1, ... , B_N，作为修过的路中每个路段的高度。由于将每一段路垫高或挖低一个单位的花费相同，修路的总支出可以表示为：

         |A_1 - B_1| + |A_2 - B_2| + ... + |A_N - B_N|

    请你计算一下，Mr. Yao在这项工程上的最小支出是多少。Mr. Yao向你保证，这个支出不会超过2^31-1。

【输入数据】

* 第1行: 输入1个整数：N

* 第2..N+1行: 第i+1行为1个整数：A_i

【输入样例】

7

1

3

2

4

5

3

9

【输出数据】

* 第1行: 输出1个正整数，表示把路修成高度不上升或高度不下降的最小花费

【输出样例】

3

【样例说明】

    将第一个高度为3的路段的高度减少为2，将第二个高度为3的路段的高度增加到5，总花费为|2-3|+|5-3| = 3，并且各路段的高度为一个不下降序列1,2,2,4,5,5,9。、

【分析】

    由于要将路面修整为单调上升的数列，因此我们可以知道，修整后的路面一定是原路面的某个值。

创建新数组c等于原数组，然后将c排序，f[i][j]表示原数组前i个数变为单调上升,且i是c中的第j个数。

f[i][j]=min(f[i-1][k]+abs(a[i]-c[j]));

但是还可以优化,用一个g[i][j]表示f[i][k]中的最优值

g[i][j]=min(g[i][j-1],f[i][j])

f[i][j]=g[i-1][j]+abs(a[i]-a[j]);

边界f[0][i]=g[0][i]=0;

g[i][0]=maxlongint;

结果为g[n][n]；

对于单调上升与递减，求出某个后，将原数组倒过来在dp一次即可。

【代码】

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
int n,a[2010],d[2010],i,j,len,k,re[2010][2],ans=0,ans1=0;
int main()
{
    freopen("grading.in","r",stdin);
    freopen("grading.out","w",stdout);
	cin>>n;
	for ( i=1;i<=n;i++) scanf("%d",&a[i]); 
    d[0]=-1;  
    d[1]=a[0];
    len=0;
	 for ( i=1;i<=n;i++) 
    {
           j=lower_bound(d+1,d+len+1,a[i])-d;
            d[j]=a[i];
			re[len+1][0]=a[i]; 
			re[len+1][1]=i; 
            if(j>len) 
			len=j; 
	}				
	re[len+1][0]=-50000;re[len+1][1]=n;	
	re[0][0]=-50000;re[0][1]=1;					
for (int i=0;i<=len;i++)
  {
   for (int j=re[i][1]+1;j<re[i+1][1];j++)
    {
      int g=min(fabs(re[i][0]-a[j]),fabs(re[i+1][0]-a[j])) ; 
	  ans+=g;	
	}
  }
  
  memset(re,0,sizeof(re));
   memset(d,0,sizeof(re));
  for (int  i=n;i>=1;i--) 
    {
           j=lower_bound(d+1,d+len+1,a[i])-d;
            d[j]=a[i];
			re[len+1][0]=a[i]; 
			re[len+1][1]=i; 
            if(j>len) 
			len=j; 
	}				
	re[len+1][0]=-50000;re[len+1][1]=1;	
	re[0][0]=-50000;re[0][1]=n;					
 for (int i=0;i<=len;i++)
  {
   for (int j=re[i][1]-1;j>re[i+1][1];j--)
    {
      int g=min(fabs(re[i][0]-a[j]),fabs(re[i+1][0]-a[j])) ; 
	  ans1+=g;	
	}
  }  
  cout<<min(ans,ans1);
return 0;
}