BZOJ 3489: A simple rmq problem

之前应该是做过这题的弱化版，没有强制在线随便就艹过去了

正解应该是和可持久化堆相关的东西，但是我不会啊怎么办QAQ

祭出暴力大法KD-Tree，话说我第二次写KD-Tree胡完什么都没看随便写了一发就艹过去了，这东西真JB好用

首先那个区间内的限制我们很套路地拆掉，把一个位置看作四元组\((pos,pre,nxt,val)\)，分别表示位置，前驱（和这个数相同的前一个位置，后面同理），后继，权值

然后对于询问的区间\([l,r]\)，一个位置可以被统计当且仅当满足\(pos\in [l,r];pre<l;nxt>r\)

容易发现这就是一个高维数点的板子，应该是可以通过CDQ分治再带个可持久化堆撤销来两个\(\log\)做的

emmm，所以我写KD-Tree

#include<cstdio>
#include<iostream>
#include<algorithm>
#define RI register int
#define CI const int&
#define Tp template <typename T>
using namespace std;
const int N=100005;
int D;
struct point
{
    int w[3],val;
    friend inline bool operator < (const point& A,const point& B)
    {
        for (RI i=0;i<3;++i) if (A.w[(D+i)%3]!=B.w[(D+i)%3])
        return A.w[(D+i)%3]<B.w[(D+i)%3]; return 0;
    }
}a[N]; int n,m,x,y,ans,rt,rst[N],pre[N],nxt[N];
class FileInputOutput
{
    private:
        static const int S=1<<21;
        #define tc() (A==B&&(B=(A=Fin)+fread(Fin,1,S,stdin),A==B)?EOF:*A++)
        #define pc(ch) (Ftop!=Fend?*Ftop++=ch:(fwrite(Fout,1,S,stdout),*(Ftop=Fout)++=ch))
        char Fin[S],Fout[S],*A,*B,*Ftop,*Fend; int pt[15];
    public:
        inline FileInputOutput(void) { Ftop=Fout; Fend=Fout+S; }
        Tp inline void read(T& x)
        {
            x=0; char ch; while (!isdigit(ch=tc()));
            while (x=(x<<3)+(x<<1)+(ch&15),isdigit(ch=tc()));
        }
        Tp inline void write(T x)
        {
            RI ptop=0; while (pt[++ptop]=x%10,x/=10);
            while (ptop) pc(pt[ptop--]+48); pc('\n');
        }
        inline void flush(void)
        {
            fwrite(Fout,1,Ftop-Fout,stdout);
        }
        #undef tc
        #undef pc
}F;
class KD_Tree
{
    private:
        struct interval
        {
            int ch[2]; point p,mi,mx; 
        }node[N]; int tot;
        #define lc(x) node[x].ch[0]
        #define rc(x) node[x].ch[1]
        #define P(x) node[x].p
        #define Mi(x) node[x].mi
        #define Mx(x) node[x].mx
        inline void pushup(CI x,CI y)
        {
            RI i; Mx(x).val=max(Mx(x).val,Mx(y).val);
            for (i=0;i<3;++i) Mi(x).w[i]=min(Mi(x).w[i],Mi(y).w[i]);
            for (i=0;i<3;++i) Mx(x).w[i]=max(Mx(x).w[i],Mx(y).w[i]);
        }
        inline bool judge(CI now)
        {
            return Mx(now).val>ans&&Mx(now).w[0]>=x&&Mi(now).w[0]<=y&&Mi(now).w[1]<x&&Mx(now).w[2]>y;
        }
    public:
        inline void build(int& now,CI l=1,CI r=n,CI d=0)
        {
            now=++tot; int mid=l+r>>1; D=d;   nth_element(a+l+1,a+mid+1,a+r+1);
            P(now)=Mi(now)=Mx(now)=a[mid];
            if (l!=mid) build(lc(now),l,mid-1,(d+1)%3),pushup(now,lc(now));
            if (r!=mid) build(rc(now),mid+1,r,(d+1)%3),pushup(now,rc(now));
        }
        inline void query(CI now)
        {
            if (!now||!judge(now)) return;
            if (P(now).w[0]>=x&&P(now).w[0]<=y&&P(now).w[1]<x&&P(now).w[2]>y) ans=max(ans,P(now).val);
            query(lc(now)); query(rc(now));
        }
        #undef lc
        #undef rc
        #undef P
        #undef Mi
        #undef Mx
}KD;
int main()
{
    //freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
    RI i; for (F.read(n),F.read(m),i=1;i<=n;++i)
    F.read(x),pre[i]=rst[x],nxt[rst[x]]=i,rst[x]=i,a[i].val=x;
    for (i=1;i<=n;++i) a[i].w[0]=i,a[i].w[1]=pre[i],a[i].w[2]=nxt[i]?nxt[i]:n+1;
    for (KD.build(rt),i=1;i<=m;++i)
    {
        F.read(x); F.read(y); x=(x+ans)%n+1; y=(y+ans)%n+1;
        if (x>y) swap(x,y); ans=0; KD.query(rt); F.write(ans);
    }
    return F.flush(),0; 
}

标签：3489,CI,undef,KD,simple,mid,int,problem,now
来源： https://www.cnblogs.com/cjjsb/p/12241114.html