B4172 学习计划 题解

B4172 学习计划 题解

思路

可以将收益式子换一下,设 c i c_i ci 为 a i a_i ai 被分到的段的编号,那收益式子变成 ∑ i = 1 n a i × b c i \sum_{i=1}^n a_i \times b_{c_i} ∑i=1nai×bci。

很显然的 dp, 设 f i , j f_{i,j} fi,j 为将 a a a 的前 i i i 个数分成 j j j 段的最大收益。

那现在有两种选择。

  1. a i − 1 a_{i-1} ai−1 选的是第 j j j 段,这样的收益是 f i − 1 , j + a i × b j f_{i-1,j}+a_i\times b_j fi−1,j+ai×bj。
  2. a i − 1 a_{i-1} ai−1 选的是第 j − 1 j-1 j−1 段,这样的收益是 f i − 1 , j − 1 + a i × b j f_{i-1,j-1}+a_i\times b_j fi−1,j−1+ai×bj。

整理一下就能得到转移方程 f i , j = max ⁡ ( f i − 1 , j , f i − 1 , j − 1 ) + a i × b j f_{i,j}= \max (f_{i-1,j},f_{i-1,j-1})+a_i\times b_j fi,j=max(fi−1,j,fi−1,j−1)+ai×bj。

由于有负数,初始值要设为 − inf ⁡ -\inf −inf, f 0 , 0 f_{0,0} f0,0 设为 0 0 0。

代码

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int N = 2010;
int a[N], b[N];
int f[N][N];
void run() {
	memset(f, -0x3f, sizeof(f));
	int n, m;
	cin >> n >> m;
	for (int i = 1; i <= n; i++) cin >> a[i];
	for (int i = 1; i <= m; i++) cin >> b[i];
	f[0][0] = 0;
	for (int i = 1; i <= n; i++)
		for (int j = 1; j <= min(i, m); j++)
			f[i][j] = max(f[i - 1][j - 1], f[i - 1][j]) + a[i] * b[j];
	cout << f[n][m] << endl;
}
int main() {
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	int _ = 1;
	cin >> _;
	while (_--) run();
	return 0;
}
相关推荐
智者知已应修善业5 小时前
【求中位数】2024-1-23
c语言·c++·经验分享·笔记·算法
地平线开发者6 小时前
PTQ 量化数值范围与优化
算法·自动驾驶
sali-tec6 小时前
C# 基于halcon的视觉工作流-章68 深度学习-对象检测
开发语言·算法·计算机视觉·重构·c#
测试人社区-小明6 小时前
智能弹性伸缩算法在测试环境中的实践与验证
人工智能·测试工具·算法·机器学习·金融·机器人·量子计算
罗西的思考7 小时前
【Agent】MemOS 源码笔记---(5)---记忆分类
人工智能·深度学习·算法
qq_4335545410 小时前
C++数位DP
c++·算法·图论
AshinGau10 小时前
Softmax 与 交叉熵损失
神经网络·算法
似水এ᭄往昔10 小时前
【C++】--AVL树的认识和实现
开发语言·数据结构·c++·算法·stl
栀秋66610 小时前
“无重复字符的最长子串”:从O(n²)哈希优化到滑动窗口封神,再到DP降维打击!
前端·javascript·算法
xhxxx10 小时前
不用 Set,只用两个布尔值:如何用标志位将矩阵置零的空间复杂度压到 O(1)
javascript·算法·面试