bzoj 4032 [HEOI2015]最短不公共子串——后缀自动机

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4032

不是 b 的子串的话就对 b 建后缀自动机，在 a 上枚举从每个位置开始的子串或者找子序列（子序列就是记录 a 的前 i 个，走到 b 的 j 状态用的最短长度），对应到自动机上看看能不能走下去就行了。

不是 b 的子序列的话就对 b 建子序列自动机？就是那个知道每个位置再填一个字符会走到哪个位置的数组。然后在 a 上枚举，看看自动机上能不能走下去就行了。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=,M=,K=;
char a[N],b[N];
int n,m,cnt=,lst=,go[M][K],fa[M],l[M],nxt[N][K],lt[K];
int dp[N][M];//M//int
int Mn(int a,int b){return a<b?a:b;}
void add(int w)
{
  int p=lst,np=++cnt;lst=np;l[np]=l[p]+;
  for(;p&&!go[p][w];p=fa[p])go[p][w]=np;
  if(!p)fa[np]=;
  else
    {
      int q=go[p][w];
      if(l[q]==l[p]+)fa[np]=q;
      else
    {
      int nq=++cnt;l[nq]=l[p]+;
      fa[nq]=fa[q];fa[q]=nq;fa[np]=nq;
      memcpy(go[nq],go[q],sizeof go[q]);
      for(;go[p][w]==q;p=fa[p])go[p][w]=nq;
    }
    }
}
void solve1()
{
  int ans=n+;
  for(int i=;i<=n;i++)
    {
      int cr=;
      for(int j=i;j<=n;j++)
    {
      if(!go[cr][a[j]-'a'+])
        {ans=Mn(ans,j-i+);break;}
      cr=go[cr][a[j]-'a'+];
    }
    }
  printf("%d\n",ans>n?-:ans);
}
void solve2()
{
  int ans=n+;
  for(int t=;t<=n;t++)
    {
      int cr=;
      for(int i=t;i<=n;i++)
    {
      int d=a[i]-'a'+;
      if(nxt[cr][d])cr=nxt[cr][d];
      else {ans=Mn(ans,i-t+);break;}
    }
    }
  printf("%d\n",ans>n?-:ans);
}
void solve3()
{
  memset(dp,0x3f,sizeof dp);
  dp[][]=; int ans=n+;
  for(int i=;i<=n;i++)
    {
      int d=a[i]-'a'+;
      for(int j=;j<=cnt;j++)
    if(dp[i-][j]<=n)
    {
      dp[i][j]=Mn(dp[i][j],dp[i-][j]);
      if(!go[j][d])ans=Mn(ans,dp[i-][j]+);
      else dp[i][go[j][d]]=Mn(dp[i][go[j][d]],dp[i-][j]+);
    }
    }
  printf("%d\n",ans>n?-:ans);
}
void solve4()
{
  memset(dp,0x3f,sizeof dp);
  dp[][]=; int ans=n+;
  for(int i=;i<=n;i++)
    {
      int d=a[i]-'a'+;
      for(int j=;j<=m;j++)
    if(dp[i-][j]<=n)
      {
        dp[i][j]=Mn(dp[i][j],dp[i-][j]);
        if(!nxt[j][d])ans=Mn(ans,dp[i-][j]+);
        else dp[i][nxt[j][d]]=Mn(dp[i][nxt[j][d]],dp[i-][j]+);
      }
    }
  printf("%d\n",ans>n?-:ans);
}
int main()
{
  scanf("%s",a+);scanf("%s",b+);
  n=strlen(a+);m=strlen(b+);
  for(int i=;i<=m;i++)add(b[i]-'a'+);
  for(int i=m;i>=;i--)
    {
      for(int j=;j<=;j++)
    nxt[i][j]=lt[j];
      lt[b[i]-'a'+]=i;
    }
  solve1();solve2();
  solve3();solve4();
  return ;
}