Making the Grade [POJ3666] [DP]

题意：

给定一个序列，以最小代价将其变成单调不增或单调不减序列，代价为Σabs(i变化后-i变化前)，序列长度<=2000，单个数字<=1e9

输入：（第一行表示序列长度，之后一行一个表示序列第i的大小）

7
1
3
2
4
5
3
9

输出：（代价）

  3

分析：

这道题有bug，只要求单调不减序列

首先，对于这种问题，我们容易想到DP

记dp[i][j]为处理到i个，最高的为j

那我们的dp[i][j]=min(dp[i-1][k])+abs(j-h[i]) (k<=j)

但是显然数字太大了，我们需要离散化

先把序列从小到大排序，存在另一个数组b[]里

把j记成第j大的，那么状态转移方程为dp[i][j]=min(dp[i-1][k])+abs(b[j]-h[i]) (k<=j)

这样空间复杂度就够了

但时间复杂度还不够

我们仔细观察可以发现，dp[i-1][k]的每次从1开始循环找最小值是浪费的

我们在j从小到大循环上来的时候，就记录下最小值mn，那么转移方程就优化成dp[i][j]=mn+abs(b[j]-h[i])

至此，本题解决。

（提示：开long long，inf要开大！）

Code：

 #include<set>
 #include<map>
 #include<queue>
 #include<stack>
 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define RG register ll
 #define rep(i,a,b)    for(RG i=a;i<=b;++i)
 #define per(i,a,b)    for(RG i=a;i>=b;--i)
 #define ll long long
 #define inf (1<<30)
 #define maxn 2005
 using namespace std;
 ll n;
 ll a[maxn],b[maxn],dp[maxn][maxn];
 inline ll read()
 {
     ll x=,f=;char c=getchar();
     while(c<''||c>''){if(c=='-')f=-;c=getchar();}
     while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
     return x*f;
 }

 void work()
 {
     rep(i,,n)
     {
         ll mn=inf;
         rep(j,,n)
         {
             mn=min(mn,dp[i-][j]);
             dp[i][j]=(a[i]-b[j]>=?a[i]-b[j]:b[j]-a[i])+mn;
         }
     }
     ll ans=inf;
     rep(i,,n) ans=min(ans,dp[n][i]);
     cout<<ans;
 }

 int main()
 {
     n=read();
     rep(i,,n) a[i]=b[i]=read();
     sort(b+,b++n);
     work();
     return ;
 }