【poj 2185】Milking Grid（字符串--KMP+问题分解）

题意：给定一个由字符组成的矩阵，求出它的面积最小的覆盖矩阵。（感觉应该是可重叠的......* (・ω・)っ）

解法：KMP。行列互不影响，可以问题分解。先求出每一行的最小重复串，利用kmp中的next[ ]数组，然后求其最小公倍数。若变成最大值就可以break了，这也是我代码本来应该O(2*m*n)，但却32ms很快的原因......(๑•ᴗ•๑) 于是同样地求每一列的最小重复串，再求最小公倍数。最终2个答案相乘就行了。

P.S.这题我十分纠结，也有人说poj数据水......但我还是AC了~屏蔽的都是我纠结的。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;

const int N=10010,M=80;
char s[N][M];
int n,m;
int next[N];//,tend[M];

int gcd(int a,int b) {return b?gcd(b,a%b):a;}
int lcm(int x,int y) {return x*y/gcd(x,y);}
int kmp()
{
    int p,xx=1,yy=1;
    for (int i=1;i<=n;i++)
    {
      memset(next,0,sizeof(next));
      next[1]=0; p=0;
      for (int j=2;j<=m;j++)
      {
        while (s[i][j]!=s[i][p+1] && p) p=next[p];
        if (s[i][j]==s[i][p+1]) p++;
        next[j]=p;
      }
      /*int t=m-p,tt;
      if (t%(t-next[t])!=0) tt=t;
      else tt=t-next[t];*/
      //if (m%(m-next[m])!=0) {xx=m;break;}
      //else 
      xx=lcm(xx,m-next[m]);
      if (xx>=m) {xx=m;break;}
    }
    for (int i=1;i<=m;i++)
    {
      memset(next,0,sizeof(next));
      next[1]=0; p=0;
      for (int j=2;j<=n;j++)
      {
        while (s[j][i]!=s[p+1][i] && p) p=next[p];
        if (s[j][i]==s[p+1][i]) p++;
        next[j]=p;
      }
      //if (n%(n-next[n])!=0) {yy=n;break;}
      //else 
      yy=lcm(yy,n-next[n]);
      if (yy>=n) {yy=n;break;}
    }
    return xx*yy;
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
      scanf("%s",s[i]+1);
    printf("%d\n",kmp());
    return 0;
}