/* This Submission is to determine how many 120/240 min const. delivery point there are. //info 120 req. steps <= 5 */ #define _CRT_SECURE_NO_WARNINGS #include #include #include #include #include #include #include #include #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; typedef string::const_iterator State; #define Ma_PI 3.141592653589793 #define eps 1e-5 #define LONG_INF 1e18 #define GOLD 1.61803398874989484820458 #define MAX_MOD 1000000007LL #define MOD 998244353LL #define seg_size 262144 #define REP(a,b) for(long long a = 0;a < b;++a) unsigned long xor128() { static unsigned long x = time(NULL), y = 362436069, z = 521288629, w = 88675123; unsigned long t = (x ^ (x << 11)); x = y; y = z; z = w; return (w = (w ^ (w >> 19)) ^ (t ^ (t >> 8))); } double dot(complex a, complex b) { return a.real() * b.real() + a.imag() * b.imag(); } double gyaku_dot(complex a, complex b) { return a.real() * b.imag() - a.imag() * b.real(); } double leng(complex a) { return sqrt(a.real() * a.real() + a.imag() * a.imag()); } double angles(complex a, complex b) { double cosine = dot(a, b) / (leng(a) * leng(b)); double sine = gyaku_dot(a, b) / (leng(a) * leng(b)); double kaku = acos(min((double)1.0, max((double)-1.0, cosine))); if (sine <= 0) { kaku = 2 * Ma_PI - kaku; } return kaku; } vector convex_hull(vector> a) { vector ans; double now_minnest = a[0].real(); int now_itr = 0; REP(i, a.size()) { if (now_minnest > a[i].real()) { now_minnest = a[i].real(); now_itr = i; } } ans.push_back(now_itr); complex ba(0, 1); while (true) { int now_go = 0; double now_min = 0; double now_length = 0; int starter = ans[ans.size() - 1]; for (int i = 0; i < a.size(); ++i) { if (i != starter) { double goa = angles(ba, a[i] - a[starter]); if (goa - now_min >= eps || (abs(goa - now_min) <= eps && (abs(a[i] - a[starter]) - now_length) >= eps)) { now_min = goa; now_go = i; now_length = abs(a[i] - a[starter]); } } } if (now_go == ans[0]) break; ans.push_back(now_go); ba = complex(a[now_go] - a[starter]); } return ans; } pair cnt_seg[seg_size * 2]; pair plu(pair A, pair B) { return make_pair(A.first + B.first, A.second + B.second); } pair ans_find(int now, int n_l, int n_r, int w_l, int w_r) { if (w_l <= n_l && n_r <= w_r) return cnt_seg[now]; if (w_r <= n_l || n_r <= w_l) return make_pair(0,0); return plu(ans_find(now * 2, n_l, (n_l + n_r) / 2, w_l, w_r) ,ans_find(now * 2 + 1, (n_l + n_r) / 2, n_r, w_l, w_r)); } void ans_set(int now) { cnt_seg[now] = plu(cnt_seg[now * 2], cnt_seg[now * 2 + 1]); if (now != 1) ans_set(now / 2); } long long gogo[200000]; long long ans[200000]; int main() { #define int long long iostream::sync_with_stdio(false); cin.tie(0); int n, query; cin >> n >> query; vector> gogo; REP(i, n) { long long a; cin >> a; gogo.push_back(make_pair(a, i)); } gogo.push_back(make_pair(-1, -100)); sort(gogo.begin(), gogo.end()); vector> que; REP(i, query) { long long a, b, c, d; cin >> a >> b >> c >> d; que.push_back(make_tuple(d, b - 1, c - 1, i)); } sort(que.begin(), que.end()); reverse(que.begin(), que.end()); REP(i, query) { while (get<0>(que[i]) <= gogo.back().first) { int hoge = gogo.back().second; cnt_seg[(seg_size + hoge)].first += gogo.back().first; cnt_seg[(seg_size + hoge)].second++; ans_set((seg_size + hoge) / 2); gogo.pop_back(); } pair tmp = ans_find(1, 0, seg_size, get<1>(que[i]), get<2>(que[i]) + 1); ans[get<3>(que[i])] = tmp.first - tmp.second * get<0>(que[i]); } REP(i, query) { cout << ans[i] << endl; } }