ural 1468

写了好久，不知道为什么不过，也不清楚到底卡在哪里。。。

只好看别人的代码，感觉除了HASH不一样外，倒没什么特别之处。同时参考那论文写的。。

http://blog.csdn.net/jyysc2010/article/details/9964513

 #include<stdio.h>
 #include<cstring>
 #include<cstdlib>
 #define Min(a,b) (a<b?a:b)
 struct node
 {
     int x,y;
 } hash[],a1,a2;
 int  p=;
 int n,m,f[][],ans,x1,x2,y1,y2,a[][],g[][],b[],next[],start[],tot;
 long long s[];

 void  insert(int x,node t)//将HASH值x插入链表
 {
     hash[++tot]=t;
     next[tot]=start[x];
     start[x]=tot;
 }

 bool  cmp(int x,int l)//比较HASH值为X的相同矩阵是否出现过
 {
     int i,j;
     for  (int g=start[x]; g; g=next[g])
     {
         a2=hash[g];
         j=l;
         for  (i=; i<l&&j>=l; i++)
             if  (f[a1.x][a1.y-i]!=f[a2.x][a2.y-i]) break;
         if  (i>=l)
         {
             if  (l>ans) ans=l,x1=a1.x,y1=a1.y,x2=a2.x,y2=a2.y;
             return true;
         }
     }
     return false;
 }

 int  main()
 {
     int i,j,k,l,r;
     char c;
     scanf("%d%d%*c",&n,&m);
     for  (i=; i<=n; i++)
     {
         for  (j=; j<=m; j++)
         {
             scanf("%c",&a[i][j]);
             a[i][j]-=;
         }
         scanf("%*c");
     }
     s[]=;
     for  (i=; i<=n; i++) s[i]=(p*s[i-])&;

     r=Min(n,m);
     l=;
     while  (l<=r)
     {
         k=l+r>>;
         for  (i=; i<=n; i++)
             for  (j=; j<=m; j++)
             {
                 if  (i<k) f[i][j]=f[i-][j]*p+a[i][j]&;
                 else  f[i][j]=((f[i-][j]-a[i-k][j]*s[k-]&)+&)*p+a[i][j]&;//减去最高位，加上最低位

             }//f[i][j]存下第j列第i行往上长度为k的串的HASH值
         memset(start,,sizeof(start));
         tot=;
         for  (j=; j<=m&&i>n; j++)
             for  (i=k; i<=n; i++)
             {
                 if  (j<k) g[i][j]=g[i][j-]*p+f[i][j]&;
                 else
                 {
                     g[i][j]=((g[i][j-]-f[i][j-k]*s[k-]&)+&)*p+f[i][j]&;
                     if  (cmp(g[a1.x=i][a1.y=j],k))   break;
                     insert(g[i][j],a1);
                 }
             }//g[i][j]存下右下角为第i行第j列长度为k的矩阵的HASH值
         if  (j<=m) l=k+;
         else r=k-;

     }
     if  (ans)  printf("%d\n%d %d\n%d %d\n",ans,x1-ans+,y1-ans+,x2-ans+,y2-ans+);
     else  printf("0\n");
     return ;
 }