【模板】KMP字符串匹配

题目描述

如题，给出两个字符串s1和s2，其中s2为s1的子串，求出s2在s1中所有出现的位置。

为了减少骗分的情况，接下来还要输出子串的前缀数组next。如果你不知道这是什么意思也不要问，去百度搜[kmp算法]学习一下就知道了。

输入输出格式

输入格式：

 

第一行为一个字符串，即为s1（仅包含大写字母）

第二行为一个字符串，即为s2（仅包含大写字母）

 

输出格式：

 

若干行，每行包含一个整数，表示s2在s1中出现的位置

接下来1行，包括length(s2)个整数，表示前缀数组next[i]的值。

 

输入输出样例

输入样例#1：

ABABABC
ABA

输出样例#1：

1
3
0 0 1 

说明

时空限制：1000ms,128M

数据规模：

设s1长度为N，s2长度为M

对于30%的数据：N<=15，M<=5

对于70%的数据：N<=10000，M<=100

对于100%的数据：N<=1000000，M<=1000

样例说明：

所以两个匹配位置为1和3，输出1、3

题解：
虽说是裸的KMP，但是这个算法还是蛮难理解的。

 

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
char s1[1000001],s2[100001];
int Next[100001];
void getnext(char *p,int *Next)
{
    int i,len=strlen(p)-1;
    Next[0]=0;
    for(i=1;i<=len;i++)
    {
        int j=Next[i-1];
        while(p[i]!=p[j]&&j>0)j=Next[j-1];
        if(p[i]==p[j])Next[i]=j+1;
        else Next[i]=0;
    }
}
void KMP(char *s,char *p)
{
    int i,j,lens=strlen(s)-1,lenp=strlen(p)-1;
    getnext(p,Next);
    i=0;j=0;
    while(i<lens&&j<lenp)
    {
        while(i<lens&&j<lenp&&s[i]==p[j]){i++;j++;}
        if(j==0)i++;
        if(j==-1){i++;j++;}
        while(j>0&&s[i]!=p[j])j=Next[j-1];
        if(j==lenp)
        {
            printf("%d\n",i-j+1);
            j=Next[j]-1;
        }
    }
}
int main()
{
    int i;
    scanf("%s%s",s1,s2);
    KMP(s1,s2);
    int len=strlen(s2);
    for(i=0;i<len;i++)
    printf("%d ",Next[i]);
    return 0;
}