苏塞克岛是一个有着30001个小岛的群岛,这些小岛沿着一条直线均匀间隔分布,从西到东编号为0到30000。众所周知,这些岛上有很多宝石,在苏塞克岛上总共有n颗宝石,并且第i颗宝石位于岛 pi上。
小法正好到达0号小岛上,他拥有卓越的跳跃能力,能根据以下规则在小岛之间向东重复跳跃:
· 首先,他会从0号岛跳到d号岛
· 此后,他会根据以下规则继续跳跃,L是上一次跳跃的长度,即,如果他上一次跳跃是从岛prev岛cur,L= cur-prev。他可以向东做一次长度为L-1,L或L+1的跳跃。即,他将会跳到岛 (cur + L - 1), (cur + L) 或 (cur + L + 1)(如果这些岛存在)。一次跳跃的长度必须是正数,即,当L=1时,他不能做一次长度为0的跳跃,如果没有有效的目的地,他将会停止跳跃。
小法将会在跳跃的过程中收集到过的岛上的宝石。我们要找到小法能收集的宝石的最大数。
样例解释:在第一个样例中,最优路径是
0
→
10 (+1宝石
)
→
19
→
27 (+2宝石
)
→...
Input
输入的第一行是两个以空格隔开的整数n和d (1 ≤ n, d ≤ 30000),分别表示苏塞克岛上的宝石数量和小法第一次跳跃的长度。 接下来n行表示这些宝石的位置,第i行(1 ≤ i ≤ n)包含一个整pi(d ≤ p1 ≤ p2 ≤ ... ≤ pn ≤ 30000),表示包含第i颗宝石的小岛的编号。
Output
输出小法能收集的宝石的最大数
Input示例
4 10 10 21 27 27
Output示例
3 思路: 动态规划, 若每一步都加一,且d = 1, 则1 +2 + 3 + ... + n = 30000,可以算出n最大值超不过300,即最大步长超不过300. 用dp[i][j]来记录调到i用j步。#include <algorithm> #include <cstring> #include <iostream> using namespace std; const int MAXN = 30100; int val[MAXN]; int dp[MAXN][600]; int main() { int n, d; cin >> n >> d; memset(val, 0, sizeof(val)); int pos; for (int i = 0; i < n; i++) { cin >> pos; val[pos]++; } memset(dp, -1, sizeof(dp)); int result = val[d]; dp[d][300] = val[d]; for (int i = d; i <= 30000; i++) { for (int j = 0; j < 590; j++) { int next = i + j + d - 300; if (dp[i][j] != -1) { if (next <= 30000) { dp[next][j] = max(dp[next][j], dp[i][j] + val[next]); result = max(result, dp[next][j]); } if (next + 1 <= 30000) { dp[next + 1][j + 1] = max(dp[next + 1][j + 1], dp[i][j] + val[next + 1]); result = max(result, dp[next + 1][j + 1]); } if (next - 1 >= 1 && j + d > 0 && next - i > 1 && next - 1 <= 30000) { dp[next - 1][j - 1] = max(dp[next - 1][j - 1], dp[i][j] + val[next - 1]); result = max(result, dp[next - 1][j - 1]); } } } } cout << result << endl; return 0; }