这题调了好久:
问题1:Trie图的节点数给小了,应该是150*70,我给的是150*26,少脑子了;
问题2:就是数组开小了,老是提交超时,导致开始怀疑算法的复杂度有问题;
另外,LA上这题的数据非常弱,非常非常弱。
AC自动机的本质:先把模板串建成Trie树,然后根据KMP的思想用BFS建Trie图,说穿了就是维护两个数组,f[u]往回找,last[u]指向其他较短模板串的节点编号,其实并不复杂,关键是理解清楚递推式。
#include <iostream>
#include <fstream>
#include <algorithm>
#include <cstring>
#include <queue>
using namespace std;
#define SUM 155
#define ML 75
#define N 11000 //155*26
#define SIG 26
#define MAXL 1000005
char p[SUM][ML], text[MAXL];
int num[SUM], cnt[SUM];
int n, T;
struct AC_auto {
int ch[N][SIG];
int val[N], f[N], last[N];
int sz;
void init() {
sz=1;
memset(ch[0],0,sizeof(ch[0]));
memset(val,0,sizeof(val));
//memset(f,0,sizeof(f));
// memset(last,0,sizeof(last));
}
void Insert(char *s,int v) {
int u=0, n=strlen(s);
for (int i=0;i<n;i++) {
int c= s[i]-'a';
if (!ch[u][c]) {
memset(ch[sz],0,sizeof(ch[sz]));
//val[sz]=0;
ch[u][c]=sz++;
}
u= ch[u][c];
}
val[u]=v;
// correct
num[v]=u;
}
void getFail(){
queue<int> q;
f[0]=0;
//init
for (int c=0;c<SIG;c++) {
int u=ch[0][c];
if (u) {
f[u]=0;
q.push(u);
last[u]=0;
}
}
//bfs
while (!q.empty()) {
int r=q.front();
for (int c=0;c<SIG;c++) {
int u=ch[r][c];
if (!u) {
//key
ch[r][c]= ch[f[r]][c];
continue;
}
q.push(u);
int v=f[r];
while (v && !ch[v][c]) v=f[v];
f[u]= ch[v][c];
last[u] = val[f[u]] ? f[u] : last[f[u]];
}
q.pop();
}
}
void Find(char* T){
int n = strlen(T);
int j = 0;
for (int i=0;i<n;i++) {
int c=T[i]-'a';
while (j && !ch[j][c]) j=f[j];
j= ch[j][c];
//correct
int temp = j;
if (val[temp]) {
while (temp) {
cnt[val[temp]]++;
temp=last[temp];
}
}
else if (last[temp]) {
temp=last[temp];
while (temp) {
cnt[val[temp]]++;
temp=last[temp];
}
}
}
}
};
AC_auto ac;
int main()
{
freopen("1.in","r",stdin);
while (scanf("%d",&n)!=EOF && n) {
ac.init();
memset(cnt,0,sizeof(cnt));
memset(num,0,sizeof(num));
for (int i=1;i<=n;i++){
scanf("%s",p[i]);
ac.Insert(p[i],i);
}
ac.getFail();
scanf("%s",text);
ac.Find(text);
int best= -1;
for (int i=1;i<=n;i++){
if (cnt[i] >best) best= cnt[i];
}
printf("%d\n",best);
for (int i=1;i<=n;i++)
if (cnt[ac.val[num[i]]]==best)
printf("%s\n",p[i]);
}
return 0;
}