最大得分（三维DP）

最新推荐文章于 2025-11-27 19:19:08 发布

原创最新推荐文章于 2025-11-27 19:19:08 发布 · 724 阅读

2 ·

本内容遵循CC 4.0 BY-SA版权协议

动态规划专栏收录该内容

3 篇文章

订阅专栏

这篇博客介绍了一个利用三维动态规划解决的数学问题，具体是关于如何将一个长度为n的序列分成k个非空且连续下标的部分，以使所有部分的GCD（最大公约数）之和最大化。作者给出了详细的思路解析，包括状态转移方程，并提供了C++代码实现。该算法适用于处理序列分割问题，涉及数论和动态规划的知识。

题目链接
Description
你有一个长度为n的a序列，要求你将这个序列分成 k 个部分。
要求：

每一个部分都非空，
部分与部分不能重叠
每一个部分都是一串连续下标的数

每一个部分的分数等于这个部分的 gcd ，输出一个sum表示这个大小为n的序列，分成k 部分所获得的最大分数。
Input
第一行包含两个整数，分别是n，k (k <=n <=10⁴,1<= k <= 50),n 表示a 序列的大小，k 表示这个序列应该要被分成的部分数。第二行n 个整数 a_i(1<=a_i<=100)

思路：
三维DP
dp[i] [j] [l]表示1~i个数，已经分成了l 个部分，当前部分的gcd为j的gcd之和。
状态转移：
第i+1个数继续放在第l部分：
dp[i+1] [g] [l]=max(dp[i+1] [g] [l]，dp[i] [j] [l]-j+g) 其中g=gcd(j,s[i+1])
第i+1个数放在新的部分：
dp[i+1] [s[i+1]] [l+1]=max(dp[i+1] [s[i+1]] [l+1]，dp[i] [j] [l]+s[i+1])
code：

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <string>
#define debug(x) cout << #x << ":" << x << endl;
#define bug cout << "********" << endl;
#define fi first
#define se second
using namespace std;
typedef long long ll;
const double eps = 1e-8;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 + 7;
const int maxn = 2e5 + 10;

int n, k, s[maxn];
int dp[10020][110][60];

int main()
{
    scanf("%d%d", &n, &k);
    for (int i = 1; i <= n;i++)
        scanf("%d", &s[i]);
    dp[1][s[1]][1] = s[1];
    for (int i = 2; i <= n;i++){
        for (int j = 1; j <= 100;j++){
            int g = __gcd(j, s[i]);
            for (int l = 1; l <= k;l++){
                if(dp[i-1][j][l-1]){
                    dp[i][s[i]][l] = max(dp[i][s[i]][l], dp[i - 1][j][l-1] + s[i]);
                }
                if(dp[i-1][j][l]){
                    dp[i][g][l] = max(dp[i][g][l], dp[i - 1][j][l] - j + g);
                }
            }
        }
    }
    int ans = 0;
    for (int i = 1; i <= 100;i++)
        ans = max(ans, dp[n][i][k]);
    printf("%d\n", ans);
    return 0;
}