#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using Int = long long; template ostream &operator<<(ostream &os, const pair &a) { return os << "(" << a.first << ", " << a.second << ")"; }; template ostream &operator<<(ostream &os, const vector &as) { const int sz = as.size(); os << "["; for (int i = 0; i < sz; ++i) { if (i >= 256) { os << ", ..."; break; } if (i > 0) { os << ", "; } os << as[i]; } return os << "]"; } template void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; } template bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; } template bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; } // T: monoid representing information of an interval. // T() should return the identity. // T(S s) should represent a single element of the array. // T::pull(const T &l, const T &r) should pull two intervals. template struct SegmentTreePoint { int logN, n; vector ts; SegmentTreePoint() {} explicit SegmentTreePoint(int n_) { for (logN = 0, n = 1; n < n_; ++logN, n <<= 1) {} ts.resize(n << 1); } template explicit SegmentTreePoint(const vector &ss) { const int n_ = ss.size(); for (logN = 0, n = 1; n < n_; ++logN, n <<= 1) {} ts.resize(n << 1); for (int i = 0; i < n_; ++i) at(i) = T(ss[i]); build(); } T &at(int i) { return ts[n + i]; } void build() { for (int u = n; --u; ) pull(u); } inline void pull(int u) { ts[u].pull(ts[u << 1], ts[u << 1 | 1]); } // Changes the value of point a to s. template void change(int a, const S &s) { assert(0 <= a); assert(a < n); ts[a += n] = T(s); for (; a >>= 1; ) pull(a); } // Applies T::f(args...) to point a. template void ch(int a, F f, Args &&... args) { assert(0 <= a); assert(a < n); (ts[a += n].*f)(args...); for (; a >>= 1; ) pull(a); } // Calculates the product for [a, b). T get(int a, int b) { assert(0 <= a); assert(a <= b); assert(b <= n); if (a == b) return T(); a += n; b += n; T prodL, prodR, t; for (int aa = a, bb = b; aa < bb; aa >>= 1, bb >>= 1) { if (aa & 1) { t.pull(prodL, ts[aa++]); prodL = t; } if (bb & 1) { t.pull(ts[--bb], prodR); prodR = t; } } t.pull(prodL, prodR); return t; } // Calculates T::f(args...) of a monoid type for [a, b). // op(-, -) should calculate the product. // e() should return the identity. template #if __cplusplus >= 201402L auto #else decltype((std::declval().*F())()) #endif get(int a, int b, Op op, E e, F f, Args &&... args) { assert(0 <= a); assert(a <= b); assert(b <= n); if (a == b) return e(); a += n; b += n; auto prodL = e(), prodR = e(); for (int aa = a, bb = b; aa < bb; aa >>= 1, bb >>= 1) { if (aa & 1) prodL = op(prodL, (ts[aa++].*f)(args...)); if (bb & 1) prodR = op((ts[--bb].*f)(args...), prodR); } return op(prodL, prodR); } // Find min b s.t. T::f(args...) returns true, // when called for the partition of [a, b) from left to right. // Returns n + 1 if there is no such b. template int findRight(int a, F f, Args &&... args) { assert(0 <= a); assert(a <= n); if ((T().*f)(args...)) return a; if (a == n) return n + 1; a += n; for (; ; a >>= 1) if (a & 1) { if ((ts[a].*f)(args...)) { for (; a < n; ) { if (!(ts[a <<= 1].*f)(args...)) ++a; } return a - n + 1; } ++a; if (!(a & (a - 1))) return n + 1; } } // Find max a s.t. T::f(args...) returns true, // when called for the partition of [a, b) from right to left. // Returns -1 if there is no such a. template int findLeft(int b, F f, Args &&... args) { assert(0 <= b); assert(b <= n); if ((T().*f)(args...)) return b; if (b == 0) return -1; b += n; for (; ; b >>= 1) if ((b & 1) || b == 2) { if ((ts[b - 1].*f)(args...)) { for (; b <= n; ) { if (!(ts[(b <<= 1) - 1].*f)(args...)) --b; } return b - n - 1; } --b; if (!(b & (b - 1))) return -1; } } }; //////////////////////////////////////////////////////////////////////////////// constexpr int INF = 1001001001; struct Node { int s, l, r, m; Node() : s(0), l(0), r(0), m(0) {} Node(int val) : s(val) { l = r = m = max(val, 0); } void pull(const Node &L, const Node &R) { s = L.s + R.s; l = max(L.l, L.s + R.l); r = max(L.r + R.s, R.r); m = max({L.m, L.r + R.l, R.m}); } }; int N, Q; vector A; vector L, R; int main() { for (; ~scanf("%d%d", &N, &Q); ) { A.resize(N); for (int i = 0; i < N; ++i) { scanf("%d", &A[i]); } L.resize(Q); R.resize(Q); for (int q = 0; q < Q; ++q) { scanf("%d%d", &L[q], &R[q]); --L[q]; } vector> as(N); for (int i = 0; i < N; ++i) { as[i] = make_pair(A[i], i); } sort(as.begin(), as.end()); vector los(Q, 0), his(Q, N); for (; ; ) { vector> mids; for (int q = 0; q < Q; ++q) { if (los[q] + 1 < his[q]) { mids.emplace_back((los[q] + his[q]) / 2, q); } } if (mids.empty()) { break; } sort(mids.begin(), mids.end()); SegmentTreePoint seg(vector(N, +1)); int pos = 0; for (const auto &mid_q : mids) { const int mid = mid_q.first; const int q = mid_q.second; for (; pos < N && as[pos].first <= mid; ++pos) { seg.change(as[pos].second, -1); } const Node res = seg.get(L[q], R[q]); int score = res.s; if (res.m > 1) { score -= res.m; score += 1; } ((score < 0) ? his[q] : los[q]) = mid; } } for (int q = 0; q < Q; ++q) { printf("%d\n", his[q]); } #ifdef LOCAL cerr<<"========"<