Submission #650365

#	Time	Username	Problem	Language	Result	Execution time	Memory
650365		600Mihnea	Fire (JOI20_ho_t5)	C++17	1 / 100	1086 ms	30056 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
	const int INF = (int) 1e9 + 7;
	vector<int> s(n, INF), f(n, -INF);
	{
		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);
	vector<int> indsfisi(n); iota(indsfisi.begin(), indsfisi.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;
	});	
	sort(indsfisi.begin(), indsfisi.end(), [&] (int i, int j) {
		return s[i] - f[i] > s[j] - f[j];
	});
	vector<ll> foo(n, 0), bar(n, 0);
	for (int i = 0; i < n; i++) {
		foo[i] += a[i] * 1;
		bar[i] += a[i];
	}
	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]--;
				foo[i] += a[i] * (i - f[i]);
				bar[i] -= a[i];
			} else {
				break;
			}
		}
				
		while (!indssi.empty()) {
			int i = indssi.back();
			if (s[i] - i <= k) {
				indssi.pop_back();
				c2[i]--;
				foo[i] += -a[i] * (i - s[i]);
				bar[i] -= a[i];
			} else {
				break;
			}
		}

		while (!indsfisi.empty()) {
			int i = indsfisi.back();
			if (f[i] > s[i] - k) {
				indsfisi.pop_back();
				foo[i] = bar[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++) {
			if (f[i] == -INF) {
				sol += (foo[i] + bar[i] * k);
			}	
		}
	
		for (int i = 0; i < n; i++) {
 			if (0 <= f[i] && f[i] <= limit - k && i <= limit) {		
				sol += (foo[i] + bar[i] * k);
 			}
		}
	
		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