bzoj 1212: [HNOI2004]L语言 AC自动机+状压

时间:2023-03-09 01:12:58
bzoj 1212: [HNOI2004]L语言 AC自动机+状压

  为什么这道题网上所有题解写的都是N*Len的trie树的暴力啊,4E的复杂度。。。

  为什么暴力还跑这么快啊TAT。。

  有一个O(Len)的做法就是先把AC自动机建出来,因为每个字典串的长度很小,所以我们可以用fail树状压一发,每个节点记录一个值ss,表示这个点向前(1~10)的长度的串是不是一个字典串,这个东西延fail树递推就行了。

  然后每次把每个串放AC自动机上走,同时记录一个值T表示(s(1,i-11)~s(1-(i-1)))这些串能否成为答案,如果如果T&ss!=0那么i这个点就可以成为答案。

  

 #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
int n,m;
char s[];
int ch[][];int cnt;int v[];
int fail[];int len[];
void insert()
{
int ln=strlen(s);int now=;
for(int i=;i<ln;i++)
{
int c=s[i]-'a';
if(ch[now][c])
{
now=ch[now][c];
}
else
{
ch[now][c]=++cnt;
len[cnt]=len[now]+;
now=ch[now][c];
}
}
v[now]=;
}
queue<int>q;
unsigned int ss[];
void build()
{
for(int i=;i<;i++)if(ch[][i])q.push(ch[][i]);
while(!q.empty())
{
int tmp=q.front();q.pop();
for(int i=;i<;i++)
{
if(ch[tmp][i])
{
q.push(ch[tmp][i]);
fail[ch[tmp][i]]=ch[fail[tmp]][i];
}
else
{
ch[tmp][i]=ch[fail[tmp]][i];
}
}
if(v[tmp])ss[tmp]=(ss[fail[tmp]]^(<<(len[tmp]-)));
else ss[tmp]=ss[fail[tmp]];
}
return ;
}
unsigned int tmp;
void solve()
{
int ans=;
unsigned int ssr=;
int ln=strlen(s+);
int now=;
for(int i=;i<=ln;i++)
{
int c=s[i]-'a';
now=ch[now][c];
tmp=ssr;
ssr<<=;
if(tmp&ss[now])ans=max(ans,i),ssr^=;
}
printf("%d\n",ans);
return ;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++)
{
scanf("%s",s);
insert();
}
build();
for(int i=;i<=m;i++)
{
scanf("%s",s+);
solve();
}
return ;
}