版权声明:本文为博主原创文章,未经博主允许不得转载,除非先点了赞。 https://blog.csdn.net/A_Bright_CH/article/details/83006946
题目
求模式串对主串的最大匹配长度。
题解
DP+AC自动机
设f[i]=true表示前i位可以匹配出来,那么转移方案就是背包的,(要求:主串的后缀与模式串j完全一致),其中len[j]是一个模式串的长度。
如果要是大(bao)力DP的话,显然会很慢,因为我们要配对每一个模式串。
显然这种后缀配前缀的问题应当交由AC自动机来处理。
把匹配串逐一insert到字典树中,同时标记一下结尾,求一个fail指针,然后我们就可以来一波DP了。DP时,对于长度为i的前缀,我们通过跳fail枚举所有 模式串的前缀 与当前主串后缀相同的树节点,一不小心碰到标记过的,说明满足这个主串以i结尾的后缀与这个模式串匹配,此时符合DP方程的转移条件,考虑更新f[i]。
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=30;
const int maxp=10*maxn;
char s[(const int)1e6+10];int len;
struct Tr{int x,fail,dep,son[25];}tr[maxp];int cnt=1,root=1;
bool flag[maxp];
void insert()
{
int x=root;
for(int i=1;i<=len;i++)
{
int k=s[i]-'a';
if(tr[x].son[k]==0)
{
tr[x].son[k]=++cnt;
tr[cnt].dep=tr[x].dep+1;
}
x=tr[x].son[k];
}
flag[x]=true;
}
int head,tail,q[maxp];
void getfail()
{
head=0,tail=1;q[0]=root;
while(head<tail)
{
int x=q[head++];
for(int k=0;k<26;k++)//debug 一共26个字母!
{
int y=tr[x].son[k];
if(y==0) continue;
else if(x==root) tr[y].fail=root;
else
{
int p=tr[x].fail;//debug int p=x;
while(p!=root && !tr[p].son[k]) p=tr[p].fail;
//debug tr[y].fail=p;
if(tr[p].son[k]) tr[y].fail=tr[p].son[k];
else tr[y].fail=root;
}
q[tail++]=y;
}
}
}
int f[(const int)1e6+10];//f[i]表示长度为i的前缀能否被翻译
void solve(int id)
{
int x=root,ans=0;
f[0]=id;
for(int i=1;i<=len;i++)
{
int k=s[i]-'a';
while(x!=root && !tr[x].son[k]) x=tr[x].fail;
x= !tr[x].son[k]?root:tr[x].son[k] ;//debug x=tr[x].son[k];
for(int p=x;p!=root;p=tr[p].fail)
{
if(flag[p] && i-tr[p].dep>=0 && f[i-tr[p].dep]==id){f[i]=id;break;}//DP转移
}
if(f[i]==id) ans=i;
}
printf("%d\n",ans);
}
int main()
{
int n,m;
scanf("%d%d",&n,&m);
tr[root].dep=0;
for(int i=1;i<=n;i++)
{
scanf("%s",s+1);
len=strlen(s+1);
insert();
}
getfail();
// for(int i=1;i<=cnt;i++)
// {
// printf("fail %d : %d\n",i,tr[i].fail);
// }
for(int i=1;i<=m;i++)
{
scanf("%s",s+1);
len=strlen(s+1);
solve(i);
}
return 0;
}