本来计划下一步去搞AC自动机的……忽发现它以KMP为基础……虽然假期练过,但基本内涵都不太清楚了,都是直接套模板Orz……是本次来再补一下(⊙﹏⊙) ,总算把自己的板给写出来了,真是可喜可贺XD
以下先是KMP的重新理解部分:用parent[]表示母串,son[]表示子串
当比较母串和子串时,以abcdefg和bcdf为例
母串: a b c d e f g
子串: b c d f
序号: 0 1 2 3 4
当比较到上面这一步时,可发现son[2]与parent[2]失配。如果是朴素算法,接下来会把子串右移,重新比较,而当比较到第x位时,子串的前x位与对应母串的的前x位相同,所以可以预处理子串的情况,先把子串与子串本身比较,当第x位失配时,子串要右移到第几位为止,这样就不需要每次都右移一位,一个个地试了。把需要右移到第几位的值保存在一个next数组中。我见过的next数组求法都不知道几组了QAQ,而且还难理解,果然要自己来写才行……
以下是今天写的一个模板:
//该模板用来找子串在母串中是否存在,若存在,返回起始点,否则,返回-1
const int maxsize=100010;
struct KMP
{
char parent[maxsize];int parentlen;
char son[maxsize];int sonlen;
int next[maxsize];
void init()
{
memset(next, 0, sizeof(next));
return ;
}
void insert()
{
scanf("%s", parent);
parentlen=strlen(parent);
scanf("%s", son);
sonlen=strlen(son);
return ;
}
void buildnext()
{
int pos=2, cnt=0;
next[0]=-1;next[1]=0;
while(pos<=sonlen)
{
if(son[pos-1]==son[cnt])
next[pos++]=++cnt;
else if(cnt>0)
cnt=next[cnt];
else
pos++;
}
return ;
}
int compare()
{
int parentcnt=0, soncnt=0;
while(parentcnt<=parentlen)
{
if(parent[parentcnt]==son[soncnt])
{
parentcnt++;soncnt++;
}
else if(soncnt>0)
soncnt=next[soncnt];
else
parentcnt++;
if(soncnt==sonlen)
return parentcnt-soncnt+1;
}
return -1;
}
}kmp;
然后找了道水题试了一下……稍微改下模板就能过了真是太赞!XD
题号是poj
#include<cstdio>
#include<cstring>
const int maxsize=1000010;
struct KMP
{
char parent[maxsize];int parentlen;
char son[maxsize];int sonlen;
int next[maxsize];
void init()
{
memset(next, 0, sizeof(next));
return ;
}
void insert()
{
scanf("%s", son);
sonlen=strlen(son);
scanf("%s", parent);
parentlen=strlen(parent);
return ;
}
void buildnext()
{
int pos=2, cnt=0;
next[0]=-1;next[1]=0;next[sonlen]=0;
while(pos<=sonlen)
{
if(son[pos-1]==son[cnt])
next[pos++]=++cnt;
else if(cnt>0)
cnt=next[cnt];
else
pos++;
}
return ;
}
int compare()
{
int parentcnt=0, soncnt=0;
int ret=0;
while(parentcnt<=parentlen)
{
if(parent[parentcnt]==son[soncnt])
{
parentcnt++;soncnt++;
}
else if(soncnt>0)
soncnt=next[soncnt];
else
parentcnt++;
if(soncnt==sonlen)
ret++;
}
return ret;
}
}kmp;
int main()
{
int ncase;
int ans;
scanf("%d", &ncase);
while(ncase--)
{
kmp.init();
kmp.insert();
kmp.buildnext();
ans=kmp.compare();
printf("%d\n", ans);
}
return 0;
}
自从trie后就喜欢上用struct了……看起来比较方便QvQ,但是总归得用class比较好吧……下次再想,明天可准备自动机了=w=