#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<iomanip>
using namespace std;

bool intersects(int x1, int y1, int x2, int y2, int x3, int y3, int x4, int y4)
{
    int a=(x2 - x1) * (y3 - y1) - (y2 - y1) * (x3 - x1);
    int b=(x2 - x1) * (y4 - y1) - (y2 - y1) * (x4 - x1);
    if(a<0 && b<0 || a>0 && b>0)
    {
        return 0;
    }
    a=(x4 - x3) * (y1 - y3) - (y4 - y3) * (x1 - x3);
    b=(x4 - x3) * (y2 - y3) - (y4 - y3) * (x2 - x3);
    if(a<0 && b<0 || a>0 && b>0)
    {
        return 0;
    }
    return 1;
}

int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    int n,m;
    int x1,y1,x2,y2;
    cin>>n>>m>>x1>>y1>>x2>>y2;
    int** tab=new int*[n/2];
    for(int i=0;i<n/2;i++)
    {
        tab[i]=new int[m/2];
    }
    for(int i=0;i<m/2;i++)
    {
        for(int j=0;j<n/2;j++)
        {
            cin>>tab[j][i];
        }
    }
    float odp=0;
    int a,b,c,d;
    a=min(x1-1,x2-1);
    b=max(x1+1,x2+1);
    c=min(y1-1,y2-1);
    d=max(y1+1,y2+1);
    for(int i=a+2;i<b;i+=2)
    {
        for(int j=c;j<d;j+=2)
        {
            if(intersects(x1,y1,x2,y2,i,j,i,j+2))
            {
                odp+=abs(tab[(i/2)-1][j/2]-tab[i/2][j/2]);
            }
        }
    }
    for(int i=a;i<b;i+=2)
    {
        for(int j=c+2;j<d;j+=2)
        {
            if(intersects(x1,y1,x2,y2,i,j,i+2,j))
            {
                odp+=abs(tab[i/2][(j/2)-1]-tab[i/2][j/2]);
            }
        }
    }
    for(int i=a+2;i<b;i+=2)
    {
        for(int j=c+2;j<d;j+=2)
        {
            if(intersects(x1,y1,x2,y2,i,j,i,j))
            {
                int maks=max(max(tab[(i/2)-1][(j/2)-1],tab[(i/2)][(j/2)-1]),max(tab[(i/2)-1][(j/2)],tab[(i/2)][(j/2)]));
                odp-=abs(tab[(i/2)-1][(j/2)-1]-tab[(i/2)-1][(j/2)]) + abs(tab[(i/2)-1][(j/2)-1]-tab[(i/2)][(j/2)-1]) + abs(tab[(i/2)][(j/2)]-tab[(i/2)-1][(j/2)]) + abs(tab[(i/2)][(j/2)]-tab[(i/2)][(j/2)-1]);
                if(intersects(x1,y1,x2,y2,i-1,j-2,i-1,j) || intersects(x1,y1,x2,y2,i-2,j-1,i,j-1) )
                {
                    odp+=abs(tab[(i/2)-1][(j/2)-1]-maks)+abs(tab[(i/2)][(j/2)]-maks);
                }else
                {
                    odp+=abs(tab[(i/2)][(j/2)-1]-maks)+abs(tab[(i/2)-1][(j/2)]-maks);
                }
            }
        }
    }
    odp+=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
    cout<<fixed<<setprecision(6)<<odp<<endl;
    return 0;
}