[BZOJ2223][Coci 2009]PATULJCI（整体二分+bit）

题目描述

传送门

题解

以前好像做过一样的题，可持久化权值线段树裸题对吧…
整体二分的话也不难想
首先需要找单调性，乍一看是没有单调性的
不过，如果当前二分的答案为mid，这个区间里所有小于等于mid的数都不够的话，mid就一定变大，反之相反
那么每一次就把小于等于mid的数都加到bit里，然后判断二分的方向就可以了
需要注意的一点是，最后当l=r时就不能再加小于等于mid的数了，必须要加入严格等于mid的数，这样才能判断是否有解

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 300005

int n,lim,m;
struct hp{int val,loc,l,r,goal,id,num;}ins[N],q[N];
int C[N],ans[N];

int read()
{
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x;
}
void add(int loc,int val)
{
    for (int i=loc;i<=n;i+=i&(-i))
        C[i]+=val;
}
int query(int loc)
{
    int ans=0;
    for (int i=loc;i>=1;i-=i&(-i))
        ans+=C[i];
    return ans;
}
int cmpq(hp a,hp b)
{
    return a.num<b.num;
}
void dvd(int l,int r,int x,int y,int a,int b)
{
    if (a>b) return;
    if (l==r)
    {
        for (int i=x;i<=y;++i)
            if (ins[i].val==l) add(ins[i].loc,1);
        for (int i=a;i<=b;++i)
        {
            int t=query(q[i].r)-query(q[i].l-1);
            if (t>=(q[i].r-q[i].l+1)/2+1) ans[q[i].id]=l;
            else ans[q[i].id]=-1;
        }
        for (int i=x;i<=y;++i)
            if (ins[i].val==l) add(ins[i].loc,-1);
        return;
    }
    int mid=(l+r)>>1;
    int pmid=0,pa=0,pb=b-a+1;
    for (int i=x;i<=y;++i)
        if (ins[i].val<=mid)
        {
            pmid=i;
            add(ins[i].loc,1);
        }
        else break;
    for (int i=a;i<=b;++i)
    {
        int t=query(q[i].r)-query(q[i].l-1);
        if (t>=q[i].goal) q[i].num=++pa;
        else q[i].goal-=t,q[i].num=++pb;
    }
    for (int i=x;i<=pmid;++i)
        add(ins[i].loc,-1);
    sort(q+a,q+b+1,cmpq);
    dvd(l,mid,x,pmid,a,a+pa-1);
    dvd(mid+1,r,pmid+1,y,a+pa,b);
}
int cmpins(hp a,hp b)
{
    return a.val<b.val;
}
int main()
{
    n=read();lim=read();
    for (int i=1;i<=n;++i)
    {
        ins[i].val=read();
        ins[i].loc=i;
    }
    sort(ins+1,ins+n+1,cmpins);
    m=read();
    for (int i=1;i<=m;++i)
    {
        q[i].l=read();q[i].r=read();
        q[i].id=i;q[i].goal=(q[i].r-q[i].l+1)/2+1;
    }
    dvd(1,lim,1,n,1,m);
    for (int i=1;i<=m;++i)
        if (ans[i]==-1) puts("no");
        else printf("yes %d\n",ans[i]);
}