#include<bits/stdc++.h>
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) ((m)*(x-1)+(y))
#define R(n) (rnd()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define fi first
#define se second
using namespace std;using ll=long long;using db=double;using lb=long db;using ui=unsigned;using ull=unsigned long long;using pii=pair<int,int>;using LL=__int128;
const int N=500+5,M=3e5+5,K=(1<<15)+5,mod=998244353,Mod=mod-1;const db eps=1e-9;const int INF=1e9+7;mt19937 rnd(time(0));
int n,q,A[N],B[N],w[N][N],ans[M];
struct ques{int x,y,l,r,id;};
int g1[N][N][2][2],g2[N][N][2][2];
void chx(int &x,int y){x=max(x,y);}
void Solve(int l,int r,int x,int y,vector<ques> q){
	if(l>r||x>y||q.empty()) return;
	if(r-l+1>=y-x+1){
		int m=l+r>>1;
		vector<ques> q1,q2,q3;
		for(ques i:q){
			if(i.r<m) q1.emplace_back(i);
			else if(i.l>m) q2.emplace_back(i);
			else q3.emplace_back(i);
		}
		Solve(l,m-1,x,y,q1);Solve(m+1,r,x,y,q2);
		for(int i=x;i<=y;i++){
			for(int o1=0;o1<2;o1++){
				for(int j=l;j<=m;j++) for(int h=x;h<=i;h++) Me(g1[j][h],-0x3f);
				g1[m][i][o1][1]=-A[m]*o1-B[i];
				for(int j=m;j>=l;j--)for(int h=i;h>=x;h--){
					for(int o3=0;o3<2;o3++) for(int o4=0;o4<2;o4++){
						if(o3&&o4) g1[j][h][1][1]+=w[j][h];
						if(!o3&&j^m) chx(g1[j][h][1][o4],g1[j][h][0][o4]-A[j]);
						if(!o4&&h^i) chx(g1[j][h][o3][1],g1[j][h][o3][0]-B[h]);
						if(j>l) chx(g1[j-1][h][0][o4],g1[j][h][o3][o4]);
						if(h>x) chx(g1[j][h-1][o3][0],g1[j][h][o3][o4]);
					}
				}
				for(int j=m;j<=r;j++) for(int h=i;h<=y;h++) Me(g2[j][h],-0x3f);
				g2[m][i][o1][1]=-A[m]*o1-B[i];
				for(int j=m;j<=r;j++) for(int h=i;h<=y;h++){
					for(int o3=0;o3<2;o3++) for(int o4=0;o4<2;o4++){
						if(o3&&o4) g2[j][h][1][1]+=w[j][h];
						if(!o3&&j^m) chx(g2[j][h][1][o4],g2[j][h][0][o4]-A[j]);
						if(!o4&&h^i) chx(g2[j][h][o3][1],g2[j][h][o3][0]-B[h]);
						if(j<r) chx(g2[j+1][h][0][o4],g2[j][h][o3][o4]);
						if(h<y) chx(g2[j][h+1][o3][0],g2[j][h][o3][o4]);
					}
				}
				for(ques j:q3) if(j.x<=i&&j.y>=i) {
					int m1=-INF,m2=-INF;
					for(int p1=0;p1<2;p1++) for(int p2=0;p2<2;p2++) m1=max(m1,g1[j.l][j.x][p1][p2]);
					for(int p1=0;p1<2;p1++) for(int p2=0;p2<2;p2++) m2=max(m2,g2[j.r][j.y][p1][p2]);
					ans[j.id]=max(ans[j.id],m1+m2+o1*A[m]+B[i]-o1*w[m][i]);
				}
			}
		}
	}else{
		int m=x+y>>1;
		vector<ques> q1,q2,q3;
		for(ques i:q){
			if(i.y<m) q1.emplace_back(i);
			else if(i.x>m) q2.emplace_back(i);
			else q3.emplace_back(i);
		}
		Solve(l,r,x,m-1,q1);Solve(l,r,m+1,y,q2);
		for(int i=l;i<=r;i++){
			for(int o2=0;o2<2;o2++){
				for(int j=l;j<=i;j++) for(int h=x;h<=m;h++) Me(g1[j][h],-0x3f);
				g1[i][m][1][o2]=-A[i]-B[m]*o2;
				for(int j=i;j>=l;j--)for(int h=m;h>=x;h--){
					for(int o3=0;o3<2;o3++) for(int o4=0;o4<2;o4++){
						if(o3&&o4) g1[j][h][1][1]+=w[j][h];
						if(!o3&&j^i) chx(g1[j][h][1][o4],g1[j][h][0][o4]-A[j]);
						if(!o4&&h^m) chx(g1[j][h][o3][1],g1[j][h][o3][0]-B[h]);
						if(j>l) chx(g1[j-1][h][0][o4],g1[j][h][o3][o4]);
						if(h>x) chx(g1[j][h-1][o3][0],g1[j][h][o3][o4]);
					}
				}
				for(int j=i;j<=r;j++) for(int h=m;h<=y;h++) Me(g2[j][h],-0x3f);
				g2[i][m][1][o2]=-A[i]-B[m]*o2;
				for(int j=i;j<=r;j++) for(int h=m;h<=y;h++){
					for(int o3=0;o3<2;o3++) for(int o4=0;o4<2;o4++){
						if(o3&&o4) g2[j][h][1][1]+=w[j][h];
						if(!o3&&j^i) chx(g2[j][h][1][o4],g2[j][h][0][o4]-A[j]);
						if(!o4&&h^m) chx(g2[j][h][o3][1],g2[j][h][o3][0]-B[h]);
						if(j<r) chx(g2[j+1][h][0][o4],g2[j][h][o3][o4]);
						if(h<y) chx(g2[j][h+1][o3][0],g2[j][h][o3][o4]);
					}
				}
				for(ques j:q3) if(j.l<=i&&j.r>=i) {
					int m1=-INF,m2=-INF;
					for(int p1=0;p1<2;p1++) for(int p2=0;p2<2;p2++) m1=max(m1,g1[j.l][j.x][p1][p2]);
					for(int p1=0;p1<2;p1++) for(int p2=0;p2<2;p2++) m2=max(m2,g2[j.r][j.y][p1][p2]);
					ans[j.id]=max(ans[j.id],m1+m2+A[i]+o2*B[m]-o2*w[i][m]);
				}
			}
		}
	}
}
void Solve(){
	int i,j;scanf("%d%d",&n,&q);
	for(i=1;i<=n;i++) scanf("%d",&A[i]);
	for(i=1;i<=n;i++) scanf("%d",&B[i]);
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++) scanf("%d",&w[i][j]);
	}
	vector<ques> Q;
	for(i=1;i<=q;i++){
		ques p;scanf("%d%d%d%d",&p.l,&p.r,&p.x,&p.y);p.id=i;
		Q.emplace_back(p);
	}
	Solve(1,n,1,n,Q);
	for(i=1;i<=q;i++) printf("%d\n",ans[i]);
}
int main(){
	int t;
	// scanf("%d",&t);
	t=1;
	while(t--) Solve();
	cerr<<clock()*1.0/CLOCKS_PER_SEC<<'\n';
}