POJ 2815 Milking Grid（KMP）

一个字符串的最小覆盖子串长度等于len-next[len]

证明：http://blog.csdn.net/fjsd155/article/details/6866991 




先算出每一行的最小覆盖长度，然后所有行的最小覆盖长度就是这些数的最小公倍数，设为r。

因为列和行相互独立，所以同样处理。


最后结果就是l*r


注意：如果r>R 那么将r设为R，l同样。


代码：
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

int next[10005][80];
int R,C;
char mp[10005][80];

void getnextR(int k){
    for(int i=1;i<C;i++){
        int j=next[k][i];
        while(j&&mp[k][i]!=mp[k][j]) j=next[k][j];
        next[k][i+1]=(mp[k][i]==mp[k][j]? j+1:0);
    }
}

void getnextC(int k){
    for(int i=1;i<R;i++){
        int j=next[i][k];
        while(j&&mp[i][k]!=mp[j][k]) j=next[j][k];
        next[i+1][k]=(mp[i][k]==mp[j][k]? j+1:0);
    }
}


int gcd(int a,int b){
    if(b==0) return a;
    return gcd(b,a%b);
}

int lcm(int a,int b){
    return a*b/gcd(a,b);
}


int main(){
    scanf("%d%d",&R,&C);
    for(int i=0;i<R;i++){
        scanf("%s",mp[i]);
    }
    int l,w;
    memset(next,0,sizeof(next));
    l=1;
    for(int i=0;i<R;i++){
        getnextR(i);
        l=lcm(l,C-next[i][C]);
    }
    if(l>C) l=C;
    memset(next,0,sizeof(next));
    w=1;
    for(int i=0;i<C;i++){
        getnextC(i);
        w=lcm(w,R-next[R][i]);
    }
    if(w>R) w=R;
    printf("%d\n",l*w);
    return 0;
}