[POJ2185]Milking Grid（KMP求解最小覆盖问题）

传送门

题意：
求一个字符串矩阵的最小覆盖矩阵。覆盖矩阵是指，字符串矩阵的子矩阵通过复制平移的方式可以形成一个新的大矩阵使得原字符串矩阵成为他的子矩阵。如:
　　　ABABA
　　　CDCDC
他的最小覆盖矩阵为
　　　AB
　　　CD
因为将上面的矩阵复制平移后可以得到这样一个新矩阵
　　　ABABAB
　　　CDCDCD
则原矩阵是这个新矩阵的子矩阵。

---

　　最小覆盖矩阵可以转化为行和列的最小覆盖问题解决。 把每一行的最小覆盖子串的长度求出来，这一个矩阵最小覆盖矩阵的长就是他们（每一行的最小覆盖子串的长度）的最小公倍数。 对于列做相同的操作，就可以得到这个最小覆盖矩阵的长和宽了。
　　那么新的问题是，如何高效的求每一行的最小覆盖子串长度，我们可以用kmp算法间接求解。 对于一个串我们把它自己kmp一边，然后用n-next[n]就可以的得到最小覆盖子串长度（n为字符串长度）
　　具体的证明可以看这两篇博客：
　　《字符串 - KMP》
　　《KMP与最小覆盖子串》
　　

　　

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int R=10010;
const int C=80;
char s[R][C];
int next1[C],next2[R];
int r,c;
int Gcd(int a,int b)
{
    if(b==0) return a;
    else return Gcd(b,a%b); 
}
int Lcm(int a,int b)
{
    if(b>a) swap(a,b);
    return (a*b)/Gcd(a,b);
}
void get_row_next(int r)
{
    memset(next1,0,sizeof(next1));
    for(int i=2,j=0;i<=c;i++)
    {
        while(j>0 && s[r][i]!=s[r][j+1]) j=next1[j];
        if(s[r][i]==s[r][j+1]) j++;
        next1[i]=j; 
    }
}
void get_col_next(int c)
{
    memset(next2,0,sizeof(next2));
    for(int i=2,j=0;i<=r;i++)
    {
        while(j>0 && s[i][c]!=s[j+1][c]) j=next2[j];
        if(s[i][c]==s[j+1][c]) j++;
        next2[i]=j; 
    }   
}
int main()
{
    scanf("%d%d",&r,&c);
    int row=1,col=1;
    for(int i=1;i<=r;i++)
    {
        scanf("%s",s[i]+1);
    }
    for(int i=1;i<=r;i++)
    {
        get_row_next(i);
        row=Lcm(row,c-next1[c]);
        if(row>=c){row=c; break;}
    }
    for(int i=1;i<=c;i++)
    {
        get_col_next(i);
        col=Lcm(col,r-next2[r]);
        if(col>=r){col=r; break;}
    }
    printf("%d\n",row*col);
}