(POJ1961)Period <KMP算法求最小循环节>

来源:互联网 发布:九品网络电视 智能电视 编辑:程序博客网 时间:2024/05/21 07:02

Period
Description

For each prefix of a given string S with N characters (each character has an ASCII code between 97 and 126, inclusive), we want to know whether the prefix is a periodic string. That is, for each i (2 <= i <= N) we want to know the largest K > 1 (if there is one) such that the prefix of S with length i can be written as AK ,that is A concatenated K times, for some string A. Of course, we also want to know the period K.
Input

The input consists of several test cases. Each test case consists of two lines. The first one contains N (2 <= N <= 1 000 000) – the size of the string S.The second line contains the string S. The input file ends with a line, having the
number zero on it.
Output

For each test case, output “Test case #” and the consecutive test case number on a single line; then, for each prefix with length i that has a period K > 1, output the prefix size i and the period K separated by a single space; the prefix sizes must be in increasing order. Print a blank line after each test case.
Sample Input

3
aaa
12
aabaabaabaab
0
Sample Output

Test case #1
2 2
3 3

Test case #2
2 2
6 2
9 3
12 4
Source

Southeastern Europe 2004

题意:
给你一个长度为n的字符串s,求他的每一个前缀的最短循环节。当循环次数大于1时,输出i,k。i表示前i个字符组成的前缀,k为循环次数。

分析:
KMP算法中,f[i]数组存的是[0,i)的最长真前缀等于最长真后缀的长度。那么若前i个字符组成的前缀是循环的,那么i%(i-f[i])==0 一定成立。 错位部分就是循环节。

AC代码:

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>using namespace std;const int maxn = 1000010;char s[maxn];int n,f[maxn];void getfail(char* p){    f[0] = 0; f[1] = 0;    for(int i=1;i<n;i++)    {        int j = f[i];        while(j && p[i] != p[j]) j = f[j];        f[i+1] = p[i]==p[j] ? j+1 : 0;    }}int main(){    scanf("%d",&n);    int kase = 1;    while(scanf("%d",&n)!=EOF && n)    {        scanf("%s",s);        getfail(s);        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]));        }        printf("\n");    }    return 0;}
0 0
原创粉丝点击