#pragma GCC optimize("Ofast")
#include "bits/stdc++.h"
#define FOR(i,a,b)for(int i=(a);i<(b);i++)
#define F0R(i,a)FOR(i,0,a)
#define ROF(i,a,b)for(int i=(b)-1;i>=(a);i--)
#define R0F(i,a)ROF(i,0,a)
#define REP(a)F0R(_,a)
using namespace std;
const int mxn = 2501;
string a[mxn];
int dp[2][mxn][mxn];
int w[mxn];
int att[mxn];
int main() {
cin.tie(nullptr)->sync_with_stdio(false);
int n, m; cin >> n >> m;
F0R(i,n)cin>>a[i];
int mw = INT_MAX;
F0R(j,m){
F0R(i,n)if(a[i][j]=='1'){
dp[0][i][j]=1;
if(i)dp[0][i][j] += dp[0][i-1][j];
}
R0F(i,n)if(a[i][j]=='1'){
dp[1][i][j]=1;
if(i<n-1)dp[1][i][j] += dp[1][i+1][j];
}
F0R(i,n){
if(a[i][j]=='1')att[i]++;
else {
if(att[i])mw=min(mw,att[i]);
att[i]=0;
}
}
}
F0R(j,m)if(att[j])mw=min(mw,att[j]);
int mh=INT_MAX;
F0R(i,n)F0R(j,m)if(a[i][j]=='1')mh=min(mh,dp[0][i][j]+dp[1][i][j]-1);
/*F0R(i,n){
F0R(j,m)cout<<dp[1][i][j]<<' ';cout<<endl;
}
cout << mw << ' ' << mh << endl;*/
fill_n(w,mw+1,mh);
F0R(i,n){
int ln=0,ma=INT_MAX,mb=INT_MAX;
F0R(j,m){
if (a[i][j]=='1'){
ln++,ma=min(ma,dp[0][i][j]),mb=min(mb,dp[1][i][j]);
w[ln]=min(w[ln],ma+mb-1);
}else ln=0;
}
}
int ans=0;
F0R(i,mw+1)ans=max(ans,i*w[i]);
cout << ans << endl;
}
