[后缀自动机 线段树] 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;
}