#include<bits/stdc++.h>
using namespace std;
int F(int x,int y){
    if(!x&&!y)return 0;
    if((x&1)&&(y&1)){
        if(x>y)swap(x,y);
        return 1+F((x+1)/2,y/2);
    }
    if(!(y&1))swap(x,y);
    return 1+F((x+1)/2,y/2);
}
int jx(int l,int r,int n){
    int o=r-l+1;l--,r=n-r;
    l=(l+o-1)/o,r=(r+o-1)/o;
    return F(l,r);
}
int n,m;
int s[1010000];
int id(int x,int y){if(x<=0||y<=0)return 0;return (x-1)*m+y;}
char c[1010000];
int pn[1010000][30],pm[1001000][30];
bool ck(int xl,int yl,int xr,int yr){
    if(!xl||!yl)return 0;xl--,yl--;
    return s[id(xr,yr)]+s[id(xl,yl)]==s[id(xl,yr)]+s[id(xr,yl)];
}
int checkn(int l,int r,int layer)
{
	if(!layer) return (l==1&&r==n);
	return pn[r][layer-1]>=max(1,l+l-r-1)||pn[min(n,r+r-l+1)][layer-1]>=l;
}
int checkm(int l,int r,int layer)
{
	if(!layer) return (l==1&&r==m);
	return pm[r][layer-1]>=max(1,l+l-r-1)||pm[min(m,r+r-l+1)][layer-1]>=l;
}
const int Q=22;
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%s",c+1);
        for(int j=1;j<=m;j++)
        s[id(i,j)]=1^(c[j]-'0');
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    s[id(i,j)]+=s[id(i,j-1)]+s[id(i-1,j)]-s[id(i-1,j-1)];
    for(int o=0;o<=Q;o++){
        for(int i=1,j=0;i<=n;i++){
            while(j<i&&checkn(j+1,i,o))j++;
            pn[i][o]=j;
        }
    }
    for(int o=0;o<=Q;o++){
        for(int i=1,j=0;i<=m;i++){
            while(j<i&&checkm(j+1,i,o))j++;
            pm[i][o]=j;
        }
    }
    int ans=1e9;
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++){
        for(int o1=0,o2=Q+1;o1<=Q;o1++){
            if(!pn[i][o1])continue;
            while(o2>0&&ck(pn[i][o1],pm[j][o2-1],i,j))o2--;
            if(o2!=Q+1)ans=min(ans,o1+o2);
        }
    }
    printf("%d\n",ans);
    return 0;
}