bzoj 4032 [ HEOI 2015 ] 最短不公共子串 —— 后缀自动机+序列自动机

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

序列自动机其实就是每个位置记录一下某字母后面第一个出现位置，为了子序列能尽量长。

对字符串B建一个后缀自动机，一个序列自动机，然后让A在上面找即可；

1.枚举A每个位置开始的子串，在SAM上走，失配就找到了；

2.枚举A子串，用B的序列自动机找失配处；

3.设 f[i][j] 表示A串前 i 个字符形成的子序列对应SAM上节点 j 时的最短子序列长度，可以转移；

4.f[i][j] 表示A串前 i 个字符的子序列对应B串第 j 个位置时的最短子序列长度，用序列自动机转移；

所以注意 f[i][j] 的第二维要开 4000！否则不是 MLE 而是 WA 。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int const xn=,inf=0x3f3f3f3f;
int fa[xn<<],l[xn<<],go[xn<<][],lst=,cnt=,f[xn][xn<<],g[xn][],nxt[];//<<1!!
char s1[xn],s2[xn];
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]+;
      memcpy(go[nq],go[q],sizeof go[q]);
      fa[nq]=fa[q]; fa[q]=fa[np]=nq;
      for(;go[p][w]==q;p=fa[p])go[p][w]=nq;
    }
    }
}
int main()
{
  scanf("%s",s1+); int l1=strlen(s1+);
  scanf("%s",s2+); int l2=strlen(s2+);
  for(int i=;i<=l1;i++)s1[i]-='a';
  for(int i=;i<=l2;i++)s2[i]-='a';
  for(int i=;i<=l2;i++)add(s2[i]);
  for(int i=;i<;i++)nxt[i]=l2+;
  for(int i=l2;i>=;i--)
    {
      for(int j=;j<;j++)g[i][j]=nxt[j];
      nxt[s2[i]]=i;
    }

  int ans=inf;
  for(int i=;i<=l1;i++)
    {
      int p=,j;
      for(j=i;j<=l1;j++)
    {
      if(!go[p][s1[j]])break;
      p=go[p][s1[j]];
    }
      if(j<=l1)ans=min(ans,j-i+);
    }
  if(ans==inf)ans=-; printf("%d\n",ans);

  ans=inf;
  for(int i=;i<=l1;i++)
    {
      if(g[][s1[i]]>l2){ans=min(ans,); continue;}
      int p=g[][s1[i]],j;
      for(j=i+;j<=l1;j++)
    {
      if(g[p][s1[j]]==l2+)break;
      p=g[p][s1[j]];
    }
      if(j<=l1)ans=min(ans,j-i+);
    }
  if(ans==inf)ans=-; printf("%d\n",ans);

  memset(f,0x3f,sizeof f); f[][]=; ans=inf;
  for(int i=;i<l1;i++)
    for(int j=,p;j<=cnt;j++)
      {
    if(f[i][j]==inf)continue;
    f[i+][j]=min(f[i+][j],f[i][j]);
    if((p=go[j][s1[i+]]))f[i+][p]=min(f[i+][p],f[i][j]+);//,printf("%d %d p=%d\n",j,s1[i+1],p);
    else ans=min(ans,f[i][j]+);
      }
  if(ans==inf)ans=-; printf("%d\n",ans);

  memset(f,0x3f,sizeof f); f[][]=; ans=inf;
  for(int i=;i<l1;i++)
    for(int j=,p;j<=l2;j++)
      {
    if(f[i][j]==inf)continue;
    f[i+][j]=min(f[i+][j],f[i][j]);
    if((p=g[j][s1[i+]])<l2+)f[i+][p]=min(f[i+][p],f[i][j]+);
    else ans=min(ans,f[i][j]+);
      }
  if(ans==inf)ans=-; printf("%d\n",ans);
  return ;
}