Making the Grade [POJ3666] [DP]

时间:2024-01-11 16:24:14

题意:

给定一个序列,以最小代价将其变成单调不增或单调不减序列,代价为Σ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 ;
}