bzoj 1592: [Usaco2008 Feb]Making the Grade 路面修整【dp】

因为是单调不降或单调不升，所以所有的bi如果都是ai中出现过的一定不会变差
以递增为例，设f[i][j]为第j段选第i大的高度，预处理出s[i][j]表示选第i大的时，前j个 a与第i大的值的差的绝对值 的和。
转移显然是
\[ f[i][j]=min{f[i-1][k]+s[i][j]-s[i][k]} \]
这样看起来是\( O(n^3) \)的，但是注意到s[i][j]固定
\[ f[i][j]=min{f[i-1][k]-s[i][k]}+s[i][j] \]
这样就可以在处理i-1的时候求出mn[i-1][k]为前k个中最小的f[i-1][k]-s[i][k]，所以时间复杂度变成了\( O(n^2) \)，空间是可以用滚动数组压到\( O(n) \)的，但是方便起见（懒）就只写了\( O(n^2) \)的

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=2005,inf=(1<<30)-1+(1<<30);
int n,a[N],b[N],f[N][N],s[N][N],mn[N][N],ans=inf;
int read()
{
    int r=0,f=1;
    char p=getchar();
    while(p>'9'||p<'0')
    {
        if(p=='-')
            f=-1;
        p=getchar();
    }
    while(p>='0'&&p<='9')
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r*f;
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
        a[i]=b[i]=read();
    sort(b+1,b+1+n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            s[i][j]=s[i][j-1]+abs(a[j]-b[i]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            f[i][j]=mn[i-1][j]+s[i][j],mn[i][j]=min(f[i][j]-s[i+1][j],mn[i][j-1]);
    for(int i=1;i<=n;i++)
        ans=min(ans,f[i][n]);
    for(int i=n;i>=1;i--)
        for(int j=1;j<=n;j++)
            f[i][j]=mn[i+1][j]+s[i][j],mn[i][j]=min(f[i][j]-s[i-1][j],mn[i][j-1]);
    for(int i=1;i<=n;i++)
        ans=min(ans,f[i][n]);
    printf("%d\n",ans);
    return 0;
}