石子归并:传送门
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堆石子的数量,计算最小合并代价。
Input
第1行:N(2 <= N <= 100) 第2 - N + 1:N堆石子的数量(1 <= Ai
<= 10000)
Output
输出最小合并代价
Sample Input
4 1 2 3 4
Sample Output
19
思路:
区间dp模板题
注意dp数组初值,对于每个石子dp[i][i]=0, 其他设为无限大
模板:
//w[i][j]表示从i 到 j的花费
for(i = 1;i <= n;i++)
dp[i][i] = 初始值;
for(len = 2;len <= n;len++){//len选择区间长度
for(i = 1;i <= n;i++){//枚举起点
j = i + len - 1;//合并终点
if(j > n)break;//不可越界
for(k = i;k < j;k++)//枚举分割点,寻找最优分割
dp[i][j] = max(dp[i][j], dp[i][k] + dp[k + 1][j] + w[i][j]);//状态转移
}
}
没有优化:
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#define MAXN 105
#define INF 0x3f3f3f3f
using namespace std;
int main()
{
int n;
scanf("%d", &n);
int i, j, k, len, a[MAXN], dp[MAXN][MAXN], sum[MAXN];
memset(dp, INF, sizeof(dp));
memset(sum, 0, sizeof(sum));
for(i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
dp[i][i] = 0;
sum[i] = sum[i-1] + a[i];
}
for(len = 2; len <= n; len++)
{
for(i = 1; i <= n; i++)
{
int j = i+len-1;
if(j > n) break;
for(k = i; k < j; k++)
{
dp[i][j] = min(dp[i][j], dp[i][k] + dp[k+1][j] + sum[j] - sum[i-1]);
}
}
}
printf("%d\n", dp[1][n]);
return 0;
}
四边形优化:
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#define MAXN 105
#define INF 0x3f3f3f3f
using namespace std;
int main()
{
int n;
scanf("%d", &n);
int i, j, k, len, a[MAXN], dp[MAXN][MAXN], sum[MAXN], s[MAXN][MAXN];
memset(dp, INF, sizeof(dp));
memset(sum, 0, sizeof(sum));
for(i = 1; i <= n; i++)
{
scanf("%d", &a[i]);
s[i][i] = i;
dp[i][i] = 0;
sum[i] = sum[i-1] + a[i];
}
for(len = 2; len <= n; len++)
{
for(i = 1; i <= n; i++)
{
int j = i+len-1;
if(j > n) break;
for(k = s[i][j-1]; k <= s[i+1][j]; k++)
{
if(dp[i][j] > dp[i][k] + dp[k+1][j] + sum[j] - sum[i-1])
{
dp[i][j] = dp[i][k] + dp[k+1][j] + sum[j] - sum[i-1];
s[i][j] = k;
}
}
}
}
printf("%d\n", dp[1][n]);
return 0;
}