kmp算法

我首先看了b站上的一个博主，虽然他现在好像停更了，蛮可惜的。

https://www.bilibili.com/video/BV1hW411a7ys?spm_id_from=333.999.0.0

里面比较详细。

模板是洛谷上的题，可能之前用find是不是比较多，嘻嘻

洛谷p3375

给出两个字符串 s_1s1​ 和 s_2s2​，若 s_1s1​ 的区间 [l, r][l,r] 子串与 s_2s2​ 完全相同，则称 s_2s2​ 在 s_1s1​ 中出现了，其出现位置为 ll。
现在请你求出 s_2s2​ 在 s_1s1​ 中所有出现的位置。

定义一个字符串 ss 的 border 为 ss 的一个非 ss 本身的子串 tt，满足 tt 既是 ss 的前缀，又是 ss 的后缀。
对于 s_2s2​，你还需要求出对于其每个前缀 s's′ 的最长 border t't′ 的长度。

输入格式

第一行为一个字符串，即为 s_1s1​。
第二行为一个字符串，即为 s_2s2​。

输出格式

首先输出若干行，每行一个整数，按从小到大的顺序输出 s_2s2​ 在 s_1s1​ 中出现的位置。
最后一行输出 |s_2|∣s2​∣ 个整数，第 ii 个整数表示 s_2s2​ 的长度为 ii 的前缀的最长 border 长度。

输入输出样例

输入 #1复制

ABABABC
ABA

输出 #1复制

1
3
0 0 1 

说明/提示

样例 1 解释

。

对于 s_2s2​ 长度为 33 的前缀 ABA，字符串 A 既是其后缀也是其前缀，且是最长的，因此最长 border 长度为 11。

数据规模与约定

本题采用多测试点捆绑测试，共有 3 个子任务。

• Subtask 1（30 points）：|s_1| \leq 15∣s1​∣≤15，|s_2| \leq 5∣s2​∣≤5。
• Subtask 2（40 points）：|s_1| \leq 10^4∣s1​∣≤104，|s_2| \leq 10^2∣s2​∣≤102。
• Subtask 3（30 points）：无特殊约定。

对于全部的测试点，保证 1 \leq |s_1|,|s_2| \leq 10^61≤∣s1​∣,∣s2​∣≤106，s_1, s_2s1​,s2​ 中均只含大写英文字母。

 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+5;
inline ll read()
{
    ll x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9')
    {
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9')
    {
        x=x*10+ch-48;
        ch=getchar();
    }
    return x*f;
}
void write(ll x)
{
    if (x < 0)
    {
        putchar('-');
        x = -x;
    }
    if (x > 9)write(x / 10);
    putchar(x % 10 + '0');
}
void prefix_table(char pattern[],int prefix[],int n)
{
    prefix[0]=0;
    int len=0;
    int i=1;
    while(i<n)
    {
        if(pattern[i]==pattern[len])
        {
            len++;
            prefix[i]=len;
            i++;
        }
        else
        {
            if(len>0)
            {
                len=prefix[len-1];
            }
            else
            {
                prefix[i]=len;
                i++;
            }
        }
    }
}
void move_prefix_table(int prefix[],int n)
{
    int i;
    for(i=n-1;i>0;i--)
    {
        prefix[i]=prefix[i-1];
    }
    prefix[0]=-1;
}
void kmp_search(char text[],char pattern[])
{
    int n=strlen(pattern);
    int m=strlen(text);
    int i=0;
    int j=0;
    int* prefix=(int*)malloc(sizeof(int)*n);
    prefix_table(pattern,prefix,n);
    move_prefix_table(prefix,n);
    while(i<m)
    {
        if(j==n-1&&text[i]==pattern[j])
        {
            printf("%d\n",i-j+1);
            j=prefix[j];
        }
        if(text[i]==pattern[j])
        {
            i++;
            j++;
        }
        else
        {
            j=prefix[j];
            if(j==-1)
            {
                i++;
                j++;
            }
        }
    }
}
int main()
{
    char pattern[maxn];
    char text[maxn];
    cin>>text;
    cin>>pattern;
    int n=strlen(pattern);
    kmp_search(text,pattern);
   int prefix[n+5];
    prefix_table(pattern,prefix,n);
    int i;
    for(i=0;i<n;i++)
    {
        printf("%d ",prefix[i]);
    }
    return 0;
}