3
aaa
12
aabaabaabaab
0

Test case #1
2 2
3 3

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

//一开始完全没有看懂题目的意思，就是说从字符串的第二个开始，看前面是循环串吗，是的话就输出此时的位置，和循环串的周期，周期必须大于1

/*
对next的理解更深入了点儿。
    字符编号从0开始，那么if(i%(i-next[i])==0)，则i前面的
串为一个轮回串，其中轮回子串出现i/(i-next[i])次。
*/

#include <string.h>
#include <iostream>
#include <stdio.h>

char str[1000005];
int next[1000005];

void getnext()
{
    int i = 0,j = -1;
    memset(next,0,sizeof(next));
    next[0] = -1;
    while (str[i])
    {
        if(j == -1 || str[i] == str[j])
        {
            i++;
            j++;
            next[i] = j;
        }
        else
        j = next[j];
    }
}

void kmp()
{
    int i,t;
    for(i = 2;str[i-1];i++)
    {
        t = i-next[i];
        if(i%t == 0 && i/t>1)
        printf("%d %d\n",i,i/t);
    }
}

int main()
{
    int n,cnt = 1;

    while(scanf("%d",&n)!=EOF && n)
    {
        scanf("%s",str);
        printf("Test case #%d\n",cnt++);
        getnext();
        kmp();
        putchar(10);
    }

    return 0;
}

%关于相位误差
clear
close all
clc

fs=100; %采样率
fc=10;  %载波
fpe=1;  %相位噪声

t=0:1/fs:100;
c_pe=cos(2*pi*fpe*t);
c_c=cos(2*pi*fc*t);
c_cpe=cos(2*pi*fc*t+c_pe);

C_C=fft(c_c);
C_C=abs(C_C)/length(C_C);
C_CPE=fft(c_cpe);
C_CPE=abs(C_CPE)/length(C_CPE);

F=linspace(0,fs,length(t));
N=floor(length(F)/2);

plot(F(1:N),C_C(1:N),'r');
hold on
plot(F(1:N),C_CPE(1:N));
set(gca,'xtick',[0:1:F(N)]);
grid on;