UESTC 1501 - Defense Lines 离散化+线段树优化DP

来源:互联网 发布:php网络硬盘源码 编辑:程序博客网 时间:2024/05/18 01:44

              题意:

                       给了一列N个数(N<=200000)..现在可以删除其中连续的一段..问最长上升子串为多少...

              题解:

                       很容易想到用dp[p][state]..当state =0表示以p结尾的连续串最长上升子串为多少..state=1表示以p结尾有一段被去掉的最长上升子串为多少..

                       state=0的时候很好维护..只和它前一个位置的状态有关..

                       state=1的时候..分为两种情况来看..一个前面已经有一截被拿去了..所以接下来要一直连续..这个更新和上面是一样的..而另一种是当前这个位置前面一截是空的(其实也可以都是连续的..不影响结果)...这里就要靠线段树来快速找比它小的最长连续上升子串为多少...

                       值得注意的是数的范围远远大于N...所以用离散化来优化空间...

 

Program:

#include<iostream>#include<stdio.h>#include<string.h>#include<set>#include <stack>#include<queue>#include<algorithm>#include<cmath>#define oo 1000000007#define ll long long#define pi acos(-1.0)#define MAXN 300005 using namespace std;   int dp[MAXN][2],P[MAXN<<2];void update(int t,int l,int r,int x,int now){      if (l==r)       {            P[now]=max(P[now],x);            return;      }      int mid=l+r>>1;      if (t<=mid) update(t,l,mid,x,now<<1);      if (t>mid)  update(t,mid+1,r,x,now<<1|1);      P[now]=max(P[now<<1],P[now<<1|1]);}int query(int L,int R,int l,int r,int now){      if (L<=l && R>=r) return P[now];      int ans=0,mid=l+r>>1;      if (L<=mid) ans=max(ans,query(L,R,l,mid,now<<1));      if (R>mid)  ans=max(ans,query(L,R,mid+1,r,now<<1|1));      return ans;}int D[MAXN],A[MAXN];int Bsearch(int x,int n){      int l=0,r=n+1,mid;      while (r-l>1)      {             mid=l+r>>1;             if (D[mid]>x) r=mid;                      else l=mid;      }      return l;}int main(){             int cases,i,n,x,ans;       scanf("%d",&cases);      while (cases--)      {               scanf("%d",&n);              for (i=1;i<=n;i++) scanf("%d",&A[i]),D[i]=A[i];              sort(D+1,D+1+n);              memset(dp,0,sizeof(dp));              memset(P,0,sizeof(P));              ans=0;              for (i=1;i<=n;i++)              {                      x=Bsearch(A[i],n);                      if (A[i]>A[i-1])                               dp[i][0]=dp[i-1][0]+1,dp[i][1]=dp[i-1][1]+1;                         else  dp[i][0]=dp[i][1]=1;                      update(x,0,n,dp[i][0],1);                      dp[i][1]=max(dp[i][1],query(0,x-1,0,n,1)+1);                      ans=max(ans,max(dp[i][0],dp[i][1]));                                    }               printf("%d\n",ans);      }       return 0;}