洛谷 P8816 [CSP-J 2022] 上升点列（T4）

P.Space

已于 2024-03-08 18:06:24 修改

阅读量1.5k

点赞数 36

文章标签： c++ CSP-J 2022

于 2024-03-08 18:05:59 首次发布

本文链接：https://blog.csdn.net/pyc0705/article/details/136566280

版权

文章详细解析了CSP-J2022上升点列问题的动态规划算法，涉及状态转移方程、记忆化搜索策略以及如何计算两点之间的曼哈顿距离。作者提供了最终的C++代码实现，并强调了记忆化搜索在处理复杂求值顺序的重要性。

摘要由CSDN通过智能技术生成

题目传送门

[CSP-J 2022] 上升点列 - 洛谷https://www.luogu.com.cn/problem/P8816

算法解析

k = 0 且 xi, yi 值域不大时，这题是非常简单的 DP，类似「数字三角形」。

记 dp(x,y) 为「以 (x,y) 为终点，最长合法序列的长度」。

则对于所有（已经存在的）整点，有：

dp(x,y) = max {dp(x − 1, y), dp(x, y − 1)} + 1

xi, yi 值域比较大时：

可以考虑记 dp(n) 表示「以 n 号点结尾的合法序列，最长能有多长」。

dp(n) = max {dp(i) + 1}

i → n ✓

不会存在环状结构——因为合法序列必须向右、上方发展。

把刚刚的DP改造一下，就是本题正解：

记 dp(n, k) 表示「以 n 号点结尾，已经使用掉了 k 个自由点，获得的收益」。

dp(n,k) = max {dp(i, k − cost) + cost + 1}

i → n ✓

实现细节：本题的求值顺序值得注意，合法路径可能形如 P1 → P3 → P2。

有两种解决方法：

记忆化搜索（记忆化搜索最擅长解决求值顺序混乱的 DP）
预先按 x, y 排序，使得编号大的点一定是从编号小的点转移过来

这里记忆化搜索比较好写一些，我这里就只讲记忆化搜索了

先写一下求 a 到 b 需要补多少个点的函数，即两点曼哈顿距离再减一（a 在左下，b 在右上，否则返回无穷）

代码中 x[u] 表示 u 点的横坐标，y[u] 表示 u 点的纵坐标

int dis(int a, int b) {
	if(x[a] > x[b])
		return inf;
	
	if(y[a] > y[b])
		return inf;
	
	return x[b] - x[a]  + y[b] - y[a] - 1;
}

然后是 dp 函数，定义上面已经说过了

int dp(int now, int k)

首先判断如果自由点已经用完了，即 k < 0，那么返回负无穷（因为最后是取最大值）

int dp(int now, int k) {
	if(k < 0)
		return -inf;
}

既然是记忆化，那么就需要记忆

用 vis[n][k] 数组记录 dp(n, k) 是否访问过，val[n][k] 数组记录如果访问过的 dp(n, k) 的值

这样如果 vis[now][k] == true（访问过），则返回 val[now][k]

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
}

然后就该枚举它的前驱（代码中的 to），然后取里面最大的收益

这个记录最大收益的变量（代码中的 res）的初值一定要是 1，因为如果哪也去不了，那么就只能走到现在这一个点，也就是 now

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
	
	int res = 1;
	
	for(int to = 1; to <= n; ++to)
	
	return res;
}

接下来需要判断 to != now，然后计算出 to 到 now 需要补多少个点（代码中的 cost）

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
	
	int res = 1;
	
	for(int to = 1; to <= n; ++to)
		if(to != now) {
			int cost = dis(to, now);
		}
	
	return res;
}

再判断费用超出运算，就 contunue（如果走不到，dis 就会返回无穷，一定大于 k，所以不用特判走不到）

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
	
	int res = 1;
	
	for(int to = 1; to <= n; ++to)
		if(to != now) {
			int cost = dis(to, now);
			
			if(cost > k)
				continue;
		}
	
	return res;
}

接着就是往下递归了，now 变成了 to，预算费用还剩 k - cost，所以传进去是

dp(to, k - cost)

然后长度还需要加上 to 到 now 的距离，即 cost + 1，然后更新最大值（代码里的 res）

代码中的 bemax 函数是把第一个参数赋成两个参数的最大值用的，具体实现方法就是用一个三目运算符

void bemax(int &a, int b) {
	a = a > b ? a : b;
}

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
	
	int res = 1;
	
	for(int to = 1; to <= n; ++to)
		if(to != now) {
			int cost = dis(to, now);
			
			if(cost > k)
				continue;
			
			bemax(res, dp(to, k - cost) + cost + 1);
		}
	
	return res;
}

最后再将 vis[now][k] 设成 true，val[now][k] 设成 res

最后 return res 就行了

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
	
	int res = 1;
	
	for(int to = 1; to <= n; ++to)
		if(to != now) {
			int cost = dis(to, now);
			
			if(cost > k)
				continue;
			
			bemax(res, dp(to, k - cost) + cost + 1);
		}
	
	vis[now][k] = true;
	val[now][k] = res;
	
	return res;
}

主函数里需要枚举 i = 1 ~ n，j = 0 ~ k，然后传进去（n 为点的个数，k 为自由点的个数）

注意长度还需要加上没用的 k - j 个点，然后更新答案（代码中的 ans）

for(int i = 1; i <= n; ++i)
	for(int j = 0; j <= k; ++j)
		bemax(ans, dp(i, j) + k - j);

最后输出 ans 即可

最终代码

#include <cstdio>
#define N 1005
using namespace std;

const int inf = 0x7fffffff;

int n, k;
int x[N], y[N];
bool vis[N][N];
int val[N][N];
int ans;

void bemax(int &a, int b) {
	a = a > b ? a : b;
}

int dis(int a, int b) {
	if(x[a] > x[b])
		return inf;
	
	if(y[a] > y[b])
		return inf;
	
	return x[b] - x[a]  + y[b] - y[a] - 1;
}

int dp(int now, int k) {
	if(k < 0)
		return -inf;
	
	if(vis[now][k])
		return val[now][k];
	
	int res = 1;
	
	for(int to = 1; to <= n; ++to)
		if(to != now) {
			int cost = dis(to, now);
			
			if(cost > k)
				continue;
			
			bemax(res, dp(to, k - cost) + cost + 1);
		}
	
	vis[now][k] = true;
	val[now][k] = res;
	
	return res;
}

int main() {
	scanf("%d%d", &n, &k);
	
	for(int i = 1; i <= n; ++i)
		scanf("%d%d", &x[i], &y[i]);
	
	for(int i = 1; i <= n; ++i)
		for(int j = 0; j <= k; ++j)
			bemax(ans, dp(i, j) + k - j);
	
	printf("%d\n", ans);
	
	return 0;
}