NOJ 1526 最长等差子序列

题目链接：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;
}