题目链接:http://ac.nbutoj.com/Problem/view.xhtml?id=1526
思路:
把所有公差离散出来(理论上公差个数应该有(n+1)*n/2 , 不过实际最多2000个)
dp[i][j] 表示 公差为j ,以i下标为结尾的子序列 最优解
因为dp只需要考虑这次和上一次,所以用滚动数组优化内存(即公差是N*N,也可以开出来)
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<set>
#include<map>
using namespace std;
#define N 1005
int abs(int x){return x>0?x:-x;}
int a[N], dp[2][N*2], n; //理论上应该是dp[2][N*N];
map<int, int>mymap;
set<int>myset;
set<int>::iterator p;
int main(){
int i, j, fir = 0;
while(~scanf("%d",&n)){
myset.clear();
mymap.clear();
for(i = 1; i <= n; i++)scanf("%d",&a[i]);
for(i = 1; i < n; i++)
for(j = i+1; j <= n; j++)
myset.insert(abs(a[i]-a[j]));
int top = 0;
for(p = myset.begin(); p != myset.end(); p++)
mymap.insert(make_pair(*p, top++));
int cur = 0, ans = 1;
for(j=0;j<top;j++)dp[cur][j] = 1;
for(i = 2; i <= n; i++)
{
cur ^= 1;
for(j=0;j<top;j++)dp[cur][j] = dp[cur^1][j];
for(j = i-1; j ;j--)
{
int cha = mymap.find(abs(a[i]-a[j])) ->second;
dp[cur][cha] = max(dp[cur][cha], dp[cur^1][cha]+1);
/* 错误剪枝,此题数据水可以这样操作
if(dp[cur][cha]>ans) ans = dp[cur][cha];
else if(ans >= j+3)break; */
}
}
if(fir++)puts("");
printf("%d",ans);
}
return 0;
}