#include<cstdio>

 #include<iostream>

 #include<cstring>

 using namespace std;

 const int maxn=;

 const int inf=;

 int f[][][maxn][maxn];

 int sm[maxn][maxn],mx1[maxn][maxn],mx2[maxn][maxn];

 int n,m,pre,now,ans,tmpmx;

 int ra;char rx;

 inline int read(){

     rx=getchar(),ra=;

     while(rx<''||rx>'')rx=getchar();

     while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;

 }

 inline int max(int a,int b){return a>b?a:b;}

 inline int min(int a,int b){return a<b?a:b;}

 int main(){

     register int i,j,k;

     n=read(),m=read();

     for(i=;i<=n;i++)

         for(j=;j<=m;j++)sm[i][j]=sm[i][j-]+read();

     for(i=;i<=m;i++)sm[][i]=i;

     now=,pre=;

     for(i=;i<=;i++)memset(f[][i],,sizeof(f[][i]));

     for(i=;i<=n;i++,swap(now,pre)){

         for(j=;j<m;j++)for(k=j+;k<=m;k++)

             if(sm[i][k]==sm[i][j-])

                 f[now][][j][k]=(sm[i-][k]==sm[i-][j-]?f[pre][][j][k]:)+k-j+;

             else f[now][][j][k]=-inf;

         for(j=;j<m;j++)for(k=m,tmpmx=-inf;k>j;k--)

             tmpmx=max(tmpmx,f[pre][][j][k]),

             mx1[j][k]=max(mx1[j-][k],tmpmx);

         for(j=;j<m;j++)for(k=j;k<m;k++)

             if(sm[i][k]==sm[i][j-]){

                 f[now][][j][k]=f[pre][][j][k]+k-j+;

                 if(sm[i-][k+]==sm[i-][j-])

                     f[now][][j][k]=max(f[now][][j][k],mx1[j-][k+]+k-j+);

             }

             else f[now][][j][k]=-inf;

         for(j=;j<m;j++)mx2[j][j]=f[pre][][j][j];

         for(j=;j<m;j++)

             for(k=;k<m-j;k++)

             mx2[k][k+j]=max(mx2[k][k+j-],max(mx2[k+][k+j],f[pre][][k][k+j]));

         for(j=;j<m;j++)for(k=j+;k<=m;k++){

             if(sm[i][k]==sm[i][j-])

                 f[now][][j][k]=max(f[pre][][j][k],mx2[j+][k-])+k-j+;

             else f[now][][j][k]=-inf;

             if(f[now][][j][k]>ans)ans=f[now][][j][k];

         }

     }

     printf("%d\n",ans);

     return ;

 }