// #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") // #define INTERACTIVE // #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") // #define INTERACTIVE #include using namespace std; namespace templates { // type using ll = long long; using ull = unsigned long long; using Pii = pair; using Pil = pair; using Pli = pair; using Pll = pair; template using pq = priority_queue; template using qp = priority_queue, greater>; // clang-format off #define vec(T, A, ...) vector A(__VA_ARGS__); #define vvec(T, A, h, ...) vector> A(h, vector(__VA_ARGS__)); #define vvvec(T, A, h1, h2, ...) vector>> A(h1, vector>(h2, vector(__VA_ARGS__))); // clang-format on // for loop #define fori1(a) for (ll _ = 0; _ < (a); _++) #define fori2(i, a) for (ll i = 0; i < (a); i++) #define fori3(i, a, b) for (ll i = (a); i < (b); i++) #define fori4(i, a, b, c) for (ll i = (a); ((c) > 0 || i > (b)) && ((c) < 0 || i < (b)); i += (c)) #define overload4(a, b, c, d, e, ...) e #define fori(...) overload4(__VA_ARGS__, fori4, fori3, fori2, fori1)(__VA_ARGS__) // declare and input // clang-format off #define INT(...) int __VA_ARGS__; inp(__VA_ARGS__); #define LL(...) ll __VA_ARGS__; inp(__VA_ARGS__); #define STRING(...) string __VA_ARGS__; inp(__VA_ARGS__); #define CHAR(...) char __VA_ARGS__; inp(__VA_ARGS__); #define DOUBLE(...) double __VA_ARGS__; STRING(str___); __VA_ARGS__ = stod(str___); #define VEC(T, A, n) vector A(n); inp(A); #define VVEC(T, A, n, m) vector> A(n, vector(m)); inp(A); // clang-format on // const value const ll MOD1 = 1000000007; const ll MOD9 = 998244353; const double PI = acos(-1); // other macro #if !defined(RIN__LOCAL) && !defined(INTERACTIVE) #define endl "\n" #endif #define spa ' ' #define len(A) ll(A.size()) #define all(A) begin(A), end(A) // function vector stoc(string &S) { int n = S.size(); vector ret(n); for (int i = 0; i < n; i++) ret[i] = S[i]; return ret; } string ctos(vector &S) { int n = S.size(); string ret = ""; for (int i = 0; i < n; i++) ret += S[i]; return ret; } template auto min(const T &a) { return *min_element(all(a)); } template auto max(const T &a) { return *max_element(all(a)); } template auto clamp(T &a, const S &l, const S &r) { return (a > r ? r : a < l ? l : a); } template inline bool chmax(T &a, const S &b) { return (a < b ? a = b, 1 : 0); } template inline bool chmin(T &a, const S &b) { return (a > b ? a = b, 1 : 0); } template inline bool chclamp(T &a, const S &l, const S &r) { auto b = clamp(a, l, r); return (a != b ? a = b, 1 : 0); } template T sum(vector &A) { T tot = 0; for (auto a : A) tot += a; return tot; } template vector compression(vector X) { sort(all(X)); X.erase(unique(all(X)), X.end()); return X; } // input and output namespace io { // __int128_t std::ostream &operator<<(std::ostream &dest, __int128_t value) { std::ostream::sentry s(dest); if (s) { __uint128_t tmp = value < 0 ? -value : value; char buffer[128]; char *d = std::end(buffer); do { --d; *d = "0123456789"[tmp % 10]; tmp /= 10; } while (tmp != 0); if (value < 0) { --d; *d = '-'; } int len = std::end(buffer) - d; if (dest.rdbuf()->sputn(d, len) != len) { dest.setstate(std::ios_base::badbit); } } return dest; } // vector template istream &operator>>(istream &is, vector &A) { for (auto &a : A) is >> a; return is; } template ostream &operator<<(ostream &os, vector &A) { for (size_t i = 0; i < A.size(); i++) { os << A[i]; if (i != A.size() - 1) os << ' '; } return os; } // vector> template istream &operator>>(istream &is, vector> &A) { for (auto &a : A) is >> a; return is; } template ostream &operator<<(ostream &os, vector> &A) { for (size_t i = 0; i < A.size(); i++) { os << A[i]; if (i != A.size() - 1) os << endl; } return os; } // pair template istream &operator>>(istream &is, pair &A) { is >> A.first >> A.second; return is; } template ostream &operator<<(ostream &os, pair &A) { os << A.first << ' ' << A.second; return os; } // vector> template istream &operator>>(istream &is, vector> &A) { for (size_t i = 0; i < A.size(); i++) { is >> A[i]; } return is; } template ostream &operator<<(ostream &os, vector> &A) { for (size_t i = 0; i < A.size(); i++) { os << A[i]; if (i != A.size() - 1) os << endl; } return os; } // tuple template struct TuplePrint { static ostream &print(ostream &os, const T &t) { TuplePrint::print(os, t); os << ' ' << get(t); return os; } }; template struct TuplePrint { static ostream &print(ostream &os, const T &t) { os << get<0>(t); return os; } }; template ostream &operator<<(ostream &os, const tuple &t) { TuplePrint::print(os, t); return os; } // io functions void FLUSH() { cout << flush; } void print() { cout << endl; } template void print(Head &&head, Tail &&...tail) { cout << head; if (sizeof...(Tail)) cout << spa; print(std::forward(tail)...); } template void prisep(vector &A, S sep) { int n = A.size(); for (int i = 0; i < n; i++) { cout << A[i]; if (i != n - 1) cout << sep; } cout << endl; } template void priend(T A, S end) { cout << A << end; } template void prispa(T A) { priend(A, spa); } template bool printif(bool f, T A, S B) { if (f) print(A); else print(B); return f; } template void inp(T &...a) { (cin >> ... >> a); } } // namespace io using namespace io; // read graph vector> read_edges(int n, int m, bool direct = false, int indexed = 1) { vector> edges(n, vector()); for (int i = 0; i < m; i++) { INT(u, v); u -= indexed; v -= indexed; edges[u].push_back(v); if (!direct) edges[v].push_back(u); } return edges; } vector> read_tree(int n, int indexed = 1) { return read_edges(n, n - 1, false, indexed); } template vector>> read_wedges(int n, int m, bool direct = false, int indexed = 1) { vector>> edges(n, vector>()); for (int i = 0; i < m; i++) { INT(u, v); T w; inp(w); u -= indexed; v -= indexed; edges[u].push_back({v, w}); if (!direct) edges[v].push_back({u, w}); } return edges; } template vector>> read_wtree(int n, int indexed = 1) { return read_wedges(n, n - 1, false, indexed); } // yes / no namespace yesno { // yes inline bool yes(bool f = true) { cout << (f ? "yes" : "no") << endl; return f; } inline bool Yes(bool f = true) { cout << (f ? "Yes" : "No") << endl; return f; } inline bool YES(bool f = true) { cout << (f ? "YES" : "NO") << endl; return f; } // no inline bool no(bool f = true) { cout << (!f ? "yes" : "no") << endl; return f; } inline bool No(bool f = true) { cout << (!f ? "Yes" : "No") << endl; return f; } inline bool NO(bool f = true) { cout << (!f ? "YES" : "NO") << endl; return f; } // possible inline bool possible(bool f = true) { cout << (f ? "possible" : "impossible") << endl; return f; } inline bool Possible(bool f = true) { cout << (f ? "Possible" : "Impossible") << endl; return f; } inline bool POSSIBLE(bool f = true) { cout << (f ? "POSSIBLE" : "IMPOSSIBLE") << endl; return f; } // impossible inline bool impossible(bool f = true) { cout << (!f ? "possible" : "impossible") << endl; return f; } inline bool Impossible(bool f = true) { cout << (!f ? "Possible" : "Impossible") << endl; return f; } inline bool IMPOSSIBLE(bool f = true) { cout << (!f ? "POSSIBLE" : "IMPOSSIBLE") << endl; return f; } // Alice Bob inline bool Alice(bool f = true) { cout << (f ? "Alice" : "Bob") << endl; return f; } inline bool Bob(bool f = true) { cout << (f ? "Bob" : "Alice") << endl; return f; } // Takahashi Aoki inline bool Takahashi(bool f = true) { cout << (f ? "Takahashi" : "Aoki") << endl; return f; } inline bool Aoki(bool f = true) { cout << (f ? "Aoki" : "Takahashi") << endl; return f; } } // namespace yesno using namespace yesno; } // namespace templates using namespace templates; template std::vector> get_segment_range(T l, T r) { std::vector> left; std::vector> right; int i = 0; while (l < r) { if ((l >> i) & 1) { left.emplace_back(l, l + (1 << i)); l += 1 << i; } if ((r >> i) & 1) { r -= 1 << i; right.emplace_back(r, r + (1 << i)); } i++; } std::reverse(right.begin(), right.end()); left.insert(left.end(), right.begin(), right.end()); return left; } struct Query { int p; ll x; int k; }; void solve() { INT(n, Q); VEC(ll, A, n); VEC(ll, B, n); vec(Query, queries, Q); vvec(Pii, T, n); fori(i, n) { T[i].push_back({0, A[i]}); } fori(i, Q) { LL(p, x, k); p--; k -= 2; queries[i] = {int(p), x, int(k)}; T[p].push_back({i, x}); } vvvec(Pii, px, 1, 1); vvvec(Pii, ki, 1, 1); while (px.back().size() < Q) { px.push_back(px.back()); px.back().resize(2 * px.back().size()); ki.push_back(ki.back()); ki.back().resize(2 * ki.back().size()); } reverse(all(px)); reverse(all(ki)); fori(i, n) { T[i].push_back({Q, -1}); fori(j, len(T[i]) - 1) { auto seg = get_segment_range(T[i][j].first, T[i][j + 1].first); for (auto [l, r] : seg) { int p = __builtin_ctz(l ^ r); int x = l >> p; px[p][x].push_back({i, T[i][j].second}); } } } fori(i, Q) { ki[0][i].push_back({queries[i].k, i}); } fori(i, 1, len(ki)) { fori(j, len(ki[i - 1])) { for (auto &p : ki[i - 1][j]) { ki[i][j >> 1].push_back(p); } } } const ll inf = 1LL << 60; vec(ll, ans, Q, inf); auto sol = [&](vector &px, vector &K) -> vector { int n = len(px); int m = len(K); auto get = [&](int i, int j) -> ll { int bi = K[i] - px[j].first; if (0 <= bi and bi < len(B)) { return B[bi] + px[j].second; } else if (bi < 0) { return inf - bi; } else { return inf + bi; } }; vec(ll, res, m, inf); auto rec = [&](auto &&self, int l, int r, int li, int ri) -> void { if (l == r) return; int mid = (l + r) / 2; int bi = -1; ll min_ = inf; fori(i, li, ri) { auto res = get(mid, i); if (chmin(min_, res)) { bi = i; } } res[mid] = min_; if (bi == -1) { self(self, l, mid, li, ri); self(self, mid + 1, r, li, ri); } else { self(self, l, mid, li, bi + 1); self(self, mid + 1, r, bi, ri); } }; rec(rec, 0, m, 0, n); return res; }; reverse(all(px)); reverse(all(ki)); fori(i, len(px)) { fori(j, len(px[i])) { if (px[i][j].empty()) continue; sort(all(px[i][j])); vec(int, K, 0); for (auto [k, _] : ki[i][j]) { K.push_back(k); } K = compression(K); auto res = sol(px[i][j], K); sort(all(ki[i][j])); int p = 0; for (auto [k, idx] : ki[i][j]) { while (K[p] != k) p++; chmin(ans[idx], res[p]); } } } prisep(ans, "\n"); } int main() { #ifndef INTERACTIVE std::cin.tie(0)->sync_with_stdio(0); #endif // std::cout << std::fixed << std::setprecision(12); int t; t = 1; // std::cin >> t; while (t--) solve(); return 0; } // // #pragma GCC target("avx2") // // #pragma GCC optimize("O3") // // #pragma GCC optimize("unroll-loops") // // #define INTERACTIVE // // // #pragma GCC target("avx2") // // #pragma GCC optimize("O3") // // #pragma GCC optimize("unroll-loops") // // #define INTERACTIVE // // #include "kyopro-cpp/template.hpp" // // #include "misc/get_segment_range.hpp" // // struct Query { // int p; // ll x; // int k; // }; // // void solve() { // INT(n, Q); // VEC(ll, A, n); // VEC(ll, B, n); // // vec(Query, queries, Q); // vvec(Pii, T, n); // fori(i, n) { // T[i].push_back({0, A[i]}); // } // // fori(i, Q) { // LL(p, x, k); // p--; // k -= 2; // queries[i] = {int(p), x, int(k)}; // T[p].push_back({i, x}); // } // // vvvec(Pii, px, 1, 1); // vvvec(Pii, ki, 1, 1); // while (px.back().size() < Q) { // px.push_back(px.back()); // px.back().resize(2 * px.back().size()); // ki.push_back(ki.back()); // ki.back().resize(2 * ki.back().size()); // } // // reverse(all(px)); // reverse(all(ki)); // // fori(i, n) { // T[i].push_back({Q, -1}); // fori(j, len(T[i]) - 1) { // auto seg = get_segment_range(T[i][j].first, T[i][j + 1].first); // for (auto [l, r] : seg) { // int p = __builtin_ctz(l ^ r); // int x = l >> p; // px[p][x].push_back({i, T[i][j].second}); // } // } // } // // fori(i, Q) { // ki[0][i].push_back({queries[i].k, i}); // } // fori(i, 1, len(ki)) { // fori(j, len(ki[i - 1])) { // for (auto &p : ki[i - 1][j]) { // ki[i][j >> 1].push_back(p); // } // } // } // // const ll inf = 1LL << 60; // vec(ll, ans, Q, inf); // // auto sol = [&](vector &px, vector &K) -> vector { // int n = len(px); // int m = len(K); // // auto get = [&](int i, int j) -> ll { // int bi = K[i] - px[j].first; // if (0 <= bi and bi < len(B)) { // return B[bi] + px[j].second; // } else if (bi < 0) { // return inf - bi; // } else { // return inf + bi; // } // }; // // vec(ll, res, m, inf); // auto rec = [&](auto &&self, int l, int r, int li, int ri) -> void { // if (l == r) return; // int mid = (l + r) / 2; // int bi = -1; // ll min_ = inf; // fori(i, li, ri) { // auto res = get(mid, i); // if (chmin(min_, res)) { // bi = i; // } // } // // res[mid] = min_; // if (bi == -1) { // self(self, l, mid, li, ri); // self(self, mid + 1, r, li, ri); // } else { // self(self, l, mid, li, bi + 1); // self(self, mid + 1, r, bi, ri); // } // }; // rec(rec, 0, m, 0, n); // // return res; // }; // // reverse(all(px)); // reverse(all(ki)); // // fori(i, len(px)) { // fori(j, len(px[i])) { // if (px[i][j].empty()) continue; // sort(all(px[i][j])); // vec(int, K, 0); // for (auto [k, _] : ki[i][j]) { // K.push_back(k); // } // K = compression(K); // auto res = sol(px[i][j], K); // sort(all(ki[i][j])); // int p = 0; // for (auto [k, idx] : ki[i][j]) { // while (K[p] != k) p++; // chmin(ans[idx], res[p]); // } // } // } // prisep(ans, "\n"); // } // // int main() { // #ifndef INTERACTIVE // std::cin.tie(0)->sync_with_stdio(0); // #endif // // std::cout << std::fixed << std::setprecision(12); // int t; // t = 1; // // std::cin >> t; // while (t--) solve(); // return 0; // }