codeforces #271(div2) E

题意：给出n,d，和一个数列，求一个最长子序列，相邻的两个数绝对值>=d;，输出长度和位置。

思路：dp[i]表示以当前这个数字为终点形成的最大长度，这样的话找位置的时候就是往前找dp[i]-j(1<=j<dp[i])，这里用到线段树来优化

　　对于一个数字A,他的前一个数字的区间在(1,A-d),(A+d,10^15),所以我们按照数字的大小排序，hsah一下，求个区间最大值，

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int N = 100005;

ll h[N],hh[N];
struct node{
    int l,r,Max;
}e[N*4];
map<ll ,int >a;
map<ll ,int >::iterator it;
ll dp[N];

void build(int L,int R,int x){
    e[x].l=L;
    e[x].r=R;
    e[x].Max=0;
    if(L==R) return ;
    int mid=(L+R)>>1;
    build(L,mid,2*x);
    build(mid+1,R,2*x+1);
}

void update(int y,int val,int x){
    if(e[x].l==y&&e[x].r==y) {
        e[x].Max=val;return ;
    }
    int mid=(e[x].l+e[x].r)>>1;
    if(y<=mid) update(y,val,2*x);
    else update(y,val,2*x+1);
    e[x].Max=max(e[x*2].Max,e[2*x+1].Max);
}

ll query(int L,int R,int x){
    if(e[x].l>=L&&e[x].r<=R) {
        return e[x].Max;
    }

    int mid=(e[x].l+e[x].r)>>1;

    if(mid>=R) return query(L,R,2*x);
    else if(L>mid) return query(L,R,2*x+1);
    else return max(query(L,mid,2*x),query(mid+1,R,2*x+1));
    e[x].Max=max(e[x*2].Max,e[2*x+1].Max);
}
int wei[N];
int main(){
    int n,d;
    scanf("%d%d",&n,&d);
    for(int i=1;i<=n;i++){
        scanf("%lld",&h[i]);
        hh[i]=h[i];
    }
    sort(h+1,h+1+n);
    int x=0;
    for(int i=1;i<=n;i++){
        if(a[h[i]]==0) a[h[i]]=++x;
    }
    build(1,x,1);
    for(int i=1;i<=n;i++){
        dp[i]=1;
        it=a.upper_bound(hh[i]-d);//得找hh[i]-d及以前的数造成的最大值,upper指向第一个>该值的位置//
        if(it!=a.begin()){
            it--;
            int xx=it->second;
            ll Max=query(1,xx,1);
            dp[i]=max(dp[i],Max+1);
        }
        it=a.lower_bound(hh[i]+d);//得找hh[i]+d及后面的数字,lower指向第一个>=该值的位置//
        if(it!=a.end()){
            int xx=it->second;
            ll Max=query(xx,x,1);
            dp[i]=max(dp[i],Max+1);
        }
        update(a[hh[i]],dp[i],1);
    }
    ll ss=1;
    int z;
    for(int i=1;i<=n;i++){
        if(dp[i]>ss){
            ss=dp[i];z=i;
        }
    }
    printf("%d
",ss);
    int t=0;
    wei[++t]=z;
    for(int i=n;i>=1;i--){
     // cout<<dp[i]<<" "<<hh[i]<<" "<<
        if(dp[i]==ss-1&&abs(hh[z]-hh[i])>=d){
            wei[++t]=i;
            z=i;
            ss--;
        }
    }
    for(int i=t;i>=1;i--){
        printf("%d ",wei[i]);
    }
    printf("
");
    return 0;
}