#include <fstream>
#include <climits>
using namespace std;
ifstream f("bmatrix.in");
ofstream g("bmatrix.out");
int a[205][205],st[205],i,j,n,m,dr[205],sus[205],jos[205],v[205],p,l,sol,y,v1[205];
char c;
int main()
{
    f>>n>>m;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
    {
        f>>c;
        a[i][j]=c-'0';
    }
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
    {
        if(a[i][j]==0) v[j]++;
        else v[j]=0;
        p=INT_MAX;
        y=0;
        for(l=j;l>=1;l--)
        {
            p=min(p,v[l]);
            y=max(y,p*(j-l+1));
        }
        st[i]=max(st[i],max(st[i-1],y));
        sus[j]=max(sus[j],max(sus[j-1],y));

    }
    for(i=n;i>=1;i--)
        for(j=m;j>=1;j--)
    {
        if(a[i][j]==0) v1[j]++;
        else v1[j]=0;
        p=INT_MAX;
        y=0;
        for(l=j;l<=m;l++)
        {
            p=min(p,v1[l]);
            y=max(y,p*(l-j+1));
        }
       dr[i]=max(dr[i],max(dr[i+1],y));
        jos[j]=max(jos[j],max(jos[j+1],y));

    }
    for(i=1;i<n;i++)
    {
        sol=max(sol,st[i]+dr[i+1]);
    }
    for(i=1;i<m;i++)
    {
        sol=max(sol,sus[i]+jos[i+1]);
    }
    g<<sol;
    return 0;
}