POJ3666——Making the Grade(动态规划)

题目链接

       这题反正我是觉得用DP做有点恶心，看别人的题解做的，转移方程是:dp[i][j]表示前i个元素的最后一个元素为全部元素第j小时的最小代价，则dp[i][j]=min(dp[i-1][k])+abs(s[i]-cop[j]),     ( 0≤i<n , 0≤j<n , 0≤k<j )从这个式子来看应该是有三重循环的，但有为大牛很巧妙的省掉了k的循环，从而降低了时间复杂度，然后由于方程的递推是顺序的，所以可以省略掉第一维，从而使dp方程一维化，这样降低了空间复杂度(详见代码)，但是自己写的时候又wrong了一次，原因也是在discuss里看见的，那就是无穷大不够大，起先把ans赋值成了99999999然后由于数据过大，这个值就不是无穷大了，改正后就好了，值得注意的是此题数据较弱只写非递减的代码就能过了，但其实非递增的只需改一下sort函数就能实现了。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>

using namespace std;

typedef long long LL;

LL s[2005]= {0};
LL cop[2005]= {0};
LL dp[2005]= {0};

bool cmp(LL a,LL b)
{
    return a>b;
}

LL Abs(LL a)
{
    if(a>=0)return a;
    else return -a;
}

int main()
{
    //freopen("in.in","r",stdin);
    int n;
    scanf("%d",&n);
    for(int i=0; i<n; i++)
    {
        scanf("%lld",&s[i]);
        cop[i]=s[i];
    }

//
    sort(cop,cop+n);//非递减
    for(int i=0; i<n; i++) //初始化
        dp[i]=Abs(s[0]-cop[i]);
    for(int i=1; i<n; i++) //dp
    {
        LL temp=dp[0];
        for(int j=0; j<n; j++)
        {
            temp=min(dp[j],temp);
            dp[j]=temp+Abs(s[i]-cop[j]);
        }
    }
    LL ans=dp[0];
    for(int i=0; i<n; i++)
        ans=min(dp[i],ans);

///
    sort(cop,cop+n,cmp);//非递增
    for(int i=0; i<n; i++) //初始化
        dp[i]=Abs(s[0]-cop[i]);
    for(int i=0; i<n; i++) //dp
    {
        LL temp=dp[0];
        for(int j=0; j<n; j++)
        {
            temp=min(dp[j],temp);
            dp[j]=temp+Abs(s[i+1]-cop[j]);
        }
    }
    LL ANS=dp[0];
    for(int i=0; i<n; i++)
        ANS=min(dp[i],ANS);

    ans=min(ans,ANS);
    printf("%lld\n",ans);
    return 0;
}