AcWing 1069 凸多边形的划分 题解 （动态规划—DP—区间DP）

原题传送门

#include<bits/stdc++.h>

using namespace std;

const int N = 55, M = 35;

typedef long long ll;

int n;
int w[N];
ll f[N][N][M];//f[i][j]表示两个底角为i、j，顶角在这个区间内的三角形的三端顶点之积的最小值，最后一维M表示这个数用高精度储存，M是这个数的最大位数

void add(ll a[], ll b[]){//高精度加法板子 
	static ll c[M];
	memset(c, 0, sizeof c);
	ll t = 0;
	for(int i = 0; i < M; i ++ ){//高精度数组的第0位表示个位 
		t += a[i] + b[i];
		c[i] = t % 10;
		t /= 10;
	}
	memcpy(a, c, sizeof c);
} 

void mul(ll a[], int b){//高精度乘法板子 
	static ll c[M];
	memset(c, 0, sizeof c);
	ll t = 0;
	for(int i = 0; i < M; i ++ ){
		t += a[i] * b;
		c[i] = t % 10;
		t /= 10;
	}
	memcpy(a, c, sizeof c);
}

int cmp(ll a[], ll b[]){//高精度比大小 
	for(int i = M - 1; i >= 0; i -- ){
		if(a[i] > b[i]) return 1;
		else if(a[i] < b[i]) return -1;
	}
	return 0;
}

void Print(ll a[]){//高精度输出 
	int k = M - 1;
	while(k && !a[k]) k -- ;
	while(k >= 0) cout<<a[k -- ];
	cout<<endl;
}

int main()
{
	cin>>n;
	for(int i = 1; i <= n; i ++ ){
		cin>>w[i];
	}
	
	ll temp[M];
	for(int len = 3; len <= n; len ++ ){//枚举三角形两个固定角之间的区间长度 
		for(int l = 1; l + len - 1 <= n; l ++ ){//枚举其中一个固定角 
			int r = l + len - 1;//确定另一个固定角 
			f[l][r][M - 1] = 1;//f数组此时用高精度表示，相当于0000000000001(右侧表示高精度的最高位)
			for(int k = l + 1; k < r; k ++ ){//枚举哪个活动的角 
				memset(temp, 0, sizeof temp);
				temp[0] = w[l];
				mul(temp, w[k]);
				mul(temp, w[r]);
				add(temp, f[l][k]);
				add(temp, f[k][r]);
				if(cmp(f[l][r], temp) > 0){
					memcpy(f[l][r], temp, sizeof temp);
				}
			}
		}
	}
	
	Print(f[1][n]);
	return 0;
}