#include <bits/stdc++.h>
using namespace std;
using i32 = int;
#define vi vector<i32>
#define int long long
const int maxn = 1e5+5;
const int B = 330;
int N,M,Q,tree[4*maxn],Max[maxn];


vector<int> build_teams(vi A1,vi B1,vi A2,vi B2,vi L1,vi R1,vi L2,vi R2) {

    Q=(int)L1.size();
    vector<int> res(Q);
    auto C = [&](int i,int j){
        return 1LL*(A1[i]+A2[j])*(B1[i]+B2[j]);
    };
    auto solve = [&]{
        N=(int)A1.size();
        M=(int)A2.size();

        function<void(int,int,int,int)> dnc = [&](int l,int r,int optl,int optr){
            int mid=(l+r)>>1,opt=-1,val=0;
            for(int i=optl;i<=optr;i++){
                if(C(mid,i)>val) opt=i,val=C(mid,i);
            }
            Max[mid]=val;
            if(l<mid) dnc(l,mid-1,opt,optr);
            if(mid<r) dnc(mid+1,r,optl,opt);
        };
        for(int l=0;l<M;l+=B){
            int r=min(l+B,M)-1;
            dnc(0,N-1,l,r);
            vector<int> tree(4*N);
            function<void(int,int,int)> build = [&](int l,int r,int id){
                if(l==r){
                    tree[id]=Max[l];
                    return;
                }
                int mid=(l+r)>>1;
                build(l,mid,id<<1);build(mid+1,r,id<<1|1);
                tree[id]=max(tree[id<<1],tree[id<<1|1]);
            };
            build(0,N-1,1);
            function<int(int,int,int,int,int)> query = [&](int l,int r,int id,int tl,int tr){
                if(tr<l || r<tl) return 0LL;
                if(tl<=l && r<=tr) return tree[id];
                int mid=(l+r)>>1;
                return max(query(l,mid,id<<1,tl,tr),query(mid+1,r,id<<1|1,tl,tr));
            };
            for(int i=0;i<Q;i++) if(L2[i]<=l && r<=R2[i]) res[i]=max(res[i],query(0,N-1,1,L1[i],R1[i]));
        }
    };
    solve();
    swap(A1,A2);
    swap(B1,B2);
    swap(L1,L2);
    swap(R1,R2);
    solve();

    int N=(int)A1.size();
    int M=(int)A2.size();
    for(int t=0;t<Q;t++){
        vector<int> a,b;
        for(int j=L1[t];j<=R1[t];j++){
            if(j%B==0 && j+B-1<=R1[t]) j+=B-1;
            else a.push_back(j);
        }
        for(int j=L2[t];j<=R2[t];j++){
            if(j%B==0 && j+B-1<=R2[t]) j+=B-1;
            else b.push_back(j);
        }
        function<void(int,int,int,int)> dnc = [&](int l,int r,int optl,int optr){
            int mid=(l+r)>>1,opt=-1,val=0;
            for(int i=optl;i<=optr;i++){
                if(C(a[mid],b[i])>val) opt=i,val=C(a[mid],b[i]);
            }
            res[t]=max(res[t],val);
            if(l<mid) dnc(l,mid-1,opt,optr);
            if(mid<r) dnc(mid+1,r,optl,opt);
        };
        dnc(0,(int)a.size()-1,0,(int)b.size()-1);
    }
	return res;
}

#undef int