題目大意：給定字符串A和B，求A最短的子串/子序列S滿足S不是B的子串/子序列
這題真TM有毒*2
搞法類似這道題
然後子串是後綴自動機 子序列自然就是序列自動機了= =
每更新一個x節點時所有沒有x的後繼的節點都連向這個節點
每個節點的parent是這個字母上一次出現的位置
每個字母記錄最後一次出現的位置 更新指針時沿著parent指針撸一遍就行了
#include
#include
#include
#include
#define M 4040
using namespace std;
char A[M],B[M];
struct Suffix_Automaton{
struct Sam{
Sam *son[26],*parent;
int max_dpt;
Sam() {}
Sam(int _):max_dpt(_) {}
}mempool[M],*C;
Sam *root,*last;
Suffix_Automaton()
{
C=mempool;
last=root=new (C++)Sam(0);
}
void Extend(int x)
{
Sam *p=last;
Sam *np=new (C++)Sam(p->max_dpt+1);
for(;p&&!p->son[x];p=p->parent)
p->son[x]=np;
if(!p)
np->parent=root;
else
{
Sam *q=p->son[x];
if(p->max_dpt+1==q->max_dpt)
np->parent=q;
else
{
Sam *nq=new (C++)Sam(p->max_dpt+1);
memcpy(nq->son,q->son,sizeof nq->son);
nq->parent=q->parent;
q->parent=nq;np->parent=nq;
for(;p&&p->son[x]==q;p=p->parent)
p->son[x]=nq;
}
}
last=np;
}
}substr_A,substr_B;
struct Sequence_Automaton{
struct Sqa{
Sqa *son[26],*parent;
}mempool[M],*C;
Sqa *root,*last[27];
Sequence_Automaton()
{
C=mempool;
last[26]=root=new (C++)Sqa;
}
void Extend(int x)
{
Sqa *p,*np=new (C++)Sqa;
int i;
for(i=0;i<=26;i++)
for(p=last[i];p&&!p->son[x];p=p->parent)
p->son[x]=np;
np->parent=last[x];
last[x]=np;
}
}subseq_A,subseq_B;
struct abcd{
int x,y,dpt;
abcd() {}
abcd(int _,int __,int ___):
x(_),y(__),dpt(___) {}
}q[M];
int v[M][M];
int BFS1()//str-str
{
int i,r=0,h=0;
q[++r]=abcd(0,0,0);
v[0][0]=true;
while(r!=h)
{
abcd sta=q[++h];
for(i=0;i<26;i++)
{
if(!substr_A.mempool[sta.x].son[i]) continue;
if(!substr_B.mempool[sta.y].son[i]) return sta.dpt+1;
int xx=substr_A.mempool[sta.x].son[i]-substr_A.root;
int yy=substr_B.mempool[sta.y].son[i]-substr_B.root;
if(v[xx][yy]==1) continue;
v[xx][yy]=1;
q[++r]=abcd(xx,yy,sta.dpt+1);
}
}
return -1;
}
int BFS2()//str-seq
{
int i,r=0,h=0;
q[++r]=abcd(0,0,0);
v[0][0]=true;
while(r!=h)
{
abcd sta=q[++h];
for(i=0;i<26;i++)
{
if(!substr_A.mempool[sta.x].son[i]) continue;
if(!subseq_B.mempool[sta.y].son[i]) return sta.dpt+1;
int xx=substr_A.mempool[sta.x].son[i]-substr_A.root;
int yy=subseq_B.mempool[sta.y].son[i]-subseq_B.root;
if(v[xx][yy]==2) continue;
v[xx][yy]=2;
q[++r]=abcd(xx,yy,sta.dpt+1);
}
}
return -1;
}
int BFS3()//seq-str
{
int i,r=0,h=0;
q[++r]=abcd(0,0,0);
v[0][0]=true;
while(r!=h)
{
abcd sta=q[++h];
for(i=0;i<26;i++)
{
if(!subseq_A.mempool[sta.x].son[i]) continue;
if(!substr_B.mempool[sta.y].son[i]) return sta.dpt+1;
int xx=subseq_A.mempool[sta.x].son[i]-subseq_A.root;
int yy=substr_B.mempool[sta.y].son[i]-substr_B.root;
if(v[xx][yy]==3) continue;
v[xx][yy]=3;
q[++r]=abcd(xx,yy,sta.dpt+1);
}
}
return -1;
}
int BFS4()//seq-seq
{
int i,r=0,h=0;
q[++r]=abcd(0,0,0);
v[0][0]=true;
while(r!=h)
{
abcd sta=q[++h];
for(i=0;i<26;i++)
{
if(!subseq_A.mempool[sta.x].son[i]) continue;
if(!subseq_B.mempool[sta.y].son[i]) return sta.dpt+1;
int xx=subseq_A.mempool[sta.x].son[i]-subseq_A.root;
int yy=subseq_B.mempool[sta.y].son[i]-subseq_B.root;
if(v[xx][yy]==4) continue;
v[xx][yy]=4;
q[++r]=abcd(xx,yy,sta.dpt+1);
}
}
return -1;
}
int main()
{
//freopen("sus.in","r",stdin);
//freopen("sus.out","w",stdout);
int i;
scanf("%s%s",A+1,B+1);
for(i=1;A[i];i++)
{
substr_A.Extend(A[i]-'a');
subseq_A.Extend(A[i]-'a');
}
for(i=1;B[i];i++)
{
substr_B.Extend(B[i]-'a');
subseq_B.Extend(B[i]-'a');
}
cout<
