Submission #650349

#	Time	Username	Problem	Language	Result	Execution time	Memory
650349		600Mihnea	Fire (JOI20_ho_t5)	C++17	1 / 100	1096 ms	26288 KiB

ho_t5

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

struct Op {
	int limit;
	int k;
	int id;
};

bool operator < (Op a, Op b) {
	return a.k < b.k;
}

vector<pair<ll, int>> treemxa;
vector<ll> a;

void buildtreemxa(int v, int tl, int tr) {
	if (tl == tr) {
		treemxa[v] = {a[tl], tl};
	} else {
		int tm = (tl + tr) / 2;
		buildtreemxa(2 * v, tl, tm);
		buildtreemxa(2 * v + 1, tm + 1, tr);
		treemxa[v] = max(treemxa[2 * v], treemxa[2 * v + 1]);
	}
}

pair<ll, int> getmxa(int v, int tl, int tr, int l, int r) {
	if (tr < l || r < tl) {
		return {-1LL, 0};
	}
	if (l <= tl && tr <= r) {
		return treemxa[v];
	}
	int tm = (tl + tr) / 2;
	return max(getmxa(2 * v, tl, tm, l, r), getmxa(2 * v + 1, tm + 1, tr, l, r));
}

int main() {
#ifndef ONPC
  ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#endif
	
	int n, q;
	cin >> n >> q;
	a.resize(n);
	for (int i = 0; i < n; i++) {
		cin >> a[i];
	}
	treemxa.resize(4 * (n + 7));
	buildtreemxa(1, 0, n - 1);
	// l, r, t
	vector<ll> prn(2 * q, 0);
	vector<Op> ops;
	for (int i = 0; i < q; i++) {
		int l, r, t;
		cin >> t >> l >> r;
		l--, r--;
		ops.push_back({l - 1, t, i});
		ops.push_back({r, t, i + q});
	}
	// next bigger
	vector<int> s(n, (int) 1e9 + 7), f(n, -((int) 1e9 + 7));
	{
		vector<int> stk;
		for (int i = n - 1; i >= 0; i--) {
			while (!stk.empty() && a[stk.back()] < a[i]) {
				stk.pop_back();
			}
			if (!stk.empty()) {
				s[i] = stk.back() - 1;
			}
			stk.push_back(i);
		}	
	}
	{
		vector<int> stk;
		for (int i = 0; i < n; i++) {
			while (!stk.empty() && a[stk.back()] <= a[i]) {
				stk.pop_back();
			}
			if (!stk.empty()) {
				f[i] = stk.back() + 1;
			}
			stk.push_back(i);
		}	
	}
	sort(ops.begin(), ops.end());
	ll h = 0;
	vector<int> c1(n, 1), c2(n, 1);
	vector<int> indsfi(n); iota(indsfi.begin(), indsfi.end(), 0);
	vector<int> indssi(n); iota(indssi.begin(), indssi.end(), 0);
	sort(indsfi.begin(), indsfi.end(), [&] (int i, int j) {
		return i - f[i] > j - f[j];
	});
	sort(indssi.begin(), indssi.end(), [&] (int i, int j) {
		return s[i] - i > s[j] - j;
	});
	for (auto &it : ops) {
		int limit = it.limit;
		if (limit < 0) {
			//cout << "skip\n";
			continue;
		}
		int k = it.k;
		int id = it.id;
			
		while (!indsfi.empty()) {
			int i = indsfi.back();
			if (i - f[i] <= k) {
				indsfi.pop_back();
				c1[i] = 0;
			} else {
				break;
			}
		}
				
		while (!indssi.empty()) {
			int i = indssi.back();
			if (s[i] - i <= k) {
				indssi.pop_back();
				c2[i] = 0;
			} else {
				break;
			}
		}


		ll sol = 0;		
		{
			int i = getmxa(1, 0, n - 1, max(0, limit - k), limit).second; 
			sol += 1LL * a[i] * (limit - min(s[i], i + k));	
		}
		
		for (int i = 0; i <= limit; i++) {
			int q = max(i, f[i] + k);
			int w = min(s[i], i + k);
 			assert(q == c1[i] * i + (1 - c1[i]) * (f[i] + k));
 			assert(w == c2[i] * (i + k) + (1 - c2[i]) * s[i]);
 			if (f[i] <= min(limit, s[i]) - k) {		
 					
				sol += 1LL * a[i] * (w - q + 1);
 			}
		}
		prn[id] = sol;
	}
	for (int i = 0; i < q; i++) {
		h = h * 777777 + prn[i];
		h = h * 777777 + prn[i + q];
		cout << prn[i + q] - prn[i] << "\n";
	}
#ifdef ONPC
	assert(h == 5799234416310137250);
	cout << "h = " << h << "\n";
#endif
	return 0;
}

• Subtask 1: 1 / 1
• Subtask 2: 0 / 6
• Subtask 3: 0 / 7
• Subtask 4: 0 / 6
• Subtask 5: 0 / 80