【能量项链】

Kent_J_Truman

已于 2024-08-26 21:10:29 修改

阅读量272

点赞数 6

分类专栏： dp 算法文章标签：算法 c++

于 2024-08-25 23:10:30 首次发布

本文链接：https://blog.csdn.net/m0_73669127/article/details/141536602

版权

算法同时被 2 个专栏收录

123 篇文章 0 订阅

订阅专栏

18 篇文章 0 订阅

订阅专栏

题目

错误代码（10过7）

#include <bits/stdc++.h>
using namespace std;
const int N = 110;
int a[N];
int n;
int f[N][N];
int cal(int i, int k, int j)
{
    if(i <= k && k <= j)
    {
        return a[i] * a[(k+1)%n] * a[(j+1)%n];
    }
    else if(k < j && j <= i)
    {
        return a[i] * a[(k+1)%n] * a[(j+1)%n];
    }
    else if(j <= i && i <= k)
    {
        return a[i] * a[(k+1)%n] * a[(j+1)%n];
    }
    
    return 0;
}
int main()
{
    memset(f, -0x3f, sizeof f);
    cin >> n;
    for(int i = 0; i < n; i++)
    {
        cin >> a[i];
        f[i][i] = 0;
    }
    
    for(int len = 2; len <= n; len++)
    {
        for(int i = 0; i < n; i++)
        { 
            int j = (i + len - 1) % n;
            bool cir = false;
            for(int t = i; ; t++)
            {
                int k = t % n;
                if(k == i)
                {
                    if(cir) break;
                    else cir = true;
                }
                if(cal(i, k, j)) 
                {
                    f[i][j] = max(f[i][j], f[i][k] + f[(k+1)%n][j] + cal(i, k, j));
                }
                
            }
        }
    }

    int res = -0x3f3f3f3f;
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
        {
            if(abs(i-j) == 1) res = max(res, f[i][j]);
        }
    cout << res;
    
    return 0;
}

思路

定义状态

$f[i][j]$ 为 $i \rightarrow j$ $(i <= j)$ 或者 $j \rightarrow n-1 \rightarrow i$ $(j < i)$ 的方法集合。

例如 $i \rightarrow j$ 指的是第 $i$ 个珠子作为聚合珠串的头，第 $j$ 个珠子为尾的情况，珠串包含序列中所有珠子。

$f[i][j]$ 的值为方法集合中最大的释放能量值。

状态转移

$f[i][j] = f[i][k] + f[k+1][j] + cal(i, k, j)$

其中 $cal(i,k,j)$ 为将两个珠串头尾聚合成 $f[i][j]$ 所释放的能量值。

目标状态

将所有珠子纳入考虑的状态： $f[0][n-1]$ 、 $f[i][j]$ $(i > j) \wedge (abs(i-j) = 1)$

正确代码

#include <bits/stdc++.h>
using namespace std;
const int N = 110, M = N << 1;
int a[M];
int n;
int f[M][M];
int cal(int i, int k, int j)
{
    return a[i] * a[k] * a[j];
}
int main()
{
    memset(f, -0x3f, sizeof f);
    cin >> n;
    for(int i = 1; i <= n; i++)
    {
        cin >> a[i];
        a[n+i] = a[i];
        f[i][i] = 0;
        f[n+i][n+i] = 0;
    }

    for(int len = 2; len <= n; len++)
    {
        for (int i = 1, j; (j = i + len - 1) < n << 1; i++)
        { 
            for(int k = i; k+1 <= j; k++)
            {
                f[i][j] = max(f[i][j], f[i][k] + f[k+1][j] + cal(i,k+1,j+1));
            }
        }
    }

    int res = -0x3f3f3f3f;
    for(int i = 1; i <= n; i++)
    {
        res = max(res, f[i][i+n-1]);
    }
    cout << res;
    
    return 0;
}