LibreOJ 10148 能量项链 ------- 区间dp

题目链接 https://loj.ac/problem/10148

题解：经典区间dp题合并石子（圆形版）变形题， 核心思路是化曲为直。

因为最后一定是2个珠子合并，每个珠子是由连续若干个珠子合并的，所以如果项链拆开成直线，不允许圆形合并的话，就枚举端点k，则dp[i][k]和dp[k+1][j]分别表示两边珠子合并最大值，然后手工合并最后一次即可，因为是圆形的，所以可以化曲为直，枚举所有的项链端点将项链拆成直线，这样解决n次直线上的该问题求最大值即可。这样做时间复杂度大了，可以优化成用dp[i][k]表示以i为起点连续k个珠子合并的最大值，则 dp[i][j] = max(dp[i][k]+dp[(i+k-1)%n][j-k]+h1*t1*t2)(h1表示最后的第一颗珠子头标记....)  

参考 https://oi.men.ci/noip2006-energy/

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int n;
int flag[110];
int dp[110][110];

// dp[i][j] = max(dp[i][k]+dp[(i+k-1)%n][j-k]+)
// 2x3 3x2  =>>  2x3x2, 3x2x3
int main() {
    scanf("%d", &n);
    for(int i = 0;i < n;i++) scanf("%d", &flag[i]);
    for(int len = 2;len <= n;len++)
        for(int l = 0;l < n;l++) {
            int r = (l+len-1) % n;
            for(int k = 1;k < len;k++) {
                int h1 = flag[l];
                int h2 = flag[(l+k)%n];
                int t1 = h2;
                int t2 = flag[(r+1)%n];
                dp[l][len] = max(dp[l][len], h1*t1*t2+dp[l][k]+dp[(l+k)%n][len-k]);               
            }
        }

    int res = 0;
    for(int i = 0;i < n;i++) res = max(res, dp[i][n]);

    printf("%d\n", res);

  //  scanf("%d", &res);

    return 0;
}