【刷题】动态规划——区间DP：凸多边形的划分【高精度】

题目链接

注意题目中的三角形是互不相交的。假设选取如下三角形：

由于三角形之间互不相交，所以剩下的三角形只能在[1, 3]、[3, 6]和[6, 1]区间选取。如果我们有[1, 3]和[3, 6]区间三角形顶点积的和的最小值f[1, 3]、f[3, 6]，那我们就能知道区间[1, 6]的最小值f[1, 6] = f[1, 3] + f[3, 6] + w[1] * w[3] * w[6]。

接下来假设选取蓝色方框的三角形最优：
那f[1,8] = f[1,6] + f[6,8] + w[1] * w[6] * w[8]
其中f[6,8] = f[6,7] + f[7,8] + w[6] * w[7] * w[8] = 0 + 0 + w[6] * w[7] * w[8]。

最后答案就是f[1, 8]

这符合环形区间DP的特性，区间和最小值由子区间和最小值相加得来。

---

f[i, j]表示区间[i, j]的所有三角形（以(i, i+1), (i+1, i+2), …, (j-1, j), (j, i) 为边的三角形）顶点权值积的和 最小。（注意包括(j, i)这条边）

f[i, j] = min(f[i, j], f[i, k] + f[k, j] + w[i] * w[k] * w[j] ) 其中i < k < j
此处k不能等于i、j，因为三角形不能有两端点在同一个点上。

---

最后递推式和能量项链十分类似。

---

答案很大，需要用高精度存储

#include <iostream>
#include <cstring>
#include <vector>
using namespace std;

typedef long long LL;

int n;
LL w[55];
vector<LL> f[55][55];
vector<int>INF(50, 9);
vector<LL> add(vector<LL> A, vector<LL> B) {
	vector<LL> C;
	LL t = 0;
	for (int i = 0; i < A.size() || i < B.size(); i ++ ) {
		if (i < A.size()) t += A[i];
		if (i < B.size()) t += B[i];
		C.push_back(t % 10);
		t /= 10;
	}
	if (t) C.push_back(1);
	return C;
}

vector<LL> mul(vector<LL> A, LL b) {
	vector<LL> C;
	LL t = 0;
	for (int i = 0; i < A.size() || t; i ++ ) {
		if(i < A.size()) t += A[i] * b;
		C.push_back(t % 10);
		t /= 10;
	}
	return C;
}

// A < B 
bool cmp(vector<LL> A, vector<LL> B) {
	if (A.size() != B.size()) return A.size() < B.size();
	else {
		for (int i = A.size() - 1; i >= 0; i -- ) {
			if (A[i] != B[i]) {
				return A[i] < B[i];
			}
		}
	}
	return false; 
}

void print(vector<LL> A) {
	for (int i = A.size() - 1; i >= 0; i -- ) {
		printf("%d", A[i]);
	}
	printf("\n");
}

int main() {
	scanf("%d", &n);
	for (int i = 1; i <= n; i ++ ) {
		scanf("%d", &w[i]);
	}
	
	for (int l = 2; l <= n; l ++ ) {
		for (int i = 1; i + l  <= n; i ++ ) {
			int j = i + l;
			f[i][j].assign(INF.begin(), INF.end());
			vector<LL> tmp1(1, 1);
			for (int k = i + 1; k < j; k ++ ) {
				vector<LL> tmp2( add(add(f[i][k], f[k][j]), mul(mul(mul(tmp1, w[i]), w[k]), w[j])) );
				if(cmp(tmp2, f[i][j])) f[i][j].assign(tmp2.begin(), tmp2.end());
			}
		}
	}
	print(f[1][n]);
	return 0;
}