BZOJ 2653 middle | 主席树

题目:

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*],rs[N*],sum[N*],lm[N*],rm[N*],n,m,b[],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>>;
    if (l==r)
    {
        sum[rt]=lm[rt]=rm[rt]=;
        return;
    }
    build(ls[rt],l,mid);
    build(rs[rt],mid+,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)>>;
    if (p<=mid) insert(ls[x],l,mid,ls[y],p,v);
    else insert(rs[x],mid+,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)>>;
    if (rr<=mid) return que1(ls[k],l,mid,ll,rr);
    else if (ll>mid) return que1(rs[k],mid+,r,ll,rr);
    else return que1(ls[k],l,mid,ll,mid)+que1(rs[k],mid+,r,mid+,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)>>;
    if (rr<=mid) return que2(ls[k],l,mid,ll,rr);
    else if (ll>mid) return que2(rs[k],mid+,r,ll,rr);
    else return max(que2(rs[k],mid+,r,mid+,rr),que1(rs[k],mid+,r,mid+,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)>>;
    if (rr<=mid) return que3(ls[k],l,mid,ll,rr);
    else if (ll>mid) return que3(rs[k],mid+,r,ll,rr);
    else return max(que3(ls[k],l,mid,ll,mid),que1(ls[k],l,mid,ll,mid)+que3(rs[k],mid+,r,mid+,rr));
}
bool judge(int k,int a,int b,int c,int d)
{
    int temp();
    if (c->b) temp+=que1(root[k],,n-,b+,c-);
    temp+=que2(root[k],,n-,a,b);
    temp+=que3(root[k],,n-,c,d);
    return temp>=;
}
int main()
{
    scanf("%d",&n);
    for (int i=; i<n; i++)
    {
        scanf("%d",&a[i].v);
        a[i].p=i;
    }
    sort(a,a+n,cmp);
    build(root[],,n-);
    for (int i=; i<n; i++)insert(root[i-],,n-,root[i],a[i-].p,-);
    scanf("%d",&q);
    for (int i=; i<=q; i++)
    {
        scanf("%d%d%d%d",&b[],&b[],&b[],&b[]);
        for (int j=; j<; j++) b[j]=(b[j]+lastans)%n;
        sort(b,b+);
        int l=,r=n-,x;
        while (l<=r)
        {
            int mid=(l+r)>>;
            if (judge(mid,b[],b[],b[],b[])) x=mid,l=mid+;
            else r=mid-;
        }
        lastans=a[x].v;
        printf("%d\n",lastans);
    }
}