P4093-[HEOI2016/TJOI2016]序列【CDQ分治,树状数组】
生活随笔
收集整理的這篇文章主要介紹了
P4093-[HEOI2016/TJOI2016]序列【CDQ分治,树状数组】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/P4093
題目大意
nnn個數字,每次有一個數字可能和原序列不同,但最多只有一個不同。
求所有情況下都滿足的最長不降子序列
解題思路
定義maximax_imaxi?表示位置iii的最大數,minimin_imini?表示位置iii的最小數。
然后有轉移方程fi=max{fj}+1(ai≥maxj,mini≥aj)f_i=max\{f_j\}+1(a_i\geq max_j,min_i\geq a_j)fi?=max{fj?}+1(ai?≥maxj?,mini?≥aj?)
可以發現這是一個三維偏序,用cdq分治+樹狀數組即可。
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #define lowbit(x) (x&-x) using namespace std; const int N=1e5+10; int n,m,a[N],b[N],s[N]; int t[N],p[N],f[N],ans; bool cmp1(int x,int y) {return b[x]<b[y];} bool cmp2(int x,int y) {return a[x]<a[y];} void Change(int x,int val){while(x<=n){t[x]=max(t[x],val);x+=lowbit(x);}return; } void Clear(int x){while(x<=n){t[x]=0;x+=lowbit(x);}return; } int Ask(int x){int ans=0;while(x){ans=max(ans,t[x]);x-=lowbit(x);}return ans; } void cdq(int l,int r){if(l==r){f[l]=max(f[l],1);return;}int mid=(l+r)>>1;cdq(l,mid);for(int i=l;i<=r;i++)p[i]=i;sort(p+l,p+mid+1,cmp1);sort(p+mid+1,p+r+1,cmp2);int k=l;for(int i=mid+1;i<=r;i++){while(k<=mid&&b[p[k]]<=a[p[i]])Change(a[p[k]],f[p[k]]),k++;f[p[i]]=max(f[p[i]],Ask(s[p[i]])+1);}for(int i=l;i<=mid;i++)Clear(a[i]);cdq(mid+1,r);return; } int main() {scanf("%d%d",&n,&m);for(int i=1;i<=n;i++)scanf("%d",&a[i]),b[i]=s[i]=a[i];for(int i=1;i<=m;i++){int x,y;scanf("%d%d",&x,&y);b[x]=max(b[x],y);s[x]=min(s[x],y);}cdq(1,n);for(int i=1;i<=n;i++)ans=max(ans,f[i]);printf("%d",ans); }總結
以上是生活随笔為你收集整理的P4093-[HEOI2016/TJOI2016]序列【CDQ分治,树状数组】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 勒布朗詹姆斯个人详细资料
- 下一篇: 行行重行行什么意思