#include #define whlie while #define pb push_back #define eb emplace_back #define fi first #define se second #define rep(i,N) for(int i = 0; i < (N); i++) #define repr(i,N) for(int i = (N) - 1; i >= 0; i--) #define rep1(i,N) for(int i = 1; i <= (N) ; i++) #define repr1(i,N) for(int i = (N) ; i > 0 ; i--) #define each(x,v) for(auto& x : v) #define all(v) (v).begin(),(v).end() #define sz(v) ((int)(v).size()) #define vrep(v,it) for(auto it = v.begin(); it != v.end(); it++) #define vrepr(v,it) for(auto it = v.rbegin(); it != v.rend(); it++) #define ini(...) int __VA_ARGS__; in(__VA_ARGS__) #define inl(...) ll __VA_ARGS__; in(__VA_ARGS__) #define ins(...) string __VA_ARGS__; in(__VA_ARGS__) using namespace std; void solve(); using ll = long long; using vl = vector; using vi = vector; using vvi = vector< vector >; constexpr int inf = 1001001001; constexpr ll infLL = (1LL << 61) - 1; struct IoSetupNya {IoSetupNya() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); cerr << fixed << setprecision(7);} } iosetupnya; template inline bool amin(T &x, U y) { return (y < x) ? (x = y, true) : false; } template inline bool amax(T &x, U y) { return (x < y) ? (x = y, true) : false; } template ostream& operator <<(ostream& os, const pair &p) { os << p.first << " " << p.second; return os; } template istream& operator >>(istream& is, pair &p) { is >> p.first >> p.second; return is; } template ostream& operator <<(ostream& os, const vector &v) { int s = (int)v.size(); rep(i,s) os << (i ? " " : "") << v[i]; return os; } template istream& operator >>(istream& is, vector &v) { for(auto &x : v) is >> x; return is; } void in(){} template void in(T &t,U &...u){ cin >> t; in(u...);} void out(){cout << "\n";} template void out(const T &t,const U &...u){ cout << t; if(sizeof...(u)) cout << " "; out(u...);} templatevoid die(T x){out(x); exit(0);} #ifdef NyaanDebug #include "NyaanDebug.h" #define trc(...) do { cerr << #__VA_ARGS__ << " = "; dbg_out(__VA_ARGS__);} while(0) #define trca(v,...) do { cerr << #v << " = "; array_out(v , __VA_ARGS__ );} while(0) #else #define trc(...) #define trca(...) int main(){solve();} #endif using P = pair; using vp = vector

