POJ - 3666 Making the Grade DP

传送门：POJ3666

题意：给定一个序列Ai，要求将其变成一个非递减或者非递增序列Bi，花费为∑ |Ai - Bi|，问最小花费是多少。

思路：首先我们明确如果需要将Ai改变，那么Ai要么变成Ai+1要么变成Ai-1，否则就不变，这样一定是最优的。我证明不了，不过可以自行举几个例子模拟一下改变过程就知道了。

因此，Bi中的数必定都属于Ai，既然有了Bi的范围集合，那么我们就可以dp转移了，

dp[i][j]代表将Ai的前i个数变成满足题意的序列并且最大值为A[j]的最小花费。

转移方程： dp[i][j] = min( dp[i - 1][k] )(Ak <=Aj) + abs(A[j] - A[i])

如果直接转移复杂度是O（n^3）,不过如果我们将Ai排序后再进行转移，求min的部分就可以边转移边动态求出来了。这样复杂度降为O（n^2）.

代码：

#include<iostream>
#include<algorithm>
#include<string.h>
#define ll long long
#define pb push_back
#define fi first
#define se second
#define pi acos(-1)
#define inf 0x3f3f3f3f
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define rep(i,x,n) for(int i=x;i<n;i++)
#define per(i,n,x) for(int i=n;i>=x;i--)
using namespace std;
typedef pair<int,int>P;
const int MAXN=100010;
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int a[2100], b[2100], dp[2][2100];
int main()
{
	int n;
	cin >> n;
	for(int i = 1; i <= n; i++)
	cin >> a[i], b[i] = a[i];
	sort(b + 1, b + n + 1);
	//memset(dp, inf, sizeof(dp));
	int e = 1, tmp = inf;
	for(int i = 1; i <= n; i++)
	{
		tmp = inf;
		for(int j = 1; j <= n; j++)
		{
			tmp = min(tmp, dp[e ^ 1][j]);
			dp[e][j] = tmp + abs(b[j] - a[i]);
		}
		e ^= 1;
	}
	cout << *min_element(dp[e ^ 1] + 1, dp[e ^ 1] + n + 1) << endl;
 	return 0;
}