// (◕ᴗ◕✿) // #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include #define rep(i, n) for (ll i = 0; i < (n); i++) #define srep(i, s, n) for (ll i = s; i < (n); i++) #define len(x) ((int)(x).size()) #define all(x) (x).begin(), (x).end() using namespace std; template using vc = vector; template using vv = vc>; template using vvv = vv>; using vi = vc;using vvi = vv; using vvvi = vv; using ll = long long;using vl = vc;using vvl = vv; using vvvl = vv; using ld = long double; using vld = vc; using vvld = vc; using vvvld = vc; using uint = unsigned int; using ull = unsigned long long; const ld pi = acos(-1.0); const int inf = 0x3f3f3f3f; const ll INF = 0x3f3f3f3f3f3f3f3f; // const ll mod = 1000000007; const ll mod = 998244353; inline bool inside(ll y, ll x, ll H, ll W) {return 0 <= (y) and (y) < (H) and 0 <= (x) and (x) < (W); } #define debug(var) do { cerr << #var << " :\n"; view(var); } while(0) templatevoid view(const T& e) {cerr << e;} templatevoid view(const pair& p) {cerr << "{" << p.first << ", " << p.second << "}";} templatevoid view(const vc& v) {for (const auto& e : v) {view(e);cerr << " ";} cerr << endl;} templatevoid view(const vv& vv) {for (const auto& v : vv) {view(v);} cerr << endl;} templatevoid view(const set& s) {for (const auto& e : s) {view(e);cerr << " ";} cerr << endl;} templatevoid view(const multiset& s) {for (const auto& e : s) {view(e);cerr << " ";} cerr << endl;} templatevoid view(const unordered_set& s) {for (const auto& e : s) {view(e);cerr << " ";} cerr << endl;} templatevoid view(const map& mp){for (const auto& e : mp) {view(e);cerr << " ";} cerr << endl;} void solve(){ int N, Q; cin >> N >> Q; vl A(N); rep(i, N) cin >> A[i]; vl B(N); rep(i, N) cin >> B[i]; vi L(Q), D(Q), R(Q), U(Q); rep(i, Q) cin >> L[i] >> D[i] >> R[i] >> U[i], L[i]--, D[i]--; vl ans(Q, 0); int x = N / 3, y = (N - x) / 2, z = N - y - x; // cout << x << ' ' << y << ' ' << z << endl; vvl cmsm(max({x, y, z}) + 1, vl(N + 1, 0)); { rep(i, x) rep(j, N){ cmsm[i + 1][j + 1] = max(A[i], B[j]) + cmsm[i + 1][j] + cmsm[i][j + 1] - cmsm[i][j]; } rep(i, Q){ int l = min(x, L[i]), r = min(x, R[i]); ans[i] += cmsm[r][U[i]] - cmsm[r][D[i]] - cmsm[l][U[i]] + cmsm[l][D[i]]; } } { rep(i, y) rep(j, N){ cmsm[i + 1][j + 1] = max(A[i + x], B[j]) + cmsm[i + 1][j] + cmsm[i][j + 1] - cmsm[i][j]; } rep(i, Q){ int l = min(y, max(x, L[i]) - x), r = min(y, max(x, R[i]) - x); ans[i] += cmsm[r][U[i]] - cmsm[r][D[i]] - cmsm[l][U[i]] + cmsm[l][D[i]]; } } { rep(i, z) rep(j, N){ cmsm[i + 1][j + 1] = max(A[i + y + x], B[j]) + cmsm[i + 1][j] + cmsm[i][j + 1] - cmsm[i][j]; } rep(i, Q){ int l = max(y + x, L[i]) - y - x, r = max(y + x, R[i]) - y - x; ans[i] += cmsm[r][U[i]] - cmsm[r][D[i]] - cmsm[l][U[i]] + cmsm[l][D[i]]; } } rep(i, Q) cout << ans[i] << endl; } int main(){ ios_base::sync_with_stdio(false); cin.tie(nullptr); int T = 1; // cin >> T; while (T--) solve(); }