#include "fish.h"
#include <algorithm>
#include <bitset>
#include <cstdint>
#include <cstring>
#include <iostream>
#include <limits.h>
#include <math.h>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>
#define ll long long
#define loop(X, N) for(int X = 0; X < (N); X++)
#define all(V) V.begin(), V.end()
#define rall(V) V.rbegin(), V.rend()
using namespace std;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int, int> ii;
typedef vector<ii> vii;
typedef vector<vector<ii>> vvii;
typedef vector<bool> vb;
typedef vector<vector<bool>> vvb;
vector<vector<pair<int, ll>>> fishes;
vector<map<int, ll>> prefixWeight;
vector<ll> totals;
inline ll weightSum(int x, int y1, int y2) {
if (x >= fishes.size() || y2 <= y1)
return 0;
auto it1 = prefixWeight[x].upper_bound(y2 - 1);
auto it2 = prefixWeight[x].upper_bound(y1 - 1);
ll res = (--it1)->second - (--it2)->second;
//cout << x << ' ' << y1 << ' ' << y2 << " = " << res << endl;
return res;
}
long long max_weights(signed N, signed M, vector<signed> X, std::vector<signed> Y, std::vector<signed> W) {
fishes = vector<vector<pair<int, ll>>>(N);
loop(i, M) {
fishes[X[i]].push_back({Y[i], W[i]});
}
fishes.push_back({});
prefixWeight = vector<map<int, ll>>(N + 1);
totals = vector<ll>(N + 1);
loop(x, N + 1) {
sort(all(fishes[x]));
ll total = 0;
prefixWeight[x][-2] = 0;
for (auto [y, w] : fishes[x]) {
total += w;
prefixWeight[x][y] = total;
}
totals[x] = total;
}
vector<map<int, ll>> beforeDp(N);
vector<map<int, ll>> bothDp(N);
vector<ll> maxBefore(N);
vector<ll> maxBoth(N);
vector<ll> maxBeforeMinCur(N), maxBeforeMinNext(N);
//bestBefore, bestBeforeAndAfter
int counter = 0;
loop(x, N) {
if (x == 0) {
for (const auto& [y, w] : fishes[x + 1]) {
int height = y + 1;
bothDp[0][height] = weightSum(x + 1, 0, height);
maxBoth[x] = max(maxBoth[x], bothDp[0][height]);
maxBeforeMinCur[x] = max(maxBeforeMinCur[x], 0 - weightSum(x, 0, height));
maxBeforeMinNext[x] = max(maxBeforeMinNext[x], 0 - weightSum(x + 1, 0, height));
}
beforeDp[0][0] = 0;
maxBefore[0] = 0;
continue;
}
ll bestBefore2 = 0;
if (x >= 2) {
bestBefore2 = maxBefore[x - 2];
}
if (x >= 3) {
bestBefore2 = max(bestBefore2, maxBoth[x - 3]);
}
for (const auto& [y, w] : fishes.at(x - 1)) {
int height = y + 1;
ll& both = bothDp[x][height];
ll& before = beforeDp[x][height];
both = maxBeforeMinCur[x - 1] + weightSum(x - 1, 0, height);
before = maxBeforeMinCur[x - 1] + weightSum(x - 1, 0, height);
// for (const auto& [prevHeight, prevScore] : beforeDp.at(x - 1)) {
// counter++;
// //ll prevCur = weightSum(x - 1, prevHeight, height);
// ll prevCur = -weightSum(x - 1, 0, prevHeight);
// both = max(both, prevScore + prevCur);
// before = max(before, prevScore + prevCur);
// }
// both += weightSum(x - 1, 0, height);
// before += weightSum(x - 1, 0, height);
both = max(both, maxBeforeMinNext[x]);
// for (const auto& [prevY, prevScore] : bothDp.at(x - 1)) {
// counter++;
// if (prevY <= height)
// continue;
// ll prevCur = -weightSum(x, 0, height);
// both = max(both, prevScore + prevCur);
// }
ll prevWeight = weightSum(x - 1, 0, height);
both = max(both, bestBefore2 + prevWeight);
before = max(before, bestBefore2 + prevWeight);
maxBefore[x] = max(maxBefore[x], before);
maxBeforeMinCur[x] = max(maxBeforeMinCur[x], before - weightSum(x, 0, height));
maxBeforeMinNext[x] = max(maxBeforeMinNext[x], before - weightSum(x + 1, 0, height));
}
for (const auto& [y, w] : fishes.at(x + 1)) {
int height = y + 1;
ll& both = bothDp[x][height];
ll& before = beforeDp[x][height];
both = maxBeforeMinCur[x - 1] + weightSum(x - 1, 0, height);
before = maxBeforeMinCur[x - 1] + weightSum(x - 1, 0, height);
// for (const auto& [prevHeight, prevScore] : beforeDp.at(x - 1)) {
// counter++;
// //ll prevCur = weightSum(x - 1, prevHeight, height);
// ll prevCur = -weightSum(x - 1, 0, prevHeight);
// both = max(both, prevScore + prevCur);
// before = max(before, prevScore + prevCur);
// }
// both += weightSum(x - 1, 0, height);
// before += weightSum(x - 1, 0, height);
both = max(both, maxBeforeMinNext[x]);
// for (const auto& [prevY, prevScore] : bothDp.at(x - 1)) {
// counter++;
// if (prevY <= height)
// continue;
// ll prevCur = -weightSum(x, 0, height);
// both = max(both, prevScore + prevCur);
// }
ll prevWeight = weightSum(x - 1, 0, height);
both = max(both, bestBefore2 + prevWeight);
before = max(before, bestBefore2 + prevWeight);
maxBefore[x] = max(maxBefore[x], before);
maxBeforeMinCur[x] = max(maxBeforeMinCur[x], before - weightSum(x, 0, height));
maxBeforeMinNext[x] = max(maxBeforeMinNext[x], before - weightSum(x + 1, 0, height));
}
for (const auto& [y, w] : fishes.at(x + 1)) {
int height = y + 1;
ll& both = bothDp[x][height];
both = max(both, maxBeforeMinCur[x - 1] + weightSum(x - 1, 0, height));
both = max(both, maxBeforeMinNext[x]);
for (const auto& [prevHeight, prevScore] : beforeDp.at(x - 1)) {
counter++;
ll prevCur = weightSum(x - 1, prevHeight, height);
both = max(both, prevScore + prevCur);
}
for (const auto& [prevY, prevScore] : bothDp.at(x - 1)) {
counter++;
if (prevY <= height)
continue;
ll prevCur = -weightSum(x, 0, height);
both = max(both, prevScore + prevCur);
}
ll prevWeight = weightSum(x - 1, 0, height);
both = max(both, bestBefore2 + prevWeight);
}
for (const auto& [height, score] : bothDp.at(x)) {
counter++;
ll nextWeight = weightSum(x + 1, 0, height);
bothDp[x][height] += nextWeight;
maxBoth[x] = max(maxBoth[x], bothDp[x][height]);
}
ll& both = bothDp[x][0];
ll& before = beforeDp[x][0];
both = max({ both, maxBefore[x - 1], bestBefore2 });
before = max(before, both);
maxBefore[x] = max(maxBefore[x], before);
maxBoth[x] = max(maxBoth[x], both);
maxBeforeMinCur[x] = max(maxBeforeMinCur[x], before);
maxBeforeMinNext[x] = max(maxBeforeMinNext[x], before);
}
//cout << "Counter: " << counter << endl;
ll res = 0;
loop(x, N) {
for (auto [y, score] : beforeDp[x]) {
res = max(res, score);
}
for (auto [y, score] : bothDp[x]) {
res = max(res, score);
}
}
return res;
}
Compilation message
fish.cpp: In function 'long long int weightSum(int, int, int)':
fish.cpp:38:11: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::vector<std::pair<int, long long int> > >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
38 | if (x >= fishes.size() || y2 <= y1)
| ~~^~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
137 ms |
54852 KB |
Output is correct |
2 |
Correct |
172 ms |
63128 KB |
Output is correct |
3 |
Correct |
27 ms |
39472 KB |
Output is correct |
4 |
Correct |
27 ms |
39476 KB |
Output is correct |
5 |
Correct |
759 ms |
148604 KB |
Output is correct |
6 |
Correct |
431 ms |
151592 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
347 ms |
78888 KB |
Output is correct |
3 |
Correct |
406 ms |
93232 KB |
Output is correct |
4 |
Correct |
149 ms |
54932 KB |
Output is correct |
5 |
Correct |
160 ms |
63176 KB |
Output is correct |
6 |
Correct |
0 ms |
344 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
27 ms |
39368 KB |
Output is correct |
11 |
Correct |
26 ms |
39472 KB |
Output is correct |
12 |
Correct |
180 ms |
60028 KB |
Output is correct |
13 |
Correct |
223 ms |
69524 KB |
Output is correct |
14 |
Correct |
163 ms |
57504 KB |
Output is correct |
15 |
Correct |
185 ms |
63880 KB |
Output is correct |
16 |
Correct |
160 ms |
57472 KB |
Output is correct |
17 |
Correct |
181 ms |
63784 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
25 ms |
39476 KB |
Output is correct |
2 |
Correct |
27 ms |
39408 KB |
Output is correct |
3 |
Correct |
60 ms |
49168 KB |
Output is correct |
4 |
Correct |
50 ms |
48684 KB |
Output is correct |
5 |
Correct |
118 ms |
63784 KB |
Output is correct |
6 |
Correct |
95 ms |
64440 KB |
Output is correct |
7 |
Correct |
94 ms |
64812 KB |
Output is correct |
8 |
Correct |
88 ms |
64812 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 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 |
444 KB |
Output is correct |
10 |
Correct |
2 ms |
1116 KB |
Output is correct |
11 |
Correct |
1 ms |
600 KB |
Output is correct |
12 |
Correct |
2 ms |
856 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
604 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 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 |
444 KB |
Output is correct |
10 |
Correct |
2 ms |
1116 KB |
Output is correct |
11 |
Correct |
1 ms |
600 KB |
Output is correct |
12 |
Correct |
2 ms |
856 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
604 KB |
Output is correct |
15 |
Correct |
1 ms |
604 KB |
Output is correct |
16 |
Correct |
9 ms |
860 KB |
Output is correct |
17 |
Correct |
558 ms |
13364 KB |
Output is correct |
18 |
Correct |
606 ms |
12788 KB |
Output is correct |
19 |
Correct |
336 ms |
13028 KB |
Output is correct |
20 |
Correct |
323 ms |
12116 KB |
Output is correct |
21 |
Correct |
321 ms |
11896 KB |
Output is correct |
22 |
Execution timed out |
1075 ms |
22312 KB |
Time limit exceeded |
23 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 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 |
444 KB |
Output is correct |
10 |
Correct |
2 ms |
1116 KB |
Output is correct |
11 |
Correct |
1 ms |
600 KB |
Output is correct |
12 |
Correct |
2 ms |
856 KB |
Output is correct |
13 |
Correct |
1 ms |
348 KB |
Output is correct |
14 |
Correct |
1 ms |
604 KB |
Output is correct |
15 |
Correct |
1 ms |
604 KB |
Output is correct |
16 |
Correct |
9 ms |
860 KB |
Output is correct |
17 |
Correct |
558 ms |
13364 KB |
Output is correct |
18 |
Correct |
606 ms |
12788 KB |
Output is correct |
19 |
Correct |
336 ms |
13028 KB |
Output is correct |
20 |
Correct |
323 ms |
12116 KB |
Output is correct |
21 |
Correct |
321 ms |
11896 KB |
Output is correct |
22 |
Execution timed out |
1075 ms |
22312 KB |
Time limit exceeded |
23 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
25 ms |
39476 KB |
Output is correct |
2 |
Correct |
27 ms |
39408 KB |
Output is correct |
3 |
Correct |
60 ms |
49168 KB |
Output is correct |
4 |
Correct |
50 ms |
48684 KB |
Output is correct |
5 |
Correct |
118 ms |
63784 KB |
Output is correct |
6 |
Correct |
95 ms |
64440 KB |
Output is correct |
7 |
Correct |
94 ms |
64812 KB |
Output is correct |
8 |
Correct |
88 ms |
64812 KB |
Output is correct |
9 |
Correct |
120 ms |
77356 KB |
Output is correct |
10 |
Correct |
79 ms |
44736 KB |
Output is correct |
11 |
Correct |
163 ms |
90408 KB |
Output is correct |
12 |
Correct |
0 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 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 |
28 ms |
39572 KB |
Output is correct |
19 |
Correct |
26 ms |
39468 KB |
Output is correct |
20 |
Correct |
28 ms |
39476 KB |
Output is correct |
21 |
Correct |
26 ms |
39476 KB |
Output is correct |
22 |
Correct |
124 ms |
77808 KB |
Output is correct |
23 |
Correct |
204 ms |
108400 KB |
Output is correct |
24 |
Correct |
215 ms |
112940 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
137 ms |
54852 KB |
Output is correct |
2 |
Correct |
172 ms |
63128 KB |
Output is correct |
3 |
Correct |
27 ms |
39472 KB |
Output is correct |
4 |
Correct |
27 ms |
39476 KB |
Output is correct |
5 |
Correct |
759 ms |
148604 KB |
Output is correct |
6 |
Correct |
431 ms |
151592 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
347 ms |
78888 KB |
Output is correct |
9 |
Correct |
406 ms |
93232 KB |
Output is correct |
10 |
Correct |
149 ms |
54932 KB |
Output is correct |
11 |
Correct |
160 ms |
63176 KB |
Output is correct |
12 |
Correct |
0 ms |
344 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
0 ms |
348 KB |
Output is correct |
16 |
Correct |
27 ms |
39368 KB |
Output is correct |
17 |
Correct |
26 ms |
39472 KB |
Output is correct |
18 |
Correct |
180 ms |
60028 KB |
Output is correct |
19 |
Correct |
223 ms |
69524 KB |
Output is correct |
20 |
Correct |
163 ms |
57504 KB |
Output is correct |
21 |
Correct |
185 ms |
63880 KB |
Output is correct |
22 |
Correct |
160 ms |
57472 KB |
Output is correct |
23 |
Correct |
181 ms |
63784 KB |
Output is correct |
24 |
Correct |
25 ms |
39476 KB |
Output is correct |
25 |
Correct |
27 ms |
39408 KB |
Output is correct |
26 |
Correct |
60 ms |
49168 KB |
Output is correct |
27 |
Correct |
50 ms |
48684 KB |
Output is correct |
28 |
Correct |
118 ms |
63784 KB |
Output is correct |
29 |
Correct |
95 ms |
64440 KB |
Output is correct |
30 |
Correct |
94 ms |
64812 KB |
Output is correct |
31 |
Correct |
88 ms |
64812 KB |
Output is correct |
32 |
Correct |
0 ms |
348 KB |
Output is correct |
33 |
Correct |
0 ms |
348 KB |
Output is correct |
34 |
Correct |
0 ms |
348 KB |
Output is correct |
35 |
Correct |
0 ms |
348 KB |
Output is correct |
36 |
Correct |
0 ms |
344 KB |
Output is correct |
37 |
Correct |
0 ms |
348 KB |
Output is correct |
38 |
Correct |
0 ms |
348 KB |
Output is correct |
39 |
Correct |
0 ms |
348 KB |
Output is correct |
40 |
Correct |
1 ms |
444 KB |
Output is correct |
41 |
Correct |
2 ms |
1116 KB |
Output is correct |
42 |
Correct |
1 ms |
600 KB |
Output is correct |
43 |
Correct |
2 ms |
856 KB |
Output is correct |
44 |
Correct |
1 ms |
348 KB |
Output is correct |
45 |
Correct |
1 ms |
604 KB |
Output is correct |
46 |
Correct |
1 ms |
604 KB |
Output is correct |
47 |
Correct |
9 ms |
860 KB |
Output is correct |
48 |
Correct |
558 ms |
13364 KB |
Output is correct |
49 |
Correct |
606 ms |
12788 KB |
Output is correct |
50 |
Correct |
336 ms |
13028 KB |
Output is correct |
51 |
Correct |
323 ms |
12116 KB |
Output is correct |
52 |
Correct |
321 ms |
11896 KB |
Output is correct |
53 |
Execution timed out |
1075 ms |
22312 KB |
Time limit exceeded |
54 |
Halted |
0 ms |
0 KB |
- |