BZOJ3682 Phorni

Description

Phorni 是一个音之妖精,喜欢在你的打字机上跳舞。
一天,阳光映射到刚刚淋浴过小雨的城市上时,Phorni 用魔法分裂出了许多个幻影,从 1 到 n 编号。
她的每一个幻影都站在打出的字符串的一个位置上,多个幻影可以站在同一个位置上。
每一个幻影代表的字符串即为从它站立位置开始的后缀,注意站立位置是从右往左数的。
让我们形式化地描述一下,若第 i 个幻影站在 Pi 上,那么它所代表的字符串就是 S[L-Pi+1…L],其中 L 是字符串 S 的长度。
每一次,她会选一段编号区间 [l..r],而编号在这个区间中的幻影中,字典序最小的一个将跳一支舞,若有多个幻影字典序相同,选编号最小的。
当然由于 Phorni 还会在打字机上跳动,所以有时字符串的前面会加入一个字符。
当然这个打字机是带加密功能的。
字典序的比较:
将两个字符串逐位比较,长度不足的向后补 0 ( 0 小于任何字符) 。直到比出大小或判定相等。
比如 “pho” > “ph” , “pb” > “pab” 。
下标从 1 开始,保证涉及到的所有字符都为小写字母。

 

Input

输入共 m + 3 行,
第一行为三个整数 n, m, len, type, 分别代表幻影个数,操作次数,初始字符串长度。type = 1 时表示所有的字符都经过了加密。
第二行为初始字符串 S 。
第三行,n 个数 Pi ,意义如题面所示。
接下来 m 行,每行表示一个操作。
1. I c 若 type = 0 表示在字符串前面加入第 c + 1 个小写字母,若 type = 1 则表示加入第 (c xor lastans) + 1 个小写字母,lastans 表示上一次的答案,初始为 0 。
2. C x pos 表示第 x 个幻影跳到了从右向左数 pos 的位置上。
3. Q l r 表示询问[l..r] 。

 

Output

对于每个询问操作输出一行,表示去跳舞的幻影编号。

 

Sample Input

3 3 5 0
horni
3 2 5
I 15
C 1 6
Q 1 3

Sample Output

3

HINT

 

样例解释

前2个操作后，三个幻影代表的后缀分别为phorni,ni,horni。

数据范围与约定

数据范围与约定

对于 100% 的数据, 1 ≤ n ≤ 500000, 1 ≤ m ≤ 800000, 1 ≤ Pi ≤ len ≤ 100000。

若 type = 0,保证 I 操作中 0 ≤ c ≤ 25;否则 0 ≤ c xor lastans ≤ 25 。

C 操作中 1 ≤ x ≤ n, 1 ≤ pos ≤ 当前字符串长度。

Q 操作中 1 ≤ l, r ≤ n,I 操作数量约占总操作数量的 1/5,C,Q 操作数量约各占总操作数量的 2/5。

后缀平衡树来在线得到rk，线段树查询即可。

minv[]忘开3倍maxn调哭了233

bzoj不能用srand(time(0))哦！

#include<cstdio>
#include<cctype>
#include<queue>
#include<ctime>
#include<cstring>
#include<algorithm>
#define mid l+r>>1
#define rep(s,t) for(int i=s;i<=t;i++)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=x*10+c-'0';
    return x*f;
}
const int maxn=500010;
typedef long long ll;
int n,m,len,type,lastans;
int ch[maxn][2],h[maxn],c[maxn],cnt,rt;
ll rk[maxn];
void rebuild(int& o,ll l,ll r) {
    if(!o) return;rk[o]=l+r;
    rebuild(ch[o][0],l,mid);rebuild(ch[o][1],mid,r);
}
void rotate(int& o,int d,ll l,ll r) {
    int k=ch[o][d^1];ch[o][d^1]=ch[k][d];ch[k][d]=o;
    o=k;rebuild(o,l,r);
}
int cmp(int x,int y) {return c[x]<c[y]||(c[x]==c[y]&&rk[x-1]<rk[y-1]);}
void insert(int& o,ll l,ll r) {
    if(!o) {o=cnt;rk[o]=l+r;return;}
    if(cmp(cnt,o)) {insert(ch[o][0],l,mid);if(h[ch[o][0]]>h[o]) rotate(o,1,l,r);}
    else {insert(ch[o][1],mid,r);if(h[ch[o][1]]>h[o]) rotate(o,0,l,r);}
}
void insert(int v) {
    c[++cnt]=v;h[cnt]=rand()*rand();
    insert(rt,1,1ll<<61);
}
int minv[maxn*3],p[maxn];
void build(int o,int l,int r) {
    if(l==r) p[l]=read(),minv[o]=l;
    else {
        int M=l+r>>1,lc=o<<1,rc=lc|1;
        build(lc,l,M);build(rc,M+1,r);
        minv[o]=rk[p[minv[lc]]]<=rk[p[minv[rc]]]?minv[lc]:minv[rc];
    }
}
void update(int o,int l,int r,int pos) {
    if(l==r) return;
    int M=l+r>>1,lc=o<<1,rc=lc|1;
    if(pos<=M) update(lc,l,M,pos);
    else update(rc,M+1,r,pos);
    minv[o]=rk[p[minv[lc]]]<=rk[p[minv[rc]]]?minv[lc]:minv[rc];
}
int query(int o,int l,int r,int ql,int qr) {
    if(ql<=l&&r<=qr) return minv[o];
    int M=l+r>>1,lc=o<<1,rc=lc|1;
    if(qr<=M) return query(lc,l,M,ql,qr);
    if(ql>M) return query(rc,M+1,r,ql,qr);
    int v1=query(lc,l,M,ql,qr),v2=query(rc,M+1,r,ql,qr);
    return rk[p[v1]]<=rk[p[v2]]?v1:v2;
}
char s[maxn],cmd[5];
int main() {
    n=read();m=read();len=read();type=read();
    scanf("%s",s);
    rep(1,len) insert(s[len-i]-'a');
    build(1,1,n);
    while(m--) {
        scanf("%s",cmd);
        if(!type) lastans=0;
        if(cmd[0]=='I') insert(read()^lastans),len++;
        else if(cmd[0]=='C') {
            int x=read();p[x]=read();
            update(1,1,n,x);
        }
        else {
            int x=read(),y=read();
            printf("%d\n",lastans=query(1,1,n,x,y));
        }
    }
    return 0;
}