hdu 小明系列问题——小明序列

Problem Description

　　大家都知道小明最喜欢研究跟序列有关的问题了，可是也就因为这样，小明几乎已经玩遍各种序列问题了。可怜的小明苦苦地在各大网站上寻找着新的序列问题，可是找来找去都是自己早已研究过的序列。小明想既然找不到，那就自己来发明一个新的序列问题吧！小明想啊想，终于想出了一个新的序列问题，他欣喜若狂，因为是自己想出来的，于是将其新序列问题命名为“小明序列”。

　　提起小明序列，他给出的定义是这样的：

　　①首先定义S为一个有序序列，S={ A1 , A2 , A3 , ... , An }，n为元素个数 ；

　　②然后定义Sub为S中取出的一个子序列，Sub={ Ai1 , Ai2 , Ai3 , ... , Aim }，m为元素个数 ；

　　③其中Sub满足 Ai1 < Ai2 < Ai3 < ... < Aij-1 < Aij < Aij+1 < ... < Aim ；

　　④同时Sub满足对于任意相连的两个Aij-1与Aij都有 ij - ij-1 > d (1 < j <= m， d为给定的整数)；

　　⑤显然满足这样的Sub子序列会有许许多多，而在取出的这些子序列Sub中，元素个数最多的称为“小明序列”(即m最大的一个Sub子序列)。

　　例如：序列S={2,1,3,4} ，其中d=1；

　　可得“小明序列”的m=2。即Sub={2,3}或者{2,4}或者{1,4}都是“小明序列”。

　　当小明发明了“小明序列”那一刻,情绪非常激动，以至于头脑凌乱，于是他想请你来帮他算算在给定的S序列以及整数d的情况下，“小明序列”中的元素需要多少个呢？

 

Input

　　输入数据多组，处理到文件结束;

　　输入的第一行为两个正整数 n 和 d；(1<=n<=10^5 , 0<=d<=10^5)

　　输入的第二行为n个整数A1 , A2 , A3 , ... , An，表示S序列的n个元素。(0<=Ai<=10^5)

 

Output

　　请对每组数据输出“小明序列”中的元素需要多少个，每组测试数据输出一行。

 

Sample Input

2 0

1 2

5 1

3 4 5 1 2

5 2

3 4 5 1 2

 

Sample Output

2

2

1

 

       四省完了，成绩不堪入目。哎都几乎没勇气再继续搞算法了，真的丢脸呀。一年的努力那一天全都白费了。真心感觉对不起自己。哎我就当那一年我白费了吧，现在我在从头再来。明年我在卷土从来。现在来说说这道题吧！

其实这道是一个线段树，但是我看见网上的大神说这道题可用lis，于是我花了一些时间看看别人的代码。但是那几天比较忙根本就没有什么时间静心看，所以就看了四，五天。今天干敲我Cocos2d的代码，感觉还不错就看了看。哎只能说做什么事情都得用心。我就看了几分钟就看懂了。

其实就是一个更新时间问题，比方说吧：

         正常的最长上升子序列：  1 4 2 3 5 

  第 1  次： 1

  第2 次： 1  4

  第 3 次： 1 2

  第4 次： 1  2  3

。。。。

      现在有一个限制条件就是d，但是我可以发现dp公式是dp(i)=max(dp(j))  j<i-d-1;很简单的一个公式吧！其实我们就会发现 i到i-d其实并不用再i时去更新。其实我只要去看i-d以前的吧！就是这个思想：

                                                

序列：  1 4 2 3  5  6 d =2

  第 1  次： 0

  第2 次： 0  

  第 3 次： 1 

  第4 次：  min(1,4)

         ......

基本上很lis算法是一样的了。大家可以都担心会不有没有距离的插入里面，其实那完全不用担心。你在纸上画就知道了。我就不多说了。

#include<stdio.h>
#include<string.h>
int arr[100005],ca[100005],n,d,dp[100005];
int find(int a)
{
	int l=0,r=n,m;
	while(l<=r)
	{
		m=(l+r+1)/2;
		if(ca[m]<a)
			l=m+1;
		else
			r=m-1;
	}
	return r;
}
int main()
{
	int i,j,mi,maxn;
	while(scanf("%d %d",&n,&d)!=EOF)
	{
		memset(ca,0x3f,sizeof ca);
		maxn=0;
		ca[0]=-1;
		for(i=1;i<=n;i++)
		{
			scanf("%d",arr+i);
			j=find(arr[i]);
			dp[i]=j+1;
			if(i>d)
				ca[dp[i-d]]=ca[dp[i-d]]<arr[i-d]?ca[dp[i-d]]:arr[i-d];
			if(maxn<dp[i])
				maxn=dp[i];
		}
		printf("%d\n",maxn);
	}

	return 0;
}