思路:KMP很好的一道题。首先易证:最小覆盖子矩阵一定靠左上角。那么,我们考虑求出每一行的最小重复串长度,所有行的最小重复串的长度的lcm就是最小重复子矩阵的宽。然后我们对列也做相同的操作。于是我们就可以求得最小重复子矩阵的大小了。(这里要注意一点:当所得的宽大于原来的宽时,就让等于原来的宽,长也如此)。算法实现:算法的核心在于高效的求出每一行和每一列的最小重复串,这个可以最原串做一次KMP中的get_next()。(注:以上部分为转载)最小覆盖就等i-next[i];
代码:
#include<stdio.h>
#define N 10005
int r,l,next[N];
char s[N][N];
int getr(int n)
{
int i,j;
i=0;j=-1;
next[0]=-1;
while(i<l)
{
if(j==-1||s[n][i]==s[n][j])
{
i++;j++;
next[i]=j;
}
else
j=next[j];
}
return i-next[i];
}
int getl(int n)
{
int i,j;
i=0;j=-1;
next[0]=-1;
while(i<r)
{
if(j==-1||s[i][n]==s[j][n])
{
i++;j++;
next[i]=j;
}
else
j=next[j];
}
return i-next[i];
}
int gcd(int a,int b)
{
return b==0 ? a : gcd(b,a%b);
}
int lcm(int a,int b)
{
return a/gcd(a,b)*b;
}
int main()
{
int i,j,ans1,ans2;
while(scanf("%d%d",&r,&l)!=EOF)
{
for(i=0;i<r;i++)
scanf("%s",s[i]);
ans1=1;
for(i=0;i<r;i++)
ans1=lcm(ans1,getr(i));
if(ans1>l)//注意判断
ans1=l;
ans2=1;
for(j=0;j<l;j++)
ans2=lcm(ans2,getl(j));
if(ans2>r)
ans2=r;
printf("%d\n",ans1*ans2);
}
return 0;
}