KMP

最近效率不高不低吧，在写usaco的题，然而真的很水==也没什么可以放上来的东西

边复习一下其他模板

KMP

http://www.matrix67.com/blog/archives/115

利用i,j两个指针检索两个字符串:s1[i-j+1..i]==s2[1..j]

s1[i+1]!=s2[j+1]，失配

next[j]应该是所有满足s2[1..next[j]]=s2[j-next[j]+1..j]的最大值

对s2预处理出next数组，在失配时直接跳到next[i]位

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1001000;
const int maxm=1001000;
char s1[maxn],s2[maxm];
int next[maxm],n,m;
int main(){
    scanf("%s%s",s1+1,s2+1);
    n=strlen(s1+1),m=strlen(s2+1);
    int j=0;
    memset(next,0,sizeof(next));
    for (int i=2;i<=m;i++){
        while (j>0&&s2[j+1]!=s2[i]) j=next[j];
        if (s2[j+1]==s2[i]) j++;
        next[i]=j;
    }
    j=0;
    for (int i=1;i<=n;i++){
        while (j>0&&s2[j+1]!=s1[i]) j=next[j];
        if (s2[j+1]==s1[i]) j++;
        if (j==m){
            printf("%d\n",i-j+1);
            j=next[j];
        }
    }
    for (int i=1;i<=m;i++)
      printf("%d ",next[i]);
    return 0;
}

 next数组的应用：n-next[n]为最小循环节长度

http://www.cnblogs.com/jackge/archive/2013/01/05/2846006.html

n/(n-next[n])表示循环节重复次数，不是整数时最后一个循环节残缺。

特别注意next[i]=0的情况，此时循环节默认判定为自身，要特判为否。

循环长度：while(next[i]>0&&i%(i-f[i])==0)  printf("%d %d\n",i,i/(i-next[i]));

#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1000100;
int n,next[maxn];
char s[maxn];
int main(){
    scanf("%d%s",&n,s+1);
    int j=0;
    memset(next,0,sizeof(next));
    for (int i=2;i<=n;i++){
        while (j>0&&s[j+1]!=s[i]) j=next[j];
        if (s[j+1]==s[i]) j++;
        next[i]=j; 
    }
    printf("%d\n",n-next[n]);
    return 0;
}