#pragma GCC optimize("O2")
#include<bits/stdc++.h>
#define ll   long long
#define pb   emplace_back
#define mp   make_pair
#define orz  1000000007
#define fi   first
#define se   second
using namespace std;
int n,m;
string s;
int a[1000005],ans=orz;
int S,dp[1000005],bac[1048666];
pair<int,int> p[1000005];
int sol(int l,int m,int r){
	if(l>r)swap(l,r);
	int x=(l+m-1)/m;
	int k=bac[x];
	int len=l+((1<<k)-x)*m;
	k+=bac[(l+m+r-1)/len];
	return k;
}
int main(){
	for(int i=1;i<=20;++i){
		for(int j=(1<<(i-1));j<(1<<i);++j)bac[j]=i;
	}
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;++i){
		cin>>s;
		s=" "+s;
		for(int j=1;j<=m;++j){
			if(s[j]=='1')++a[j];
			else a[j]=0;
		}
		S=0;
		for(int j=1;j<=m;++j){
			if(a[j]==0)S=0;
			else{
				int l=j;
				while(p[S].fi>=a[j]){
					l=p[S].se;
					--S;
				}
				p[++S]=mp(a[j],l);
				dp[j]=l;
			}
		}
		S=0;
		for(int j=m;j;--j){
			if(a[j]==0)S=0;
			else{
				int r=j;
				while(p[S].fi>=a[j]){
					r=p[S].se;
					--S;
				}
				p[++S]=mp(a[j],r);
				int x,y,l=dp[j];
				x=sol(l-1,r-l+1,m-r);
				y=sol(i-a[j],a[j],n-i);
				ans=min(ans,x+y);
			}
		}
	}
	cout<<ans<<'\n';
    return 0;
}