#pragma GCC optimize("Ofast", "unroll-loops")
#include "bits/extc++.h"
using namespace std;
template <typename T>
void dbgh(const T& t) {
cerr << t << endl;
}
template <typename T, typename... U>
void dbgh(const T& t, const U&... u) {
cerr << t << " | ";
dbgh(u...);
}
#ifdef DEBUG
#define dbg(...) \
cerr << "L" << __LINE__ << " [" << #__VA_ARGS__ << "]: "; \
dbgh(__VA_ARGS__);
#else
#define dbg(...)
#define cerr \
if (false) \
cerr
#endif
#define endl "\n"
#define long int64_t
#define sz(x) int((x).size())
#ifdef __cplusplus
extern "C" {
#endif
void init(int R, int C, int H[5000][200], int V[5000][200]);
void changeH(int P, int Q, int W);
void changeV(int P, int Q, int W);
int escape(int V1, int V2);
#ifdef __cplusplus
}
#endif
template <typename T, size_t N>
ostream& operator<<(ostream& out, const array<T, N>& arr) {
out << "{";
for (size_t i = 0; i < N; i++) {
out << arr[i];
if (i + 1 < N) {
out << ", ";
}
}
return out << "}";
}
const int maxn = 5005, maxm = 205, bsize = 320, maxt = 16;
// const int maxn = 32, maxm = 10, bsize = 2, maxt = 16;
static_assert((maxn - 1) / bsize < maxt);
template <size_t N, size_t M>
using mat = array<array<int, M>, N>;
template <size_t A, size_t B, size_t C>
mat<A, C> mul(const mat<A, B>& x, const mat<B, C>& y) {
mat<A, C> ans;
for (auto& a : ans) {
a.fill(1e9);
}
for (size_t i = 0; i < A; i++) {
for (size_t j = 0; j < B; j++) {
for (size_t k = 0; k < C; k++) {
ans[i][k] = min(ans[i][k], x[i][j] + y[j][k]);
}
}
}
return ans;
}
int n, m, harr[5000][200], varr[5000][200];
mat<maxm, maxm> v[maxt * 2];
mat<maxm, maxm> solve(int l, int r) {
mat<maxm, maxm> ans {};
for (int i = 0; i < maxm; i++) {
ans[i].fill(1e9);
ans[i][i] = 0;
}
for (int i = 0; i < m; i++) {
auto& dp = ans[i];
for (int j = l; j < r; j++) {
array<int, maxm> ndp;
ndp.fill(1e9);
int psum[m + 1];
psum[0] = 0;
partial_sum(harr[j], harr[j] + m, psum + 1);
int ans = 1e9;
for (int k = 0; k < m; k++) {
ans = min(ans, dp[k] - psum[k]);
ndp[k] = min(ndp[k], ans + varr[j][k] + psum[k]);
}
ans = 1e9;
for (int k = m - 1; k >= 0; k--) {
ans = min(ans, dp[k] + psum[k]);
ndp[k] = min(ndp[k], ans + varr[j][k] - psum[k]);
}
swap(dp, ndp);
}
}
return ans;
}
void update_leaf(int ind) {
int x = ind - maxt;
v[ind] = solve(x * bsize, min(n, (x + 1) * bsize));
}
void maintain(int o) {
v[o] = mul(v[o * 2], v[o * 2 + 1]);
}
void init(int _n, int _m, int _harr[5000][200], int _varr[5000][200]) {
n = _n;
m = _m;
for (int i = 0; i < n; i++) {
copy(_harr[i], _harr[i] + m, harr[i]);
copy(_varr[i], _varr[i] + m, varr[i]);
}
for (int i = maxt; i < 2 * maxt; i++) {
update_leaf(i);
}
for (int i = maxt - 1; i; i--) {
maintain(i);
}
}
void update(int ind) {
static int cnt = 0;
assert(++cnt <= 500);
ind /= bsize;
ind += maxt;
update_leaf(ind);
while (true) {
ind >>= 1;
if (!ind) {
break;
}
maintain(ind);
}
}
void changeH(int x, int y, int w) {
harr[x][y] = w;
update(x);
}
void changeV(int x, int y, int w) {
varr[x][y] = w;
update(x);
}
int escape(int a, int b) {
return v[1][a][b];
}
Compilation message
grader.c: In function 'int main()':
grader.c:15:6: warning: variable 'res' set but not used [-Wunused-but-set-variable]
15 | int res;
| ^~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3375 ms |
17464 KB |
Output is correct |
2 |
Correct |
3382 ms |
17460 KB |
Output is correct |
3 |
Correct |
3405 ms |
18996 KB |
Output is correct |
4 |
Correct |
3458 ms |
17456 KB |
Output is correct |
5 |
Correct |
3534 ms |
17456 KB |
Output is correct |
6 |
Correct |
33 ms |
5716 KB |
Output is correct |
7 |
Correct |
29 ms |
5656 KB |
Output is correct |
8 |
Correct |
28 ms |
5676 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
29 ms |
5648 KB |
Output is correct |
2 |
Correct |
38 ms |
5716 KB |
Output is correct |
3 |
Correct |
29 ms |
5680 KB |
Output is correct |
4 |
Correct |
28 ms |
5700 KB |
Output is correct |
5 |
Correct |
28 ms |
5768 KB |
Output is correct |
6 |
Correct |
30 ms |
5676 KB |
Output is correct |
7 |
Correct |
27 ms |
5732 KB |
Output is correct |
8 |
Correct |
31 ms |
5712 KB |
Output is correct |
9 |
Correct |
27 ms |
5716 KB |
Output is correct |
10 |
Correct |
31 ms |
5724 KB |
Output is correct |
11 |
Correct |
93 ms |
6692 KB |
Output is correct |
12 |
Correct |
27 ms |
5796 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1122 ms |
6052 KB |
Output is correct |
2 |
Correct |
1119 ms |
6052 KB |
Output is correct |
3 |
Correct |
1140 ms |
6048 KB |
Output is correct |
4 |
Correct |
1088 ms |
6056 KB |
Output is correct |
5 |
Correct |
1133 ms |
6056 KB |
Output is correct |
6 |
Correct |
28 ms |
5708 KB |
Output is correct |
7 |
Correct |
27 ms |
5708 KB |
Output is correct |
8 |
Correct |
27 ms |
5708 KB |
Output is correct |
9 |
Correct |
5370 ms |
6056 KB |
Output is correct |
10 |
Correct |
55 ms |
5648 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3516 ms |
21492 KB |
Output is correct |
2 |
Correct |
3379 ms |
21456 KB |
Output is correct |
3 |
Correct |
3416 ms |
21392 KB |
Output is correct |
4 |
Correct |
3449 ms |
22140 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1087 ms |
6052 KB |
Output is correct |
2 |
Correct |
1109 ms |
6048 KB |
Output is correct |
3 |
Correct |
1113 ms |
6048 KB |
Output is correct |
4 |
Correct |
1116 ms |
6172 KB |
Output is correct |
5 |
Correct |
1075 ms |
6052 KB |
Output is correct |
6 |
Correct |
3238 ms |
21376 KB |
Output is correct |
7 |
Correct |
3297 ms |
21424 KB |
Output is correct |
8 |
Correct |
3269 ms |
21372 KB |
Output is correct |
9 |
Correct |
3329 ms |
22184 KB |
Output is correct |
10 |
Correct |
3229 ms |
17464 KB |
Output is correct |
11 |
Correct |
3238 ms |
17460 KB |
Output is correct |
12 |
Correct |
3279 ms |
18968 KB |
Output is correct |
13 |
Correct |
3250 ms |
17460 KB |
Output is correct |
14 |
Correct |
3233 ms |
17464 KB |
Output is correct |
15 |
Correct |
28 ms |
5708 KB |
Output is correct |
16 |
Correct |
27 ms |
5708 KB |
Output is correct |
17 |
Correct |
27 ms |
5672 KB |
Output is correct |
18 |
Correct |
27 ms |
5772 KB |
Output is correct |
19 |
Correct |
27 ms |
5716 KB |
Output is correct |
20 |
Correct |
27 ms |
5784 KB |
Output is correct |
21 |
Correct |
28 ms |
5776 KB |
Output is correct |
22 |
Correct |
27 ms |
5688 KB |
Output is correct |
23 |
Correct |
28 ms |
5732 KB |
Output is correct |
24 |
Correct |
29 ms |
5784 KB |
Output is correct |
25 |
Correct |
91 ms |
6836 KB |
Output is correct |
26 |
Correct |
28 ms |
5780 KB |
Output is correct |
27 |
Correct |
5217 ms |
6052 KB |
Output is correct |
28 |
Correct |
9537 ms |
21964 KB |
Output is correct |
29 |
Correct |
9840 ms |
19252 KB |
Output is correct |
30 |
Correct |
9446 ms |
19228 KB |
Output is correct |
31 |
Correct |
9904 ms |
22844 KB |
Output is correct |
32 |
Correct |
41 ms |
5716 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1080 ms |
6060 KB |
Output is correct |
2 |
Correct |
1059 ms |
6052 KB |
Output is correct |
3 |
Correct |
1061 ms |
6092 KB |
Output is correct |
4 |
Correct |
1074 ms |
6048 KB |
Output is correct |
5 |
Correct |
1053 ms |
6060 KB |
Output is correct |
6 |
Correct |
3231 ms |
21452 KB |
Output is correct |
7 |
Correct |
3262 ms |
21580 KB |
Output is correct |
8 |
Correct |
3236 ms |
21372 KB |
Output is correct |
9 |
Correct |
3291 ms |
22332 KB |
Output is correct |
10 |
Correct |
3309 ms |
17456 KB |
Output is correct |
11 |
Correct |
3276 ms |
17460 KB |
Output is correct |
12 |
Correct |
3291 ms |
19116 KB |
Output is correct |
13 |
Correct |
3276 ms |
17460 KB |
Output is correct |
14 |
Correct |
3293 ms |
17460 KB |
Output is correct |
15 |
Correct |
1108 ms |
21372 KB |
Output is correct |
16 |
Execution timed out |
20102 ms |
23112 KB |
Time limit exceeded |
17 |
Halted |
0 ms |
0 KB |
- |