Submission #561766

# Submission time Handle Problem Language Result Execution time Memory
561766 2022-05-13T13:00:36 Z DanShaders Fish 2 (JOI22_fish2) C++17
48 / 100
4000 ms 9704 KB
//bs:sanitizers
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
using namespace std;

namespace x = __gnu_pbds;
template <typename T>
using ordered_set = x::tree<T, x::null_type, less<T>, x::rb_tree_tag, x::tree_order_statistics_node_update>;

template <typename T>
using normal_queue = priority_queue<T, vector<T>, greater<>>;

#define all(x) begin(x), end(x)
#define sz(x) ((int) (x).size())
#define x first
#define y second
using ll = long long;
using ld = long double;

const int N = 1 << 17, INF = 0x3f3f3f3f;

pair<int, int> t[2 * N];
int mod[2 * N];

pair<int, int> merge(pair<int, int> a, pair<int, int> b) {
	if (a.x == b.x) {
		return {a.x, a.y + b.y};
	}
	return min(a, b);
}

void tapply(int i, int x) {
	t[i].x += x;
	mod[i] += x;
}

void tpush(int i) {
	if (mod[i]) {
		tapply(2 * i + 1, mod[i]);
		tapply(2 * i + 2, mod[i]);
		mod[i] = 0;
	}
}

void tpull(int i) {
	t[i] = merge(t[2 * i + 1], t[2 * i + 2]);
}

void tadd(int ql, int qr, int x, int ci = 0, int cl = 0, int cr = N) {
	if (ql <= cl && cr <= qr) {
		tapply(ci, x);
		return;
	}
	if (qr <= cl || cr <= ql) {
		return;
	}
	int mid = (cl + cr) / 2;
	tpush(ci);
	tadd(ql, qr, x, 2 * ci + 1, cl, mid);
	tadd(ql, qr, x, 2 * ci + 2, mid, cr);
	tpull(ci);
}

pair<int, int> tget(int ql, int qr, int ci = 0, int cl = 0, int cr = N) {
	if (ql <= cl && cr <= qr) {
		return t[ci];
	}
	if (qr <= cl || cr <= ql) {
		return {+INF, 0};
	}
	int mid = (cl + cr) / 2;
	tpush(ci);
	return merge(
		tget(ql, qr, 2 * ci + 1, cl, mid),
		tget(ql, qr, 2 * ci + 2, mid, cr)
	);
}

pair<int, int> g[2 * N];
ll gs[2 * N];

void gset(int i, int x) {
	g[i + N] = {x, i};
	i += N;
	gs[i] = x;
	i /= 2;
	while (i) {
		g[i] = max(g[2 * i], g[2 * i + 1]);
		gs[i] = gs[2 * i] + gs[2 * i + 1];
		i /= 2;
	}
}

pair<int, int> gget(int l, int r) {
	l += N, r += N - 1;
	pair<int, int> ans = {-1, -1};
	while (l <= r) {
		if (l & 1)	ans = max(ans, g[l]);
		if (!(r & 1))	ans = max(ans, g[r]);
		l = (l + 1) / 2;
		r = (r - 1) / 2;
	}
	return ans;
}

ll gsum(int l, int r) {
	l += N, r += N - 1;
	ll ans = 0;
	while (l <= r) {
		if (l & 1)	ans += gs[l];
		if (!(r & 1))	ans += gs[r];
		l = (l + 1) / 2;
		r = (r - 1) / 2;
	}
	return ans;
}

int gfirstgreater(int ql, int qr, ll x, int ci = 1, int cl = 0, int cr = N) {
	if (qr <= cl || cr <= ql || g[ci].x <= x) {
		return qr;
	}
	if (cr - cl == 1) {
		return cl;
	}
	int mid = (cl + cr) / 2;
	int res = gfirstgreater(ql, qr, x, 2 * ci, cl, mid);
	if (res >= qr) {
		return gfirstgreater(ql, qr, x, 2 * ci + 1, mid, cr);
	} else {
		return res;
	}
}

