• [BZOJ3682] Phorni


    Description

    有一个长度为 (len) 的字符串和 (n) 个幻影,每个幻影站在字符串的某个字符上,代表从该位置到字符串末尾的一个后缀。要求支持以下三个操作:在字符串开头添加一个字符 改变幻影位置 询问编号在 ([l,r]) 的幻影代表字符串字典序最小的幻影编号是多少 强制在线

    Solution

    后缀平衡树。

    注意后缀平衡树如果要判断两个后缀字典序关系的话并不需要去(log)的查rk数组,直接判断Treap上对应节点的 (tag) 哪个更小就好了。因为根据后缀平衡树的性质保证了 (tag) 互不相同且中序遍历是递增的。

    然后这题再拿线段树维护一下每个幻影就好了

    拿替罪羊树写会跑的飞快吧.... 真的懒得写

    洛谷博客太垃圾了

    计划把前一阵子在洛谷写的博客搬过来

    Code

    #pragma GCC optimize(2)
    #include<bits/stdc++.h>
    using std::min;
    using std::max;
    using std::swap;
    using std::vector;
    typedef double db;
    typedef long long ll;
    #define pb(A) push_back(A)
    #define pii std::pair<int,int>
    #define all(A) A.begin(),A.end()
    #define mp(A,B) std::make_pair(A,B)
    const int N=2e6+5;
    const db inf=1e18;
    const db alph=0.9;
    
    char s[N];
    int n,m,cnt;
    int len,type,p[N];
    
    int getint(){
        int X=0,w=0;char ch=getchar();
        while(!isdigit(ch))w|=ch=='-',ch=getchar();
        while( isdigit(ch))X=X*10+ch-48,ch=getchar();
        if(w) return -X;return X;
    }
    
    namespace treap{
        #define ls ch[x][0]
        #define rs ch[x][1]
    
        int tot,root,cnt,ch[N][2],d[N];
        int sze[N],prio[N];db tag[N];
    
        void pushup(int x){
            sze[x]=sze[ls]+sze[rs]+1;
            d[x]=max(d[ls],d[rs])+1;
        }
    
        bool cmp(int x,int y){
            return s[x]<s[y] or s[x]==s[y] and tag[x-1]<tag[y-1];
        }
    
        void dfs(int now,db l,db r){
            if(!now) return;
            db mid=(l+r)/2;tag[now]=mid;
            dfs(ch[now][0],l,mid);dfs(ch[now][1],mid,r);
            pushup(now);
        }
    
        void rotate(int &x,int d,db l,db r){
            int y=ch[x][d],z=ch[y][d^1];
            ch[x][d]=z;ch[y][d^1]=x;
            x=y; dfs(x,l,r);
        }
    
        void insert(int &x,db l,db r){
            if(!x){
                x=++tot;tag[x]=(l+r)/2;
                sze[x]=1;prio[x]=rand();
                return;
            } db mid=(l+r)/2;int d=cmp(x,len);
            if(!d) insert(ls,l,mid);
            else insert(rs,mid,r);
            if(prio[ch[x][d]]<prio[x]) rotate(x,d,l,r);
            pushup(x);
        }
        #undef ls
        #undef rs
    }
    
    namespace seg{
        #define ls x<<1
        #define rs x<<1|1
        #define lss ls,l,mid,ql
        #define rss rs,mid+1,r,ql
    
        int mn[N];
    
        void pushup(int x){
            mn[x]=treap::tag[p[mn[ls]]]<=treap::tag[p[mn[rs]]]?mn[ls]:mn[rs];
        }
    
        void build(int x,int l,int r){
            if(l==r) return mn[x]=l,void();
            int mid=l+r>>1;
            build(ls,l,mid),build(rs,mid+1,r);
            pushup(x);
        }
    
        void modify(int x,int l,int r,int ql){
            if(l==r) return;
            int mid=l+r>>1;
            ql<=mid?modify(lss):modify(rss);
            pushup(x);
        }
    
        int query(int x,int l,int r,int ql,int qr){
            if(ql<=l and r<=qr) return mn[x];
            int mid=l+r>>1;
            if(qr<=mid) return query(lss,qr);
            if(ql>mid) return query(rss,qr);
            int a=query(lss,qr),b=query(rss,qr);
            return treap::tag[p[a]]<=treap::tag[p[b]]?a:b;
        }
        #undef ls
        #undef rs
    }
    
    signed main(){
        srand(20020619);
        n=getint(),m=getint(),len=getint(),type=getint();
        scanf("%s",s+1);
        int mxlen=len,lasans=0;
        std::reverse(s+1,s+1+len);
        for(len=1;len<=mxlen;len++) treap::insert(treap::root,1,inf);len--;
        for(int i=1;i<=n;i++) p[i]=getint();
        seg::build(1,1,n);
        while(m--){
            char ch[5];scanf("%s",ch);
            if(ch[0]=='I'){
                int x=(getint()^lasans*type)+1;
                s[++len]=x+'a'-1; treap::insert(treap::root,1,inf);
            } else if(ch[0]=='C'){
                int x=getint(),pos=getint();
                p[x]=pos;seg::modify(1,1,n,x);
            } else{
                int l=getint(),r=getint();
                printf("%d
    ",lasans=seg::query(1,1,n,l,r));
            }
        } 
        return 0;
    }
    
    
  • 相关阅读:
    ETL开发
    ETL 运行维护
    软件开发的核心
    索引的本质
    也谈谈学习
    oracle内部结构
    [kuangbin带你飞]专题十六 KMP & 扩展KMP & Manacher H
    [kuangbin带你飞]专题十六 KMP & 扩展KMP & Manacher :G
    [kuangbin带你飞]专题十六 KMP & 扩展KMP & Manacher F
    [kuangbin带你飞]专题十六 KMP & 扩展KMP & Manacher E
  • 原文地址:https://www.cnblogs.com/YoungNeal/p/10300520.html
Copyright © 2020-2023  润新知