codevs 3289 花匠

  1 #include<cstdio>
  2 #include<algorithm>
  3 #include<iostream>
  4 using namespace std;
  5 const int maxn=100010;
  6 
  7 int n,maxfj[maxn*4],maxfo[maxn*4],f[maxn],p[maxn],maxpj[maxn*4],maxpo[maxn*4];
  8 
  9 struct u
 10 {
 11     int v,r;
 12     bool operator <(const u &rhs) const
 13     {
 14         return v<rhs.v;
 15     }
 16 }A[maxn];
 17 
 18 bool cmp(u a,u b)
 19 {
 20     return a.r<b.r;
 21 }
 22 
 23 int w,v;
 24 void updatefj(int o,int L,int R)
 25 {
 26     if(L==R) maxfj[o]=max(maxfj[o],v);
 27     else
 28     {
 29         int M=(L+R)/2;
 30         if(w<=M) updatefj(o*2,L,M); else updatefj(o*2+1,M+1,R);
 31         maxfj[o]=max(maxfj[o*2],maxfj[o*2+1]);
 32     }
 33 }
 34 
 35 int y1,y2,ans;
 36 void queryfj(int o,int L,int R)
 37 {
 38     if(y1<=L && R<=y2) ans=max(ans,maxfj[o]);
 39     else
 40     {
 41         int M=(L+R)/2;
 42         if(y1<=M) queryfj(o*2,L,M);
 43         if(y2>M) queryfj(o*2+1,M+1,R);
 44     }
 45 }
 46 void updatefo(int o,int L,int R)
 47 {
 48     if(L==R) maxfo[o]=max(maxfo[o],v);
 49     else
 50     {
 51         int M=(L+R)/2;
 52         if(w<=M) updatefo(o*2,L,M); else updatefo(o*2+1,M+1,R);
 53         maxfo[o]=max(maxfo[o*2],maxfo[o*2+1]);
 54     }
 55 }
 56 
 57 void queryfo(int o,int L,int R)
 58 {
 59     if(y1<=L && R<=y2) ans=max(ans,maxfo[o]);
 60     else
 61     {
 62         int M=(L+R)/2;
 63         if(y1<=M) queryfo(o*2,L,M);
 64         if(y2>M) queryfo(o*2+1,M+1,R);
 65     }
 66 }
 67 void updatepj(int o,int L,int R)
 68 {
 69     if(L==R) maxpj[o]=max(maxpj[o],v);
 70     else
 71     {
 72         int M=(L+R)/2;
 73         if(w<=M) updatepj(o*2,L,M); else updatepj(o*2+1,M+1,R);
 74         maxpj[o]=max(maxpj[o*2],maxpj[o*2+1]);
 75     }
 76 }
 77 
 78 void querypj(int o,int L,int R)
 79 {
 80     if(y1<=L && R<=y2) ans=max(ans,maxpj[o]);
 81     else
 82     {
 83         int M=(L+R)/2;
 84         if(y1<=M) querypj(o*2,L,M);
 85         if(y2>M) querypj(o*2+1,M+1,R);
 86     }
 87 }
 88 
 89 void updatepo(int o,int L,int R)
 90 {
 91     if(L==R) maxpo[o]=max(maxpo[o],v);
 92     else
 93     {
 94         int M=(L+R)/2;
 95         if(w<=M) updatepo(o*2,L,M); else updatepo(o*2+1,M+1,R);
 96         maxpo[o]=max(maxpo[o*2],maxpo[o*2+1]);
 97     }
 98 }
 99 
100 void querypo(int o,int L,int R)
101 {
102     if(y1<=L && R<=y2) ans=max(ans,maxpo[o]);
103     else
104     {
105         int M=(L+R)/2;
106         if(y1<=M) querypo(o*2,L,M);
107         if(y2>M) querypo(o*2+1,M+1,R);
108     }
109 }
110 
111 int main()
112 {
113     cin>>n;
114     for(int i=1;i<=n;i++) cin>>A[i].v,A[i].r=i;
115     //离散化 
116     sort(A+1,A+n+1);
117     for(int i=1,nw=1;i<=n;i++)
118     {
119         f[i]=1,p[i]=1;//顺便做的f[],p[]初始化 
120         if(A[i+1].v!=A[i].v) nw++,A[i].v=nw-1;
121         else A[i].v=nw;
122     } 
123     sort(A+1,A+n+1,cmp);
124      
125     w=A[1].v,v=f[1];
126     updatefj(1,1,n);
127     updatepj(1,1,n);
128 
129     for(int i=2;i<=n;i++)
130     {
131         y1=1,y2=A[i].v-1,ans=0;
132         if(y2>=y1) queryfj(1,1,n);//注意，因为有A[i]-1，所以要判断区间的存在 
133         y1=A[i].v+1,y2=n;
134         if(y2>=y1) queryfo(1,1,n);
135         f[i]=ans+1;
136         v=f[i],w=A[i].v;
137         if(f[i]%2==1) updatefj(1,1,n);
138         else updatefo(1,1,n);
139         
140         y1=1,y2=A[i].v-1,ans=0;
141         if(y2>=y1) querypo(1,1,n);
142         y1=A[i].v+1,y2=n;
143         if(y2>=y1) querypj(1,1,n);
144         p[i]=ans+1;
145         v=p[i],w=A[i].v;
146         if(p[i]%2==1) updatepj(1,1,n);
147         else updatepo(1,1,n);
148     }
149     
150     cout<<max(f[n],p[n]);
151     return 0;
152 }