HDU - 5324（分治+树状数组）

本题目的原来意思是，给定两个长度为n的数组，L，R，要求一个子序列（可以不连续）使的L递减，R递增。

分析：

加上下标递增，两个维度增，一个减，那么考虑  d[ i ]代表以i为起点的串往后找能得到的最大长度，用分治方法更新最有值。那么下面说一下，怎么分治维护最优。

首先，数组原顺序保持不变，就是下标递增，对于l - > r区间的每个点的最优值，先求出右半区间每个点的最优值，然后把左右区间分别以R值增序排序，那么就可以用树状数组，来更新左半区间的d值，这样消除了右半区间对左半区间的影响，然后再求左半区间的最优值即可。

#include <set>
#include <cstdio>
#include <vector>
#include <iostream>
#include <cstring>
#include <queue>
#include <map>
#include <string>
#include<iomanip>
#include <algorithm>
using namespace std;
#define rep(i,n) for(int i=0;i<(int)n;i++)
#define rep1(i,x,y) for(int i=x;i<=(int)y;i++)
#define lowbit(x) (x&(-x))
#define pii(x,y) make_pair(x,y)
typedef pair<int,int> pii;
const int N = 51010;
const int inf = 0x3f3f3f3f;

struct node{
    int l,r,id;
    bool operator <(const node& rhs)const{
      return r<rhs.r;
    }
}a[N],tem[N];
struct San :vector<int>{
    void prepare(){
         sort(begin(),end());
         erase(unique(begin(),end()),end());
    }
    int get(int x){
        return (int)(lower_bound(begin(),end(),x)-begin())+1;
    }
}sl,sr;
struct Bit{
     pii c[N];
     void modify(int x,pii v){
          while(x<N){
              c[x]=v; x+=lowbit(x);
          }
     }
     void update(int x,pii v){
          while(x<N){
              c[x]=max(v,c[x]); x+=lowbit(x);
          }
     }
     pii query(int x){
         pii res = pii(0,-N);
          while(x>0){
              res=max(res,c[x]); x-=lowbit(x);
          }
         return res;
     }
}bit;

pii d[N];
void solve(int l,int r){
   if(l == r) return ;
   int m=(l+r)>>1;
   solve(m+1,r);
   rep1(i,l,r) tem[i]=a[i];
   sort(tem+l,tem+m+1);
   sort(tem+m+1,tem+r+1);
   int pr=r;
   for(int i=m;i>=l;i--){
        while(pr>m && tem[pr].r>=tem[i].r){
            bit.update(tem[pr].l,pii(d[tem[pr].id].first+1,-tem[pr].id));
            pr--;
        }
        pii can=bit.query(tem[i].l);
        d[tem[i].id]=max(d[tem[i].id],can);
   }
   for(int i=m+1;i<=r;i++) bit.modify(tem[i].l,pii(0,-N));
   solve(l,m);
}
int n;
void read(){
     sl.clear(); sr.clear();
     rep1(i,1,n){
         scanf("%d",&a[i].l);
         sl.push_back(a[i].l);
     }
     rep1(i,1,n){
         scanf("%d",&a[i].r);
         sr.push_back(a[i].r);
     }
     sl.prepare(); sr.prepare();
     rep1(i,1,n){
         a[i].l=sl.get(a[i].l);
         a[i].r=sr.get(a[i].r);
         a[i].id=i;
     }
     rep(i,N) bit.c[i]=pii(0,-N);
     for(int i=1;i<=n;i++){
        d[i]=pii(1,-N);
     }
     solve(1,n);
     int ans=-1,cnt=0,p;
     rep1(i,1,n){
       if(d[i].first>ans) {
          ans=d[i].first; p=i;
       }
     }
     printf("%d\n",ans);
     for(int i=p;i!=(-N);i=d[i].second){
          i=abs(i);
          if(cnt++) printf(" ");
          printf("%d",i);
     }
     printf("\n");
}
int main()
{
    while(scanf("%d",&n)==1){
       read();
    }
    return 0;
}
/*
4
3 2 1 1
2 3 1 1
*/