kmp算法巩固 以poj3461为试验

本来计划下一步去搞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=