HDU：4521 小明系列问题——小明序列（LIS-n*logn解法变形+思维+技巧）

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

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 2951    Accepted Submission(s): 924

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
  
  

 

Source

2013腾讯编程马拉松初赛第四场（3月24日）

 

Recommend

liuyiding   |   We have carefully selected several similar problems for you:   5831  5830  5829  5828  5827 

 

题目大意：求一个最长的小明序列长度。

解题思路：这道题是LIS的变形题，因为所选的序列中有间隔要求，所以扩展一下经典LIS算法，原算法的间距是0，这里间距是d，那么我们可以开一个dp数组记录每个数没插入辅助数组时能查到第几位，且这个辅助数组肯定有一个最值是最长小明序列长度（因为每扫描一个新元素，要么插入c序列末尾，要么替换c序列原有的元素，当插入末尾时，dp[i]就更新了此时的小明序列长度），c辅助数组，记录更新小明序列。

代码如下：

#include <cstdio>
#include <algorithm>
using namespace std;
#define INF 0x3f3f3f3f
int dp[100010];//记录之前那个数能插到哪，或者说当前小明序列的元素个数 
int a[100010];
int c[100010];//辅助数组，记录插入的数 
int n,d;
int erfen(int x)
{
	int l=1,r=n;
	int mid,pos;
	while(l<=r)
	{
		mid=(l+r)/2;
		if(c[mid]>=x)//大于等于 1 2  2 3查找2时，最终返回的角标是2因为所找的小明序列是严格单调递增的 
		{
			pos=mid;
			r=mid-1;
		}
		else
		{
			l=mid+1;
		}
	}
	return pos;
}
int main()
{
	while(scanf("%d%d",&n,&d)!=EOF)
	{
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&a[i]);
			c[i]=INF;//初始化辅助数组 
		}
		int ans=0;
		for(int i=1;i<=n;i++)
		{
			dp[i]=erfen(a[i]);//找到当前能插入c小明序列的位置 
			ans=max(ans,dp[i]);//更新ans 
			int id=i-d;
			if(id>=1)//不能越界 
			{
				if(c[dp[id]]>a[id])//更新c数组 
				{
					c[dp[id]]=a[id];
				}
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}