给定n个字符串(S1,S2,„,Sn),要求找到一个最短的字符串T,使得这n个字符串(S1,S2,„,Sn)都是T的子串。
BZOJ_1195_[HNOI2006]最短母串_AC自动机+BFS+分层图
给定n个字符串(S1,S2,„,Sn),要求找到一个最短的字符串T,使得这n个字符串(S1,S2,„,Sn)都是T的子串。
然后因为n很小,每个节点维护fail到根的路径出现了哪些字符串(出现字符串的状态)
分1<<n个层,f[i][j]表示到第i个节点,当前已经有了的字符串状态为j,最短需要几个字符。
因为边权为1可以直接BFS一遍。
记录一下从哪来的方便输出路径。
代码:
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 650
int ch[N][26],fail[N],Q[N*4250],l,r,n,sta[N],cnt=1,flg[N],path[N][1<<12],zifu[N][1<<12];
char w[N];
void insert(int x) {
int p=1,i;
for(i=1;w[i];i++) {
int &k=ch[p][w[i]-‘A‘];
if(!k) k=++cnt;
p=k;
}
sta[p]|=(1<<(x-1));
}
void build_ac() {
int i,p;
for(i=0;i<26;i++) ch[0][i]=1;
Q[r++]=1;
while(l<r) {
p=Q[l++];
for(i=0;i<26;i++) {
if(ch[p][i]) fail[ch[p][i]]=ch[fail[p]][i],Q[r++]=ch[p][i];
else ch[p][i]=ch[fail[p]][i];
}
sta[p]|=sta[fail[p]];
}
}
void print(int x,int y) {
if(path[x][y]) {
int d=path[x][y];
print(Q[d],Q[d+1]);
}
printf("%c",zifu[x][y]+‘A‘);
}
int main() {
scanf("%d",&n);
int mask=(1<<n)-1;
int i;
for(i=1;i<=n;i++) {
scanf("%s",w+1); insert(i);
}
build_ac();
l=r=0;
memset(path,-1,sizeof(path));
Q[r++]=1; Q[r++]=0; path[1][0]=0;
while(l<r) {
// puts("FUCK");
int x=Q[l++],y=Q[l++];
for(i=0;i<26;i++) {
int c=ch[x][i];
if(path[c][y|sta[c]]==-1) {
path[c][y|sta[c]]=l-2;
zifu[c][y|sta[c]]=i;
if((y|sta[c])==mask) {
print(c,y|sta[c]); return 0;
}
Q[r++]=c; Q[r++]=y|sta[c];
}
}
}
}
BZOJ_1195_[HNOI2006]最短母串_AC自动机+BFS+分层图
原文:https://www.cnblogs.com/suika/p/9186948.html