poj 2406 && poj 1961 (KMP)

poj2406 http://poj.org/problem?id=2406

给定一个字符串，问最多是多少个相同子串不重叠连接构成。

KMP的next数组应用。这里主要是如何判断是否有这样的子串，和子串的个数。

next数组的作用是在匹配无法进行下去时，可继续匹配的位置。这就要求next[j]必须为满足str[1..next[j]] = str[j-next[j]+1..j]的最大值。对于串的每一个前缀都取满足条件的最大值，则next[len]的值即为前一个子串的尾位置。

例如ababab，next[4]的值为2，即指向第一个b，因为str[1..2] = str[3..4]，同理next[6]的值为4，这样就使得str[1..2] = str[3..4] = str[5..6]，显然ab为满足条件的子串。

若为abababa，显然除其本身外，没有子串满足条件。而分析其next数组，next[7] = 5，next[5] = 3，next[3] = 1，即str[2..7]可由ba子串连接构成，那怎么否定这样的情况呢？很简单，若该子串满足条件，则len%sublen必为0。sunlen可由上面的分析得到为len-next[len]。

（这里可以加一个优化，长度为素数的输出为１，不过效率没见提高0 0）

因为子串是首尾相接，len/sublen即为substr的个数。

#include<cstdio>
#include<cstring>
char str[1000001] ;
int next[1000001] ;
int len ;
void get_next(){//根据已知next推出
    next[1] = 0 ;
    int j = 0 ;
    for(int i=2; i<=len; i++){
        while(j>0&&str[j+1]!=str[i])
            j = next[j] ;
        if(str[j+1]==str[i])    j ++ ;
        next[i] = j ;
    }
}
int main(){
    while(scanf("%s", str+1)!=EOF){
        if(str[1]=='.') break ;
        len = strlen(str+1) ;
        get_next() ;
        if(len%(len-next[len])==0){
            printf("%d\n", len/(len-next[len])) ;
            //for(int i=0; i<=len; i++)
                //printf("%d ", next[i]) ;
            continue ;
        }
        printf("1\n") ;
        //for(int i=0; i<=len; i++)
            //printf("%d ", next[i]) ;
    }
    return 0 ;

} 

poj1961 http://poj.org/problem?id=1961

1961与 2406只是多了枚举的一步，而且要求sublen>1。

code: 

#include<cstdio>
char str[1000001] ;
int next[1000001] ;
int len ;
void get_next(){//根据已知next推出
    next[1] = 0 ;
    int j = 0 ;
    for(int i=2; i<=len; i++){
        while(j>0&&str[j+1]!=str[i])
            j = next[j] ;
        if(str[j+1]==str[i])    j ++ ;
        next[i] = j ;
    }
}
int main(){
    int t = 0 ;
    while(~scanf("%d", &len)&&len){
        scanf("%s", str+1) ;
        t ++ ;
        get_next() ;
        printf("Test case #%d\n", t) ;
        for(int i=2; i<=len; i++){
            if(next[i]&&i%(i-next[i])==0)//K>1,所以加上next[i]!=0的判断
                printf("%d %d\n", i, i/(i-next[i])) ;
        }
        printf("\n") ;
    }

}