# |
Submission time |
Handle |
Problem |
Language |
Result |
Execution time |
Memory |
753237 |
2023-06-04T22:32:00 Z |
Olympia |
K blocks (IZhO14_blocks) |
C++17 |
|
1000 ms |
59624 KB |
#include <vector>
#include <iostream>
#include <cassert>
#include <cmath>
#pragma GCC target ("avx2")
#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")
#include <map>
#include <set>
using namespace std;
vector<int64_t> p;
struct SegmentTree {
int64_t orig[262144];
int64_t vec[262144];
int64_t lazy[262144];
const int64_t INF = 1e18;
void push (int dum) {
if (lazy[dum] == -1) {
return;
}
lazy[2 * dum + 1] = lazy[dum];
lazy[2 * dum + 2] = lazy[dum];
vec[2 * dum + 1] = lazy[dum] + orig[2 * dum + 1];
vec[2 * dum + 2] = lazy[dum] + orig[2 * dum + 2];
lazy[dum] = -1;
}
void upd (int dum, int tl, int tr, int l, int r, int64_t x) {
if (tl >= l and tr <= r) {
lazy[dum] = x;
vec[dum] = orig[dum] + x;
return;
}
if (tl > r or l > tr) {
return;
}
push(dum);
upd(2 * dum + 1, tl, (tl + tr)/2, l, r, x);
upd(2 * dum + 2, (tl + tr)/2 + 1, tr, l, r, x);
vec[dum] = min(vec[2 * dum + 1], vec[2 * dum + 2]);
}
int64_t query (int dum, int tl, int tr, int l, int r) {
if (tl >= l and tr <= r) {
return vec[dum];
}
if (tl > r or l > tr) {
return INF;
}
push(dum);
return min(query(2 * dum + 1, tl, (tl + tr)/2, l, r), query(2 * dum + 2, (tl + tr)/2 + 1, tr, l, r));
}
void upd (int l, int r, int64_t x) {
upd(0, 0, (int)p.size() - 1, l, r, x);
}
int64_t query (int l, int r) {
return query(0, 0, (int)p.size() - 1, l, r);
}
void build (int dum, int tl, int tr) {
if (tl == tr) {
orig[dum] = p[tl];
return;
}
build(2 * dum + 1, tl, (tl + tr)/2);
build(2 * dum + 2, (tl + tr)/2 + 1, tr);
orig[dum] = min(orig[2 * dum + 1], orig[2 * dum + 2]);
}
void reset () {
for (int i = 0; i < 2 * p.size(); i++) {
orig[i] = 0, lazy[i] = -1;
}
build(0, 0, (int)p.size() - 1);
for (int i = 0; i < 2 * p.size(); i++) {
vec[i] = orig[i];
}
}
SegmentTree () {
int n = p.size();
for (int i = 0; i < 2 * n; i++) {
orig[i] = 0;
lazy[i] = -1;
}
build(0, 0, n - 1);
for (int i = 0; i < 2 * n; i++) {
vec[i] = orig[i];
}
}
};
struct SparseTable {
int64_t dp[131072][25];
int64_t query (int l, int r) {
int sz = log2(r - l + 1);
return min(dp[l][sz], dp[r - (1 << sz) + 1][sz]);
}
void upd (vector<int64_t> &a) {
int n = a.size();
for (int j = 0; j < 25; j++) {
for (int i = 0; i < n; i++) {
if (j == 0) {
dp[i][j] = a[i];
} else {
dp[i][j] = min(dp[i][j - 1], dp[min(i + (1 << (j - 1)), 131071)][j - 1]);
}
}
}
}
};
int main () {
ios_base::sync_with_stdio(false);
cin.tie(NULL);
int N, K;
cin >> N >> K;
vector<int64_t> arr(N);
for (int i = 0; i < N; i++) {
cin >> arr[i];
}
int64_t INF = 1e17;
vector<int64_t> cur(N + 1);
p.assign(N + 1, INF), cur.assign(N + 1, INF);
for (int j = 1; j <= N; j++) {
p[j] = ((j == 1) ? arr[0] : max(p[j - 1], arr[j - 1]));
}
for (int i = 0; i < arr.size(); i++) {
arr[i] = -arr[i];
}
SparseTable st_a;
st_a.upd(arr);
for (int i = 0; i < arr.size(); i++) {
arr[i] = -arr[i];
}
int nxt[N];
for (int j = 0; j < N; j++) {
int l = 0;
int r = j;
while (l != r) {
int m = (l + r)/2;
if (st_a.query(m, j) == -arr[j]) {
r = m;
} else {
l = m + 1;
}
}
nxt[j] = l;
}
while (__builtin_popcount(p.size()) != 1) {
p.push_back(0);
cur.push_back(0);
}
SparseTable st_p;
SegmentTree val;
for (int i = 2; i <= K; i++) {
val.reset();
st_p.upd(p);
for (int j = 1; j <= N; j++) {
if (j == 1) {
cur[j] = p[0] + arr[0];
} else {
cur[j] = INF;
int ind = nxt[j - 1];
val.upd(ind, j - 1, arr[j - 1]);
cur[j] = min(val.query(0, ind - 1), arr[j - 1] + st_p.query(ind, j - 1));
}
}
swap(cur, p);
}
cout << p[N];
}
Compilation message
blocks.cpp:6: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
6 | #pragma GCC optimization ("O3")
|
blocks.cpp:7: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
7 | #pragma GCC optimization ("unroll-loops")
|
blocks.cpp: In member function 'void SegmentTree::reset()':
blocks.cpp:67:27: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
67 | for (int i = 0; i < 2 * p.size(); i++) {
| ~~^~~~~~~~~~~~~~
blocks.cpp:71:27: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
71 | for (int i = 0; i < 2 * p.size(); i++) {
| ~~^~~~~~~~~~~~~~
blocks.cpp: In function 'int main()':
blocks.cpp:121:23: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
121 | for (int i = 0; i < arr.size(); i++) {
| ~~^~~~~~~~~~~~
blocks.cpp:126:23: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
126 | for (int i = 0; i < arr.size(); i++) {
| ~~^~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
33 ms |
57732 KB |
Output is correct |
2 |
Correct |
33 ms |
57724 KB |
Output is correct |
3 |
Correct |
34 ms |
57684 KB |
Output is correct |
4 |
Correct |
31 ms |
57648 KB |
Output is correct |
5 |
Correct |
30 ms |
57772 KB |
Output is correct |
6 |
Correct |
35 ms |
57696 KB |
Output is correct |
7 |
Correct |
31 ms |
57684 KB |
Output is correct |
8 |
Correct |
29 ms |
57680 KB |
Output is correct |
9 |
Correct |
35 ms |
57644 KB |
Output is correct |
10 |
Correct |
35 ms |
57700 KB |
Output is correct |
11 |
Correct |
31 ms |
57736 KB |
Output is correct |
12 |
Correct |
28 ms |
57684 KB |
Output is correct |
13 |
Correct |
35 ms |
57728 KB |
Output is correct |
14 |
Correct |
30 ms |
57656 KB |
Output is correct |
15 |
Correct |
34 ms |
57756 KB |
Output is correct |
16 |
Correct |
33 ms |
57672 KB |
Output is correct |
17 |
Correct |
32 ms |
57652 KB |
Output is correct |
18 |
Correct |
32 ms |
57744 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
35 ms |
57748 KB |
Output is correct |
2 |
Correct |
35 ms |
57672 KB |
Output is correct |
3 |
Correct |
34 ms |
57704 KB |
Output is correct |
4 |
Correct |
31 ms |
57712 KB |
Output is correct |
5 |
Correct |
32 ms |
57680 KB |
Output is correct |
6 |
Correct |
30 ms |
57700 KB |
Output is correct |
7 |
Correct |
36 ms |
57708 KB |
Output is correct |
8 |
Correct |
28 ms |
57752 KB |
Output is correct |
9 |
Correct |
28 ms |
57684 KB |
Output is correct |
10 |
Correct |
32 ms |
57660 KB |
Output is correct |
11 |
Correct |
33 ms |
57688 KB |
Output is correct |
12 |
Correct |
27 ms |
57684 KB |
Output is correct |
13 |
Correct |
30 ms |
57736 KB |
Output is correct |
14 |
Correct |
34 ms |
57748 KB |
Output is correct |
15 |
Correct |
30 ms |
57736 KB |
Output is correct |
16 |
Correct |
32 ms |
57840 KB |
Output is correct |
17 |
Correct |
36 ms |
57752 KB |
Output is correct |
18 |
Correct |
31 ms |
57756 KB |
Output is correct |
19 |
Correct |
30 ms |
57688 KB |
Output is correct |
20 |
Correct |
31 ms |
57660 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
33 ms |
57732 KB |
Output is correct |
2 |
Correct |
33 ms |
57724 KB |
Output is correct |
3 |
Correct |
34 ms |
57684 KB |
Output is correct |
4 |
Correct |
31 ms |
57648 KB |
Output is correct |
5 |
Correct |
30 ms |
57772 KB |
Output is correct |
6 |
Correct |
35 ms |
57696 KB |
Output is correct |
7 |
Correct |
31 ms |
57684 KB |
Output is correct |
8 |
Correct |
29 ms |
57680 KB |
Output is correct |
9 |
Correct |
35 ms |
57644 KB |
Output is correct |
10 |
Correct |
35 ms |
57700 KB |
Output is correct |
11 |
Correct |
31 ms |
57736 KB |
Output is correct |
12 |
Correct |
28 ms |
57684 KB |
Output is correct |
13 |
Correct |
35 ms |
57728 KB |
Output is correct |
14 |
Correct |
30 ms |
57656 KB |
Output is correct |
15 |
Correct |
34 ms |
57756 KB |
Output is correct |
16 |
Correct |
33 ms |
57672 KB |
Output is correct |
17 |
Correct |
32 ms |
57652 KB |
Output is correct |
18 |
Correct |
32 ms |
57744 KB |
Output is correct |
19 |
Correct |
35 ms |
57748 KB |
Output is correct |
20 |
Correct |
35 ms |
57672 KB |
Output is correct |
21 |
Correct |
34 ms |
57704 KB |
Output is correct |
22 |
Correct |
31 ms |
57712 KB |
Output is correct |
23 |
Correct |
32 ms |
57680 KB |
Output is correct |
24 |
Correct |
30 ms |
57700 KB |
Output is correct |
25 |
Correct |
36 ms |
57708 KB |
Output is correct |
26 |
Correct |
28 ms |
57752 KB |
Output is correct |
27 |
Correct |
28 ms |
57684 KB |
Output is correct |
28 |
Correct |
32 ms |
57660 KB |
Output is correct |
29 |
Correct |
33 ms |
57688 KB |
Output is correct |
30 |
Correct |
27 ms |
57684 KB |
Output is correct |
31 |
Correct |
30 ms |
57736 KB |
Output is correct |
32 |
Correct |
34 ms |
57748 KB |
Output is correct |
33 |
Correct |
30 ms |
57736 KB |
Output is correct |
34 |
Correct |
32 ms |
57840 KB |
Output is correct |
35 |
Correct |
36 ms |
57752 KB |
Output is correct |
36 |
Correct |
31 ms |
57756 KB |
Output is correct |
37 |
Correct |
30 ms |
57688 KB |
Output is correct |
38 |
Correct |
31 ms |
57660 KB |
Output is correct |
39 |
Correct |
41 ms |
57756 KB |
Output is correct |
40 |
Correct |
43 ms |
57728 KB |
Output is correct |
41 |
Correct |
43 ms |
57648 KB |
Output is correct |
42 |
Correct |
36 ms |
57644 KB |
Output is correct |
43 |
Correct |
26 ms |
57720 KB |
Output is correct |
44 |
Correct |
62 ms |
57744 KB |
Output is correct |
45 |
Correct |
38 ms |
57660 KB |
Output is correct |
46 |
Correct |
38 ms |
57752 KB |
Output is correct |
47 |
Correct |
37 ms |
57756 KB |
Output is correct |
48 |
Correct |
39 ms |
57748 KB |
Output is correct |
49 |
Correct |
30 ms |
57716 KB |
Output is correct |
50 |
Correct |
27 ms |
57676 KB |
Output is correct |
51 |
Correct |
33 ms |
57696 KB |
Output is correct |
52 |
Correct |
28 ms |
57720 KB |
Output is correct |
53 |
Correct |
32 ms |
57772 KB |
Output is correct |
54 |
Correct |
35 ms |
57752 KB |
Output is correct |
55 |
Correct |
30 ms |
57684 KB |
Output is correct |
56 |
Correct |
41 ms |
57744 KB |
Output is correct |
57 |
Correct |
31 ms |
57680 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
33 ms |
57732 KB |
Output is correct |
2 |
Correct |
33 ms |
57724 KB |
Output is correct |
3 |
Correct |
34 ms |
57684 KB |
Output is correct |
4 |
Correct |
31 ms |
57648 KB |
Output is correct |
5 |
Correct |
30 ms |
57772 KB |
Output is correct |
6 |
Correct |
35 ms |
57696 KB |
Output is correct |
7 |
Correct |
31 ms |
57684 KB |
Output is correct |
8 |
Correct |
29 ms |
57680 KB |
Output is correct |
9 |
Correct |
35 ms |
57644 KB |
Output is correct |
10 |
Correct |
35 ms |
57700 KB |
Output is correct |
11 |
Correct |
31 ms |
57736 KB |
Output is correct |
12 |
Correct |
28 ms |
57684 KB |
Output is correct |
13 |
Correct |
35 ms |
57728 KB |
Output is correct |
14 |
Correct |
30 ms |
57656 KB |
Output is correct |
15 |
Correct |
34 ms |
57756 KB |
Output is correct |
16 |
Correct |
33 ms |
57672 KB |
Output is correct |
17 |
Correct |
32 ms |
57652 KB |
Output is correct |
18 |
Correct |
32 ms |
57744 KB |
Output is correct |
19 |
Correct |
35 ms |
57748 KB |
Output is correct |
20 |
Correct |
35 ms |
57672 KB |
Output is correct |
21 |
Correct |
34 ms |
57704 KB |
Output is correct |
22 |
Correct |
31 ms |
57712 KB |
Output is correct |
23 |
Correct |
32 ms |
57680 KB |
Output is correct |
24 |
Correct |
30 ms |
57700 KB |
Output is correct |
25 |
Correct |
36 ms |
57708 KB |
Output is correct |
26 |
Correct |
28 ms |
57752 KB |
Output is correct |
27 |
Correct |
28 ms |
57684 KB |
Output is correct |
28 |
Correct |
32 ms |
57660 KB |
Output is correct |
29 |
Correct |
33 ms |
57688 KB |
Output is correct |
30 |
Correct |
27 ms |
57684 KB |
Output is correct |
31 |
Correct |
30 ms |
57736 KB |
Output is correct |
32 |
Correct |
34 ms |
57748 KB |
Output is correct |
33 |
Correct |
30 ms |
57736 KB |
Output is correct |
34 |
Correct |
32 ms |
57840 KB |
Output is correct |
35 |
Correct |
36 ms |
57752 KB |
Output is correct |
36 |
Correct |
31 ms |
57756 KB |
Output is correct |
37 |
Correct |
30 ms |
57688 KB |
Output is correct |
38 |
Correct |
31 ms |
57660 KB |
Output is correct |
39 |
Correct |
41 ms |
57756 KB |
Output is correct |
40 |
Correct |
43 ms |
57728 KB |
Output is correct |
41 |
Correct |
43 ms |
57648 KB |
Output is correct |
42 |
Correct |
36 ms |
57644 KB |
Output is correct |
43 |
Correct |
26 ms |
57720 KB |
Output is correct |
44 |
Correct |
62 ms |
57744 KB |
Output is correct |
45 |
Correct |
38 ms |
57660 KB |
Output is correct |
46 |
Correct |
38 ms |
57752 KB |
Output is correct |
47 |
Correct |
37 ms |
57756 KB |
Output is correct |
48 |
Correct |
39 ms |
57748 KB |
Output is correct |
49 |
Correct |
30 ms |
57716 KB |
Output is correct |
50 |
Correct |
27 ms |
57676 KB |
Output is correct |
51 |
Correct |
33 ms |
57696 KB |
Output is correct |
52 |
Correct |
28 ms |
57720 KB |
Output is correct |
53 |
Correct |
32 ms |
57772 KB |
Output is correct |
54 |
Correct |
35 ms |
57752 KB |
Output is correct |
55 |
Correct |
30 ms |
57684 KB |
Output is correct |
56 |
Correct |
41 ms |
57744 KB |
Output is correct |
57 |
Correct |
31 ms |
57680 KB |
Output is correct |
58 |
Correct |
131 ms |
58348 KB |
Output is correct |
59 |
Correct |
54 ms |
59624 KB |
Output is correct |
60 |
Correct |
233 ms |
59592 KB |
Output is correct |
61 |
Execution timed out |
1077 ms |
59588 KB |
Time limit exceeded |
62 |
Halted |
0 ms |
0 KB |
- |