【BZOJ3489】A simple rmq problem

Description

因为是OJ上的题，就简单点好了。给出一个长度为n的序列，给出M个询问：在[l,r]之间找到一个在这个区间里只出现过一次的数，并且要求找的这个数尽可能大。如果找不到这样的数，则直接输出0。我会采取一些措施强制在线。

Input

第一行为两个整数N,M。M是询问数，N是序列的长度（N<=100000，M<=200000)

第二行为N个整数，描述这个序列{ai}，其中所有1<=ai<=N

再下面M行，每行两个整数x，y，

询问区间[l,r]由下列规则产生（OIER都知道是怎样的吧>_<)：

l=min(（x+lastans)mod n+1,(y+lastans）mod n+1);

r=max(（x+lastans)mod n+1,(y+lastans）mod n+1);

Lastans表示上一个询问的答案，一开始lastans为0

Output

一共M行，每行给出每个询问的答案。

Sample Input

10 10

6 4 9 10 9 10 9 4 10 4

3 8

10 1

3 4

9 4

8 1

7 8

2 9

1 1

7 3

9 9

Sample Output

4

10

10

0

0

10

0

4

0

4

HINT

注意出题人为了方便，input的第二行最后多了个空格。

2015.6.24新加数据一组，但未重测

Source

by zhzqkkk

因为在集中刷kdtree,正好从qiancl的kdtree系列整合里发现了这个题,然后就发现这题原来也可以不用树套树用kdtree做..

搞成三维kdtree,坐标分别代表 序列中位置,前一个和后一个相同数出现的位置,然后维护最大权值就行了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 100010
#define GET (ch>='0'&&ch<='9')
#define MAXINT 0x3f3f3f3f
#define Dnum 3
using namespace std;
inline void in(int &x)
{
char ch=getchar();x=0;
while (!GET)    ch=getchar();
while (GET) x=x*10+ch-'0',ch=getchar();
}
int n,m,cmp_d,ans,root;
int a,b,l,r;
int pre[MAXN],nxt[MAXN],lst[MAXN],w[MAXN];
struct KDtree
{
int ch[2],d[Dnum],minn[Dnum],maxn[Dnum],val,maxv;
inline void init()  {   for (int i=0;i<Dnum;++i)    minn[i]=maxn[i]=d[i];   }
inline bool operator < (const KDtree& a)const   {   return d[cmp_d]<a.d[cmp_d]; }
}tree[MAXN];
inline void push_up(int rt)
{
tree[rt].maxv=tree[rt].val;
for (int i=0,x=0;i<2;++i)
if ((x=tree[rt].ch[i]))
{
for (int j=0;j<Dnum;++j)
tree[rt].minn[j]=min(tree[rt].minn[j],tree[x].minn[j]),
tree[rt].maxn[j]=max(tree[rt].maxn[j],tree[x].maxn[j]);
tree[rt].maxv=max(tree[rt].maxv,tree[x].maxv);
}
}
int rebuild(int l=1,int r=n,int d=0)
{
cmp_d=d;int mid=(l+r)>>1,nxtd=(d+1)%3;nth_element(tree+l,tree+mid,tree+r+1);
tree[mid].init();
if (l!=mid) tree[mid].ch[0]=rebuild(l,mid-1,nxtd);
if (r!=mid) tree[mid].ch[1]=rebuild(mid+1,r,nxtd);
return push_up(mid),mid;
}
inline bool check(int rt)
{
if (tree[rt].maxn[1]<=r||tree[rt].minn[2]>=l||tree[rt].minn[0]>r||tree[rt].maxn[0]<l)   return 0;
return 1;
}
void query(int rt=root)
{
if (tree[rt].minn[0]>=l&&tree[rt].maxn[0]<=r&&tree[rt].minn[1]>r&&tree[rt].maxn[2]<l)   {   ans=max(ans,tree[rt].maxv);return;  }
if (tree[rt].d[0]>=l&&tree[rt].d[0]<=r&&tree[rt].d[1]>r&&tree[rt].d[2]<l)   ans=max(ans,tree[rt].val);
int ls=tree[rt].ch[0],rs=tree[rt].ch[1];
if (tree[ls].maxv>tree[rs].maxv)
{
if (ls&&tree[ls].maxv>ans&&check(ls))   query(ls);
if (rs&&tree[rs].maxv>ans&&check(rs))   query(rs);
}
else
{
if (rs&&tree[rs].maxv>ans&&check(rs))   query(rs);
if (ls&&tree[ls].maxv>ans&&check(ls))   query(ls);
}
}
int main()
{
in(n);in(m);
for (int i=1;i<=n;i++)  in(w[i]),pre[i]=lst[w[i]],lst[w[i]]=i;
for (int i=1;i<=n;i++)  nxt[pre[i]]=i;nxt[0]=0;
for (int i=1;i<=n;i++)
{
if (!nxt[i])    nxt[i]=n+1;
tree[i].d[0]=i;tree[i].d[1]=nxt[i];tree[i].d[2]=pre[i];tree[i].val=w[i];
}
for (root=rebuild();m;m--)
{
in(a);in(b);l=min((a+ans)%n+1,(b+ans)%n+1);r=max((a+ans)%n+1,(b+ans)%n+1);
ans=0;query();printf("%d\n",ans);
}
}