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

using T = long long; //o double o long double
typedef complex<T> pt;
const T EPS = 1e-10; //per i double o long double
#define X real()
#define Y imag()
#define M_PI acos(-1)
// *** Punti ***
T dot(pt v, pt w) { return (conj(v) * w).X; }
T cross(pt v, pt w) { return (conj(v) * w).Y; }
T orient(pt a, pt b, pt c) { return cross(b-a, c-a); }
pt translate(pt v, pt p) { return p+v; }
pt perp(pt p) { return {-p.Y, p.X}; }
bool isPerp(pt v, pt w) { return dot(v,w) == 0; }
pt scale(pt c, T factor, pt p) { return c+(p-c)*factor; }

struct line {
    pt v; T c;
    // From direction vector v and offset c
    line(pt v, T c) : v(v), c(c) {}
    // From equation ax+by=c
    line(T a, T b, T c) : v({b,-a}), c(c) {}
    // From points P and Q
    line(pt p, pt q) : v(q-p), c(cross(v,p)) {}
    T side(pt p) { return cross(v,p)-c; }
    T dist(pt p) {return abs(side(p)) / abs(v);}
};



void solve(int t){
    long long h,w,sx,sy,ex,ey;
    cin>>h>>w>>sx>>sy>>ex>>ey;
    
    line l({sx,sy},{ex,ey});
    
    vector mat(h,vector(w,0ll));
    for(int i=0;i<h/2;i++){
        for(int j=0;j<w/2;j++){
            cin>>mat[i][j];
        }
    }
    
    long double ans=0;
    for(int i=2;i<w;i+=2){
        for(int j=0;j<h;j+=2){
            
            pt p1(j,i),p2(j+2,i);
            
            if(l.side(p1)*l.side(p2)<0){
                ans+=abs(mat[j/2][i/2-1]-mat[j/2][i/2]);
            }
        }
    }
    
    for(int i=2;i<h;i+=2){
        for(int j=0;j<w;j+=2){
            
            pt p1(i,j),p2(i,j+2);
            
            if(l.side(p1)*l.side(p2)<0){
                ans+=abs(mat[i/2-1][j/2]-mat[i/2][j/2]);
            }
        }
    }
    
    for(int i=2;i<h;i+=2){
        for(int j=2;j<w;j+=2){
            pt p(i,j);
            
            if(l.side(p)==0){
                long long maxi = max({mat[i/2][j/2],mat[i/2-1][j/2],mat[i/2][j/2-1],mat[i/2-1][j/2-1]});
                long long part,dest;
                
                if((sx<i && sy<j) || (ex<i && ey<j)){
                    part = mat[i/2][j/2];
                    dest = mat[i/2-1][j/2-1];
                }else{
                    part = mat[i/2-1][j/2];
                    dest = mat[i/2][j/2-1];
                }
                    
                ans+=abs(maxi-part)+abs(maxi-dest);
            }
        }
    }
    
    ans+=sqrtl((sx-ex)*(sx-ex)+(sy-ey)*(sy-ey));
    cout<<setprecision(20)<<ans<<"\n";
}

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    
    int t=1;
    //cin>>t;
    for(int i=1;i<=t;i++)solve(i);
    
    return 0;
}