Period MP算法

         本来该题，不会难但是在学习的过程中看到很多人把此题归入了KMP求Next[]数组当中，为了提醒后来者的注意才写了该题。

       此题题意明显没有什么好讲的，就是用我们一般很少见到的MP算法当中的状态转移图，从而来求得周期。

       而其中在求得了f[]数组，即配对失败的数组后，我们可以根据周期性质得到 如果存在周期则 I%（I-f[I] ） ==  0   必定成立。

 具体的解题过程自己看代码吧！！

#include <stdio.h>

const int MAX = 1000000+5;
char pattern[MAX];
int f[MAX];
int n;

void GetFail()
{
    f[0] = f[1] = 0;
    for(int i = 1;i < n;i++)
    {
        int j = f[i];
        while(j&&pattern[i]!=pattern[j]) j = f[j];  //一直往下查找配对的
        f[i+1] = (pattern[i] == pattern[j] ? j+1 : 0); //配对成功+1,否则归零在找
    }
}

int main()
{
    int kase = 1;
    while(scanf("%d",&n),n)
    {
        scanf("%s",pattern);
        GetFail();
        printf("Test case #%d\n",kase++);
        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("");
    }
    return 0;
}

 

另一种方法即使很多博客上所说的求Next[]数组。其实也就是KMP的未优化版，即上面所说的MP算法。

代码如下：

#include <stdio.h>

const int MAX = 1000000+5;
char pattern[MAX];
int n,next[MAX];

void GetNext()
{
    int j = 0,k = -1;
    next[0] = -1;
    while(j<n)
    {
        if(k==-1||pattern[j]==pattern[k]){
            ++j; ++k;
            //if(pattern[j]!=pattern[k])
              next[j] = k;
            //else
            //  next[j] = next[k]; // 以上的注释的三处即为KMP相对于MP优化之处。
        }else
        {
            k = next[k];
        }
    }
}

int main()
{
    int kase = 1;
    while(scanf("%d",&n),n)
    {
        scanf("%s",pattern);
        GetNext();
        printf("Test case #%d\n",kase++);
        for(int i = 1;i < n;i++)if(next[i+1]>0&&(i+1)%(i+1-next[i+1])==0)
        {
            printf("%d %d\n",i+1,(i+1)/(i+1-next[i+1]));
        }
        puts("");
    }
    return 0;
}