HDU 4270 SAM 后缀自动机

题目：给出一个串，有3种操作，
操作1：在原串后面，添加一个串
操作2：查询长度为len的子串或者长度小于等于len的后缀中字典序最小的
操作3：删除最后的len个字符



关于删除

首先明确一点： 新建的点只与模板点有关系

由新加入的点7产生的点8是根据点3为模板进行复制的（复制所有信息，包括位置）。那么点8是有点3控制的。

可以从两个图看出，删除点7，不删除点8也是可以达到第一个图的效果。

只是图不是最简的SAM。即a-a-b的父亲变为了a-b。但是无关大雅。

从这个可以看出新建立的辅助点是根据某点为模板，并且为这个模板分担了一些线路的点。

模板点控制这这些辅助点。当且仅当这个模板点被删除后，以这个模板点为模板扩展出去的辅助点才会也才能删除。

声明一个bool is_del[]

结构体里加上一个 bool *d 指针

字符串里的每一个结点对应一个域，即d=is_del[x]，那么其复制出来的点也共享这个is_del[x]

当这个字符被删除后，也就是说这个域被删除后，所有的点的*d==0 ，全部删除了

删除从最后一个字符（域）开始删除，删除len个字符（域）即可

删除完之后last需要更新到删除后的第一个点

关于询问

node记录下该点的位置，如果这个状态是后缀，则ans=seq_cnt-i+1

否则为ans=p->pos-i+1

#include <iostream>
#include <cstring>
#include <cstdio>
#include <string>
using namespace std;

const int maxn=200005;
struct suffixautomaton
{
struct node
{
int len,suf;//到这个状态允许的最大长长度，即max(s)
node *f,*ch[26];  //如果suf==v_idx,则意味这个状态是后缀子串
bool *d;  //
int pos;  //字符原来的位置
node(){}
node(int l)
{
len=l;
f=NULL;
memset(ch,0,sizeof(ch));
}
};
node *root,*last;
node pool[maxn*2];  //储蓄结点用的
node *seq[maxn];  //用于储存输入进来的字符串,,其实就是储存pool数组中主心轴的点。
bool is_del[maxn*2];  //实际上is_del[i] 是代表第i个字符被删除没。node里的 *d是某一个is_del的地址，同一个域里的点共享一个is_del，同在一个域里面的点共享一个is_del[x].
int del_cnt,seq_cnt;  //is_del数组的大小，seq数组的大小
int cnt;             //结点的数量
int v_idx;        //此处的v_idx方便设置suf,避免出错
void init()
{
v_idx=0;
memset(is_del,0,sizeof(is_del));
del_cnt=seq_cnt=0;
root=last=pool;
seq[0]=root;    //此处必须先设置第0位root,避免数组越界
memset(root,0,sizeof(node));
cnt=1;
}
node * new_node(int l=0)
{
node *x=pool+cnt++;
memset(x,0,sizeof(node));
if(l!=0) x->len=l;
return x;
}
void add(char ch)
{
int c=ch-'a';
node *p=last,*np=new_node(last->len+1);
np->pos=np->len;
np->d=is_del+del_cnt++;
*np->d=0;
seq[seq_cnt=np->len]=np;      //以上四行初始化

last=np;
for(;NULL!=p&&(NULL==p->ch[c]||*p->ch[c]->d);p=p->f)  //如果这个点没有被删除且无下一个点ch
p->ch[c]=np;
if(NULL==p)
np->f=root;
else
{
if(p->ch[c]->len==p->len+1)
np->f=p->ch[c];
else
{ //新建立的点是以p->ch[c]为模板复制的点，共享一个is_del[],当且仅当模板点被删除后新建立的点才会被删除。
node *q=p->ch[c],*nq=new_node();
*nq=*q;
nq->len=p->len+1;
q->f=np->f=nq;
for(;NULL!=p&&p->ch[c]==q;p=p->f)
p->ch[c]=nq;
}
}
}
void get_suf()//如果suf==v_idx,则意味这个状态是后缀子串
{
v_idx++;
node *p=last;
while(p!=NULL)
p->suf=v_idx,p=p->f;
root->suf=0;
}
void del(int len)  //从最后一个域开始删除
{
while(len--)
*seq[seq_cnt--]->d=1;
last=seq[seq_cnt];  //last需要改变
}
void query(int l)
{
get_suf();
node *p=root;
for(int i=1;i<=l;i++)
{
for(int j=0;j<26;j++)
if(!(NULL==p->ch[j]||*p->ch[j]->d))
{
p=p->ch[j];
if(i==l)
cout<<p->pos-i+1<<endl;
else if(p->suf==v_idx)
{
cout<<seq_cnt-i+1<<endl;
return ;
}
break;
}
}
}
};
suffixautomaton sam;
char str[maxn];
int main()
{
while(scanf("%s",str)!=EOF)
{
sam.init();
for(int i=0;str[i];i++)
sam.add(str[i]);

int q,cas,len;
scanf("%d",&q);
while(q--)
{
scanf("%d",&cas);
if(cas==1)
{
scanf("%s",str);
for(int i=0;str[i];i++)
sam.add(str[i]);
}
else if(cas==2)
{
scanf("%d",&len);
sam.query(len);
}
else
{
scanf("%d",&len);
sam.del(len);
}
}
}
return 0;
}