hdu3450(线段树+DP)

/*
分析：
    囧~，俩小时弄到正确思路、俩小时敲完代码，
跟脑卡了似的，想问题想不进去，第一次想出来的
方法竟然还和题意不对。。。还好1次ac了，囧~~~
    (线段树(&&离散化)+二分)&&DP。
    DP的想法应该挺容易想到的，dp[i]表示以元素
i结尾的，符合条件的子序列，那么dp[i]=∑dp[j](
1<=j&&j<i && abs(num[j]-num[i])<=limit)。
    然后是对于每个i，怎么找出这些j了。由于limit
的限制结果，使得每个dp[i]能更新的是一段区间。那
么，将num从小到大排序，并离散化，然后建立线段树，
问题就转化为了线段树中，基础的：成段更新、按点
查询的问题。


    也可以用树状数组做，只用做一些小调整就行了。
实在是没有状态，树状数组的就暂时不敲了，囧~~~


                                          2012-10-19
*/

#include"stdio.h"
#include"string.h"
#include"stdlib.h"
#include"math.h"
int n,limit;
int dp[100011],num[100011],hash0[100011];
struct B
{
	int n,index;
}C[100011];
struct Seg
{
	int l,r,mid;
	int num;
}T[400011];
void build(int l,int r,int k)
{
	T[k].l=l;
	T[k].r=r;
	T[k].num=0;
	T[k].mid=(l+r)>>1;
	if(l==r)	return ;
	build(l,T[k].mid,2*k);
	build(T[k].mid+1,r,2*k+1);
}
void update(int l,int r,int dir,int k)
{
	if(T[k].l==l && T[k].r==r)	{T[k].num+=dir;T[k].num%=9901;return ;}
	if(r<=T[k].mid)		update(l,r,dir,2*k);
	else if(l>T[k].mid)	update(l,r,dir,2*k+1);
	else
	{
		update(l,T[k].mid,dir,2*k);
		update(T[k].mid+1,r,dir,2*k+1);
	}
}
int find(int aim,int k)
{
	int ans=0;
	if(T[k].l==T[k].r && T[k].l==aim)	return T[k].num%9901;
	if(T[k].l<=aim && aim<=T[k].r)
	{
		ans+=T[k].num;
		if(aim<=T[k].mid)	ans+=find(aim,2*k);
		else				ans+=find(aim,2*k+1);
	}
	return ans%9901;
}
int cmp(const void *a,const void *b)
{
	struct B *c,*d;
	c=(struct B *)a;
	d=(struct B *)b;
	return c->n-d->n;
}
int main()
{
	int ans;
	int i;
	int l,r;
	int low,up,mid;
	while(scanf("%d%d",&n,&limit)!=-1)
	{
		build(1,n,1);
		for(i=1;i<=n;i++)	{scanf("%d",&num[i]);C[i].n=num[i];C[i].index=i;}
		C[0].n=-111111111;
		qsort(C,n+1,sizeof(C[0]),cmp);
		for(i=1;i<=n;i++)	hash0[C[i].index]=i;

		memset(dp,0,sizeof(dp));
		ans=0;
		for(i=1;i<=n;i++)
		{
			dp[i]=find(hash0[i],1);
			ans+=dp[i];
			ans%=9901;
			low=1;up=hash0[i];mid=(low+up)>>1;
			while(low<=up)
			{
				if(C[mid].n-num[i]<-limit)	low=mid+1;
				else						up=mid-1;
				mid=(low+up)>>1;
			}
			l=low;
			low=hash0[i];up=n;mid=(low+up)>>1;
			while(low<=up)
			{
				if(C[mid].n-num[i]<=limit)	low=mid+1;
				else						up=mid-1;
				mid=(low+up)>>1;
			}
			r=low-1;
			update(l,r,(dp[i]+1)%9901,1);
		}
		printf("%d\n",ans);
	}
	return 0;
}