#include #define int long long using namespace std; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (a>b) { a=b; return 1; } return 0; } template struct wavelet_matrix { vector dat; vector> cnt; vector dic; wavelet_matrix(const vector &a) : dat(a) { const int n = dat.size(); dic.assign(a.begin(), a.end()); sort(dic.begin(), dic.end()); dic.erase(unique(dic.begin(), dic.end()), dic.end()); for (int i = 0; i < n; i++) { dat[i] = lower_bound(dic.begin(), dic.end(), a[i]) - dic.begin(); } int H = 1; while ((1 << H) < dic.size()) H++; cnt.resize(H, vector(n + 1)); auto dfs = [&](auto dfs, int l, int r, int h) -> void { if (r - l == 0) return; if (h == H) return; for (int i = l; i < r; i++) { cnt[h][i + 1] = cnt[h][i] + (~dat[i] >> (H - 1 - h) & 1); } int m = stable_partition(dat.begin() + l, dat.begin() + r, [&](T x) { return ~x >> (H - 1 - h) & 1; }) - dat.begin(); dfs(dfs, l, m, h + 1); dfs(dfs, m, r, h + 1); }; dfs(dfs, 0, n, 0); } T quantile(int l, int r, int k) { const int n = cnt[0].size() - 1; int a = 0; int b = n; for (int i = 0; i < cnt.size(); i++) { int cnt0 = cnt[i][b] - cnt[i][a]; int c = cnt[i][r] - cnt[i][l]; if (k < c) { l = a + cnt[i][l] - cnt[i][a]; r = a + cnt[i][r] - cnt[i][a]; b = a + cnt0; } else { l += cnt[i][b] - cnt[i][l]; r += cnt[i][b] - cnt[i][r]; k -= c; a += cnt0; } } return dic[dat[l]]; } }; vector merge(vector &a,vector &b) { int N=a.size(),M=b.size(); int i=0,j=0; vector res; while(i> data; vector> sum; public: SegmentTree(vector a){ n=1; while(n=0;i--){ data[i]=merge(data[2*i+1],data[2*i+2]); } build(); } void build(){ for(int i=0;i<2*n-1;i++){ sum[i].resize(data[i].size()+1); for(int j=0;j get(int a,int b,int x,int k,int l,int r) { // cout< res={0,0}; if(r<=a or b<=l) return res; if(a<=l and r<=b){ int i=upper_bound(data[k].begin(),data[k].end(),x)-data[k].begin(); res.first=sum[k][i]; res.second=i; return res; }else{ pair L=get(a,b,x,2*k+1,l,(l+r)/2); pair R=get(a,b,x,2*k+2,(l+r)/2,r); L.first+=R.first; L.second+=R.second; return L; } } pair get(int a,int b,int x){ return get(a,b,x,0,0,n); } }; signed main() { cin.tie(0); ios::sync_with_stdio(false); int N,Q; cin>>N>>Q; vector a(N); for(int i=0;i>a[i]; wavelet_matrix wm(a); SegmentTree seg(a); vector sum(N+1); for(int i=0;i>l>>r; l--; int med=wm.quantile(l,r,(r-l)/2); pair p=seg.get(l,r,med); int less=p.first,cnt=p.second; int more=(sum[r]-sum[l])-less; int ans=med*cnt-less+more-med*(r-l-cnt); // cout<