#include #include #define rep(i,n) for(int (i)=0;(i)<(n);(i)++) #define rep2(i,m,n) for(int (i)=(m);(i)<(n);(i)++) #define rep2ll(i,m,n) for(ll (i)=(m);(i)<(n);(i)++) #define ALL(obj) (obj).begin(), (obj).end() #define rALL(obj) (obj).rbegin(), (obj).rend() using namespace std; using ll = long long; using P = pair; using mint = atcoder::modint998244353; using VL = vector; using VVL = vector; using VVVL = vector; using VM = vector; using VVM = vector; using VVVM = vector; using VD = vector; using VVD = vector; using VVVD = vector; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b struct PersistentArray { struct Node { T data; Node *child[1 << LOG] = {}; Node() {} Node(const T &data) : data(data) {} }; Node *root; PersistentArray() : root(nullptr) {} T get(Node *t, int k) { if (k == 0) return t->data; return get(t->child[k & ((1 << LOG) - 1)], k >> LOG); } T get(const int &k) { return get(root, k); } pair mutable_get(Node *t, int k) { t = t ? new Node(*t) : new Node(); if (k == 0) return {t, &t->data}; auto p = mutable_get(t->child[k & ((1 << LOG) - 1)], k >> LOG); t->child[k & ((1 << LOG) - 1)] = p.first; return {t, p.second}; } T *mutable_get(const int &k) { auto ret = mutable_get(root, k); root = ret.first; return ret.second; } Node *build(Node *t, const T &data, int k) { if (!t) t = new Node(); if (k == 0) { t->data = data; return t; } auto p = build(t->child[k & ((1 << LOG) - 1)], data, k >> LOG); t->child[k & ((1 << LOG) - 1)] = p; return t; } void build(const vector &v) { root = nullptr; for (int i = 0; i < v.size(); i++) { root = build(root, v[i], i); } } }; /* * @brief Persistent-Union-Find(永続Union-Find) */ struct PersistentUnionFind { PersistentArray data; PersistentUnionFind() {} PersistentUnionFind(int sz) { data.build(vector(sz, -1)); } int find(int k) { int p = data.get(k); return p >= 0 ? find(p) : k; } int size(int k) { return (-data.get(find(k))); } bool unite(int x, int y) { x = find(x); y = find(y); if (x == y) return false; auto u = data.get(x); auto v = data.get(y); if (u < v) { auto a = data.mutable_get(x); *a += v; auto b = data.mutable_get(y); *b = x; } else { auto a = data.mutable_get(y); *a += u; auto b = data.mutable_get(x); *b = y; } return true; } }; int main(){ ios::sync_with_stdio(false); cin.tie(0); ll h, w; cin>>h>>w; VVL a(h, VL(w)); rep(x,h) rep(y,w) cin>>a[x][y], a[x][y]--; vector> b(h*w); rep(x,h) rep(y,w){ if(x){ b[max(a[x][y], a[x-1][y])].emplace_back(w*x+y, w*(x-1)+y); } if(y){ b[max(a[x][y], a[x][y-1])].emplace_back(w*x+y, w*x+(y-1)); } } vector uf(h*w+1); uf[0] = PersistentUnionFind(h*w); rep(i,h*w){ uf[i+1] = uf[i]; for(auto [u,v]:b[i]) uf[i+1].unite(u, v); } ll q; cin>>q; while(q--){ ll rs1, cs1, rt1, ct1; cin>>rs1>>cs1>>rt1>>ct1; rs1--,cs1--,rt1--,ct1--; ll u = w*rs1+cs1, v = w*rt1+ct1; ll ng = 0, ok = h*w+1; while(ng+1