（SAM）HDU4270 Dynamic Lover

数了一下，在HDU上总共提交了52次……基本上清一色的都是TLE，其中绝大多数的原因是memcpy时没有按照给定的范围，而是莫名其妙的只memcpy了一部分orz……

利用性质：在SAM上沿匹配边不断走就能得到所有子串。

故对题目中的3中操作可以对应于在SAM上进行相应操作。

1、直接继续在SAM上添加字符串即可

2、从根节点开始走最多x（x为给定长度）步，每次向可走的最小处（e.g. a不可走，b可走，就走b）转移。走x步，或者遇到包含当前串的后缀串的结点时结束该过程（某结点是否包含当前串的后缀，只需要在此之前从最后添加字符新得到的结点开始沿着slink边走到根节点并一路标记即可），即得到字典序最小的串。

3、对于删除操作，采用的策略是对每个结点额外记录一个属性，记录该结点是否被删除过。并且记录添加字符串到SAM过程中，某个位置在SAM中的结点编号。由于SAM中某个结点表示一个后缀区间，在往SAM中做添加操作时拆开的点，其对应字符串的结尾是相同的，故其是否被删除过的状态是相同的，为了减少修改操作，将其指向同一个标记数组中的位置即可（用指针、数组记录下标均可）。被删除的点当作无法转移一样处理即可。每次删除只需将连续L个结点修改状态即可，未被删除的结点必定存在某个endpos，满足对应的字符串不包含任何一个被删除的位置。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
#include <list>
#include <map>
#include <string>
#include <cstring>
#include <stack>
#include <queue>
#include <cmath>
#include <ctime>
#include <bitset>
#include <utility>
#include <complex>
#include <assert.h>
#include <limits.h>
#include <iomanip>
//#include <bits/stdc++.h>
using namespace std;
#define rank rankk
#define mp make_pair
#define pb push_back
#define xo(a,b) ((b)&1?(a):0)
#define tm tmp

//#define LL ll
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef long long ll;
typedef pair<ll,int> pli;
typedef pair<int,ll>pil;
typedef pair<ll,ll> pll;
//const double INF=1e20;
const int INF=0x3f3f3f3f;
//const int INF= 0x7fffffff;
//const ll INF=0x3f3f3f3f3f3f3f3fll;
const ll INFF=0x3f3f3f3f3f3f3fll;
//const ll INFF=1e14+5;
const int MAX=2e5+5;
//const ll MAXN=2e8;
//const int MAX_N=MAX;
//const ll MOD=1e9+7;
const ll MOD=998244353;
//const long double pi=acos(-1.0);
//const double eps=0.00000001;
int gcd(int a,int b){return b?gcd(b,a%b):a;}
template<typename T>inline T abs(T a) {return a>0?a:-a;}
//#define double long double
template<class T> inline
void read(T& num) {
    bool start=false,neg=false;
    char c;
    num=0;
    while((c=getchar())!=EOF) {
        if(c=='-') start=neg=true;
        else if(c>='0' && c<='9') {
            start=true;
            num=num*10+c-'0';
        } else if(start) break;
    }
    if(neg) num=-num;
}
inline ll powMM(ll a,ll b,ll M){
    ll ret=1;
    a%=M;
//    b%=M;
    while (b){
        if (b&1) ret=ret*a%M;
        b>>=1;
        a=a*a%M;
    }
    return ret;
}
//const long double eps=-1.0;
//clock_t t1 = clock();
//fprintf(stderr, "%ld ms
", clock() - t1);
void open()
{
    freopen("in2.txt","r",stdin);
    freopen("out2.txt","w",stdout);
}

const int MAXL=400005;
const int MAXCH=26;
char s[MAXL],q[MAXL];
const char CH='a';
int tot=0;
int cur_len;
int maxlen[MAXL],trans[MAXL][MAXCH],slink[MAXL],pos[MAXL];
int leftmost[MAXL];
int deled[MAXL];
bool D[MAXL];
int DID;
int bck[MAXL];

int pe,ID;
inline int new_state(int _maxlen,int *_trans,int _slink)
{
    maxlen[tot]=_maxlen;
    bck[tot]=0;
    if(_trans==NULL)
        memset(trans[tot],-1,sizeof(trans[tot]));
    else
        for(int j=0;j<26;j++)trans[tot][j]=_trans[j];
    slink[tot]=_slink;
    return tot++;
}
char ch;
inline void add_char()
{
    ++cur_len;
    int c=ch-CH;
    int z=new_state(cur_len,NULL,-1);
    leftmost[z]=cur_len;
    deled[z]=++DID;
    D[DID]=0;
    pos[cur_len]=z;
    int v=pe;
    while(v!=-1&&(trans[v][c]==-1||D[deled[trans[v][c]]]))
    {
        trans[v][c]=z;
        v=slink[v];
    }
    if(v==-1)
        slink[z]=0;
    else
    {
        int x=trans[v][c];
        if(maxlen[v]+1==maxlen[x])
            slink[z]=x;
        else
        {
            int y=new_state(maxlen[v]+1,trans[x],slink[x]);
            leftmost[y]=leftmost[x];
            slink[y]=slink[x];
            deled[y]=deled[x];
            slink[x]=y;
            slink[z]=y;
            int w=v;
            while(w!=-1&&trans[w][c]==x)
            {
                trans[w][c]=y;
                w=slink[w];
            }
        }
    }
    pe=z;
    return;
}
inline void Del(int len)
{
    while(len--)
        D[deled[pos[cur_len--]]]=1;
    pe=pos[cur_len];
}
inline void col()
{
    int now=pe;
    while(now>0)
    {
        bck[now]=ID;
        now=slink[now];
    }
}
int main()
{
    while(~scanf("%s",s))
    {
        DID=tot=pe=ID=0;
        cur_len=0;
        int L=strlen(s);
        pos[0]=0;
        pe=new_state(0,NULL,-1);
        for(int i=0;i<L;i++)
        {
            ch=s[i];
            add_char();
        }
        int m,k,n;
        read(m);
        while(m--)
        {
            read(k);
            if(k==1)
            {
                scanf("%s",s);
                int ls=strlen(s);
                for(int i=0;i<ls;i++)
                {
                    ch=s[i];
                    add_char();
                }
            }
            else if(k==2)
            {
                read(n);
                ++ID;
                col();
                int nowlen,len=n;
                int now=0;
                for(nowlen=0;;++nowlen)
                {
                    if(nowlen==len)
                    {
                        printf("%d
",leftmost[now]-len+1 );
                        break;
                    }
                    if(bck[now]==ID)
                    {
                        printf("%d
",cur_len-nowlen+1 );
                        break; ;
                    }
                    for(int i=0;i<26;i++)
                    {
                        int to=trans[now][i];
                        if(to!=-1&&!D[deled[to]]){now=to;break;}
                    }
                }
            }
            else
            {
                read(n);
                Del(n);
            }
        }
    }
    return 0;
}
/*
bcdacdaa
10
3 1
2 2
*/