#include "candies.h"
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int BS = 768;
vector<ll> compute_function(const vector<int>& ops) {
vector<ll> stk;
stk.push_back(4e18);
stk.push_back(2e18);
for (auto delta : ops) {
if (delta == 0) continue;
if (delta > 0) {
ll rem = delta, px = 0;
while (true) {
ll x = stk.back();
if (x - px > rem) {
stk.push_back(px + rem);
stk.push_back(0);
break;
}
rem -= (x - px);
stk.pop_back();
px = stk.back();
stk.pop_back();
}
} else {
ll rem = -delta;
while (true) {
ll px = stk.back(); stk.pop_back();
ll x = stk.back();
if (x - px > rem) {
stk.push_back(px + rem);
break;
}
rem -= (x - px);
stk.pop_back();
}
}
}
stk.push_back(0);
stk.push_back(0);
reverse(stk.begin(), stk.end());
return stk;
}
vector<int> distribute_candies(
vector<int> c, vector<int> l,
vector<int> r, vector<int> v) {
int n = c.size();
// Sort each bucket increasingly by c[i] to facilitate easy evaluation of f.
vector<int> order(n);
iota(order.begin(), order.end(), 0);
for (int i = 0; i < n; i += BS) {
int j = min(i + BS, n);
sort(order.begin() + i, order.begin() + j, [&](int a, int b) {
return c[a] < c[b];
});
}
vector<ll> stk;
vector<int> a(n, 0), na(n, 0);
vector<vector<int>> ops(n / BS + 1);
auto commit = [&](int bi) {
if (ops[bi].empty()) return;
int l = bi * BS, r = min(n, (bi + 1) * BS);
// compute f(.) and evaluate it in all c[i] of this bucket
auto stk = compute_function(ops[bi]);
ll value = 0;
int ptr = 0;
for (int i = l; i < r; ++i) {
int ci = c[order[i]];
while (stk[ptr + 1] < ci) {
if (ptr % 2 == 0)
value += stk[ptr + 1] - stk[ptr];
++ptr;
}
if (ptr % 2 == 0) {
na[order[i]] = value + (ci - stk[ptr]);
} else {
na[order[i]] = value;
}
}
// Adapt for non-zero initial values.
ll mind = 0, maxd = 0, nowd = 0;
for (auto x : ops[bi]) {
nowd += x;
mind = min(mind, nowd);
maxd = max(maxd, nowd);
}
for (int i = l; i < r; ++i) {
ll hi = min(c[i] - maxd, (ll)a[i]);
ll lo = -mind;
if (lo < hi)
na[i] += hi - lo;
}
// Housekeeping.
for (int i = l; i < r; ++i)
a[i] = na[i];
ops[bi].clear();
};
int q = l.size();
for (int i = 0; i < q; ++i) {
int bl = l[i] / BS, br = r[i] / BS;
commit(bl); commit(br);
for (int j = l[i]; j <= r[i]; ) {
if (j % BS == 0 && (j + BS - 1) <= r[i]) {
ops[j / BS].push_back(v[i]);
j += BS;
} else {
a[j] = max(0, min(a[j] + v[i], c[j]));
j += 1;
}
}
if (i % BS == 0)
for (int i = 0; i < n; i += BS)
commit(i / BS);
}
for (int i = 0; i < n; i += BS)
commit(i / BS);
return a;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
204 KB |
Output is correct |
2 |
Correct |
1 ms |
256 KB |
Output is correct |
3 |
Correct |
2 ms |
332 KB |
Output is correct |
4 |
Correct |
2 ms |
332 KB |
Output is correct |
5 |
Correct |
10 ms |
376 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
4977 ms |
9388 KB |
Output is correct |
2 |
Correct |
4870 ms |
9480 KB |
Output is correct |
3 |
Correct |
4452 ms |
9496 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
204 KB |
Output is correct |
2 |
Correct |
627 ms |
5036 KB |
Output is correct |
3 |
Correct |
99 ms |
4420 KB |
Output is correct |
4 |
Correct |
3240 ms |
9392 KB |
Output is correct |
5 |
Correct |
3239 ms |
9396 KB |
Output is correct |
6 |
Correct |
3254 ms |
9288 KB |
Output is correct |
7 |
Correct |
3163 ms |
9392 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
204 KB |
Output is correct |
2 |
Correct |
1 ms |
204 KB |
Output is correct |
3 |
Correct |
1628 ms |
5044 KB |
Output is correct |
4 |
Correct |
99 ms |
5288 KB |
Output is correct |
5 |
Correct |
2753 ms |
10012 KB |
Output is correct |
6 |
Correct |
2749 ms |
10008 KB |
Output is correct |
7 |
Correct |
2542 ms |
10004 KB |
Output is correct |
8 |
Correct |
2715 ms |
10012 KB |
Output is correct |
9 |
Correct |
2284 ms |
10016 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
204 KB |
Output is correct |
2 |
Correct |
1 ms |
256 KB |
Output is correct |
3 |
Correct |
2 ms |
332 KB |
Output is correct |
4 |
Correct |
2 ms |
332 KB |
Output is correct |
5 |
Correct |
10 ms |
376 KB |
Output is correct |
6 |
Correct |
4977 ms |
9388 KB |
Output is correct |
7 |
Correct |
4870 ms |
9480 KB |
Output is correct |
8 |
Correct |
4452 ms |
9496 KB |
Output is correct |
9 |
Correct |
1 ms |
204 KB |
Output is correct |
10 |
Correct |
627 ms |
5036 KB |
Output is correct |
11 |
Correct |
99 ms |
4420 KB |
Output is correct |
12 |
Correct |
3240 ms |
9392 KB |
Output is correct |
13 |
Correct |
3239 ms |
9396 KB |
Output is correct |
14 |
Correct |
3254 ms |
9288 KB |
Output is correct |
15 |
Correct |
3163 ms |
9392 KB |
Output is correct |
16 |
Correct |
1 ms |
204 KB |
Output is correct |
17 |
Correct |
1 ms |
204 KB |
Output is correct |
18 |
Correct |
1628 ms |
5044 KB |
Output is correct |
19 |
Correct |
99 ms |
5288 KB |
Output is correct |
20 |
Correct |
2753 ms |
10012 KB |
Output is correct |
21 |
Correct |
2749 ms |
10008 KB |
Output is correct |
22 |
Correct |
2542 ms |
10004 KB |
Output is correct |
23 |
Correct |
2715 ms |
10012 KB |
Output is correct |
24 |
Correct |
2284 ms |
10016 KB |
Output is correct |
25 |
Correct |
1 ms |
204 KB |
Output is correct |
26 |
Correct |
115 ms |
4472 KB |
Output is correct |
27 |
Correct |
671 ms |
4944 KB |
Output is correct |
28 |
Execution timed out |
5050 ms |
9488 KB |
Time limit exceeded |
29 |
Halted |
0 ms |
0 KB |
- |