luogu P3375 【模板】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

#include<cstdio>
#include<cstring>
using namespace std;
int i,j,len1,len2,ne[1001];
char s1[1000001],s2[1001];
void findnext()
{
    for (i=2,j=0; i<=len2; i++)//i不从1（开头）开始是为了防止后缀最大长度等于前缀最大长度等于整串长，至于j那是已经匹配的前缀最大长和后缀最大长 
    {    
        while (s2[i]!=s2[j+1]&&j>0) j=ne[j];//如果当前位与前缀已匹配后一位不相等，就递归找最大公共前后缀，否则最优跳出******************************label 1;
        if (s2[i]==s2[j+1]) ne[i]=++j;//如果尽可能地匹配后，还是不能使当前位与之前匹配的后一位相等，就把ne赋为0（不赋值）否则在已匹配位+1 
    }
}
void kmp()
{
    for (i=1,j=0; i<=len1; i++)//这里i指要努力匹配的，j指已经匹配的； 
    {
        while (s1[i]!=s2[j+1]&&j>0) j=ne[j];//label1（手比较懒） 
        if (s1[i]==s2[j+1]) ++j;//如果尽可能匹配后还是不匹配，理论上应该j清零，{可我为什么没写呢？因为j在上一句话就清0了}，否则增加已匹配位数 （大括号表示层次） 
        if (j==len2) {printf("%d\n",i-len2+1); j=ne[len2];}//如果全部匹配输出，然后由于我们已经知道j的最大公共前后缀了，可以直接将j设为他的最大公共前后缀用以加速 
    }
}
int main()
{
    scanf("%s%s",s1+1,s2+1); 
    len1=strlen(s1+1); len2=strlen(s2+1);
    findnext();
    kmp();
    for (i=1; i<=len2; i++) printf("%d ",ne[i]);
    return 0; 
}

 

 

#include<iostream>//***28&&16行相当于递推算法****** 
#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN=1000001;
char a[MAXN],b[MAXN];
int la,lb;
int p[MAXN];
void makep()
{
    int j=0;// 记录最长公共前后缀的长度；
    p[0]=0;// 第一个数的前缀和后缀都是空集，故为0 
    for(int i=1;i<lb;i++)// 第一个是0，所以从1开始，半闭半开区间 
    {
        while(j>0&&b[j]!=b[i])//当出现不匹配的情况的的话，就让j从p数组中取一个能继续匹配的元素 
        j=p[j-1];// 因为对j的最后一次处理是j++，而p是从0开始，故j=p[j-1]
        if(b[j]==b[i])//如果相同的话就让j++ 自己可以举例子试试 
        j++;
        p[i]=j;// 将每一个元素的最长公共前后缀的长度记录下来，方便以后处理 
    }
}
void KMP()
{
    int j=0;// 已经匹配的最长的长度 
    for(int i=0;i<la;i++)
    {
        if(j>0&&b[j]!=a[i])// 当将要匹配的元素不想同且找到相同的匹配元素的时候 
        j=p[j-1];// j==0的时候p[-1]不存在，所以j>0 
        if(b[j]==a[i])// 同理，当将要匹配的元素相同的时候，j++ 
        j++;
        if(j==lb)// 找到一个匹配 
            printf("%d\n",i+1-lb+1),j=p[j-1];
        // 因为i，lb都是从0开始计数的，so+1 
    }
}
int main()
{
    scanf("%s %s",a,b);
    la=strlen(a);lb=strlen(b);//统计两个字符串的长度 
    makep();// 构造p数组 
    KMP();//
    for(int i=0;i<lb;i++)printf("%d ",p[i]);
    return 0;
}