#pragma GCC optimize("O3,unroll-loops")
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
using namespace std;
#include "wombats.h"
#ifndef ROW
#define ROW 5000
#endif
#define COL 200
constexpr int B = 130;
constexpr int ROW_BY_B = (ROW / B) + 2;
int r, c, h[ROW][COL], v[ROW + 1][COL];
int L[9999], R[9999], By[9999], a[ROW_BY_B][COL][COL], ans[COL][COL];
void calc(int blk) {
for (int ii = 0; ii < c; ++ii) {
static int dp[2][COL];
/* find answer for (V1, V2) where V1 = ii */
memset(dp[0], 63, sizeof dp[0]);
dp[0][ii] = 0;
int I = 0;
for (int i = L[blk]; i <= R[blk]; ++i, I ^= 1) {
static int pre[COL], suf[COL], ph[COL]; /* annoying ahh math */
ph[0] = 0;
for (int j = 0; j + 1 < c; ++j)
ph[j + 1] = ph[j] + h[i][j];
pre[0] = dp[I][0];
for (int j = 1; j < c; ++j) pre[j] = min(pre[j - 1], dp[I][j] - ph[j]);
suf[c - 1] = dp[I][c - 1] + ph[c - 1];
for (int j = c - 2; j >= 0; --j) suf[j] = min(suf[j + 1], dp[I][j] + ph[j]);
for (int j = 0; j < c; ++j) dp[I][j] = min(ph[j] + pre[j], suf[j] - ph[j]);
for (int j = 0; j < c; ++j) dp[!I][j] = dp[I][j] + v[i][j];
}
memcpy(a[blk][ii], dp[I], sizeof *dp);
}
}
void merg(int up[COL][COL], int down[COL][COL], int out[COL][COL]) {
static int buf[COL][COL];
for (int i = 0; i < c; ++i) for (int j = 0; j < c; ++j) {
buf[i][j] = 1e9;
for (int k = 0; k < c; ++k)
buf[i][j] = min({buf[i][j], up[i][k]+ down[k][j]});
}
for (int i = 0; i < c; ++i) for (int j = 0; j < c; ++j) out[i][j] = buf[i][j];
}
void CalcAns() {
for (int i = 0; i < c; ++i) for (int j = 0; j < c; ++j) ans[i][j] = a[0][i][j];
for (int i = 1; i * B < r; ++i)
merg(ans, a[i], ans);
}
int outdated = 0;
void changeH(int P, int Q, int W) {
h[P][Q] = W;
calc(By[P]);
outdated = 1;
}
void changeV(int P, int Q, int W) {
v[P][Q] = W;
calc(By[P]);
outdated = 1;
}
void init(int R, int C, int H[ROW][COL], int V[ROW][COL]) {
r = R, c = C;
for (int i = 0; i + 1 < r; ++i) for (int j = 0; j < c; ++j) v[i][j] = V[i][j];
for (int i = 0 ; i < r; ++i) for (int j = 0; j + 1 < c; ++j) h[i][j] = H[i][j];
for (int i = 0; i * B < r; ++i) {
L[i] = i * B;
::R[i] = min(r - 1, i * B + B - 1);
calc(i);
for (int j = L[i]; j <= ::R[i]; ++j) By[j] = i;
}
outdated = 1;
}
int escape(int V1, int V2) { if (outdated) --outdated, CalcAns(); return ans[V1][V2]; }
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 |
4 ms |
8280 KB |
Output is correct |
2 |
Correct |
4 ms |
8284 KB |
Output is correct |
3 |
Correct |
47 ms |
9976 KB |
Output is correct |
4 |
Correct |
4 ms |
8280 KB |
Output is correct |
5 |
Correct |
4 ms |
8284 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
1 ms |
348 KB |
Output is correct |
8 |
Correct |
1 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
1 ms |
344 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
1 ms |
348 KB |
Output is correct |
10 |
Correct |
1 ms |
348 KB |
Output is correct |
11 |
Correct |
37 ms |
1364 KB |
Output is correct |
12 |
Correct |
1 ms |
344 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
796 ms |
916 KB |
Output is correct |
2 |
Correct |
697 ms |
860 KB |
Output is correct |
3 |
Correct |
799 ms |
924 KB |
Output is correct |
4 |
Correct |
815 ms |
920 KB |
Output is correct |
5 |
Correct |
789 ms |
924 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
3959 ms |
860 KB |
Output is correct |
10 |
Correct |
1 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
10 ms |
16220 KB |
Output is correct |
2 |
Correct |
10 ms |
16180 KB |
Output is correct |
3 |
Correct |
10 ms |
16228 KB |
Output is correct |
4 |
Correct |
29 ms |
16896 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
817 ms |
924 KB |
Output is correct |
2 |
Correct |
710 ms |
852 KB |
Output is correct |
3 |
Correct |
787 ms |
1100 KB |
Output is correct |
4 |
Correct |
788 ms |
924 KB |
Output is correct |
5 |
Correct |
750 ms |
924 KB |
Output is correct |
6 |
Correct |
10 ms |
16220 KB |
Output is correct |
7 |
Correct |
9 ms |
16220 KB |
Output is correct |
8 |
Correct |
10 ms |
16220 KB |
Output is correct |
9 |
Correct |
29 ms |
16976 KB |
Output is correct |
10 |
Correct |
4 ms |
8284 KB |
Output is correct |
11 |
Correct |
4 ms |
8284 KB |
Output is correct |
12 |
Correct |
38 ms |
10044 KB |
Output is correct |
13 |
Correct |
4 ms |
8284 KB |
Output is correct |
14 |
Correct |
4 ms |
8284 KB |
Output is correct |
15 |
Correct |
0 ms |
348 KB |
Output is correct |
16 |
Correct |
0 ms |
348 KB |
Output is correct |
17 |
Correct |
0 ms |
348 KB |
Output is correct |
18 |
Correct |
0 ms |
348 KB |
Output is correct |
19 |
Correct |
1 ms |
344 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
0 ms |
348 KB |
Output is correct |
23 |
Correct |
1 ms |
348 KB |
Output is correct |
24 |
Correct |
0 ms |
348 KB |
Output is correct |
25 |
Correct |
36 ms |
1364 KB |
Output is correct |
26 |
Correct |
1 ms |
344 KB |
Output is correct |
27 |
Correct |
3927 ms |
860 KB |
Output is correct |
28 |
Execution timed out |
20070 ms |
19612 KB |
Time limit exceeded |
29 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
981 ms |
860 KB |
Output is correct |
2 |
Correct |
691 ms |
860 KB |
Output is correct |
3 |
Correct |
805 ms |
916 KB |
Output is correct |
4 |
Correct |
804 ms |
860 KB |
Output is correct |
5 |
Correct |
770 ms |
1108 KB |
Output is correct |
6 |
Correct |
10 ms |
16216 KB |
Output is correct |
7 |
Correct |
10 ms |
16148 KB |
Output is correct |
8 |
Correct |
10 ms |
16216 KB |
Output is correct |
9 |
Correct |
29 ms |
17116 KB |
Output is correct |
10 |
Correct |
4 ms |
8284 KB |
Output is correct |
11 |
Correct |
5 ms |
8384 KB |
Output is correct |
12 |
Correct |
38 ms |
9948 KB |
Output is correct |
13 |
Correct |
4 ms |
8280 KB |
Output is correct |
14 |
Correct |
4 ms |
8284 KB |
Output is correct |
15 |
Correct |
3341 ms |
22464 KB |
Output is correct |
16 |
Execution timed out |
20012 ms |
22604 KB |
Time limit exceeded |
17 |
Halted |
0 ms |
0 KB |
- |