Making the Grade(POJ3666)

题目大意：

给出长度为n的整数数列,每次可以将一个数加1或者减1,最少要多少次可以将其变成单调增或者单调减(不严格).

题解：

1.一开始我有一个猜想,就是不管怎么改变,最终的所有数都是原来的某个数。然而我并不会证明,然而我属于那种不彻底弄清楚就不会去写的那种顽固分子,于是就拖了好几天。网络上有很多关于此题的题解,确实用了这个猜想来离散化,但是都是讲怎么dp,然后最后扯一句“由于数据比较大,可以离散化”之类的话,要么就是相当粗略的证明(也许已经说的够清楚了只不过我没理解...)。

2.今天早上起来洗漱的时候,感觉头脑比较清醒,再次想了一下这个问题,想到一个自认为正确的证明：

记原来的数从小到大排序后分别是$a_1\ a_2\ a_3\cdots a_n$ 修改后从左到右分别是$b_1\ b_2\ b_3\cdots b_n$. 为了方便描述,在数轴上标出这些点,称为关键点。

假设存在$a_s<b_i<=b_{i+1}<=\cdots <=b_j<a_{s+1}$

情况一：如果这些b都相等,那么把这些b都改成$a_s$或者$a_{s+1}$ 肯定会有一种更优。

情况二：如果不全相等,那么肯定存在 $b_p\ b_{p+1}\ b_{p+2}\cdots b_q$,他们的值相等,那么把他们移到左边的关键点或者右边的关键点,肯定有一种会更加优. 不断这样移动,最后就变成情况一了。

综上至少存在一种最优方案,最终的所有数都是原来的某个数。

因此可以离散化之后做dp,dp[i][j]表示把前i个数变成单调增(不严格)且第i个数变成原来第j大的数的最小代价。

dp[i][j]=min{dp[i-1][1...j]}+abs(a[i]-b[j]).

单调减(不严格)的情况也一样,更加方便的是可以把原数组倒转后做单调增的dp。

另外这题和CF Codeforces Round #371 (Div. 1)的C题类似,有个超简单的nlogn做法。可以看这篇：http://www.cnblogs.com/vb4896/p/5894578.html

代码：

 #include <iostream>
 #include <cstdio>
 #include <cmath>
 #include <cstring>
 #include <algorithm>
 #include <vector>
 #include <map>
 #include <cstdlib>
 #include <set>
 using namespace std;

 #define X first
 #define Y second
 #define Mod 1000000007
 #define N 3010
 typedef long long ll;
 typedef pair<int,int> pii;

 inline int Mul(int x,int y){return 1ll*x*y%Mod;}
 inline int Add(int x,int y){return ((x+y)%Mod+Mod)%Mod;}

 int n,m;
 int a[N],b[N];
 ll dp[N][N],Ans=1ll<<;

 void Solve()
 {
     for (int j=;j<=m;j++) dp[][j]=abs(b[j]-a[]);
     for (int i=;i<=n;i++)
     {
         ll tmp=1ll<<;
         for (int j=;j<=m;j++)
         {
             tmp=min(tmp,dp[i-][j]);
             dp[i][j]=abs(b[j]-a[i]);
             dp[i][j]+=tmp;

         }
     }
     for (int j=;j<=m;j++) Ans=min(Ans,dp[n][j]);
 }

 int main()
 {
     //freopen("in.in","r",stdin);
     //freopen("out.out","w",stdout);

     scanf("%d",&n);
     for (int i=;i<=n;i++) scanf("%d",&a[i]),b[i]=a[i]=a[i]-i;
     sort(b+,b+n+);
     for (int i=;i<=n;i++) if (i== || b[i]!=b[i-]) b[++m]=b[i];

     Solve();
     printf("%I64d\n",Ans);

     return ;
 }