区间dp的四边形优化问题

四边形优化

首先我们先要了解满足四边形不等式 所需要的条件。

当 $a1<a2⩽b1<b2$

且有 $m[a1,b1]+m[a2,b2]⩽m[a1,b2]+m[a2,b1]$时，满足四边形不等式
很显然，对于区间动态规划问题，上述不等式通常是成立的
设$s[i,j]$为$m[i,j]$中取得最小值的k值，也就是$m[i,j]=m[i,k]+m[k,j]$中的k
则有$s[i,j-1]⩽s[i,j]⩽s[i+1][j]$
这个式子很容易证明，可以自己在草稿纸上尝试一下。
有了这个式子，那么k的取值也就是缩小到了$s[i,j-1]到s[i+1,j]$之间了，从而降低了时间复杂度。

接下来上一道需要使用四边形优化的区间dp题。
题目链接 < 点我 >

石子归并问题

N堆石子摆成一个环。现要将石子有次序地合并成一堆。规定每次只能选相邻的2堆石子合并成新的一堆，并将新的一堆石子数记为该次合并的代价。计算将N堆石子合并成一堆的最小代价。

例如： 1 2 3 4，有不少合并方法 1 2 3 4 => 3 3 4(3) => 6 4(9) => 10(19) 1 2 3 4 =>
1 5 4(5) => 1 9(14) => 10(24) 1 2 3 4 => 1 2 7(7) => 3 7(10) => 10(20)

括号里面为总代价可以看出，第一种方法的代价最低，现在给出n堆石子的数量，计算最小合并代价。

这里要注意的是，环形石头的表示方法可以创建两条线性的石头序列，通过移动即可表示出环形。

#include<bits/stdc++.h>
#define bug(a) (cout<<'*'<<a<<endl)
#define bugg(a,b) (cout<<'*'<<a<<' '<<b<<endl)
#define buggg(a,b,c) (cout<<'*'<<a<<' '<<b<<' '<<c<<endl)
#define pn (cout<<endl)
typedef long long ll;
using namespace std;
int ar[5000];   
int sum[5000];  //前缀和，用来求区间和
int dp[2002][2002];
int f[2002][2002]; //k取值范围的数组
int main()
{
	int n, i, j, dis;
	memset(dp, 0x3f, sizeof(dp));
	scanf("%d", &n);
	for (i = 1; i <= n; i++)
	{
		scanf("%d", &ar[i]);
		ar[i + n] = ar[i];
		sum[i] = sum[i - 1] + ar[i];
		dp[i][i] = 0;
		dp[i + n][i + n] = 0;
		f[i][i] = i;//自己到自己的k值当然为本身
	}
	for (i = n + 1; i <= 2 * n; i++)
		sum[i] = sum[i - 1] + ar[i], f[i][i] = i;
	for (int k = 1; k < n; k++)
	{
		for (int i = 1; i + k <= 2 * n; i++)
		{
			int j = i + k;
			for (int p = f[i][j - 1]; p <= f[i + 1][j]; p++)//四边形优化
			{
				if (dp[i][j] > dp[i][p] + dp[p + 1][j] + sum[j] - sum[i - 1])
				{
					dp[i][j] = dp[i][p] + dp[p + 1][j] + sum[j] - sum[i - 1];
					f[i][j] = p;//记得更新！！！
				}
			}
		}
	}
	int ans = 0x3f3f3f3f;
	for (i = 1; i <= n; i++)
		ans = min(dp[i][i + n - 1], ans);
	printf("%d\n", ans);
	return 0;
}