#include <bits/stdc++.h>
using namespace std;

const int N=1000050;
char s[N];
int lg[N];
int n,m,ans=N;
int Solve2(int l,int r,int n){
    int sz=r-l+1;
    int cnt=lg[(l+sz-1)/sz];
    sz=r+1;
    cnt+=lg[(n-r-1+sz-1)/sz];
    return cnt;
}
int Solve(int l,int r,int n){
    return min(Solve2(l,r,n),Solve2(n-r-1,n-l-1,n));
}
void Try(int xl,int xr,int yl,int yr){
    if(xl>xr)return;
    //printf("Try (%i %i) (%i %i)\n",xl,xr,yl,yr);
    ans=min(ans,Solve(xl,xr,n)+Solve(yl,yr,m));
}
int main(){
    for(int i=1;i<N;i++)lg[i]=lg[i>>1]+1;
    scanf("%i %i",&n,&m);
    bool sw=false;
    if(n<m)sw=true,swap(n,m);
    vector<vector<int>> a(n,vector<int>(m,0));
    if(sw){
        for(int j=0;j<m;j++){
            scanf("%s",s);
            for(int i=0;i<n;i++){
                a[i][j]=s[i]-'0';
            }
        }
    }else{
        for(int i=0;i<n;i++){
            scanf("%s",s);
            for(int j=0;j<m;j++){
                a[i][j]=s[j]-'0';
            }
        }
    }
    vector<int> h(m,0);
    for(int i=0;i<n;i++){
        vector<int> stk;
        stk.push_back(-1);
        for(int j=0;j<m;j++){
            if(a[i][j]==1)h[j]++;
            else h[j]=0;
            while(stk.size()>1 && h[stk.back()]>=h[j]){
                stk.pop_back();
            }
            stk.push_back(j);
            for(int k=0;k+1<stk.size();k++){
                Try(i-h[stk[k+1]]+1,i,stk[k]+1,j);
            }
        }
    }
    printf("%i\n",ans);
    return 0;
}