int glastgreater(int ql, int qr, ll x, int ci = 1, int cl = 0, int cr = N) {
	if (qr <= cl || cr <= ql || g[ci].x <= x) {
		return ql - 1;
	}
	if (cr - cl == 1) {
		return cl;
	}
	int mid = (cl + cr) / 2;
	int res = glastgreater(ql, qr, x, 2 * ci + 1, mid, cr);
	if (res < ql) {
		return glastgreater(ql, qr, x, 2 * ci, cl, mid);
	} else {
		return res;
	}
}

int a[N], n;
int up[N];

struct S {
	int i, x, cl, cr;

	bool operator<(const S &other) const {
		return i < other.i;
	}

	bool operator==(const S &other) const {
		return i == other.i;
	}
};

void set_up_to(const S &s) {
	if (up[s.i] != s.x) {
		up[s.i] = s.x;
		// cout << "tadd " << s.cl + 1 << " " << s.cr << " " << (s.x ? 1 : -1) << endl;
		if (!s.x) {
			tadd(s.cl + 1, s.cr, -1);
		} else {
			tadd(s.cl + 1, s.cr, 1);
		}
	}
}

void loop_set_up(int u, int l, int r, int x, vector<S> &change) {
	if (!(l <= u && u < r)) {
		return;
	}
	auto ai = [&](int i) {
		if (i < l || i >= r) {
			return +INF;
		}
		return a[i];
	};

	int cl = glastgreater(l, u, a[u]),
		cr = gfirstgreater(u + 1, r, a[u] - 1);
	ll sm = gsum(cl + 1, cr);
	while (cl != l - 1 || cr != r) {
		if (sm < min(ai(cl), ai(cr))) {
			int v = gget(cl + 1, cr).y;
			// cout << "set " << cl << " " << cr << " " << v << endl;
			change.push_back({v, x, cl, cr});
			if (ai(cl) <= ai(cr)) {
				sm += a[cl];
				--cl;
			} else {
				sm += a[cr];
				++cr;
			}
		}
		cl = glastgreater(l, cl + 1, sm);
		cr = gfirstgreater(cr - 1, r, sm);
		sm = gsum(cl + 1, cr);
	}
}

int ai(int i) {
	if (i < 0 || i >= n) {
		return +INF;
	}
	return a[i];
}

vector<S> change;

int query(int l, int r) {
	change.clear();

	if (l) {
		loop_set_up(l, 0, n, 0, change);
		loop_set_up(l, l, r, 1, change);
	}
	if (r != n) {
		loop_set_up(r - 1, 0, n, 0, change);
		loop_set_up(r - 1, l, r, 1, change);
	}

	stable_sort(all(change));
	reverse(all(change));
	change.resize(unique(all(change)) - begin(change));

	for (auto &s : change) {
		// cout << "temp " << s.cl << " " << s.cr << " " << s.i << " " << s.x << endl;
		int backup = up[s.i];
		set_up_to(s);
		s.x = backup;
	}
	auto [ansx, ansy] = tget(l, r);
	// cout << ansx << " " << ansy << endl;
	for (const auto &s : change) {
		set_up_to(s);
	}
	return ansx ? 0 : ansy;
}

