HDU 4521 小明系列问题——小明序列（线段树+树状数组）

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

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)

Total Submission(s): 2499    Accepted Submission(s): 778

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
  
  

 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <vector>
#include <map>
#include <list>
#include <cmath>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#define LL long long
#define INF 0x3f3f3f3f
#define RR freopen("in.txt","r",stdin)
#define WW freopen("out.txt","w",stdout)
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define offset 1000
#define eps 1e-8
#define zero(x) (((x)>0?(x):-(x))<eps)
#define _sign(x) ((x)>eps?1:((x)<-eps?2:0))

using namespace std;
const int maxn = 100010;

int Sub[maxn * 4];
int a[maxn];
int dp[maxn];

void PushUp(int rt)
{
    Sub[rt] = max(Sub[rt<<1],Sub[rt<<1|1]);
}

void Build(int l, int r, int rt)
{
    int mid = (l + r) >> 1;
    Sub[rt] = 0;
    if(l == r)
        return ;
    Build(l, mid, rt<<1);
    Build(mid+1, r, rt<<1|1);
}

void Updata(int l, int r, int rt, int pos, int n)
{
    if(l == r)
    {
        Sub[rt] = max(Sub[rt], n);
        return ;
    }
    int mid = (l + r) >> 1;
    if(pos <= mid)
        Updata(l, mid, rt<<1, pos, n);
    else
        Updata(mid+1, r, rt<<1|1, pos, n);
    PushUp(rt);
}
int Query(int l,int r, int rt, int L, int R)
{
    if(L <= l && r <= R)
    {
        return Sub[rt];
    }
    int mid = (l + r)>>1;
    if(L > mid)
    {
        return Query(mid +1, r, rt<<1|1, L, R);
    }
    else if(R <= mid)
    {
        return Query(l, mid, rt<<1, L, R);
    }
    else
        return max(Query(l, mid, rt<<1, L, R), Query(mid + 1, r, rt<<1|1, L, R));
}

int main()
{
    int n,d;
    while(~scanf("%d%d",&n, &d))
    {
        memset(dp, 0, sizeof(dp));
        int m = 0;
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&a[i]);
            a[i]++;
            m = max(m,a[i]);
        }
        Build(1,m,1);
        int MAX = 0;
        for(int i=1; i<=n; i++)
        {
            if(i - 1 > d)
                Updata(1,m,1,a[i-d-1],dp[i-d-1]);
            if(a[i] > 1)
                dp[i] = Query(1,m,1,1,a[i] - 1) + 1;
            else
                dp[i] = 1;
            MAX = max(MAX, dp[i]);
        }
        printf("%d\n",MAX);
    }
}