这道题的数据范围为1<=n<=10^5,若用普通的O(n^2)算法肯定超时,所以此处要用到优化后的LIS。
但这题又不是单纯的LIS,因为它加了一个限制条件。
算法思想:此题要用到一个附加数组c[i],其中c[i]表示长度为i的子序列的最小尾部,但因为此题有个限制条件:相邻两个数字的距离大于d,因此在读第i个元素时,c数组其实还只更新到i-d,这也就是处理此题的技巧了。
dp[i]:以第i号元素结尾的最长子序列
c[i]:长度为i的子序列的最小尾部
a[i]:集合元素
#include<iostream>
#include<cstdio>
using namespace std;
const int N=100005;
int a[N],c[N],dp[N];
int bsearch(int low,int up,int key)
{
while(up>=low)
{
int mid=(low+up)>>1;
if(key>c[mid])
low=mid+1;
else
up=mid-1;
}
return low-1;
}
int main()
{
int n,d,ans;
while(~scanf("%d%d",&n,&d))
{
for(int i=1;i<=n;i++)
{
scanf("%d",a+i);
c[i]=10000000;
}
ans=0;
for(int i=1;i<=n;i++)
{
dp[i]=bsearch(1,n,a[i])+1;
if(dp[i]>ans)
ans=dp[i];
int j=i-d;
if(j>0&&c[dp[j]]>a[j])
c[dp[j]]=a[j];
}
printf("%d\n",ans);
}
return 0;
}