分级(LIS) 转移状态优化，三维优化为二维

分级

给定长度为N的序列A，构造一个长度为N的序列B，满足：

1、B非严格单调，即B1≤B2≤…≤BN或B1≥B2≥…≥BN。
2、最小化 S=∑Ni=1|Ai−Bi|。

只需要求出这个最小值S。

输入格式

第一行包含一个整数N。

接下来N行，每行包含一个整数Ai。

输出格式

输出一个整数，表示最小S值。

数据范围

1≤N≤2000,
1≤|Ai|≤10^9

输入样例：

7
1
3
2
4
5
3
9

输出样例：

3

分析：

很明显这个状态转移是一个三重循环，决策集合中的元素只增多不减少，所以可以优化为二维，用一个变量记录f[i-1][1]...f[i-1][j]的最小值，可以去掉第三重循环

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int N = 2010,INF=0x3f3f3f3f;
int A[N],B[N],f[N][N],n;
int work(){
    for(int i=1;i<=n;i++) B[i]=A[i];
    sort(B+1,B+1+n);
    /*
    朴素版做法
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            f[i][j]=INF;
            for(int k=1;k<=j;k++){
                f[i][j]=min(f[i][j],f[i-1][k]+abs(A[i]-B[j]));
            }
        }
    }
    */
    for(int i=1;i<=n;i++){
        int minv=INF;
        for(int j=1;j<=n;j++){
            minv=min(minv,f[i-1][j]);
            f[i][j]=minv+abs(A[i]-B[j]);
        }
    }
    int res=INF;
    for(int i=1;i<=n;i++)
        res=min(res,f[n][i]);
    return res;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&A[i]);
    int res=work();
    reverse(A+1,A+1+n);
    res=min(res,work());
    printf("%d\n",res);
    return 0;
}