#include <iostream>
#include <cmath>
#include <algorithm>
#include <iomanip>
#include <vector>
using namespace std;
int main(){
  int h,w;cin>>h>>w;
  vector<vector<char>> A(h,vector<char>(w));
  for(int i = 0; h > i; i++){
    for(int j = 0; w > j; j++){
      cin>>A[i][j];
    }
  }
  long double ans = 1e9;
  for(int i = 0; h > i; i++){
    // -1,w+1
    long double tmp = 0,tmp2 = 0;
    for(int j = 0; h > j; j++){
      for(int k = 0; w > k; k++){
        if(A[j][k] == '1'){
          tmp += sqrt((k+1)*(k+1)+(j-i)*(j-i));
          tmp2 += sqrt((w-k)*(w-k)+(j-i)*(j-i));
        }
      }
    }
    ans = min(ans,min(tmp,tmp2));
  }
  for(int i = 0; w > i; i++){
    // -1,w+1
    long double tmp = 0,tmp2 = 0;
    for(int j = 0; h > j; j++){
      for(int k = 0; w > k; k++){
        if(A[j][k] == '1'){
          tmp += sqrt((j+1)*(j+1)+(k-i)*(k-i));
          tmp2 += sqrt((h-j)*(h-j)+(k-i)*(k-i));
        }
      }
    }
    ans = min(ans,min(tmp,tmp2));
  }
  cout << fixed << setprecision(10) << ans << endl;
}