#include<bits/stdc++.h>

using namespace std;
int n,q;
int a[2][505],b[2][505];
int w[2][505][505];
struct query{
    int l,r,L,R;
}Q[300005];
long long ans[300005];
long long f[505][505][2][2],dp[505][505];
long long sav[505];
void Solve(int l,int r,int L,int R,vector<int>&vt,int op){
    if(vt.empty())return;
    // cerr<<l<<" "<<r<<" "<<L<<" "<<R<<endl;
    if(l==r){
        for(auto x:vt){
            long long sum=-a[0][l];
            for(int i=Q[x].L;i<=Q[x].R;i++)sum+=max(0,-a[1][i]+w[0][l][i]);
            ans[x]=max(sum,0ll);   
        }
        return;
    }
    if(L==R){
        for(auto x:vt){
            long long sum=-a[1][L];
            for(int i=Q[x].l;i<=Q[x].r;i++)sum+=max(0,-a[0][i]+w[0][i][L]);
            ans[x]=max(sum,0ll);   
        }
        return;
    }
    if(op){
        swap(l,L);swap(r,R);
        for(auto x:vt)
            swap(Q[x].l,Q[x].L),swap(Q[x].r,Q[x].R);
    }
    int mid=l+r>>1;
    vector<int>vm,vl,vr;
    for(auto x:vt)
        if(Q[x].l<=mid&&Q[x].r>mid)vm.push_back(x);
        else if(Q[x].r<=mid)vl.push_back(x);
        else vr.push_back(x);
    if(!vm.empty()){
        for(int m=L;m<=R;m++){
            for(int s=0;s<2;s++){
                for(int i=l;i<=r;i++)
                    for(int j=L;j<=R;j++)
                        for(int u=0;u<2;u++)
                            for(int v=0;v<2;v++)
                                f[i][j][u][v]=-1e18;
                for(auto x:vm)sav[x]=(s?-a[op^1][m]:0);
                f[mid][m][0][s]=0;
                for(int i=mid;i>=l;i--)
                    for(int j=m;j>=L;j--){
                        dp[i][j]=-1e18;
                        for(int u=0;u<2;u++)
                            for(int v=0;v<2;v++){
                                dp[i][j]=max(dp[i][j],f[i][j][u][v]);
                                if(!u||!v)f[i][j][1][1]=max(f[i][j][1][1],f[i][j][u][v]+w[op][i][j]-(u?0:a[op][i])-(v?0:a[op^1][j]));
                                f[i-1][j][0][v]=max(f[i-1][j][0][v],f[i][j][u][v]);
                                f[i][j-1][u][0]=max(f[i][j-1][u][0],f[i][j][u][v]);
                            }
                    }
                for(auto x:vm)
                    if(Q[x].L<=m&&Q[x].R>=m)sav[x]+=dp[Q[x].l][Q[x].L];
                f[mid+1][m][0][s]=0;
                for(int i=mid+1;i<=r;i++)
                    for(int j=m;j<=R;j++){
                        dp[i][j]=-1e18;
                        for(int u=0;u<2;u++)
                            for(int v=0;v<2;v++){
                                dp[i][j]=max(dp[i][j],f[i][j][u][v]);
                                if(!u||!v)f[i][j][1][1]=max(f[i][j][1][1],f[i][j][u][v]+w[op][i][j]-(u?0:a[op][i])-(v?0:a[op^1][j]));
                                f[i+1][j][0][v]=max(f[i+1][j][0][v],f[i][j][u][v]);
                                f[i][j+1][u][0]=max(f[i][j+1][u][0],f[i][j][u][v]);
                            }
                    }
                for(auto x:vm)
                    if(Q[x].L<=m&&Q[x].R>=m)sav[x]+=dp[Q[x].r][Q[x].R];
                for(auto x:vm)
                    if(Q[x].L<=m&&Q[x].R>=m)ans[x]=max(ans[x],sav[x]);
            }
        }
    }
    if(op){
        swap(l,L);swap(r,R);
        for(auto x:vt)
            swap(Q[x].l,Q[x].L),swap(Q[x].r,Q[x].R);
        Solve(l,r,L,mid,vl,op^1);
        Solve(l,r,mid+1,R,vr,op^1);
    }
    else {
        Solve(l,mid,L,R,vl,op^1);
        Solve(mid+1,r,L,R,vr,op^1);
    }
}
int main(){
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++)scanf("%d",&a[0][i]);
    for(int i=1;i<=n;i++)scanf("%d",&a[1][i]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&w[0][i][j]),w[1][j][i]=w[0][i][j];
    vector<int>vt;
    for(int i=1;i<=q;i++){
        scanf("%d%d%d%d",&Q[i].l,&Q[i].r,&Q[i].L,&Q[i].R);
        vt.push_back(i);
    }
    Solve(1,100,1,100,vt,0);
    for(int i=1;i<=q;i++)printf("%lld\n",ans[i]);
    return 0;
}