void aset(int i, int x) {
	int cl = glastgreater(0, i, a[i]),
		cr = gfirstgreater(i + 1, n, a[i] - 1),
		nl = glastgreater(0, i, x),
		nr = gfirstgreater(i + 1, n, x - 1);
	// cout << cr << endl;
	int clef = gget(cl + 1, i).y, crig = gget(i + 1, cr).y,
		nlef = gget(nl + 1, i).y, nrig = gget(i + 1, nr).y;

	// cout << "counted " << clef << " " << crig << " " << nlef << " " << nrig << endl;

	change.clear();
	loop_set_up(i, 0, n, 0, change);
	loop_set_up(cl, 0, n, 0, change);
	loop_set_up(cr, 0, n, 0, change);
	loop_set_up(clef, 0, n, 0, change);
	loop_set_up(crig, 0, n, 0, change);
	loop_set_up(nl, 0, n, 0, change);
	loop_set_up(nr, 0, n, 0, change);
	loop_set_up(nlef, 0, n, 0, change);
	loop_set_up(nrig, 0, n, 0, change);

	stable_sort(all(change));
	reverse(all(change));
	change.resize(unique(all(change)) - begin(change));
	for (auto s : change) {
		// cout << "rollback " << s.cl << " " << s.cr << " " << s.i << " " << s.x << endl;
		set_up_to(s);
	}

	a[i] = x;
	gset(i, x);

	change.clear();
	loop_set_up(i, 0, n, 1, change);
	loop_set_up(cl, 0, n, 1, change);
	loop_set_up(cr, 0, n, 1, change);
	loop_set_up(clef, 0, n, 1, change);
	loop_set_up(crig, 0, n, 1, change);
	loop_set_up(nl, 0, n, 1, change);
	loop_set_up(nr, 0, n, 1, change);
	loop_set_up(nlef, 0, n, 1, change);
	loop_set_up(nrig, 0, n, 1, change);

	stable_sort(all(change));
	reverse(all(change));
	change.resize(unique(all(change)) - begin(change));
	for (auto s : change) {
		// cout << "apply " << s.cl << " " << s.cr << " " << s.i << " " << s.x << endl;
		set_up_to(s);
	}
}

