链接:http://acm.hdu.edu.cn/showproblem.php?pid=4521
题目中文
这个题我用线段树做的,也有更快的算法,这里只介绍线段树的方式
叶子节点为输入的值
从头开始遍历,线段树天然地排好序,所以不用管大小的问题,设加入了点a[i],那么在一个子序列中这个点最近的下一个点就是a[i+d+1],相当于对于a[i+d+1]来说,当加入了点a[i]时,就可以查询出以a[i+d+1]为最后一个点的子序列的最长长度,于是,实现的时候就插入a[i],查询a[i+d+1],并保存当前结果
刚开始用#define 来定义max函数,结果一个多小时都在超时,总之长经验了,宏定义的函数每次替换的时候都要做解析,如果大量使用到max函数就会效率很低,再也不得用宏定义函数了,坑啊、、、
#include<stdio.h>
#include<string.h>
#define MAXN 100010
int max(int a,int b)//define max会超时
{
return a>b?a:b;
}
struct node
{
int l,r;
int s;
}t[MAXN*4];
int a[MAXN];
int rec[MAXN];
void construct(int l,int r,int p)
{
t[p].l=l,t[p].r=r,t[p].s=0;
if(l==r)
return ;
int m=(l+r)>>1;
construct(l,m,p*2);
construct(m+1,r,p*2+1);
}
void insert(int x,int p,int val)
{
if(t[p].l==t[p].r)
{
if(val>t[p].s) t[p].s=val;
return ;
}
int m=(t[p].l+t[p].r)>>1;
if(x<=m)
insert(x,p*2,val);
else
insert(x,p*2+1,val);
t[p].s=max(t[p*2].s,t[p*2+1].s);
}
int query(int l,int r,int p)
{
if(t[p].l==l&&t[p].r==r)
return t[p].s;
int m=(t[p].l+t[p].r)>>1;
if(r<=m)
return query(l,r,p*2);
else if(l>m)
return query(l,r,p*2+1);
else return max(query(l,m,p*2),query(m+1,r,p*2+1));
}
int main()
{
int n,d,i;
while(scanf("%d%d",&n,&d)!=EOF)
{
int up=0;
for(i=0;i<n;i++)
{
scanf("%d",&a[i]);
if(a[i]>up) up=a[i];
}
construct(0,up,1);
for(i=0;i<=d;i++) rec[i]=1;//这几个点没办法查询其前的子序列
int ans=1;
for(i=0;i<n;i++)
{
insert(a[i],1,rec[i]);
if(i+d+1<n)
{
if(a[i+d+1]!=0)
{
int t=query(0,a[i+d+1]-1,1);
rec[i+d+1]=t+1;
if(t+1>ans) ans=t+1;
}
else rec[i+d+1]=1;
}
else
break;
}
printf("%d\n",ans);
}
return 0;
}