LA 3026 Peroid

题目链接：https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1027

本题联系字符串匹配算法KMP。

KMP算法的核心在于失效函数f[]的求得。f[i]代表Pi的真后缀P的最长前缀的长度，注意如果i是从0开始的，那么f[i]实际代表了i-1之前的字符串匹配效果。

核心函数如下：

void getFail(char *p,int *f)
{
    int m = strlen(p);
    f[0] = f[1] = 0;
    int j = 0;
    for(int i=1;i<m;i++)
    {
        j = f[i];
        while(j && p[i]!=p[j]) j = f[j];
        if(p[i] == p[j]) f[i+1] = j + 1;
        else f[i+1] = 0;
    }
}
void kmp(char * T,char *p, int *f)
{
    int n = strlen(T);
    int m = strlen(p);
    getFail(p,f);
    int j = 0;
    for(int i = 0;i<n;i++)
    {
        while(j && p[j]!=T[i]) j = f[j];
        if(p[j] == T[i]) j++;
        if(j == m) printf("%d\n",i-m+1);
    }
}

本题而言，错位部分长度为i-f[i]，如果i能被错位部分整除，那么说明到目前为止的字符串是周期性质的。 

另外，如果数组开的很大，最好放在全局，因为即便是main函数，也是有栈限制的，太大会爆栈。

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <stdlib.h>
#include <time.h>
#include <vector>
#include <set>
#include <map>
#include <algorithm>

using namespace std;
#define Maxn 1000005

void getFail(char *p,int *f)
{
    int m = strlen(p);
    f[0] = f[1] = 0;
    int j = 0;
    for(int i=1;i<m;i++)
    {
        j = f[i];
        while(j && p[i]!=p[j]) j = f[j];
        if(p[i] == p[j]) f[i+1] = j + 1;
        else f[i+1] = 0;
    }
}

char p[Maxn];
//数组开的大的，声明在外部，否则会爆栈
int f[Maxn];

int main()
{

#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
#endif
    int n;
    int cas = 0;
    while(scanf(" %d",&n)!=EOF && n!=0)
    {
        cas++;
        printf("Test case #%d\n",cas);
        scanf(" %s",p);
        getFail(p,f);
        for(int i=2;i<=n;i++)
        {
            if(f[i]>0 && i%(i-f[i])==0) printf("%d %d\n",i,i/(i-f[i]));
        }
        puts("");
    }
}