#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;
}
AcWing 1069 凸多边形的划分 题解 (动态规划—DP—区间DP)
最新推荐文章于 2022-08-03 22:22:10 发布