oj.uz

Submission #757372

# Submission time Handle Problem Language Result Execution time Memory
757372 2023-06-13T06:32:03 Z SanguineChameleon Pyramid Base (IOI08_pyramid_base) C++17
80 / 100
 5000 ms 55112 KB 
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")

#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;
}

const int maxN = 1e6 + 20;
const int maxK = 4e5 + 20;
int tree[maxN << 2];
int lazy[maxN << 2];
int Lx[maxK];
int Ly[maxK];
int Rx[maxK];
int Ry[maxK];
int ly[maxK];
int ry[maxK];
int cnt_x[maxN];
int cnt_y[maxN];
int cost[maxK];
int order[maxK];
int queries[maxK << 2];
pair<int, int> events_x[maxK << 1];
pair<int, int> events_y[maxK << 1];
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 lim) {
	int n = N - sz + 1;
	int m = M - sz + 1;
	cnt_x[1] = 1;
	cnt_x[n + 1] = 1;
	for (int i = 2; i <= n; i++) {
		cnt_x[i] = 0;
	}
	cnt_y[1] = 1;
	cnt_y[m + 1] = 1;
	for (int i = 2; i <= m; i++) {
		cnt_y[i] = 0;
	}
	for (int pos = 1; pos <= min(K, lim); pos++) {
		int i = order[pos];
		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) {
			cnt_x[lx]++;
			cnt_x[rx + 1]++;
			cnt_y[ly]++;
			cnt_y[ry + 1]++;
		}
	}
	int sz_x = 0;
	for (int i = 1; i <= n + 1; i++) {
		swap(sz_x, cnt_x[i]);
		sz_x += cnt_x[i];
	}
	int sz_y = 0;
	for (int i = 1; i <= m + 1; i++) {
		swap(sz_y, cnt_y[i]);
		sz_y += cnt_y[i];
	}
	events_x[cnt_x[1]].first = 1;
	events_x[cnt_x[1]].second = 0;
	cnt_x[1]++;
	events_x[cnt_x[n + 1]].first = n + 1;
	events_x[cnt_x[n + 1]].second = 0;
	cnt_x[n + 1]++;
	events_y[cnt_y[1]].first = 1;
	events_y[cnt_y[1]].second = 0;
	cnt_y[1]++;
	events_y[cnt_y[m + 1]].first = m + 1;
	events_y[cnt_y[m + 1]].second = 0;
	cnt_y[m + 1]++;
	for (int pos = 1; pos <= min(K, lim); pos++) {
		int i = order[pos];
		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) {
			events_x[cnt_x[lx]].first = lx;
			events_x[cnt_x[lx]].second = i;
			cnt_x[lx]++;
			events_x[cnt_x[rx + 1]].first = rx + 1;
			events_x[cnt_x[rx + 1]].second = -i;
			cnt_x[rx + 1]++;
			events_y[cnt_y[ly]].first = ly;
			events_y[cnt_y[ly]].second = i;
			cnt_y[ly]++;
			events_y[cnt_y[ry + 1]].first = ry + 1;
			events_y[cnt_y[ry + 1]].second = -i;
			cnt_y[ry + 1]++;
		}
	}
	n = 0;
	int sz_q = 0;
	for (int i = 0; i < sz_x; i++) {
		if (i == 0 || events_x[i].first != events_x[i - 1].first) {
			if (n) {
				queries[sz_q++] = 0;
			}
			n++;
		}
		int id = events_x[i].second;
		if (id) {
			queries[sz_q++] = id;
		}
	}
	n--;
	m = 0;
	for (int i = 0; i < sz_y; i++) {
		if (i == 0 || events_y[i].first != events_y[i - 1].first) {
			m++;
		}
		int id = events_y[i].second;
		if (id > 0) {
			ly[id] = m;
		}
		else {
			ry[-id] = m - 1;
		}
	}
	m--;
	reset(1, 1, m);
	for (int i = 0; i < sz_q; i++) {
		int id = queries[i];
		if (id > 0) {
			update(1, 1, m, ly[id], ry[id], cost[id]);
		}
		else if (id < 0) {
			update(1, 1, m, ly[-id], ry[-id], -cost[-id]);
		}
		else if (tree[1] <= B) {
			return true;
		}
	}
	return false;
}

mt19937 gen(69420);
const int max_iters = 2;

void just_do_it() {
	cin >> N >> M >> B >> K;
	for (int i = 1; i <= K; i++) {
		order[i] = i;
	}
	int lt = 1;
	int rt = min(N, M);
	for (int i = 1; i <= K; i++) {
		cin >> Lx[i] >> Ly[i] >> Rx[i] >> Ry[i] >> cost[i];
		if (cost[i] > B) {
			int s1 = min(Lx[i] - 1, M);
			int s2 = min(N - Rx[i], M);
			int s3 = min(N, Ly[i] - 1);
			int s4 = min(N, M - Ry[i]);
			rt = min(rt, max(max(max(s1, s2), s3), s4));
		}
	}
	int res = rt;
	int lim = 10000;
	for (int iter = 0; iter < max_iters; iter++) {
		if (iter == max_iters - 1) {
			lim = K;
		}
		shuffle(order + 1, order + K + 1, gen);
		lt = 1;
		rt = res;
		res = 0;
		while (lt <= rt) {
			int mt = (lt + rt) >> 1;
			if (check(mt, lim)) {
				res = mt;
				lt = mt + 1;
			}
			else {
				rt = mt - 1;
			}
		}
	}
	cout << res;
}

Subtask #1
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 340 KB Output is correct

Subtask #2
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 340 KB Output is correct

Subtask #3
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 3 ms 468 KB Output is correct

Subtask #4
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 3 ms 468 KB Output is correct

Subtask #5
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 10 ms 1276 KB Output is correct

Subtask #6
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 16 ms 5220 KB Output is correct
2 Correct 51 ms 8352 KB Output is correct

Subtask #7
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 81 ms 7884 KB Output is correct
2 Correct 21 ms 5516 KB Output is correct

Subtask #8
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 11 ms 908 KB Output is correct
2 Correct 58 ms 1108 KB Output is correct
3 Correct 30 ms 980 KB Output is correct

Subtask #9
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 149 ms 2696 KB Output is correct
2 Correct 207 ms 2864 KB Output is correct
3 Correct 156 ms 2784 KB Output is correct

Subtask #10
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 141 ms 6776 KB Output is correct
2 Correct 38 ms 5784 KB Output is correct
3 Correct 112 ms 6744 KB Output is correct
4 Correct 306 ms 10124 KB Output is correct

Subtask #11
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 276 ms 7704 KB Output is correct
2 Correct 402 ms 11084 KB Output is correct
3 Correct 88 ms 6396 KB Output is correct

Subtask #12
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 192 ms 7088 KB Output is correct
2 Correct 490 ms 11464 KB Output is correct
3 Correct 467 ms 11560 KB Output is correct
4 Correct 520 ms 11620 KB Output is correct
5 Correct 528 ms 11592 KB Output is correct
6 Correct 69 ms 6756 KB Output is correct

Subtask #13
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 2714 ms 31772 KB Output is correct
2 Correct 982 ms 20152 KB Output is correct

Subtask #14
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5056 ms 39444 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #15
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5064 ms 55112 KB Time limit exceeded
2 Halted 0 ms 0 KB -