signed main() {
	cin.tie(0)->sync_with_stdio(0);
	cin >> n;
	for (int i = 0; i < n; ++i) {
		cin >> a[i];
		g[N + i] = {a[i], i};
		gs[N + i] = a[i];
	}
	for (int i = 0; i < n; ++i) {
		t[N - 1 + i] = {0, 1};
	}
	for (int i = N - 1; i--; ) {
		t[i] = merge(t[2 * i + 1], t[2 * i + 2]);
		mod[i] = 0;
	}
	for (int i = N; i--; ) {
		g[i] = max(g[2 * i], g[2 * i + 1]);
		gs[i] = gs[2 * i] + gs[2 * i + 1];
	}

	for (int i = 0; i < n; ++i) {
		int cl = glastgreater(0, i, a[i]);
		int cr = gfirstgreater(i + 1, n, a[i] - 1);
		if (gsum(cl + 1, cr) < min(ai(cl), ai(cr)) &&
			(cl != -1 || cr != n)) {
			set_up_to({i, 1, cl, cr});
		}
	}

	int queries;
	cin >> queries;
	while (queries--) {
		int type, l, r;
		cin >> type >> l >> r;
		if (type == 1) {
			aset(l - 1, r);
		} else {
			cout << query(l - 1, r) << "\n";
		}
	}
}
# Verdict Execution time Memory Grader output
1 Correct 3 ms 3924 KB Output is correct
2 Correct 4 ms 3888 KB Output is correct
3 Correct 3 ms 3924 KB Output is correct
4 Correct 4 ms 3924 KB Output is correct
5 Correct 23 ms 3980 KB Output is correct
6 Correct 16 ms 3924 KB Output is correct
7 Correct 24 ms 3996 KB Output is correct
8 Correct 20 ms 3988 KB Output is correct
9 Correct 15 ms 3924 KB Output is correct
10 Correct 6 ms 3924 KB Output is correct
11 Correct 5 ms 3924 KB Output is correct
12 Correct 9 ms 3916 KB Output is correct
13 Correct 10 ms 3924 KB Output is correct
14 Correct 9 ms 3924 KB Output is correct
15 Correct 10 ms 3916 KB Output is correct
16 Correct 9 ms 3980 KB Output is correct
17 Correct 10 ms 3924 KB Output is correct
18 Correct 10 ms 3924 KB Output is correct
19 Correct 7 ms 3924 KB Output is correct
20 Correct 8 ms 3972 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 3924 KB Output is correct
2 Correct 59 ms 7372 KB Output is correct
3 Correct 57 ms 7352 KB Output is correct
4 Correct 60 ms 7420 KB Output is correct
5 Correct 58 ms 7404 KB Output is correct
6 Correct 48 ms 7372 KB Output is correct
7 Correct 59 ms 7336 KB Output is correct
8 Correct 68 ms 7372 KB Output is correct
9 Correct 51 ms 7456 KB Output is correct
10 Correct 53 ms 7412 KB Output is correct
11 Correct 46 ms 7436 KB Output is correct
12 Correct 43 ms 7328 KB Output is correct
13 Correct 44 ms 7436 KB Output is correct
14 Correct 50 ms 7324 KB Output is correct
15 Correct 59 ms 7348 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 3924 KB Output is correct
2 Correct 4 ms 3888 KB Output is correct
3 Correct 3 ms 3924 KB Output is correct
4 Correct 4 ms 3924 KB Output is correct
5 Correct 23 ms 3980 KB Output is correct
6 Correct 16 ms 3924 KB Output is correct
7 Correct 24 ms 3996 KB Output is correct
8 Correct 20 ms 3988 KB Output is correct
9 Correct 15 ms 3924 KB Output is correct
10 Correct 6 ms 3924 KB Output is correct
11 Correct 5 ms 3924 KB Output is correct
12 Correct 9 ms 3916 KB Output is correct
13 Correct 10 ms 3924 KB Output is correct
14 Correct 9 ms 3924 KB Output is correct
15 Correct 10 ms 3916 KB Output is correct
16 Correct 9 ms 3980 KB Output is correct
17 Correct 10 ms 3924 KB Output is correct
18 Correct 10 ms 3924 KB Output is correct
19 Correct 7 ms 3924 KB Output is correct
20 Correct 8 ms 3972 KB Output is correct
21 Correct 3 ms 3924 KB Output is correct
22 Correct 59 ms 7372 KB Output is correct
23 Correct 57 ms 7352 KB Output is correct
24 Correct 60 ms 7420 KB Output is correct
25 Correct 58 ms 7404 KB Output is correct
26 Correct 48 ms 7372 KB Output is correct
27 Correct 59 ms 7336 KB Output is correct
28 Correct 68 ms 7372 KB Output is correct
29 Correct 51 ms 7456 KB Output is correct
30 Correct 53 ms 7412 KB Output is correct
31 Correct 46 ms 7436 KB Output is correct
32 Correct 43 ms 7328 KB Output is correct
33 Correct 44 ms 7436 KB Output is correct
34 Correct 50 ms 7324 KB Output is correct
35 Correct 59 ms 7348 KB Output is correct
36 Correct 99 ms 8064 KB Output is correct
37 Correct 83 ms 7720 KB Output is correct
38 Correct 76 ms 7952 KB Output is correct
39 Correct 92 ms 8072 KB Output is correct
40 Correct 82 ms 7824 KB Output is correct
41 Correct 57 ms 8304 KB Output is correct
42 Correct 53 ms 8276 KB Output is correct
43 Correct 65 ms 7668 KB Output is correct
44 Correct 63 ms 7644 KB Output is correct
45 Correct 97 ms 8012 KB Output is correct
46 Correct 70 ms 8012 KB Output is correct
47 Correct 77 ms 7792 KB Output is correct
48 Correct 64 ms 7704 KB Output is correct
49 Correct 57 ms 7628 KB Output is correct
50 Correct 62 ms 8000 KB Output is correct
51 Correct 75 ms 8020 KB Output is correct
52 Correct 66 ms 8008 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 3924 KB Output is correct
2 Correct 59 ms 7372 KB Output is correct
3 Correct 57 ms 7352 KB Output is correct
4 Correct 60 ms 7420 KB Output is correct
5 Correct 58 ms 7404 KB Output is correct
6 Correct 48 ms 7372 KB Output is correct
7 Correct 59 ms 7336 KB Output is correct
8 Correct 68 ms 7372 KB Output is correct
9 Correct 51 ms 7456 KB Output is correct
10 Correct 53 ms 7412 KB Output is correct
11 Correct 46 ms 7436 KB Output is correct
12 Correct 43 ms 7328 KB Output is correct
13 Correct 44 ms 7436 KB Output is correct
14 Correct 50 ms 7324 KB Output is correct
15 Correct 59 ms 7348 KB Output is correct
16 Correct 3 ms 3924 KB Output is correct
17 Correct 2028 ms 7920 KB Output is correct
18 Correct 1843 ms 7828 KB Output is correct
19 Correct 2059 ms 8112 KB Output is correct
20 Correct 2716 ms 7788 KB Output is correct
21 Correct 1932 ms 7760 KB Output is correct
22 Correct 1857 ms 7908 KB Output is correct
23 Correct 2032 ms 7828 KB Output is correct
24 Correct 2736 ms 7732 KB Output is correct
25 Correct 2059 ms 7812 KB Output is correct
26 Correct 2835 ms 7800 KB Output is correct
27 Correct 502 ms 7952 KB Output is correct
28 Correct 516 ms 8012 KB Output is correct
29 Correct 435 ms 7972 KB Output is correct
30 Correct 2450 ms 7772 KB Output is correct
31 Correct 2588 ms 7748 KB Output is correct
32 Correct 3972 ms 7972 KB Output is correct
33 Correct 1378 ms 7936 KB Output is correct
34 Correct 3979 ms 7876 KB Output is correct
35 Correct 2219 ms 7976 KB Output is correct
36 Correct 2887 ms 7956 KB Output is correct
37 Correct 871 ms 7780 KB Output is correct
38 Correct 849 ms 7784 KB Output is correct
39 Correct 648 ms 8028 KB Output is correct
40 Correct 660 ms 8020 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 3924 KB Output is correct
2 Correct 59 ms 7372 KB Output is correct
3 Correct 57 ms 7352 KB Output is correct
4 Correct 60 ms 7420 KB Output is correct
5 Correct 58 ms 7404 KB Output is correct
6 Correct 48 ms 7372 KB Output is correct
7 Correct 59 ms 7336 KB Output is correct
8 Correct 68 ms 7372 KB Output is correct
9 Correct 51 ms 7456 KB Output is correct
10 Correct 53 ms 7412 KB Output is correct
11 Correct 46 ms 7436 KB Output is correct
12 Correct 43 ms 7328 KB Output is correct
13 Correct 44 ms 7436 KB Output is correct
14 Correct 50 ms 7324 KB Output is correct
15 Correct 59 ms 7348 KB Output is correct
16 Correct 3 ms 3924 KB Output is correct
17 Correct 3719 ms 7656 KB Output is correct
18 Correct 1615 ms 9596 KB Output is correct
19 Correct 2566 ms 9132 KB Output is correct
20 Correct 1285 ms 9704 KB Output is correct
21 Execution timed out 4057 ms 9008 KB Time limit exceeded
22 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 3 ms 3924 KB Output is correct
2 Correct 4 ms 3888 KB Output is correct
3 Correct 3 ms 3924 KB Output is correct
4 Correct 4 ms 3924 KB Output is correct
5 Correct 23 ms 3980 KB Output is correct
6 Correct 16 ms 3924 KB Output is correct
7 Correct 24 ms 3996 KB Output is correct
8 Correct 20 ms 3988 KB Output is correct
9 Correct 15 ms 3924 KB Output is correct
10 Correct 6 ms 3924 KB Output is correct
11 Correct 5 ms 3924 KB Output is correct
12 Correct 9 ms 3916 KB Output is correct
13 Correct 10 ms 3924 KB Output is correct
14 Correct 9 ms 3924 KB Output is correct
15 Correct 10 ms 3916 KB Output is correct
16 Correct 9 ms 3980 KB Output is correct
17 Correct 10 ms 3924 KB Output is correct
18 Correct 10 ms 3924 KB Output is correct
19 Correct 7 ms 3924 KB Output is correct
20 Correct 8 ms 3972 KB Output is correct
21 Correct 3 ms 3924 KB Output is correct
22 Correct 59 ms 7372 KB Output is correct
23 Correct 57 ms 7352 KB Output is correct
24 Correct 60 ms 7420 KB Output is correct
25 Correct 58 ms 7404 KB Output is correct
26 Correct 48 ms 7372 KB Output is correct
27 Correct 59 ms 7336 KB Output is correct
28 Correct 68 ms 7372 KB Output is correct
29 Correct 51 ms 7456 KB Output is correct
30 Correct 53 ms 7412 KB Output is correct
31 Correct 46 ms 7436 KB Output is correct
32 Correct 43 ms 7328 KB Output is correct
33 Correct 44 ms 7436 KB Output is correct
34 Correct 50 ms 7324 KB Output is correct
35 Correct 59 ms 7348 KB Output is correct
36 Correct 99 ms 8064 KB Output is correct
37 Correct 83 ms 7720 KB Output is correct
38 Correct 76 ms 7952 KB Output is correct
39 Correct 92 ms 8072 KB Output is correct
40 Correct 82 ms 7824 KB Output is correct
41 Correct 57 ms 8304 KB Output is correct
42 Correct 53 ms 8276 KB Output is correct
43 Correct 65 ms 7668 KB Output is correct
44 Correct 63 ms 7644 KB Output is correct
45 Correct 97 ms 8012 KB Output is correct
46 Correct 70 ms 8012 KB Output is correct
47 Correct 77 ms 7792 KB Output is correct
48 Correct 64 ms 7704 KB Output is correct
49 Correct 57 ms 7628 KB Output is correct
50 Correct 62 ms 8000 KB Output is correct
51 Correct 75 ms 8020 KB Output is correct
52 Correct 66 ms 8008 KB Output is correct
53 Correct 3 ms 3924 KB Output is correct
54 Correct 2028 ms 7920 KB Output is correct
55 Correct 1843 ms 7828 KB Output is correct
56 Correct 2059 ms 8112 KB Output is correct
57 Correct 2716 ms 7788 KB Output is correct
58 Correct 1932 ms 7760 KB Output is correct
59 Correct 1857 ms 7908 KB Output is correct
60 Correct 2032 ms 7828 KB Output is correct
61 Correct 2736 ms 7732 KB Output is correct
62 Correct 2059 ms 7812 KB Output is correct
63 Correct 2835 ms 7800 KB Output is correct
64 Correct 502 ms 7952 KB Output is correct
65 Correct 516 ms 8012 KB Output is correct
66 Correct 435 ms 7972 KB Output is correct
67 Correct 2450 ms 7772 KB Output is correct
68 Correct 2588 ms 7748 KB Output is correct
69 Correct 3972 ms 7972 KB Output is correct
70 Correct 1378 ms 7936 KB Output is correct
71 Correct 3979 ms 7876 KB Output is correct
72 Correct 2219 ms 7976 KB Output is correct
73 Correct 2887 ms 7956 KB Output is correct
74 Correct 871 ms 7780 KB Output is correct
75 Correct 849 ms 7784 KB Output is correct
76 Correct 648 ms 8028 KB Output is correct
77 Correct 660 ms 8020 KB Output is correct
78 Correct 3 ms 3924 KB Output is correct
79 Correct 3719 ms 7656 KB Output is correct
80 Correct 1615 ms 9596 KB Output is correct
81 Correct 2566 ms 9132 KB Output is correct
82 Correct 1285 ms 9704 KB Output is correct
83 Execution timed out 4057 ms 9008 KB Time limit exceeded
84 Halted 0 ms 0 KB -