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=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

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

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