; constexpr int MOD = /**/ 1000000007; //*/ 998244353; ///////////////// // セグ木 template struct ST{ int size; vector seg; const F func; const T UNIT; ST(int N,F func , T UNIT): func(func) , UNIT(UNIT) { size = 1; while(size < N) size <<= 1; seg.assign(2 * size, UNIT); } ST(const vector &v,F func , T UNIT) : func(func) , UNIT(UNIT){ // initialize int N = (int)v.size(); size = 1; while(size < N) size <<= 1; seg.assign(2 * size , UNIT); // set for(int i = 0; i < N; i++){ seg[i + size] = v[i]; } build(); } void set(int k, T x){ seg[k + size] = x; } void build(){ for(int k = size-1; k > 0; k--){ seg[k] = func(seg[2 * k] , seg[2 * k + 1] ); } } void update(int k, T x){ k += size; seg[k] = x; while(k >>= 1){ seg[k] = func( seg[2 * k] , seg[2 * k + 1] ); } } void add(int k , T x){ k += size; seg[k] += x; while(k >>= 1){ seg[k] = func(seg[2 * k] , seg[2 * k + 1] ); } } // 半開区間[a,b)へのクエリ T query(int a, int b){ T L = UNIT, R = UNIT; for(a+=size,b+=size; a>=1,b>>=1){ if(a & 1) L = func(L,seg[a++]); if(b & 1) R = func(seg[--b],R); } return func(L, R); } T operator[](const int &k) const{ return seg[k + size]; } }; // 遅延伝播セグメント木 template struct LST{ int n,height; F f; G g; H h; T ti; E ei; vector dat; vector laz; LST(int n, F f,G g,H h,T ti,E ei): f(f),g(g),h(h),ti(ti),ei(ei){ init(n); } LST(const vector &v, F f,G g,H h,T ti,E ei): f(f),g(g),h(h),ti(ti),ei(ei){ init( (int)v.size() ); build(v); } void init(int n_){ n=1;height=0; while(n &v){ int n_=v.size(); init(n_); for(int i=0;i>i); } inline void recalc(int k){ while(k>>=1) dat[k]=f(reflect((k<<1)|0),reflect((k<<1)|1)); } void update(int a,int b,E x){ thrust(a+=n); thrust(b+=n-1); for(int l=a,r=b+1;l>=1,r>>=1){ if(l&1) laz[l]=h(laz[l],x),l++; if(r&1) --r,laz[r]=h(laz[r],x); } recalc(a); recalc(b); } void set_val(int a,T x){ thrust(a+=n); dat[a]=x;laz[a]=ei; recalc(a); } T query(int a,int b){ thrust(a+=n); thrust(b+=n-1); T vl=ti,vr=ti; for(int l=a,r=b+1;l>=1,r>>=1) { if(l&1) vl=f(vl,reflect(l++)); if(r&1) vr=f(reflect(--r),vr); } return f(vl,vr); } }; // BIT template< typename T > struct BIT { int N; int max_2beki; vector< T > data; // 初期化 1-indexedでデータを管理する 0で初期化 BIT(int size){ N = ++size; data.assign(N, 0); max_2beki = 1; while(max_2beki * 2 <= N) max_2beki *= 2; } // [0,k](閉区間)の総和 閉区間に注意! T sum(int k) { if(k < 0) return 0; // k<0のとき0を返す T ret = 0; for(++k; k > 0; k -= k & -k) ret += data[k]; return (ret); } // [l,r](閉区間)の総和 inline T sum(int l,int r){ return sum(r) - sum(l-1); } // 一点取得 更新はできないことに注意 inline T operator[](int k){ return sum(k) - sum(k-1); } // data[k] += x; void add(int k, T x) { for(++k; k < N; k += k & -k) data[k] += x; } // imos法 [l,r]にxを加算 void imos(int l,int r,T x){ add(l , x); add(r + 1 , -x); } // lower_bound sum(i)がval以上となる最小のi int lower_bound(T w){ if(w <= 0) return 0; int x = 0; for(int k = max_2beki; k > 0; k /= 2){ if(x+k <= N - 1 && data[x + k] < w){ w -= data[x + k]; x += k; } } return x; } // upper_bound sum(i)がvalより大きくなる最小のi int upper_bound(T w){ if(w < 0) return 0; int x = 0; for(int k = max_2beki; k > 0; k /= 2){ if(x+k <= N - 1 && data[x + k] <= w){ w -= data[x + k]; x += k; } } return x; } }; template struct BIT2D{ int H,W; vector> bit; BIT2D(int H,int W):H(H),W(W){ bit.resize(H + 1 , vector(W + 1 , 0)); } // 関数の入力のindexは0-originを想定 // (x,y)にwを足す // 範囲外の時は足さない void add(int x, int y, T w){ if(x < 0 || x >= H || y < 0 || y >= W) return; for(int a = (++y , ++x); a <= H; a += a & -a){ for(int b = y ; b <= W; b += b & -b){ bit[a][b] += w; } } } // imos法で[(x1,y1) , (x2,y2)]にwを足す void imos(int x1, int y1, int x2, int y2, T w){ add(x1 , y1 , w); add(x1 , y2 + 1 , -w); add(x2 + 1 , y1 , -w); add(x2 + 1 , y2 + 1 , w); } // [(0,0) , (x,y)]の和 閉区間に注意! // x,y<0の時は0 x>=H y>=Wのときはx=H-1,y=W-1とみなす // ( imos法の時は (x,y)の値を返す ) T sum(int x, int y){ if(x < 0 || y < 0) return 0; if(x >= H) return 0; if(y >= W) return 0; T ret = 0; for(int a = (++y , ++x) ; a > 0; a -= a & -a){ for(int b = y; b > 0; b -= b & -b){ ret += bit[a][b]; } } return ret; } // [(x1,y1) , (x2,y2)] の和 // x1 > x2, y1 > y2の時はswap T sum(int x1,int y1, int x2, int y2){ if(x1 > x2) swap(x1 , x2); if(y1 > y2) swap(y1 , y2); return sum(x2,y2)-sum(x2,y1-1)-sum(x1-1,y2)+sum(x1-1,y1-1); } }; void solve(){ /* using Monoid = ll; auto f = [](Monoid a,Monoid b){return a + b; }; using SEG = ST; */ ini(N,Q); vl a(N) ; in(a); vi ord(N); iota(all(ord) , 0); sort(all(ord),[&](int l,int r){return a[l] < a[r];}); trc(ord); vi invord(N); rep(i,N) invord[ord[i]] = i; int bitn = min(N,200); BIT2D bit(N , bitn); BIT2D bitexist(N,bitn); auto bidx = [&](int idx)->int{return int(1LL * idx * bitn / N);}; vi l(bitn) , r(bitn) , li(bitn) , ri(bitn); l[0] = a[ord[0]] , r[bitn-1] = a[ord[N-1]]; li[0] = 0 , ri[bitn-1] = N-1; rep(i,N-1){ if(bidx(i) + 1 == bidx(i+1)){ r[bidx(i)] = a[ord[i]] , l[bidx(i+1)] = a[ord[i+1]]; ri[bidx(i)] = i, li[bidx(i+1)] = i + 1; } } trc(l,r); rep(i,N){ trc(i , bidx(invord[i])); bit.add(i , bidx(invord[i]) , a[i]); bitexist.add(i , bidx(invord[i]) , 1); } rep(hoge,Q){ ini(c,ln,rn); ln--; rn--; inl(x); auto itr = lower_bound(all(r) , x); if(itr == r.end()){out(0); continue;} int nya = itr - r.begin(); ll ans = 0; trc(nya); if(nya != bitn - 1){ ans += bit.sum(ln , nya + 1 , rn , bitn-1); trc(ans); ans -= x * bitexist.sum(ln , nya + 1 , rn, bitn-1); trc(ans); } // itrが示す配列内を適当に探索 for(int idx = li[nya]; idx <= ri[nya]; idx++){ trc(idx , ord[idx]); if(ln <= ord[idx] && ord[idx] <= rn && a[ord[idx]] >= x) ans += a[ord[idx]] - x; } out(ans); } }