[bzoj4032] [HEOI2015]最短不公共子串

思路

我们可以发现1、2、3、4询问其实有相同点，都是在找A串中匹配、B串不匹配的最短xxx的长度，于是对A、B串分别建出后缀自动机和序列自动机，然后对于四个询问以相同的方式BFS找到第一个满足条件的串即可。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define maxn 2020
#define maxc 26
using namespace std;
char a[maxn],b[maxn];
inline int idx(char c){return c-'a';}
struct node{int o1,o2,d;};
queue<node>q;bool vis[maxn<<1][maxn<<1];
struct Suffix_Automation{
	int ch[maxn<<1][maxc],tot,la,f[maxn<<1],l[maxn<<1];
	void init(){la=tot=1;}
	inline void extend(int c){
		int p=la,np=++tot;la=np;l[np]=l[p]+1;
		for(;p&&(!ch[p][c]);p=f[p])ch[p][c]=np;
		if(!p){f[np]=1;return;}
		int q=ch[p][c];
		if(l[q]==l[p]+1){f[np]=q;return;}
		int nq=++tot;memcpy(ch[nq],ch[q],sizeof(ch[q]));f[nq]=f[q];l[nq]=l[p]+1;
		f[np]=f[q]=nq;for(;p&&(ch[p][c]==q);p=f[p])ch[p][c]=nq;
	}
	void insert(char *s){
		int len=strlen(s+1);
		for(int i=1;i<=len;++i)
		extend(idx(s[i]));
	}
}sam[2];
struct Sequence_AutoMaton{
	int ch[maxn][maxc],la[maxc],f[maxn],tot;
	void init(){tot=1;for(int i=0;i<maxc;++i)la[i]=1;}
	inline void extend(int c){
		int np=++tot;f[np]=la[c];
		for(int i=0;i<maxc;++i)
		for(int j=la[i];j&&(!ch[j][c]);j=f[j])ch[j][c]=np;
		la[c]=np;
	}
	void insert(char *s){
		int len=strlen(s+1);
		for(int i=1;i<=len;++i)
		extend(idx(s[i]));
	}
}sqa[2];
void init(){
	sam[0].init(),sam[1].init();
	sqa[0].init(),sqa[1].init();
	scanf("%s%s",a+1,b+1);
	sam[0].insert(a),sqa[0].insert(a);
	sam[1].insert(b),sqa[1].insert(b);
}
void solve1(){
	while(!q.empty())q.pop();
    memset(vis,0,sizeof(vis));
	q.push((node){1,1,0});vis[1][1]=1;
	while(!q.empty()){
		node u=q.front();q.pop();
		for(int c=0;c<maxc;++c){
			node v=(node){sam[0].ch[u.o1][c],sam[1].ch[u.o2][c],u.d+1};
			if(vis[v.o1][v.o2]||!v.o1)continue;
			if(!v.o2){printf("%d\n",v.d);return;}
			q.push(v);vis[v.o1][v.o2]=1;
		}
	}
	printf("-1\n");
}
void solve2(){
	while(!q.empty())q.pop();
    memset(vis,0,sizeof(vis));
	q.push((node){1,1,0});vis[1][1]=1;
	while(!q.empty()){
		node u=q.front();q.pop();
		for(int c=0;c<maxc;++c){
			node v=(node){sam[0].ch[u.o1][c],sqa[1].ch[u.o2][c],u.d+1};
			if(vis[v.o1][v.o2]||!v.o1)continue;
			if(!v.o2){printf("%d\n",v.d);return;}
			q.push(v);vis[v.o1][v.o2]=1;
		}
	}
	printf("-1\n");
}
void solve3(){
	while(!q.empty())q.pop();
    memset(vis,0,sizeof(vis));
	q.push((node){1,1,0});vis[1][1]=1;
	while(!q.empty()){
		node u=q.front();q.pop();
		for(int c=0;c<maxc;++c){
			node v=(node){sqa[0].ch[u.o1][c],sam[1].ch[u.o2][c],u.d+1};
			if(vis[v.o1][v.o2]||!v.o1)continue;
			if(!v.o2){printf("%d\n",v.d);return;}
			q.push(v);vis[v.o1][v.o2]=1;
		}
	}
	printf("-1\n");
}
void solve4(){
	while(!q.empty())q.pop();
	memset(vis,0,sizeof(vis));
	q.push((node){1,1,0});vis[1][1]=1;
	while(!q.empty()){
		node u=q.front();q.pop();
		for(int c=0;c<maxc;++c){
			node v=(node){sqa[0].ch[u.o1][c],sqa[1].ch[u.o2][c],u.d+1};
			if(vis[v.o1][v.o2]||!v.o1)continue;
			if(!v.o2){printf("%d\n",v.d);return;}
			q.push(v);vis[v.o1][v.o2]=1;
		}
	}
	printf("-1\n");
}
void solve(){
	solve1(),solve2();
	solve3(),solve4();
}
int main(){
	init();
	solve();
	return 0;
}