【问题描述】
大老板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;
}