|算法讨论|KMP 学习笔记

模板及讲解

输入两个字符串s1和s2，输出s2在s1中出现的位置(位置从0开始)

常见题型
1. 匹配问题
Q：输入两个字符串s1和s2，输出s2在s1中出现的位置。
解：直接KMP比较即可。
例题：poj 3461
2.求最短循环节
Q：输入字符串s1，求s1中最短的循环节。
解：KMP中$n-f[n]$即为最短循环节。
例题：poj 2406
3.求前缀和后缀最长公共部分
Q：输入字符串s1，求s1中前缀和后缀最长公共部分的长度。
解：$f[i]$即为前缀和后缀最长公共部分长度。
例题：poj 2752

//实质上是MP算法
#include<cstdio>    
#include<algorithm>    
#include<cstring>     
#define ms(i,j) memset(i,j, sizeof i);    
using namespace std;
char s1[1000000 + 5], s2[1000 + 5]; 
int next[1000 + 5];
int n,m; 
void getNext()//获得失配函数 
{
    next[0] = 0;
    next[1] = 0;
    for (int i=1;i<m;i++)
    {
        int j = next[i];
        while (j && s2[i]!=s2[j]) j = next[j];
        next[i+1] = ((s2[i]==s2[j]) ? (j+1) : (0));
    }
}
void kmp()//在s1里找s2 
{
    getNext();
    int j = 0;
    for (int i=0;i<n;i++)
    {
        while (j && s1[i]!=s2[j]) j = next[j];
        if (s1[i]==s2[j]) j++;
        if (j==m) printf("%d\n", i-m+1);
    }
}
int main()   
{     
    scanf("%s%s", s1, s2);
    n = strlen(s1), m = strlen(s2);
    kmp();
    return 0;    
}    

//KMP求字符串中的最短周期(循环节)
#include<cstdio>    
#include<algorithm>    
#include<cstring>     
#define ms(i,j) memset(i,j, sizeof i);    
using namespace std;
const int MAXN = 1000000 + 5;
int n;
char s[MAXN];
int f[MAXN];
int main()   
{     
    while (scanf("%s", s)&&(s[0]!='.'))
    {
        n = strlen(s);
        f[0] = f[1] = 0;
        for (int i=1;i<n;i++)
        {
            int j = f[i];
            while (j && s[i]!=s[j]) j = f[j];
            f[i+1] = (s[i]==s[j]) ? (j+1) : (0);
        }
        int ans = n-f[n];
        if (n%ans==0) ans = n/ans; else ans = 1;
        printf("%d\n", ans);
    }
    return 0;    
}