[后缀自动机 线段树] BZOJ 1396 识别子串 & BZOJ 2865 字符串识别

后缀自动机建起来发现出现一次的子串必然由|right(x)|==1的状态也就是叶子节点贡献
对于每个叶子x 在串上的位置为p 令l为maxs(x) 令r为maxs(fat(x))
那么区间[p-r+1,p]用r+1更新答案
区间[p-l+1,p-r]用一条斜率为-1的线段更新答案
这个可以暴力跑李超线段树
但是这道题很特殊啊 直接上两颗线段树分别维护两种情况就好了

上2865的代码吧 这题tm卡空间

#include<cstdio>#include<cstdlib>#include<algorithm>#include<cstring>#include<map>using namespace std;const int N=500005;char S[N]; int n;struct state{  int link,len; map<int,int> next;}st[N<<1];int ncnt,last;int cnt[N<<1];inline void Extend(int c){  int cur=++ncnt,p;  st[cur].len=st[last].len+1; cnt[cur]=1;  for (p=last;p!=-1 && !st[p].next.count(c);p=st[p].link)    st[p].next[c]=cur;  if (p==-1)    st[cur].link=0;  else{    int q=st[p].next[c];    if (st[q].len==st[p].len+1)      st[cur].link=q;    else{      int nq=++ncnt;      st[nq].len=st[p].len+1;      st[nq].next=st[q].next;      st[nq].link=st[q].link;      for (;p!=-1 && st[p].next[c]==q;p=st[p].link)    st[p].next[c]=nq;      st[cur].link=st[q].link=nq;    }  }  last=cur;}int tmp[N<<1],sum[N<<1];inline void Sort(){  for (int i=0;i<=ncnt;i++) sum[st[i].len]++;  for (int i=1;i<=st[last].len;i++) sum[i]+=sum[i-1];  for (int i=0;i<=ncnt;i++) tmp[sum[st[i].len]--]=i;}struct SEG{  int T[N<<2],F[N<<2],TH,M;  void Build(int n){    for (M=1,TH=0;M<n+2;M<<=1,TH++);    for (int i=2*M-1;i;i--) F[i]=T[i]=1<<30;  }  void pushdown(int rt){    int p;    for (int i=TH;i;i--)      if (F[p=rt>>i]!=1<<30){    F[p<<1]=min(F[p<<1],F[p]); T[p<<1]=min(T[p<<1],F[p]);    F[p<<1|1]=min(F[p<<1|1],F[p]); T[p<<1|1]=min(T[p<<1|1],F[p]);    F[p]=1<<30;      }  }  void Modify(int s,int t,int r){    for (pushdown(s+=M-1),pushdown(t+=M+1);s^t^1;){      if (~s&1) F[s^1]=min(F[s^1],r),T[s^1]=min(T[s^1],r);      if ( t&1) F[t^1]=min(F[t^1],r),T[t^1]=min(T[t^1],r);      T[s>>=1]=min(T[s<<1],T[s<<1|1]);      T[t>>=1]=min(T[t<<1],T[t<<1|1]);    }    while (s>>=1)      T[s]=min(T[s<<1],T[s<<1|1]);  }  int Query(int s,int t){    int ret=1<<30;    for (pushdown(s+=M-1),pushdown(t+=M+1);s^t^1;s>>=1,t>>=1){      if (~s&1) ret=min(ret,T[s^1]);      if ( t&1) ret=min(ret,T[t^1]);    }    return ret;  }}Seg1,Seg2;int main(){  freopen("t.in","r",stdin);  freopen("t.out","w",stdout);  scanf("%s",S+1); n=strlen(S+1);  st[0].link=-1;  for (int i=1;i<=n;i++) Extend((int)S[i]-'a');  Sort();  for (int i=ncnt+1;i>1;i--) cnt[st[tmp[i]].link]+=cnt[tmp[i]];  Seg1.Build(n); Seg2.Build(n);  for (int i=1;i<=ncnt;i++)    if (cnt[i]==1){      int p=st[i].len,l=st[i].len,r=st[st[i].link].len;      if (r) Seg1.Modify(p-r+1,p,r+1);      Seg2.Modify(p-l+1,p-r,p);    }  int ans;  for (int i=1;i<=n;i++){    ans=Seg1.Query(i,i);    ans=min(ans,Seg2.Query(i,i)-i+1);    printf("%d\n",ans);  }  return 0;}