#include <bits/stdc++.h>

using namespace std;
ifstream in("bmatrix.in");
ofstream out("bmatrix.out");
const int NMAX = 201;

int a[NMAX][NMAX], dir[NMAX][NMAX];
int directii[4][NMAX];
int n, m;
void debug( int N, int M){
    for( int l = 1; l <= N; l++ ){
        for( int c = 1; c <= M; c++ )
            out << dir[l][c];
        out << "\n";
    }
    out << "\n";
}
int skyline( deque <int> dq ){
    int x = dq.size(), aria = 0;
    for( int i = 0; i < x; i++ )
        aria = max( aria, dq.front() * (x-i) ), dq.pop_front();
    return aria;
}
void bmatrix( int N, int M, int mat[NMAX][NMAX], int tip){
    int up[NMAX];
    for( int i = 1; i <= M; i++ )
        up[i] = 0;
    for( int l = 1; l <= N; l++ ){
        int aria = 0;
        deque <int> dq;
        for( int c = 1; c <= M; c++ ){
            if( mat[l][c] == 0 ){
                up[c]++;
                aria = max( aria, up[c] );
                int lenght = 1;
                while( dq.size() && up[c] < dq.back() ){
                    dq.pop_back();
                    lenght++;
                }
                aria = max( aria, lenght * up[c]);
                dq.push_back(up[c]);
                aria = max( aria, skyline(dq) );
            }
            else{
                up[c] = 0;
                dq.clear();
            }
        }
        directii[tip][l] = max(aria, directii[tip][l-1]);
    }
}
int main()
{
    in >> n >> m;
    for( int l = 1; l <= n; l++ ){
        char ch;
        for( int c = 1; c <= m; c++ ){
            in >> ch;
            a[l][c] = ch - '0';
        }
    }
    for( int l = 1; l <= n; l++ )
        for( int c = 1; c <= m; c++ )
            dir[l][c] = a[l][c];
//    debug(n, m);
    bmatrix(n, m, dir, 0);

    for( int l = 1; l <= n/2; l++ )
        for( int c = 1; c <= m; c++ ){
            dir[l][c] = a[n-l+1][c];
            dir[n-l+1][c] = a[l][c];
        }
//    debug(n, m);
    bmatrix(n, m, dir, 1);

    for( int l = 1; l <= m; l++ )
        for( int c = 1; c <= n; c++ )
            dir[l][c] = a[n-c+1][l];
//    debug(m, n);
    bmatrix(m, n, dir, 2);

    for( int l = 1; l <= m; l++ )
        for( int c = 1; c <= n; c++ )
            dir[l][c] = a[c][m-l+1];
    bmatrix(m, n, dir, 3);


    int rez = 0;
    for( int i = 2; i < n; i++ ){
        rez = max( rez, directii[0][i] + directii[1][n-i] );
//        out << lin_sus[i] << " " << lin_jos[i] << "\n";
    }
//    out << "\n";
    for( int i = 2; i < m; i++ ){
        rez = max( rez, directii[2][i] + directii[3][m-i]);
//        out << col_dr[i] << " " << col_st[i] << "\n";
    }
    out << rez;
    return 0;
}