传送门: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;
}