Submission #757462

# Submission time Handle Problem Language Result Execution time Memory
757462 2023-06-13T08:16:44 Z SanguineChameleon Pyramid Base (IOI08_pyramid_base) C++17
70 / 100
1776 ms 140912 KB
#include <bits/stdc++.h>
using namespace std;

void just_do_it();

int main() {
	#ifdef KAMIRULEZ
		freopen("kamirulez.inp", "r", stdin);
		freopen("kamirulez.out", "w", stdout);
	#endif
	ios_base::sync_with_stdio(0);
	cin.tie(0);
	just_do_it();
	return 0;
}

namespace sub12 {
	const int maxN = 1e6 + 20;
	int tree[maxN << 2];
	int lazy[maxN << 2];
	int Lx[maxN];
	int Ly[maxN];
	int Rx[maxN];
	int Ry[maxN];
	int cost[maxN];
	vector<int> queries[maxN];
	int N, M, B, K;

	void reset(int id, int lt, int rt) {
		tree[id] = 0;
		lazy[id] = 0;
		if (lt == rt) {
			return;
		}
		int mt = (lt + rt) >> 1;
		reset(id << 1, lt, mt);
		reset(id << 1 | 1, mt + 1, rt);
	}

	void update(int id, int lt, int rt, int ql, int qr, int add) {
		if (lt == ql && rt == qr) {
			tree[id] += add;
			lazy[id] += add;
			return;
		}
		tree[id << 1] += lazy[id];
		lazy[id << 1] += lazy[id];
		tree[id << 1 | 1] += lazy[id];
		lazy[id << 1 | 1] += lazy[id];
		lazy[id] = 0;
		int mt = (lt + rt) >> 1;
		if (qr <= mt) {
			update(id << 1, lt, mt, ql, qr, add);
		}
		else if (ql >= mt + 1) {
			update(id << 1 | 1, mt + 1, rt, ql, qr, add);
		}
		else {
			update(id << 1, lt, mt, ql, mt, add);
			update(id << 1 | 1, mt + 1, rt, mt + 1, qr, add);
		}
		tree[id] = min(tree[id << 1], tree[id << 1 | 1]);
	}

	bool check(int sz) {
		int n = N - sz + 1;
		int m = M - sz + 1;
		for (int i = 1; i <= n; i++) {
			queries[i].clear();
		}
		for (int i = 1; i <= K; i++) {
			int lx = max(1, Lx[i] - sz + 1);
			int rx = min(n, Rx[i]);
			int ly = max(1, Ly[i] - sz + 1);
			int ry = min(m, Ry[i]);
			if (lx <= rx && ly <= ry) {
				queries[lx].push_back(i);
				if (rx < n) {
					queries[rx + 1].push_back(-i);
				}
			}
		}
		reset(1, 1, m);
		for (int i = 1; i <= n; i++) {
			for (auto id: queries[i]) {
				if (id > 0) {
					int lt = max(1, Ly[id] - sz + 1);
					int rt = min(m, Ry[id]);
					update(1, 1, m, lt, rt, cost[id]);
				}
				else {
					int lt = max(1, Ly[-id] - sz + 1);
					int rt = min(m, Ry[-id]);
					update(1, 1, m, lt, rt, -cost[-id]);
				}
			}
			if (tree[1] <= B) {
				return true;
			}
		}
		return false;
	}

	void solve() {
		for (int i = 1; i <= K; i++) {
			cin >> Lx[i] >> Ly[i] >> Rx[i] >> Ry[i] >> cost[i];
		}
		int res = 0;
		int lt = 1;
		int rt = min(N, M);
		while (lt <= rt) {
			int mt = (lt + rt) >> 1;
			if (check(mt)) {
				res = mt;
				lt = mt + 1;
			}
			else {
				rt = mt - 1;
			}
		}
		cout << res;
	}
}

namespace sub3 {
	struct node {
		int len, min_val, best_all, best_pref, best_suf, lazy;

		node() {};

		node(int _len): len(_len), min_val(0), best_all(_len), best_pref(_len), best_suf(_len), lazy(0) {};
	};

	const int maxN = 1e6 + 20;
	vector<pair<int, int>> rem[maxN];
	vector<pair<int, int>> add[maxN];
	node tree[maxN << 2];
	int N, M, B, K;

