// #pragma GCC target("popcnt,lzcnt,abm,bmi,bmi2") #pragma GCC optimize("Ofast") /************* This code requires C++17. ***************/ #include using namespace std; /* macro helpers */ #define __NARGS(...) std::tuple_size::value #define __DECOMPOSE_S(a, x) auto x = a; #define __DECOMPOSE_N(a, ...) auto [__VA_ARGS__] = a; constexpr void __() {} #define __AS_PROCEDURE(...) __(); __VA_ARGS__; __() #define __as_typeof(container) remove_reference::type /* type aliases */ #if LONG_LONG_MAX != INT64_MAX using ll = int64_t; using ull = uint64_t; #else using ll = long long; using ull = unsigned long long; #endif using int128 = __int128_t; using uint128 = __uint128_t; using ld = long double; using pii = pair; using pil = pair; using pid = pair; using pli = pair; using pll = pair; using pld = pair; using pdi = pair; using pdl = pair; using pdd = pair; using tiii = tuple; using tiil = tuple; using tiid = tuple; using tili = tuple; using till = tuple; using tild = tuple; using tidi = tuple; using tidl = tuple; using tidd = tuple; using tlii = tuple; using tlil = tuple; using tlid = tuple; using tlli = tuple; using tlll = tuple; using tlld = tuple; using tldi = tuple; using tldl = tuple; using tldd = tuple; using tdii = tuple; using tdil = tuple; using tdid = tuple; using tdli = tuple; using tdll = tuple; using tdld = tuple; using tddi = tuple; using tddl = tuple; using tddd = tuple; template using max_heap = priority_queue; template using min_heap = priority_queue, greater<>>; template using oi = ostream_iterator; template using ii = istream_iterator; /* constants */ constexpr int INF = 0x3f3f3f3f; constexpr ll INFLL = 0x3f3f3f3f3f3f3f3fLL; constexpr ll MDL = 1e9 + 7; constexpr ll PRIME = 998'244'353; constexpr ll MDL1 = 8784491; constexpr ll MDL2 = PRIME; constexpr int128 INT128_MAX = numeric_limits::max(); constexpr uint128 UINT128_MAX = numeric_limits::max(); constexpr int128 INT128_MIN = numeric_limits::min(); constexpr uint128 UINT128_MIN = numeric_limits::min(); /* random */ mt19937_64 rd(chrono::duration_cast(chrono::system_clock::now().time_since_epoch()).count()); /* bit-wise operations */ #define lowbit(x) ((x) & -(x)) #define popcount(x) (__builtin_popcountll(ll(x))) #define parity(x) (__builtin_parityll(ll(x))) #define msp(x) (63LL - __builtin_clzll(ll(x))) #define lsp(x) (__builtin_ctzll(ll(x))) /* arithmetic operations */ #define mod(x, y) ((((x) % (y)) + (y)) % (y)) /* fast pairs */ #define upair ull #define umake(x, y) (ull(x) << 32 | (ull(y) & ((1ULL << 32) - 1))) #define u1(p) ((p) >> 32) #define u2(p) ((p) & ((1ULL << 32) - 1)) #define ult std::less #define ugt std::greater #define ipair ull #define imake(x, y) (umake(x, y)) #define i1(p) (int(u1(ll(p)))) #define i2(p) (ll(u2(p) << 32) >> 32) struct ilt { bool operator()(const ipair& a, const ipair& b) const { if (i1(a) == i1(b)) return i2(a) < i2(b); else return i1(a) < i1(b); } }; struct igt { bool operator()(const ipair& a, const ipair& b) const { if (i1(a) == i1(b)) return i2(a) > i2(b); else return i1(a) > i1(b); } }; /* conditions */ #define loop while (1) #define if_or(var, val) if (!(var == val)) var = val; else #define continue_or(var, val) __AS_PROCEDURE(if (var == val) continue; var = val;) #define break_or(var, val) __AS_PROCEDURE(if (var == val) break; var = val;) /* hash */ struct safe_hash { // https://codeforces.com/blog/entry/62393 static uint64_t splitmix64(uint64_t x) { // http://xorshift.di.unimi.it/splitmix64.c x += 0x9e3779b97f4a7c15; x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9; x = (x ^ (x >> 27)) * 0x94d049bb133111eb; return x ^ (x >> 31); } size_t operator()(uint64_t x) const { static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count(); return splitmix64(x + FIXED_RANDOM); } }; struct pair_hash { template size_t operator()(const pair& a) const { auto hash1 = safe_hash()(a.first); auto hash2 = safe_hash()(a.second); if (hash1 != hash2) { return hash1 ^ hash2; } return hash1; } }; uniform_int_distribution dist(PRIME); const size_t __array_hash_b = 31, __array_hash_mdl1 = dist(rd), __array_hash_mdl2 = dist(rd); struct array_hash { template size_t operator()(const Sequence& arr) const { size_t pw1 = 1, pw2 = 1; size_t res1 = 0, res2 = 0; for (auto&& x : arr) { res1 = (res1 + x * pw1) % __array_hash_mdl1; res2 = (res2 + x * pw2) % __array_hash_mdl2; pw1 = (pw1 * __array_hash_b) % __array_hash_mdl1; pw2 = (pw2 * __array_hash_b) % __array_hash_mdl2; } return res1 + res2; } }; /* build data structures */ #define faster(um) __AS_PROCEDURE((um).reserve(1024); (um).max_load_factor(0.25);) #define unordered_counter(from, to) __AS_PROCEDURE(unordered_map<__as_typeof(from), size_t, safe_hash> to; for (auto&& x : from) ++to[x];) #define counter(from, to, cmp) __AS_PROCEDURE(map<__as_typeof(from), size_t, cmp> to; for (auto&& x : from) ++to[x];) #define pa(a) __AS_PROCEDURE(__typeof(a) pa; pa.push_back({}); for (auto&&x : a) pa.push_back(pa.back() + x);) #define sa(a) __AS_PROCEDURE(__typeof(a) sa(a.size() + 1); {int n = a.size(); for (int i = n - 1; i >= 0; --i) sa[i] = sa[i + 1] + a[i];};) #define adj(ch, n) __AS_PROCEDURE(vector> ch((n) + 1);) #define edge(ch, u, v) __AS_PROCEDURE(ch[u].push_back(v), ch[v].push_back(u);) #define edgew(ch, u, v, ...) __AS_PROCEDURE(ch[u].emplace_back(v, __VA_ARGS__), ch[v].emplace_back(u, __VA_ARGS__);) #define Edge(ch, u, v) __AS_PROCEDURE(ch[u].push_back(v);) #define Edgew(ch, u, v, ...) __AS_PROCEDURE(ch[u].emplace_back(v, __VA_ARGS__);) template pair> discretize(Iterator __first, Iterator __last) { set st(__first, __last); size_t N = 0; map mp; for (auto&& x : st) mp[x] = ++N; return {N, mp}; } template pair> unordered_discretize(Iterator __first, Iterator __last) { set st(__first, __last); size_t N = 0; unordered_map mp; for (auto&& x : st) mp[x] = ++N; return {N, mp}; } /* io */ #define untie __AS_PROCEDURE(ios_base::sync_with_stdio(0), cin.tie(NULL)) template istream& operator>>(istream& in, pair& p) { return in >> p.first >> p.second; } template ostream& operator<<(ostream& out, const pair& p) { out << "{" << p.first << ", " << p.second << "}"; return out; } template void print_tuple_impl(std::basic_ostream& os, const Tuple& t, std::index_sequence) { using swallow = int[]; // guaranties left to right order (void)swallow { 0, (void(os << (Index == 0 ? "" : ", ") << std::get(t)), 0)... }; } template decltype(auto) operator<<(std::basic_ostream& os, const std::tuple& t) { os << "{"; print_tuple_impl(os, t, std::index_sequence_for{}); return os << "}"; } template ostream& operator<<(ostream& out, const vector& vec) { for (auto&& i : vec) out << i << ' '; return out; } std::ostream& operator<<(std::ostream& dest, const int128& value) { // https://stackoverflow.com/a/25115163/23881100 std::ostream::sentry s( dest ); if ( s ) { uint128 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; } template void __read(T& x) { cin >> x; } template void __read(T& x, U&... args) { cin >> x; __read(args...); } #define read(type, ...) __AS_PROCEDURE(type __VA_ARGS__; __read(__VA_ARGS__);) #define readvec(type, a, n) __AS_PROCEDURE(vector a(n); for (auto& x : a) cin >> x;) #define readvec1(type, a, n) __AS_PROCEDURE(vector a((n) + 1); copy_n(ii(cin), (n), a.begin() + 1);) #define putvec(a) __AS_PROCEDURE(copy(a.begin(), a.end(), oi<__as_typeof(a)::value_type>(cout, " ")); cout << endl;) #define putvec1(a) __AS_PROCEDURE(copy(a.begin() + 1, a.end(), oi<__as_typeof(a)::value_type>(cout, " ")); cout << endl;) #define putvec_eol(a) __AS_PROCEDURE(copy(a.begin(), a.end(), oi<__as_typeof(a)::value_type>(cout, "\n"));) #define putvec1_eol(a) __AS_PROCEDURE(copy(a.begin() + 1, a.end(), oi<__as_typeof(a)::value_type>(cout, "\n"));) #define debug(x) __AS_PROCEDURE(cerr << #x" = " << (x) << endl;) #define debugvec(a) __AS_PROCEDURE(cerr << #a" = "; for (auto&& x : a) cerr << x << ' '; cerr << endl;) #define deb(...) debug(make_tuple(__VA_ARGS__)) /* pops */ #define poptop(q, ...) __AS_PROCEDURE(auto [__VA_ARGS__] = q.top(); q.pop();) #define popback(q, ...) __AS_PROCEDURE(auto [__VA_ARGS__] = q.back(); q.pop_back();) #define popfront(q, ...) __AS_PROCEDURE(auto [__VA_ARGS__] = q.front();q.pop_front();) /* math */ template return_t qpow(ll b, ll p) { if (b == 0 and p != 0) return 0; if (p == 0) return 1; return_t half = qpow(b, p / 2); if (p % 2 == 1) return half * half * b; else return half * half; } #define comb(n, k) ((n) < 0 or (k) < 0 or (n) < (k) ? 0 : fact[n] / fact[k] / fact[(n) - (k)]) #define fastcomb(n, k) ((n) < 0 or (k) < 0 or (n) < (k) ? 0 : fact[n] * factrev[k] * factrev[(n) - (k)]) constexpr inline int lg2(ll x) { return x == 0 ? -1 : sizeof(ll) * 8 - 1 - __builtin_clzll(x); } void __exgcd(ll a, ll b, ll& x, ll& y) { if (b == 0) { x = 1, y = 0; return; } __exgcd(b, a % b, y, x); y -= a / b * x; } ll inverse(ll a, ll b) { ll x, y; __exgcd(a, b, x, y); return mod(x, b); } vector> decompose(ll x) { // return (factor, count, factor ** count) vector> res; for (int i = 2; i * i <= x; i++) { if (x % i == 0) { int cnt = 0; ll pw = 1; while (x % i == 0) ++cnt, x /= i, pw *= i; res.emplace_back(i, cnt, pw); } } if (x != 1) { res.emplace_back(x, 1, x); } return res; } vector decompose_prime(int N) { // return (factor, count) vector result; for (int i = 2; i * i <= N; i++) { if (N % i == 0) { int cnt = 0; while (N % i == 0) N /= i, ++cnt; result.emplace_back(i, cnt); } } if (N != 1) { result.emplace_back(N, 1); } return result; } /* string algorithms */ vector calc_next(string t) { // pi function of t int n = (int)t.length(); vector pi(n); for (int i = 1; i < n; i++) { int j = pi[i - 1]; while (j > 0 && t[i] != t[j]) j = pi[j - 1]; if (t[i] == t[j]) j++; pi[i] = j; } return pi; } vector calc_z(string t) { // z function of t int m = t.length(); vector z; z.push_back(m); pair prev = {1, -1}; for (int i = 1; i < m; ++i) { if (z[i - prev.first] + i <= prev.second) { z.push_back(z[i - prev.first]); } else { int j = max(i, prev.second + 1); while (j < m && t[j] == t[j - i]) ++j; z.push_back(j - i); prev = {i, j - 1}; } } return z; } vector kmp(string s, string t) { // find all t in s string cur = t + '#' + s; int sz1 = s.size(), sz2 = t.size(); vector v; vector lps = calc_next(cur); for (int i = sz2 + 1; i <= sz1 + sz2; i++) { if (lps[i] == sz2) v.push_back(i - 2 * sz2); } return v; } int period(string s) { // find the length of shortest recurring period int n = s.length(); auto z = calc_z(s); for (int i = 1; i <= n / 2; ++i) { if (n % i == 0 && z[i] == n - i) { return i; } } return n; } /* modular arithmetic */ template struct MLL { ll val; MLL(ll v = 0) : val(mod(v, mdl)) {} MLL(const MLL& other) : val(other.val) {} friend MLL operator+(const MLL& lhs, const MLL& rhs) { return mod(lhs.val + rhs.val, mdl); } friend MLL operator-(const MLL& lhs, const MLL& rhs) { return mod(lhs.val - rhs.val, mdl); } friend MLL operator*(const MLL& lhs, const MLL& rhs) { return mod(lhs.val * rhs.val, mdl); } friend MLL operator/(const MLL& lhs, const MLL& rhs) { return mod(lhs.val * mod(inverse(rhs.val, mdl), mdl), mdl); } friend MLL operator%(const MLL& lhs, const MLL& rhs) { return mod(lhs.val - (lhs / rhs).val, mdl); } friend bool operator==(const MLL& lhs, const MLL& rhs) { return lhs.val == rhs.val; } friend bool operator!=(const MLL& lhs, const MLL& rhs) { return lhs.val != rhs.val; } void operator+=(const MLL& rhs) { val = (*this + rhs).val; } void operator-=(const MLL& rhs) { val = (*this - rhs).val; } void operator*=(const MLL& rhs) { val = (*this * rhs).val; } void operator/=(const MLL& rhs) { val = (*this / rhs).val; } void operator%=(const MLL& rhs) { val = (*this % rhs).val; } }; template ostream& operator<<(ostream& out, const MLL& num) { return out << num.val; } template istream& operator>>(istream& in, MLL& num) { return in >> num.val; } // miscancellous template bool chmax(T& lhs, const U& rhs) { bool ret = lhs < rhs; if (ret) { lhs = rhs; } return ret; } template bool chmin(T& lhs, const U& rhs) { bool ret = lhs > rhs; if (ret) { lhs = rhs; } return ret; } #define functor(func) ([&](auto&&... val) \ noexcept(noexcept(func(std::forward(val)...))) -> decltype(auto) \ {return func(std::forward(val)...);}) #define expr(ret, ...) ([&] (__VA_ARGS__) { return (ret); }) template void sort_by_key(RandomIt first, RandomIt last, Func extractor) { std::sort(first, last, [&] (auto&& a, auto&& b) { return std::less<>()(extractor(a), extractor(b)); }); } template void sort_by_key(RandomIt first, RandomIt last, Func extractor, Compare comp) { std::sort(first, last, [&] (auto&& a, auto&& b) { return comp(extractor(a), extractor(b)); }); } template vector> zip(Iterator_T a_first, Iterator_T a_last, Iterator_U b_first, Iterator_U b_last) { vector> res; auto a_it = a_first; auto b_it = b_first; for (; not (a_it == a_last) and not (b_it == b_last); ++a_it, ++b_it) { res.emplace_back(*a_it, *b_it); } return res; } template vector> zip_n(Iterator_T a_first, Iterator_U b_first, size_t n) { vector> res; if (n > 0) { res.emplace_back(*a_first, *b_first); for (size_t i = 1; i != n; ++i) { res.emplace_back(*++a_first, *++b_first); } } return res; } template class ArithmeticIterator : bidirectional_iterator_tag { public: using difference_type = ptrdiff_t; using value_type = T; private: value_type value; public: ArithmeticIterator(const T& value) : value(value) {} value_type operator*() const { return value; } ArithmeticIterator& operator++() { ++value; return *this; } ArithmeticIterator& operator--() { --value; return *this; } bool operator==(const ArithmeticIterator& rhs) const { return value == rhs.value; } }; template vector> enumerate(const vector& container) { return zip(ArithmeticIterator(0), ArithmeticIterator(INT_MAX), container.begin(), container.end()); } #define initarray(init, N) (__initarray::type, (N)>(init)) namespace detail { template constexpr std::array make_array(const T& value, std::index_sequence) { return {{(static_cast(Is), value)...}}; } } template constexpr std::array __initarray(const T& value) { return detail::make_array(value, std::make_index_sequence()); } /*******************************************************/ // #define SINGLE_TEST_CASE // #define DUMP_TEST_CASE 7219 // #define TOT_TEST_CASE 10000 void dump() {} void dump_ignore() {} void prep() { } template> class segtree { private: using size_type = uint64_t; using info_type = Addable_Info_t; using tag_type = Tag_t; size_type _max; vector d; vector b; void pull(size_type p) { d[p] = d[p * 2] + d[p * 2 + 1]; } void push(size_type p, size_type left_len, size_type right_len) { if (b[p].val == INFLL) return; pushtag(p * 2, left_len, b[p]); pushtag(p * 2 + 1, right_len, b[p]); b[p] = tag_type(); } void pushtag(size_type p, size_type len, const tag_type& c) { if (d[p].mx <= c.val) return; d[p].val += (c.val - d[p].mx) * d[p].cnt; d[p].mx = c.val; b[p] = c; } void set(size_type s, size_type t, size_type p, size_type x, const info_type& c) { if (s == t) { d[p] = c; return; } size_type m = s + (t - s >> 1); if (s != t) push(p, m - s + 1, t - m); if (x <= m) set(s, m, p * 2, x, c); else set(m + 1, t, p * 2 + 1, x, c); pull(p); } void range_apply(size_type s, size_type t, size_type p, size_type l, size_type r, const tag_type& c) { if (d[p].mx <= c.val) return; // deb(s, t, p, l, r, c.val); if (l <= s && t <= r && d[p].se < c.val) { pushtag(p, t - s + 1, c); return; } size_type m = s + (t - s >> 1); push(p, m - s + 1, t - m); if (l <= m) range_apply(s, m, p * 2, l, r, c); if (r > m) range_apply(m + 1, t, p * 2 + 1, l, r, c); pull(p); } info_type range_query(size_type s, size_type t, size_type p, size_type l, size_type r) { if (l <= s && t <= r) { return d[p]; } size_type m = s + (t - s >> 1); info_type res = {}; push(p, m - s + 1, t - m); if (l <= m) res = res + range_query(s, m, p * 2, l, r); if (r > m) res = res + range_query(m + 1, t, p * 2 + 1, l, r); return res; } void build(const Sequence& a, size_type s, size_type t, size_type p) { if (s == t) { d[p] = a[s]; return; } int m = s + (t - s >> 1); build(a, s, m, p * 2); build(a, m + 1, t, p * 2 + 1); pull(p); } public: segtree(size_type __max) : d(4 * __max), b(4 * __max), _max(__max - 1) {} segtree(const Sequence& a) : segtree(a.size()) { build(a, {}, _max, 1); } void set(size_type i, const info_type& c) { set({}, _max, 1, i, c); } void range_apply(size_type l, size_type r, const tag_type& c) { range_apply({}, _max, 1, l, r, c); } void apply(size_type i, const tag_type& c) { range_apply(i, i, c); } info_type range_query(size_type l, size_type r) { return range_query({}, _max, 1, l, r); } info_type query(size_type i) { return range_query(i, i); } Sequence serialize() { Sequence res = {}; for (size_type i = 0; i <= _max; ++i) { res.push_back(query(i)); } return res; } const vector& get_d() { return d; } }; struct Tag { ll val = INFLL; }; struct Info { ll val = 0, mx = 0, se = 0, cnt = 0; // void apply(const Tag& rhs, size_t len) { } }; Info operator+(const Info &a, const Info &b) { if (a.mx == b.mx) { return { .val = a.val + b.val, .mx = a.mx, .se = max(a.se, b.se), .cnt = a.cnt + b.cnt, }; } else if (a.mx > b.mx) { return { .val = a.val + b.val, .mx = a.mx, .se = max(a.se, b.mx), .cnt = a.cnt, }; } else { return { .val = a.val + b.val, .mx = b.mx, .se = max(a.mx, b.se), .cnt = b.cnt, }; } } // constexpr int N = 2e5 + 10; // int n; // int mx[N << 2], se[N << 2], cn[N << 2], tag[N << 2]; // long long sum[N << 2]; // // void pushup(int u) { // const int ls = u << 1, rs = u << 1 | 1; // sum[u] = sum[ls] + sum[rs]; // if (mx[ls] == mx[rs]) { // mx[u] = mx[rs]; // se[u] = max(se[ls], se[rs]); // cn[u] = cn[ls] + cn[rs]; // } else if (mx[ls] > mx[rs]) { // mx[u] = mx[ls]; // se[u] = max(se[ls], mx[rs]); // cn[u] = cn[ls]; // } else { // mx[u] = mx[rs]; // se[u] = max(mx[ls], se[rs]); // cn[u] = cn[rs]; // } // } // // void build(int u = 1, int l = 1, int r = n) { // 建树 // tag[u] = -1; // if (l == r) { // sum[u] = mx[u] = INF, se[u] = -1, cn[u] = 1; // return; // } // int mid = (l + r) >> 1; // build(u << 1, l, mid), build(u << 1 | 1, mid + 1, r); // pushup(u); // } // // void pushtag(int u, int tg) { // 单纯地打标记,不暴搜 // if (mx[u] <= tg) return; // sum[u] += (1ll * tg - mx[u]) * cn[u]; // mx[u] = tag[u] = tg; // } // // void pushdown(int u) { // 下传标记 // if (tag[u] == -1) return; // pushtag(u << 1, tag[u]), pushtag(u << 1 | 1, tag[u]); // tag[u] = -1; // } // // void modify_min(int L, int R, int v, int u = 1, int l = 1, int r = n) { // if (mx[u] <= v) return; // if (L <= l && r <= R && se[u] < v) return pushtag(u, v); // int mid = (l + r) >> 1; // pushdown(u); // if (L <= mid) modify_min(L, R, v, u << 1, l, mid); // if (mid < R) modify_min(L, R, v, u << 1 | 1, mid + 1, r); // pushup(u); // } // // int query_max(int L, int R, int u = 1, int l = 1, int r = n) { // 查询最值 // if (L <= l && r <= R) return mx[u]; // int mid = (l + r) >> 1, r1 = -1, r2 = -1; // pushdown(u); // if (L <= mid) r1 = query_max(L, R, u << 1, l, mid); // if (mid < R) r2 = query_max(L, R, u << 1 | 1, mid + 1, r); // return max(r1, r2); // } // // long long query_sum(int L, int R, int u = 1, int l = 1, int r = n) { // 数值 // if (L <= l && r <= R) return sum[u]; // int mid = (l + r) >> 1; // long long res = 0; // pushdown(u); // if (L <= mid) res += query_sum(L, R, u << 1, l, mid); // if (mid < R) res += query_sum(L, R, u << 1 | 1, mid + 1, r); // return res; // } // __attribute__((target("popcnt"))) void solve() { read(int, n, k, q); // cin >> n; // memset(mx, 0, (n << 2) * sizeof(int)); // memset(se, 0, (n << 2) * sizeof(int)); // memset(cn, 0, (n << 2) * sizeof(int)); // memset(tag, 0, (n << 2) * sizeof(int)); // memset(sum, 0, (n << 2) * sizeof(ll)); // // read(int, k, q); readvec(int, a, n); set> st; vector cnt(2 * n + 1); for (int i = 0; i <= 2 * n; ++i) { st.emplace(0, i); } vector> queries(n); for (int i = 0; i < q; ++i) { read(int, l, r); --l, --r; r -= k - 1; // if (l < 0 or l >= n) exit(825); queries[l].emplace_back(r, i); } auto modify = [&] (int i, int x) -> void { // if (i < 0 or i > 2 * n) exit(826); // if (not st.count({cnt[i], i})) exit(827); st.erase({cnt[i], i}); cnt[i] += x; st.emplace(cnt[i], i); }; vector b(n - k + 1); for (int i = 0; i < n; ++i) { modify(a[i] - i + n, 1); if (i >= k - 1) { // if (i - k + 1 >= n - k + 1 or i - k + 1 < 0) exit(828); // if (st.empty()) exit(829); b[i - k + 1] = k - st.begin()->first; modify(a[i - k + 1] - (i - k + 1) + n, -1); } } segtree tr(n - k + 1); for (int i = 0; i < n - k + 1; ++i) { tr.set(i, { .val = INFLL, .mx = INFLL, .se = -1, .cnt = 1, }); } // build(); vector res(q); // auto verify = [&] (int l, int r) -> void { // if (l >= 0 and r <= n - k and l <= r) ; // else exit(830); // }; for (int i = n - k; ~i; --i) { // verify(i, n - k); tr.range_apply(i, n - k, { b[i] }); // range min // modify_min(i + 1, n - k + 1, b[i]); for (auto&& [r, idx] : queries[i]) { // verify(i, r); // if (idx < 0 or idx >= q) exit(831); res[idx] = tr.range_query(i, r).val; // range sum // res[idx] = query_sum(i + 1, r + 1); } } for (auto&& x : res) { cout << x << '\n'; } // putvec_eol(res); } int main() { #if __cplusplus < 201402L or defined(_MSC_VER) and not defined(__clang__) assert(false && "incompatible compiler variant detected."); #endif untie; prep(); #ifdef SINGLE_TEST_CASE solve(); #else read(int, t); for (int i = 0; i < t; ++i) { #ifdef DUMP_TEST_CASE if (t != (TOT_TEST_CASE)) { solve(); } else if (i + 1 == (DUMP_TEST_CASE)) { dump(); } else { dump_ignore(); } #else solve(); #endif } #endif }