题目出处: http://acm.hdu.edu.cn/showproblem.php?pid=5064
题意: 给定n个数,求满足以下两个条件的子序列的最大长度:
(1)C 1 <C 2 <C 3 <......<C t ;
(2)C 2 -C 1 <C 3 -C 2 <......<C t -C t-1 .
分析: 解结构一定为为N 1 ,N 1 ,N 1 ,......,N 1 ,N 2 ,N 3 ,......,N t
设dp[i][j]表示以num[i], num[j]结尾的有效序列的长度,则有
dp[i][j] = max(dp[k][i]+ 1) ,,,,,,,k<=i且num[i]-num[k]<=num[j]-num[i]
剪枝一:
考虑非递减序列 num[s],num[k],num[i],num[j] ,num[e],
容易得到 dp[s][i]<=dp[s][i] (例如1,2,3) 即右端点固定 左端点越远越小,
同时因为序列从小到大排列,因此 左端点越远num[i]-num[k]越大,
所以,对于算某一特定的dp[i][j]时num[j]-num[i]是一个定值从k=i开始递减遍历dp[k][i]+ 1 取最大值,待条件不满足时跳出循环。
剪枝二:
因为 若num[i]-num[k]<=num[j]-num[i]成立 则num[i]-num[k]<=num[j+1]-num[i]必成立,
即 dp[i][j]<=dp[i][j+1] (例如 1(num[k]),4(num[i]),5(num[j]),8(num[j+1]))
即 左端点固定 小区间能取到的序列 大区间一定可以取到,
所以,每次计算新的dp[i][j]时k只要在上次的基础上继续向左遍历dp[k][i]+ 1,取最大值即可。
源代码:
1
#include<cstdio>
2
#include<iostream>
3
#include<map>
4
using
namespace
std;
5
int
n,M,num[
3000
],cnt[
3000
],size,dp[
3000
][
3000
];
6
7
int
maxVal(
int
a,
int
b){
return
a>b?
a:b;}
8
9
int
main()
10
{
11
int
t;
12
scanf(
"
%d
"
,&
t);
13
while
(t--
)
14
{
15
scanf(
"
%d%d
"
,&n,&
M);
16
int
i,j,k;
17
size =
0
;
//
序列中有size个不同的数
18
map<
int
,
int
>
mp;
19
map<
int
,
int
>
::iterator it;
20
//
记录序列中某个数i出现了几次,
21
for
(i=
0
;i<n;i++
)
22
{
23
int
a;
24
scanf(
"
%d
"
,&
a);
25
mp[a]++
;
26
}
27
//
用map是可以默认取数据时从小到大
28
//
PS:网络上不用map用数组预先处理的执行时间约600MS,用map约200MS
29
for
(it=mp.begin(); it!=mp.end();it++
)
30
{
31
num[size] = it->
first;
32
cnt[size] = it->
second;
33
size++
;
34
}
35
36
int
res =
1
;
37
for
(i=
0
;i<size;i++
)
38
{
39
dp[i][i] =
cnt[i];
40
res =
maxVal(res, dp[i][i]);
41
k=
i;
42
int
temp = -
1
;
43
for
(j=i+
1
;j<size;j++
)
44
{
45
for
(;k>=
0
;k--
)
46
{
47
if
(num[i]-num[k]<=num[j]-
num[i])
48
{
49
temp = maxVal(dp[k][i] +
1
,temp);
50
}
51
else
52
{
53
break
;
54
}
55
}
56
dp[i][j] =
temp;
57
res =
maxVal(res, dp[i][j]);
58
}
59
}
60
printf(
"
%d\n
"
,res);
61
62
}
63
return
0
;
64
}