	node merge(node L, node R) {
		node res;
		res.len = L.len + R.len;
		res.min_val = min(L.min_val, R.min_val);
		res.best_all = 0;
		int sum = 0;
		if (L.min_val == res.min_val) {
			res.best_all = max(res.best_all, L.best_all);
			sum += L.best_suf;
			if (L.best_pref == L.len) {
				res.best_pref = L.len + R.best_pref;
			}
			else {
				res.best_pref = L.best_pref;
			}
		}
		else {
			res.best_pref = 0;
		}
		if (R.min_val == res.min_val) {
			res.best_all = max(res.best_all, R.best_all);
			sum += R.best_pref;
			if (R.best_suf == R.len) {
				res.best_suf = R.len + L.best_suf;
			}
			else {
				res.best_suf = R.best_suf;
			}
		}
		else {
			res.best_suf = 0;
		}
		res.best_all = max(res.best_all, sum);
		res.lazy = 0;
		return res;
	}

	void build(int id, int lt, int rt) {
		tree[id] = node(rt - lt + 1);
		if (lt == rt) {
			return;
		}
		int mt = (lt + rt) >> 1;
		build(id << 1, lt, mt);
		build(id << 1 | 1, mt + 1, rt);
	}

	void update(int id, int lt, int rt, int ql, int qr, int add) {
		if (lt == ql && rt == qr) {
			tree[id].lazy += add;
			tree[id].min_val += add;
			return;
		}
		tree[id << 1].lazy += tree[id].lazy;
		tree[id << 1].min_val += tree[id].lazy;
		tree[id << 1 | 1].lazy += tree[id].lazy;
		tree[id << 1 | 1].min_val += tree[id].lazy;
		tree[id].lazy = 0;
		int mt = (lt + rt) >> 1;
		if (qr <= mt) {
			update(id << 1, lt, mt, ql, qr, add);
		}
		else if (ql >= mt + 1) {
			update(id << 1 | 1, mt + 1, rt, ql, qr, add);
		}
		else {
			update(id << 1, lt, mt, ql, mt, add);
			update(id << 1 | 1, mt + 1, rt, mt + 1, qr, add);
		}
		tree[id] = merge(tree[id << 1], tree[id << 1 | 1]);
	}

	void solve() {
		for (int i = 1; i <= K; i++) {
			int lx, ly, rx, ry, cost;
			cin >> lx >> ly >> rx >> ry >> cost;
			add[lx].emplace_back(ly, ry);
			rem[rx].emplace_back(ly, ry);
		}
		build(1, 1, M);
		int rt = 0;
		int res = 0;
		for (int i = 1; i <= N; i++) {
			while (rt < N) {
				rt++;
				for (auto p: add[rt]) {
					update(1, 1, M, p.first, p.second, 1);
				}
				int best = (tree[1].min_val == 0 ? tree[1].best_all : 0);
				res = max(res, min(rt - i + 1, best));
				if (best < rt - i + 1) {
					for (auto p: add[rt]) {
						update(1, 1, M, p.first, p.second, -1);
					}
					if (rt > i) {
						rt--;
					}
					break;
				}
			}
			for (auto p: rem[i]) {
				update(1, 1, M, p.first, p.second, -1);
			}
		}
		cout << res;
	}
}

void just_do_it() {
	int N, M, B, K;
	cin >> N >> M >> B >> K;
	if (B > 0 || K <= 1000) {
		sub12::N = N;
		sub12::M = M;
		sub12::B = B;
		sub12::K = K;
		sub12::solve();
	}
	else {
		sub3::N = N;
		sub3::M = M;
		sub3::B = B;
		sub3::K = K;
		sub3::solve();
	}
}
# Verdict Execution time Memory Grader output
1 Correct 34 ms 70732 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 37 ms 70740 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 36 ms 70804 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 40 ms 71196 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 69 ms 73276 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 76 ms 79228 KB Output is correct
2 Correct 232 ms 87600 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 224 ms 87784 KB Output is correct
2 Correct 106 ms 79152 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 38 ms 71092 KB Output is correct
2 Correct 71 ms 71460 KB Output is correct
3 Correct 68 ms 71500 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 169 ms 75316 KB Output is correct
2 Correct 276 ms 75648 KB Output is correct
3 Correct 199 ms 75772 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 219 ms 83816 KB Output is correct
2 Correct 102 ms 75620 KB Output is correct
3 Correct 168 ms 87856 KB Output is correct
4 Correct 524 ms 95976 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 574 ms 94880 KB Output is correct
2 Correct 876 ms 97460 KB Output is correct
3 Correct 146 ms 81100 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 270 ms 84444 KB Output is correct
2 Correct 1148 ms 100068 KB Output is correct
3 Correct 1157 ms 99584 KB Output is correct
4 Correct 1226 ms 100048 KB Output is correct
5 Correct 1377 ms 100092 KB Output is correct
6 Correct 87 ms 80316 KB Output is correct
# Verdict Execution time Memory Grader output
1 Incorrect 1460 ms 131436 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 1774 ms 136360 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 1776 ms 140912 KB Output isn't correct
2 Halted 0 ms 0 KB -