#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include using namespace std; using int128 = __int128_t; using int64 = long long; using int32 = int; using uint128 = __uint128_t; using uint64 = unsigned long long; using uint32 = unsigned int; #define ALL(obj) (obj).begin(),(obj).end() template using priority_queue_reverse = priority_queue,greater>; constexpr int64 MOD = 1'000'000'000LL + 7; //' constexpr int64 MOD2 = 998244353; constexpr int64 HIGHINF = 1'000'000'000'000'000'000LL; constexpr int64 LOWINF = 1'000'000'000'000'000LL; //' constexpr long double PI = 3.1415926535897932384626433L; template vector multivector(size_t N,T init){return vector(N,init);} template auto multivector(size_t N,T... t){return vector(N,multivector(t...));} template void corner(bool flg, T hoge) {if (flg) {cout << hoge << endl; exit(0);}} template ostream &operator<<(ostream &o, const map&obj) {o << "{"; for (auto &x : obj) o << " {" << x.first << " : " << x.second << "}" << ","; o << " }"; return o;} template ostream &operator<<(ostream &o, const set&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const multiset&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const vector&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const deque&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const pair&obj) {o << "{" << obj.first << ", " << obj.second << "}"; return o;} void print(void) {cout << endl;} template void print(Head&& head) {cout << head;print();} template void print(Head&& head, Tail&&... tail) {cout << head << " ";print(forward(tail)...);} template void chmax(T& a, const T b){a=max(a,b);} template void chmin(T& a, const T b){a=min(a,b);} vector split(const string &str, const char delemiter) {vector res;stringstream ss(str);string buffer; while( getline(ss, buffer, delemiter) ) res.push_back(buffer); return res;} inline constexpr int msb(int x) {return x?31-__builtin_clz(x):-1;} inline constexpr int64 ceil_div(const int64 a,const int64 b) {return (a+(b-1))/b;}// return ceil(a/b) void YN(bool flg) {cout << (flg ? "YES" : "NO") << endl;} void Yn(bool flg) {cout << (flg ? "Yes" : "No") << endl;} void yn(bool flg) {cout << (flg ? "yes" : "no") << endl;} /* * @title BitVector * @docs md/segment/BitVector.md */ class BitVector{ inline static constexpr size_t BIT_BLOCK_SIZE = 5; inline static constexpr size_t BIT_BLOCK_NUM = 1< vec, acc; bool is_builded = false; public: BitVector(uint32 N) { uint32 tmp = (N + BIT_BLOCK_NUM-1) / BIT_BLOCK_NUM; vec.assign(tmp, 0); acc.assign(tmp, 0); } void build() { for (size_t i = 0,sum = 0; i < acc.size(); ++i) acc[i] = (sum += popcount(vec[i]) ); is_builded = true; } //[0,r) count of bit uint32 rank(uint32 r, const bool bit) const { assert(is_builded); uint32 bit_r = r >> BIT_BLOCK_SIZE; uint32 sum = (bit_r ? acc[bit_r - 1] : 0) + popcount(vec[bit_r] & ((1U << (r & (BIT_BLOCK_NUM-1))) - 1)); return bit ? sum : r - sum; } //[l,l+1) = bit void update(uint32 l, const bool bit) { uint32 bit_l = l >> BIT_BLOCK_SIZE; if (bit) vec[bit_l] |= 1U << (l & (BIT_BLOCK_NUM-1)); else vec[bit_l] &= ~(1U << (l & (BIT_BLOCK_NUM-1))); } //[l,l+1) bool operator[](uint32 l) const { assert(is_builded); return ((vec[l >> BIT_BLOCK_SIZE] >> (l & (BIT_BLOCK_NUM-1))) & 1); } }; /* * @title WaveletMatrix * @docs md/segment/WaveletMatrix.md */ template class WaveletMatrix{ size_t length; size_t depth; vector multi_bit_vector; vector sum_bit_off; vector vec; vector zarts(const vector& ar) { vector ord(ar.size()),compressed(ar.size()); iota(ord.begin(),ord.end(),0); sort(ord.begin(),ord.end(),[&](size_t l, size_t r){return ar[l] &arg_vec): vec(arg_vec) { length = vec.size(); auto compressed = zarts(vec); sort(vec.begin(),vec.end()); vec.erase(unique(vec.begin(),vec.end()),vec.end()); uint32 maxi = *max_element(compressed.begin(),compressed.end()) + 1; for(depth=0; (1< prev = compressed, next = prev; for(uint32 j = 0; j < depth; ++j,swap(prev,next)) { uint32 bit = 1UL << (depth - j - 1); for(uint32 i = 0; i < length; ++i) sum_bit_off[j] += !(prev[i] & bit); uint32 idx_bit_off = 0, idx_bit_on = sum_bit_off[j]; for(uint32 i = 0; i < length; ++i) { if (prev[i] & bit) multi_bit_vector[j].update(i,1), next[idx_bit_on++] = prev[i]; else next[idx_bit_off++] = prev[i]; } multi_bit_vector[j].build(); } } //[l,l+1) element T operator[](uint32 l) const { uint32 val = 0; for (uint32 j = 0; j < depth; ++j) { const bool bit = multi_bit_vector[j][l]; val = ((val << 1) | bit); l = multi_bit_vector[j].rank(l, bit); if(bit) l += sum_bit_off[j]; } return vec[val]; } //[l,r) range count of more and less than c (c' rank_all(T c, uint32 l, uint32 r) const // { // if (c > vec.back()) return make_tuple(r - l, 0U, 0U); // array more_and_less = {0,0}; // for (uint32 j = 0; j < depth; ++j) { // const bool bit = (c >> (depth - j - 1)) & 1; // more_and_less[bit] += r - l; // l = multi_bit_vector[j].rank(l, bit); // r = multi_bit_vector[j].rank(r, bit); // if (bit) l += sum_bit_off[j], r += sum_bit_off[j]; // more_and_less[bit] -= r - l; // } // return make_tuple(more_and_less[1], r - l, more_and_less[0]); // } //[l,r) range k (0-indexed) th smallest number T quantile(uint32 l, uint32 r, uint32 k) const { uint32 val = 0; for (uint32 j = 0; j < depth; ++j) { uint32 cnt_bit_off = multi_bit_vector[j].rank(r, 0) - multi_bit_vector[j].rank(l, 0); const bool bit = (k >= cnt_bit_off); val = ((val << 1) | bit); l = multi_bit_vector[j].rank(l, bit); r = multi_bit_vector[j].rank(r, bit); if (bit) l += sum_bit_off[j], r += sum_bit_off[j], k -= cnt_bit_off; } return vec[val]; } }; /** * @url * @est */ int main() { cin.tie(0);ios::sync_with_stdio(false); int N; cin >> N; vector A(N); for(int i=0;i> A[i]; WaveletMatrix wm(A); int Q; cin >> Q; while(Q--) { int l,r; int64 x; cin >> l >> r >> x; l--; int k = r-l; int64 ans = LOWINF; chmin(ans,abs(x-A[l])); chmin(ans,abs(x-A[r-1])); if(k<=2) { cout << ans << endl; continue; } // print(l,r,x,ans); int64 mini = wm.quantile(l,r,0); int64 maxi = wm.quantile(l,r,k-1); chmin(ans,abs(x-mini)); chmin(ans,abs(x-maxi)); if(x <= mini || maxi <= x) { cout << ans << endl; continue; } // print(l,r,x,ans); int ok=0,ng=k,md; while(ng-ok>1) { md = (ok+ng)/2; int64 y = wm.quantile(l,r,md); (y