51nod 1455 宝石猎人（记忆化搜索）

1455 宝石猎人

题目来源：  CodeForces

基准时间限制：2 秒 空间限制：131072 KB 分值: 40  难度：4级算法题

 收藏

 关注

苏塞克岛是一个有着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

System Message  (题目提供者)

C++的运行时限为：2000 ms ，空间限制为：131072 KB  示例及语言说明请按这里

 允许其他 AC 的用户查看此代码，分享代码才能查看别人的代码并有机会获得勋章

题解：典型的dp,d[i][j]:走到i且最后一步的步长为j的最大价值。一开始懵逼了，一直纠结于内存，因为怎么想30000*30000都开不出来啊，后来发现其实d的变化不超过250个数，可以证明：假如n是30000，走了n次，每次是上一次的步数减一（这样能保证最小），呢总的次数是

n*(n+1)/2=30000  n<=250,因此二维开差不多300就好了，我这里开500.

#include<set>  
#include<map>         
#include<stack>                
#include<queue>  
#include<vector>        
#include<string>      
#include<math.h>       
#include<stdio.h>                
#include<iostream>                
#include<string.h>                
#include<stdlib.h>        
#include<algorithm>       
#include<functional>        
using namespace std;
typedef long long ll;
#define inf 1000000000           
#define mod 1000000007                 
#define maxn  31005    
#define PI 3.1415926  
#define lowbit(x) (x&-x)     
#define eps 1e-9   
int a[maxn], dp[maxn][500], ans, mx;
int solve(int x, int y)
{
	if (x > mx) return 0;
	if (dp[x][y] != -1) return dp[x][y];
	if (y > 1)
		dp[x][y] = max(dp[x][y], solve(x + y - 1, y - 1) + a[x]);
	dp[x][y] = max(dp[x][y], solve(x + y, y) + a[x]);
	dp[x][y] = max(dp[x][y], solve(x + y + 1, y + 1) + a[x]);
	ans = max(ans, dp[x][y]);
	return dp[x][y];
}
int main(void)
{
	int n, d, i, x, y;
	scanf("%d%d", &n, &d);
	for (i = 1;i <= n;i++)
	{
		scanf("%d", &x);
		a[x]++;mx = max(mx, x);
	}
	memset(dp, -1, sizeof(dp));
	solve(d, d);
	printf("%d\n", ans);
	return 0;
}