#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<iomanip>
#include<math.h>
#include<complex>
#include<queue>
#include<deque>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define REP(i,m,n) for(int i=(int)m ; i < (int) n ; i++ )
#define rep(i,n) REP(i,0,n)
typedef long long ll;
typedef pair<int,int> pint;
const int inf=1e9+7;
const ll longinf=1LL<<60 ;
const ll mod=998244353 ;
int dx[4]={1,0,-1,0} , dy[4]={0,1,0,-1} ;

int main(){
    int h,w;
    cin>>h>>w;
    string s[h];
    rep(i,h)cin>>s[i];
    double ans=100000;
    rep(i,h+2){
        double ret=0;
        rep(j,h)rep(k,w){
            if(s[j][k]=='1')ret+=sqrt((k+1)*(k+1)+(j-i+1)*(j-i+1));
        }
        ans=min(ans,ret);
    }
    rep(i,h+2){
        double ret=0;
        rep(j,h)rep(k,w){
            if(s[j][k]=='1')ret+=sqrt((w-k)*(w-k)+(j-i+1)*(j-i+1));
        }
        ans=min(ans,ret);
    }
    rep(i,w+2){
        double ret=0;
        rep(j,h)rep(k,w){
            if(s[j][k]=='1')ret+=sqrt((k-i+1)*(k-i+1)+(j+1)*(j+1));
        }
        ans=min(ans,ret);
    }
    rep(i,w+2){
        double ret=0;
        rep(j,h)rep(k,w){
            if(s[j][k]=='1')ret+=sqrt((k-i+1)*(k-i+1)+(h-j)*(h-j));
        }
        ans=min(ans,ret);
    }
cout<<fixed<<setprecision(10)<<ans<<endl;
return 0;
}