四边形不等式的动态规划优化算法不多说,只对它做一个简单的介绍:这个是优化动态规划算法时间复杂度的方法之一,空间复杂度只会增加几个常数,不用考虑。方法说实话,验证是极其麻烦的,所以很多时候我都是要看题解再知道,这里也就不多说了。
直接上代码:
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<ctime>
using namespace std;
#define REP(i, a, b) for(i = (a);i <= (b); ++i)
const int INF = 1000000000;
const int maxn = 500;
int n, m;
int a[maxn], d[maxn][maxn], dp[maxn][maxn];
int main() {
int i, j, k;
scanf("%d", &n);
scanf("%d", &m);
REP(i, 1, n)
scanf("%d", &a[i]);
REP(i, 1, n)
REP(j, i + 1, n) {
int mid = (i + j) >> 1;
REP(k, i, j)
d[i][j] += abs(a[k] - a[mid]);
}
REP(i, 1, m)
REP(j, 1, n)
dp[i][j] = INF;
REP(i, 1, n)
dp[1][i] = d[1][i];
REP(i, 2, m)
REP(j, i, n)
REP(k, 1, j)
if(dp[i][j] > dp[i - 1][k] + d[k + 1][j])
dp[i][j] = dp[i - 1][k] + d[k + 1][j];
printf("%d\n", dp[m][n]);
return 0;
}