#include #include #include #include #include #define _USE_MATH_DEFINES #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() const int INF = 0x3f3f3f3f; const long long LINF = 0x3f3f3f3f3f3f3f3fLL; const double EPS = 1e-8; const int MOD = 1000000007; // const int MOD = 998244353; const int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1}; // const int dy[] = {1, 1, 0, -1, -1, -1, 0, 1}, // dx[] = {0, -1, -1, -1, 0, 1, 1, 1}; struct IOSetup { IOSetup() { cin.tie(nullptr); ios_base::sync_with_stdio(false); cout << fixed << setprecision(20); cerr << fixed << setprecision(10); } } iosetup; /*-------------------------------------------------*/ template struct BIT { BIT(int n, const Abelian &UNITY = 0) : n(n), UNITY(UNITY), dat(n, UNITY) {} void add(int idx, const Abelian &value) { while (idx < n) { dat[idx] += value; idx |= idx + 1; } } Abelian sum(int idx) { Abelian res = UNITY; while (idx >= 0) { res += dat[idx]; idx = (idx & (idx + 1)) - 1; } return res; } Abelian sum(int left, int right) { if (right < left) return UNITY; return sum(right) - sum(left - 1); } Abelian operator[](const int idx) { return sum(idx, idx); } int lower_bound(Abelian value) { if (value <= UNITY) return 0; int res = 0, exponent = 1; while (exponent <= n) exponent <<= 1; for (int mask = exponent >> 1; mask > 0; mask >>= 1) { if (res + mask - 1 < n && dat[res + mask - 1] < value) { value -= dat[res + mask - 1]; res += mask; } } return res; } private: int n; const Abelian UNITY; vector dat; }; int main() { int n, q; cin >> n >> q; vector a(n), comp(n); REP(i, n) { cin >> a[i]; comp[i] = a[i]; } sort(ALL(comp)); comp.erase(unique(ALL(comp)), comp.end()); int m = comp.size(); vector > sorted(m); REP(i, n) sorted[lower_bound(ALL(comp), a[i]) - comp.begin()].emplace_back(i); vector l(q), r(q); REP(i, q) { cin >> l[i] >> r[i]; --l[i]; --r[i]; } vector lb(q, -1), ub(q, m - 1); while (true) { bool is_updated = false; vector > mids(m); REP(i, q) { if (ub[i] - lb[i] > 1) { is_updated = true; mids[(lb[i] + ub[i]) / 2].emplace_back(i); } } if (!is_updated) break; BIT bit(n); REP(i, m) { for (int e : sorted[i]) bit.add(e, 1); for (int e : mids[i]) { (bit.sum(l[e], r[e]) > (r[e] - l[e] + 1) / 2 ? ub[e] : lb[e]) = i; } } } vector > medians(m); REP(i, q) medians[ub[i]].emplace_back(i); vector ans(q, 0); REP(i, q) { if ((r[i] - l[i]) % 2 == 0) ans[i] -= comp[ub[i]]; } BIT bit(n); REP(i, m) { for (int e : medians[i]) ans[e] -= bit.sum(l[e], r[e]) * 2; for (int e : sorted[i]) bit.add(e, comp[i]); } REP(i, q) { ans[i] += bit.sum(l[i], r[i]); cout << ans[i] << '\n'; } return 0; }