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

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 2346    Accepted Submission(s): 734

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 01 25 13 4 5 1 25 23 4 5 1 2

 

Sample Output

221

 

这题可以用(nlogn)的最长上升子序列算法做，经典的nlogn算法是间隔大于0的算法，那么这题就是这个算法的扩展，取出的每个数的下标差大于d。那么我们只要算到i的时候，先记录下已经形成的c[]函数第一个大于等于a[i]的数的下标，然后用dp[j]记录下来，然后在算到i+d的时候，把它更新进去，那么对于每个数j就避免了j-d,j-d+2...j-1的影响。

#include<iostream>#include<stdio.h>#include<stdlib.h>#include<string.h>#include<math.h>#include<vector>#include<map>#include<set>#include<queue>#include<stack>#include<string>#include<algorithm>using namespace std;#define ll long long#define inf 10000050#define maxn 100060int a[maxn],dp[maxn],c[maxn];int main(){    int n,m,i,j,len,d,maxx;    while(scanf("%d%d",&n,&d)!=EOF)    {        for(i=1;i<=n;i++){            scanf("%d",&a[i]);            c[i]=inf;        }        maxx=1;        for(i=1;i<=n;i++){            dp[i]=lower_bound(c+1,c+1+n,a[i])-c;            maxx=max(maxx,dp[i]);//更新最大值，因为可能这是最后一个数            j=i-d;            if(j>=1 && c[dp[j] ]>a[j]){  //把j更新进去，dp[j]是在循环到j时算出来的                c[dp[j] ]=a[j];            }        }        printf("%d\n",maxx);    }    return 0;}

另一种写法：

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <vector>#include <queue>#include <set>#include <map>#include <string>#include <cmath>#include <cstdlib>#include <ctime>#include <stack>using namespace std;typedef long long ll;typedef long double ldb;#define inf 99999999#define pi acos(-1.0)#define eps 1e-15#define maxn 100050int dp[maxn],a[maxn],c[maxn];int main(){    int n,m,i,j,len,d;    while(scanf("%d%d",&n,&d)!=EOF)    {        for(i=1;i<=n;i++){            scanf("%d",&a[i]);        }        len=0;        int maxx=0;        for(i=1;i<=n;i++){            dp[i]=lower_bound(c+1,c+1+len,a[i])-c;            maxx=max(maxx,dp[i]);            if(i>d){                j=i-d;                if(dp[j]>len){                    len++;c[len]=a[j];                }                else if(c[dp[j] ]>a[j]) {                    c[dp[j] ]=a[j];                }            }        }        printf("%d\n",maxx);    }    return 0;}

也可以用线段树做：用线段树维护值为1~temp的最大子序列的长度，然后每次循环到i，先不要立刻更新，等k时间后再更新。附：用线段树来优化算普通最长上升子序列的方法有两种，一种是用线段树维护值为1~temp的最大子序列的长度，然后每次查找小于a[i]的最大子序列的长度，这里保证下标是递增的；另一种是用线段树维护下标为1~index的最长子序列的长度，然后先把输入的数按值从小到大排序，每次循环的时候只要找到线段树中下标为1~a[i].index的最长上升子序列的长就行。但是在这题中第二种方法不能用，因为二级排序之后，编号都打乱了，那么就不能有顺序地把i-k-1循环进来，所以只能用第一种。

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <vector>#include <queue>#include <set>#include <map>#include <string>#include <cmath>#include <cstdlib>#include <ctime>#include <stack>using namespace std;#define maxn 100050#define inf 999999999struct node1{    int num,idx;    int len;}a[maxn];bool cmp(node1 a,node1 b){    if(a.num==b.num)return a.idx>b.idx;    return a.num<b.num;}struct node{    int l,r,maxnum;}b[4*maxn];void build(int l,int r,int i){    int mid;    b[i].l=l;b[i].r=r;b[i].maxnum=0;    if(l==r){return;}    mid=(l+r)/2;    build(l,mid,i*2);    build(mid+1,r,i*2+1);}void update(int idx,int num,int i){    int mid;    if(b[i].l==idx && b[i].r==idx){        b[i].maxnum=num;        return;    }    mid=(b[i].l+b[i].r)/2;    if(idx<=mid)update(idx,num,i*2);    else{        update(idx,num,i*2+1);    }    b[i].maxnum=max(b[i*2].maxnum,b[i*2+1].maxnum);}int question(int l,int r,int i){    int mid;    if(b[i].l==l && b[i].r==r){        return b[i].maxnum;    }    mid=(b[i].l+b[i].r)/2;    if(r<=mid)return question(l,r,i*2);    else if(l>mid)return question(l,r,i*2+1);    else{        return max(question(l,mid,i*2),question(mid+1,r,i*2+1));    }}int main(){    int n,m,i,j,T,d;    while(scanf("%d%d",&n,&d)!=EOF)    {        int temp=0;        for(i=1;i<=n;i++){            scanf("%d",&a[i].num);            a[i].num++;            temp=max(temp,a[i].num);        }        build(1,temp,1);        int maxx=1;        for(i=1;i<=n;i++){            if(a[i].num==1)a[i].len=1;            else a[i].len=question(1,a[i].num-1,1)+1;            if(i-d>0){                j=i-d;                update(a[j].num,a[j].len,1);            }            maxx=max(maxx,a[i].len);        }        printf("%d\n",maxx);    }    return 0;}