BZOJ 2555 SubString LCT 后缀自动机

博客介绍了如何在线处理字符串问题,利用后缀自动机和LCT数据结构来动态维护字符串出现次数,特别是在末尾添加字符时,通过LCT更新相关状态的right集合大小。
摘要由CSDN通过智能技术生成

题目大意:给出一个初始字符串,要求支持,在末尾添加一段字符串,查询一个串出现了几次,强制在线。

构造后缀自动机,一个串出现次数等于代表这个串的状态的right集合大小。

然而强制在线,right集合在添加之后会发生改变。
考虑添加一个字母时,只有产生的新状态的所有parent的right集合大小会加1,于是用LCT维护一下就好了。

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 3000005
using namespace std;
namespace Link_Cut_Tree {
    struct Node {
        Node *ch[2],*pa;
        int dir() { return pa->ch[0]==this ? 0 : pa->ch[1]==this ? 1 : -1; }
        int val,add_mark;
        Node();
        void add(int);
        void pushdown();
        void* operator new(size_t) {
            static Node *mempool,*C;
            if(mempool==C) mempool=(C=new Node[1<<20])+(1<<20);
            return C++;
        }
    }*null=new Node();
    Node :: Node() {
        ch[0]=ch[1]=pa=null;
        val=add_mark=0;
    }
    void Node :: add(int x) {
        val+=x;
        add_mark+=x;
        return ;
    }
    void Node :: pushdown() {
        if(add_mark) {
            ch[0]->add(add_mark);
            ch[1]->add(add_mark);
            add_mark=0;
        }
        return ;
    }
    void Rotate(Node* o,int d) {
        Node* k=o->ch[d^1]; int d2;
        o->ch[d^1]=k->ch[d], k->ch[d]->pa=o;
        k->ch[d]=o;
        //o->maintain(), k->maintain();
        if(~(d2=o->dir())) o->pa->ch[d2]=k;
        k->pa=o->pa; o->pa=k;
        return ;
    }
    void To_pushdown(Node* o) {
        static Node* s[N];
        int top=0;
        while(~(o->dir())) s[++top]=o, o=o->pa;
        s[++top]=o;
        while(top) s[top--]->pushdown();
        return ;
    }
    void Splay(Node* o) {
        To_pushdown(o);
        int d;
        while(~(d=o->dir())) {
            if(o->pa->dir()==d) Rotate(o->pa->pa,d^1);
            Rotate(o->pa,d^1);
        }
        return ;
    }
    void Access(Node* o) {
        Node* p=null;
        while(o!=null) {
            Splay(o);
            o->ch[1]=p;
            p=o;
            o=o->pa;
        }
        return ;
    }
    void Cut(Node* x) {
        Access(x), Splay(x);
        x->ch[0]->pa=null;
        x->ch[0]=null;
        return ;
    }
    void Link(Node* x,Node* y) {
        Cut(x);
        x->pa=y;
        return ;
    }
}
#define LCT Link_Cut_Tree
namespace Suffix_Automaton {
    struct Node {
        Node *ch[26],*pa;
        int val;
        LCT::Node* o;
        Node(int _val=0):val(_val),pa(NULL) {
            memset(ch,0,sizeof ch);
            o=new LCT::Node();
        }
        void* operator new(size_t) {
            static Node *mempool,*C;
            if(mempool==C) mempool=(C=new Node[1<<20])+(1<<20);
            return C++;
        }
    }*root=new Node(0),*last=root;
    void Extend(char c) {
        int z=c-'A';
        Node *p=last,*np=new Node(p->val+1);
        last=np;
        while(p && !p->ch[z]) p->ch[z]=np, p=p->pa;
        if(!p) {
            np->pa=root;
            LCT::Link(np->o,root->o);
            LCT::Access(np->o), LCT::Splay(np->o);
            np->o->add(1);
            return ;
        }
        Node* q=p->ch[z];
        if(q->val==p->val+1) {
            np->pa=q;
            LCT::Link(np->o,q->o);
            LCT::Access(np->o), LCT::Splay(np->o);
            np->o->add(1);
            return ;
        }
        Node* nq=new Node(p->val+1);
        memcpy(nq->ch,q->ch,sizeof q->ch);
        LCT::Link(nq->o,q->pa->o);
        LCT::Link(q->o,nq->o);
        LCT::Link(np->o,nq->o);
        nq->pa=q->pa;
        q->pa=np->pa=nq;
        nq->o->val=q->o->val;
        while(p && p->ch[z]==q) p->ch[z]=nq, p=p->pa;
        LCT::Access(np->o), LCT::Splay(np->o);
        np->o->add(1);
        return ;
    }
    void Insert(char s[]) {
        int len=strlen(s);
        for(int i=0;i<len;i++) Extend(s[i]);
        return ;
    }
    int Query(char s[]) {
        Node* o=root;
        int len=strlen(s);
        for(int i=0;i<=len && o;i++)  {
            if(i==len) return LCT::To_pushdown(o->o), o->o->val;
            o=o->ch[s[i]-'A'];
        }
        return 0;
    }
}
#define SAM Suffix_Automaton
void decodeWithMask(char s[],int mask) {
    int len=strlen(s);
    for(int j=0;j<len;j++) {
        mask=(mask*131+j)%len;
        char t=s[j];
        s[j]=s[mask];
        s[mask]=t;
    }
    return ;
}
int main() {
    static char s[N];
    int T;
    scanf("%d%s",&T,s);
    SAM::Insert(s);
    int mask=0;
    while(T--) {
        char mode[10];
        scanf("%s%s",mode,s);
        decodeWithMask(s,mask);
        if(mode[0]=='Q') {
            int ans=SAM::Query(s);
            printf("%d\n",ans);
            mask^=ans;
        }
        else SAM::Insert(s);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值