1. 程式人生 > >BZOJ 2653 middle | 主席樹

BZOJ 2653 middle | 主席樹

scanf blank lin ret space spa 題解 www. algorithm

題目:

http://www.lydsy.com/JudgeOnline/problem.php?id=2653


題解:

設答案為ans,把大於等於ans的記為1,小於的記為-1,這樣可以知道當前ans是大了還是小了

然後二分答案,就是求最大子段和的問題,根據網上的題解:[b,c]是必選的,然後選[a,b]和[c,d]的最大字段和就行了

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 20010
using namespace std;
int sz,q,root[N],ls[N*20],rs[N*20
],sum[N*20],lm[N*20],rm[N*20],n,m,b[5],lastans; struct node { int p,v; }a[N]; bool cmp(node a,node b) { return a.v<b.v; } void update(int k) { sum[k]=sum[ls[k]]+sum[rs[k]]; lm[k]=max(lm[ls[k]],sum[ls[k]]+lm[rs[k]]); rm[k]=max(rm[rs[k]],sum[rs[k]]+rm[ls[k]]); } void build(int &rt,int
l,int r) { rt=++sz; int mid=l+r>>1; if (l==r) { sum[rt]=lm[rt]=rm[rt]=1; return; } build(ls[rt],l,mid); build(rs[rt],mid+1,r); update(rt); } void insert(int x,int l,int r,int &y,int p,int v) { y=++sz; ls[y]=ls[x]; rs[y]=rs[x];
if (l==r) { sum[y]=lm[y]=rm[y]=v; return; } int mid=(l+r)>>1; if (p<=mid) insert(ls[x],l,mid,ls[y],p,v); else insert(rs[x],mid+1,r,rs[y],p,v); update(y); } int que1(int k,int l,int r,int ll,int rr) { if (ll==l&&r==rr) return sum[k]; int mid=(l+r)>>1; if (rr<=mid) return que1(ls[k],l,mid,ll,rr); else if (ll>mid) return que1(rs[k],mid+1,r,ll,rr); else return que1(ls[k],l,mid,ll,mid)+que1(rs[k],mid+1,r,mid+1,rr); } int que2(int k,int l,int r,int ll,int rr) { if (ll==l&&r==rr) return rm[k]; int mid=(l+r)>>1; if (rr<=mid) return que2(ls[k],l,mid,ll,rr); else if (ll>mid) return que2(rs[k],mid+1,r,ll,rr); else return max(que2(rs[k],mid+1,r,mid+1,rr),que1(rs[k],mid+1,r,mid+1,rr)+que2(ls[k],l,mid,ll,mid)); } int que3(int k,int l,int r,int ll,int rr) { if (ll==l&&r==rr) return lm[k]; int mid=(l+r)>>1; if (rr<=mid) return que3(ls[k],l,mid,ll,rr); else if (ll>mid) return que3(rs[k],mid+1,r,ll,rr); else return max(que3(ls[k],l,mid,ll,mid),que1(ls[k],l,mid,ll,mid)+que3(rs[k],mid+1,r,mid+1,rr)); } bool judge(int k,int a,int b,int c,int d) { int temp(0); if (c-1>b) temp+=que1(root[k],0,n-1,b+1,c-1); temp+=que2(root[k],0,n-1,a,b); temp+=que3(root[k],0,n-1,c,d); return temp>=0; } int main() { scanf("%d",&n); for (int i=0; i<n; i++) { scanf("%d",&a[i].v); a[i].p=i; } sort(a,a+n,cmp); build(root[0],0,n-1); for (int i=1; i<n; i++)insert(root[i-1],0,n-1,root[i],a[i-1].p,-1); scanf("%d",&q); for (int i=1; i<=q; i++) { scanf("%d%d%d%d",&b[0],&b[1],&b[2],&b[3]); for (int j=0; j<4; j++) b[j]=(b[j]+lastans)%n; sort(b,b+4); int l=0,r=n-1,x; while (l<=r) { int mid=(l+r)>>1; if (judge(mid,b[0],b[1],b[2],b[3])) x=mid,l=mid+1; else r=mid-1; } lastans=a[x].v; printf("%d\n",lastans); } }

BZOJ 2653 middle | 主席樹