BZOJ4032: [HEOI2015]最短不公共子串（后缀自动机+序列自动机）

题目描述

在虐各种最长公共子串、子序列的题虐的不耐烦了之后，你决定反其道而行之。

一个串的“子串”指的是它的连续的一段，例如bcd是abcdef的子串，但bde不是。

一个串的“子序列”指的是它的可以不连续的一段，例如bde是abcdef的子串，但bdd不是。

下面，给两个小写字母串A，B，请你计算：

(1) A的一个最短的子串，它不是B的子串

(2) A的一个最短的子串，它不是B的子序列

(3) A的一个最短的子序列，它不是B的子串

(4) A的一个最短的子序列，它不是B的子序列

输入输出格式

输入格式：

有两行，每行一个小写字母组成的字符串，分别代表A和B。

输出格式：

输出4行，每行一个整数，表示以上4个问题的答案的长度。如果没有符合要求的答案，输出-1.

输入输出样例

输入样例#1： 复制

aabbcc
abcabc

输出样例#1： 复制

2
4
2
4

输入样例#2： 复制

aabbcc
aabbcc

输出样例#2： 复制

-1
-1
2
-1

说明

对于20%的数据，A和B的长度都不超过20

对于50%的数据，A和B的长度都不超过500

对于100%的数据，A和B的长度都不超过2000

题解

　　这题好烦……

　　据说还有大佬是用dp的方法做的……相当于做了四道题……先膜为敬……

　　然后看到有大佬是直接建上后缀自动机和回文自动机，然后在上面跑BFS，如果跑不动了就返回答案，可以保证最优

　　然后我因为一堆sb错误调了一个多小时……

　　为啥取队头元素不能直接整个结构体取啊！为啥我SAM都能打错啊！

 //minamoto
 #include<cstring>
 #include<algorithm>
 #include<iostream>
 #include<cstdio>
 #include<queue>
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 char sr[<<],z[];int C=-,Z;
 inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
 inline void print(int x){
     if(C><<)Ot();if(x<)sr[++C]=,x=-x;
     while(z[++Z]=x%+,x/=);
     while(sr[++C]=z[Z],--Z);sr[++C]='\n';
 }
 const int N=;
 struct SAM{
     int fa[N],ch[N][],l[N];int last,cnt;
     SAM(){last=cnt=;}
     void ins(int c){
         int p=last,np=++cnt;last=np,l[np]=l[p]+;
         for(;p&&!ch[p][c];p=fa[p]) ch[p][c]=np;
         if(!p) fa[np]=;
         else{
             int q=ch[p][c];
             if(l[q]==l[p]+) fa[np]=q;
             else{
                 int nq=++cnt;l[nq]=l[p]+;
                 memcpy(ch[nq],ch[q],sizeof(ch[q]));
                 fa[nq]=fa[q],fa[q]=fa[np]=nq;
                 for(;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;
             }
         }
     }
 }SA,SB;
 struct SqAM{
     int ch[N][],lst[],fa[N];
     int root,cnt;
     SqAM(){root=cnt=;for(int i=;i<;++i) lst[i]=;}
     void ins(int c){
         int p=lst[c],np=++cnt;
         fa[np]=p;
         for(int i=;i<;++i)
         for(int j=lst[i];j&&!ch[j][c];j=fa[j])
         ch[j][c]=np;
         lst[c]=np;
     }
 }SQA,SQB;
 struct node{
     int a,b,step;
     node(int a=,int b=,int step=):a(a),b(b),step(step){}
 };
 int vis[N][N];
 int bfs1(){
     queue<node> q;
     q.push(node(,,));
     vis[][]=;
     while(!q.empty()){
         int ra=q.front().a,rb=q.front().b,st=q.front().step;
         for(int i=;i<;i++){
             int va=SA.ch[ra][i],vb=SB.ch[rb][i];
             if(va&&vb){
                 if(vis[va][vb]==) continue;
                 q.push(node(va,vb,st+));
                 vis[va][vb]=;
             }
             if(va&&!vb) return st+;
         }
         q.pop();
     }
     return -;
 }
 int bfs2(){
     queue<node> q;
     q.push(node(,,));
     vis[][]=;
     while(!q.empty()){
         int ra=q.front().a,rb=q.front().b,st=q.front().step;
         for(int i=;i<;i++){
             int va=SA.ch[ra][i],vb=SQB.ch[rb][i];
             if(va&&vb){
                 if(vis[va][vb]==) continue;
                 q.push(node(va,vb,st+));
                 vis[va][vb]=;
             }
             if(va&&!vb) return st+;
         }
         q.pop();
     }
     return -;
 }
 int bfs3(){
     queue<node> q;
     q.push(node(,,));
     vis[][]=;
     while(!q.empty()){
         int ra=q.front().a,rb=q.front().b,st=q.front().step;
         for(int i=;i<;i++){
             int va=SQA.ch[ra][i],vb=SB.ch[rb][i];
             if(va&&vb){
                 if(vis[va][vb]==) continue;
                 q.push(node(va,vb,st+));
                 vis[va][vb]=;
             }
             if(va&&!vb) return st+;
         }
         q.pop();
     }
     return -;
 }
 int bfs4(){
     queue<node> q;
     q.push(node(,,));
     vis[][]=;
     while(!q.empty()){
         int ra=q.front().a,rb=q.front().b,st=q.front().step;
         for(int i=;i<;i++){
             int va=SQA.ch[ra][i],vb=SQB.ch[rb][i];
             if(va&&vb){
                 if(vis[va][vb]==) continue;
                 q.push(node(va,vb,st+));
                 vis[va][vb]=;
             }
             if(va&&!vb) return st+;
         }
         q.pop();
     }
     return -;
 }
 char a[N],b[N];int tot1,tot2;char ch;
 int main(){
     scanf("%s%s",a+,b+);
     tot1=strlen(a+),tot2=strlen(b+);
     for(int i=;i<=tot1;++i) SA.ins(a[i]-'a'),SQA.ins(a[i]-'a');
     for(int i=;i<=tot2;++i) SB.ins(b[i]-'a'),SQB.ins(b[i]-'a');
     printf("%d\n%d\n%d\n%d\n",bfs1(),bfs2(),bfs3(),bfs4());
     return ;
 }