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

              题意:

                       给了一列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;
}