首页 > 其他分享> 文章详细

BZOJ 3489: A simple rmq problem

2020-01-29 20:00:56 阅读：240 来源： 互联网

标签：3489 CI undef KD simple mid int problem now

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

正解应该是和可持久化堆相关的东西，但是我不会啊怎么办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

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

BZOJ 3489: A simple rmq problem