#include <iostream>
#include <complex>
#include <iomanip>
#include <vector>
#include <algorithm>
#include <functional>
#include <bitset>
#include <queue>
#include <map>
#include <stack>
#include <cmath>
#include <cstdint>
#include <cassert>
using namespace std;
const int Max = 3001, Inf = 1e9, TLE = 5e6;
#define ll unsigned long long
int lg[Max], cnt;
struct SPTMIN
{
vector <vector <int>> spt;
int query(int i, int j){
int x = lg[j-i+1];
return min(spt[i][x], spt[j - (1LL << x) + 1][x]);
}
SPTMIN(vector <int> v)
{
int n = v.size();
spt.assign(n + 1, vector <int> (lg[n] + 1, Inf));
for(int j = 0; j < lg[n]+1; j++)
{
for(int i = 0; i < n; i++){
if(j == 0){
spt[i][j] = v[i];
continue;
}
spt[i][j] = spt[i][j-1];
if(i + (1LL << (j-1)) < n) spt[i][j] = min(spt[i][j-1], spt[i + (1LL << (j-1))][j-1]);
}
}
}
};
struct SPTMAX
{
vector <vector <int>> spt;
int query(int i, int j){
int x = lg[j-i+1];
return max(spt[i][x], spt[j - (1LL << x) + 1][x]);
}
SPTMAX(vector <int> v)
{
int n = v.size();
spt.assign(n + 1, vector <int> (lg[n] + 1, 0));
for(int j = 0; j < lg[n]+1; j++)
{
for(int i = 0; i < n; i++){
if(j == 0){
spt[i][j] = v[i];
continue;
}
spt[i][j] = spt[i][j-1];
if(i + (1LL << (j-1)) < n) spt[i][j] = max(spt[i][j-1], spt[i + (1LL << (j-1))][j-1]);
}
}
}
};
int rnd(){
return abs(int(rand() ^ (rand() << 3LL) ^ (rand() << 6LL)));
}
int biggest_stadium(int n, std::vector<std::vector<int>> f)
{
vector <vector<int>> A(n, vector <int> (n, 0)),
B(n, vector <int> (n, 0)),
C(n, vector <int> (n, 0)),
D(n, vector <int> (n, 0)),
T(n, vector <int> (n, 0));
vector <SPTMIN> sp1;
vector <SPTMAX> sp2;
cnt = 0;
for(int i = 1; i <= n; i++){
lg[i] = log2(i);
}
for(int i = 0; i < n; i++)
{
int last = -1;
for(int j = 0; j < n; j++){
if(f[i][j] == 1) last = j;
A[i][j] = last;
}
last = n;
for(int j = n-1; j >= 0; j--){
T[i][j] = last;
if(f[i][j] == 1) {
if(j != n-1 && last == j+1){
T[i][j] = T[i][j+1];
}else T[i][j] = j;
last = j;
}
B[i][j] = last;
}
}
for(int i = 0; i < n; i++)
{
int last = -1;
for(int j = 0; j < n; j++){
if(f[j][i] == 1) last = j;
C[j][i] = last+1;
}
last = n;
for(int j = n-1; j >= 0; j--){
if(f[j][i] == 1) last = j;
D[j][i] = last-1;
}
}
for(int i = 0; i < n; i++)
{
sp1.push_back(SPTMIN(D[i]));
sp2.push_back(SPTMAX(C[i]));
}
//return 0;
//vector <vector<vector<int>>> dp1(n + 1, vector <vector<int>>
// (n + 1, vector <int> (n + 1, -1)));
//vector <vector<vector<int>>> dp2(n + 1, vector <vector<int>>
// (n + 1, vector <int> (n + 1, -1)));
unordered_map <ll, int> dp;
auto F = [&](auto F, int i, int j, int a, int b, int c, int d) -> int
{ cnt++;
long long x = ((ll) b) + (((ll) c) << 11LL) + (((ll) d) << 22LL) + (((ll) j) << 33LL)
+ (((ll) i) << 44LL) + (((ll) a) << 55LL);
if(dp[x] != 0) return dp[x];
dp[x] = (b - a + 1);
if(i != j) dp[x] += (d - c + 1);
int mx = 0;
if(i != 0)
{
for(int k = a; k <= b; k++) if(f[i-1][k] == 0)
{
int l1 = max(A[i-1][k]+1, a), r1 = min(B[i-1][k]-1, b);
int l2 = max(c, l1), r2 = min(d, r1);
int ni = sp2[i].query(l1, r1), d = abs(ni - i) - 1, t = (r1 - l1 + 1);
mx = max(mx, F(F, ni, j, l1, r1, l2, r2) - (r2 - l2 + 1) + d * t);
k = r1;
} else k = T[i-1][k];
}
if(j != n-1)
{
for(int k = c; k <= d; k++) if(f[j+1][k] == 0)
{
int l2 = max(A[j+1][k]+1, c), r2 = min(B[j+1][k]-1, d);
int l1 = max(a, l2), r1 = min(b, r2);
int nj = sp1[j].query(l2, r2), d = abs(nj - j) - 1, t = (r2 - l2 + 1);
mx = max(mx, F(F, i, nj, l1, r1, l2, r2) - (r1 - l1 + 1) + d * t);
k = r2;
} else k = T[j+1][k];
}
return dp[x] = (dp[x] + mx);
};
int ans = 0;
srand(time(NULL));
vector <pair<int, pair<int, int>>> fx;
for(int i = 0; i < n; i++)
{
for(int j = 0; j < n; j++)
{
if(f[i][j] == 1) continue;
int l = A[i][j]+1, r = B[i][j]-1;
int ni = sp2[i].query(j, j);
int nj = sp1[i].query(j, j);
int d = r-l+1, h = (nj - ni + 1);
fx.push_back({ d+ h, { i, j }});
}
}
sort(fx.begin(), fx.end());
while (fx.size() > 0)
{
pair<int, int> u = fx.back().second; fx.pop_back();
int i = u.first, j = u.second;
if(f[i][j] == 1) continue;
int l = A[i][j]+1, r = B[i][j]-1;
ans = max(ans, F(F, i, i, l, r, l, r));
if(n > 500 && cnt > TLE) {
break;
}
}
return ans;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
1 ms |
348 KB |
ok |
4 |
Correct |
1 ms |
348 KB |
ok |
5 |
Correct |
1 ms |
348 KB |
ok |
6 |
Correct |
1 ms |
348 KB |
ok |
7 |
Correct |
3 ms |
2396 KB |
ok |
8 |
Correct |
77 ms |
46708 KB |
ok |
9 |
Correct |
1387 ms |
850572 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
0 ms |
348 KB |
ok |
4 |
Correct |
0 ms |
348 KB |
ok |
5 |
Correct |
1 ms |
348 KB |
ok |
6 |
Correct |
1 ms |
348 KB |
ok |
7 |
Correct |
0 ms |
348 KB |
ok |
8 |
Correct |
0 ms |
348 KB |
ok |
9 |
Correct |
1 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
1 ms |
344 KB |
ok |
12 |
Correct |
1 ms |
348 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
0 ms |
348 KB |
ok |
4 |
Correct |
0 ms |
348 KB |
ok |
5 |
Correct |
0 ms |
348 KB |
ok |
6 |
Correct |
1 ms |
348 KB |
ok |
7 |
Correct |
1 ms |
348 KB |
ok |
8 |
Correct |
0 ms |
348 KB |
ok |
9 |
Correct |
0 ms |
348 KB |
ok |
10 |
Correct |
1 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
1 ms |
344 KB |
ok |
13 |
Correct |
1 ms |
348 KB |
ok |
14 |
Correct |
0 ms |
348 KB |
ok |
15 |
Correct |
1 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
348 KB |
ok |
17 |
Correct |
0 ms |
348 KB |
ok |
18 |
Correct |
1 ms |
348 KB |
ok |
19 |
Correct |
0 ms |
348 KB |
ok |
20 |
Correct |
0 ms |
348 KB |
ok |
21 |
Correct |
0 ms |
348 KB |
ok |
22 |
Correct |
1 ms |
348 KB |
ok |
23 |
Correct |
0 ms |
348 KB |
ok |
24 |
Correct |
0 ms |
348 KB |
ok |
25 |
Correct |
0 ms |
348 KB |
ok |
26 |
Correct |
0 ms |
348 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
0 ms |
348 KB |
ok |
4 |
Correct |
1 ms |
348 KB |
ok |
5 |
Correct |
1 ms |
348 KB |
ok |
6 |
Correct |
0 ms |
348 KB |
ok |
7 |
Correct |
0 ms |
348 KB |
ok |
8 |
Correct |
1 ms |
348 KB |
ok |
9 |
Correct |
1 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
1 ms |
348 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
14 |
Correct |
1 ms |
344 KB |
ok |
15 |
Correct |
1 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
348 KB |
ok |
17 |
Correct |
1 ms |
348 KB |
ok |
18 |
Correct |
0 ms |
348 KB |
ok |
19 |
Correct |
0 ms |
348 KB |
ok |
20 |
Correct |
1 ms |
348 KB |
ok |
21 |
Correct |
0 ms |
348 KB |
ok |
22 |
Correct |
0 ms |
348 KB |
ok |
23 |
Correct |
0 ms |
348 KB |
ok |
24 |
Correct |
1 ms |
348 KB |
ok |
25 |
Correct |
0 ms |
348 KB |
ok |
26 |
Correct |
0 ms |
348 KB |
ok |
27 |
Correct |
0 ms |
348 KB |
ok |
28 |
Correct |
0 ms |
348 KB |
ok |
29 |
Correct |
1 ms |
348 KB |
ok |
30 |
Correct |
1 ms |
604 KB |
ok |
31 |
Correct |
1 ms |
468 KB |
ok |
32 |
Correct |
1 ms |
348 KB |
ok |
33 |
Correct |
0 ms |
348 KB |
ok |
34 |
Correct |
0 ms |
348 KB |
ok |
35 |
Correct |
1 ms |
348 KB |
ok |
36 |
Correct |
1 ms |
348 KB |
ok |
37 |
Correct |
0 ms |
348 KB |
ok |
38 |
Correct |
1 ms |
344 KB |
ok |
39 |
Correct |
1 ms |
348 KB |
ok |
40 |
Correct |
1 ms |
344 KB |
ok |
41 |
Correct |
1 ms |
604 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
0 ms |
348 KB |
ok |
4 |
Correct |
1 ms |
348 KB |
ok |
5 |
Correct |
1 ms |
348 KB |
ok |
6 |
Correct |
0 ms |
348 KB |
ok |
7 |
Correct |
0 ms |
348 KB |
ok |
8 |
Correct |
1 ms |
348 KB |
ok |
9 |
Correct |
1 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
1 ms |
348 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
14 |
Correct |
1 ms |
344 KB |
ok |
15 |
Correct |
1 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
348 KB |
ok |
17 |
Correct |
1 ms |
348 KB |
ok |
18 |
Correct |
0 ms |
348 KB |
ok |
19 |
Correct |
0 ms |
348 KB |
ok |
20 |
Correct |
1 ms |
348 KB |
ok |
21 |
Correct |
0 ms |
348 KB |
ok |
22 |
Correct |
0 ms |
348 KB |
ok |
23 |
Correct |
0 ms |
348 KB |
ok |
24 |
Correct |
1 ms |
348 KB |
ok |
25 |
Correct |
0 ms |
348 KB |
ok |
26 |
Correct |
0 ms |
348 KB |
ok |
27 |
Correct |
0 ms |
348 KB |
ok |
28 |
Correct |
0 ms |
348 KB |
ok |
29 |
Correct |
1 ms |
348 KB |
ok |
30 |
Correct |
1 ms |
604 KB |
ok |
31 |
Correct |
1 ms |
468 KB |
ok |
32 |
Correct |
1 ms |
348 KB |
ok |
33 |
Correct |
0 ms |
348 KB |
ok |
34 |
Correct |
0 ms |
348 KB |
ok |
35 |
Correct |
1 ms |
348 KB |
ok |
36 |
Correct |
1 ms |
348 KB |
ok |
37 |
Correct |
0 ms |
348 KB |
ok |
38 |
Correct |
1 ms |
344 KB |
ok |
39 |
Correct |
1 ms |
348 KB |
ok |
40 |
Correct |
1 ms |
344 KB |
ok |
41 |
Correct |
1 ms |
604 KB |
ok |
42 |
Correct |
298 ms |
58900 KB |
ok |
43 |
Correct |
268 ms |
56396 KB |
ok |
44 |
Correct |
366 ms |
64960 KB |
ok |
45 |
Correct |
308 ms |
64796 KB |
ok |
46 |
Correct |
320 ms |
64700 KB |
ok |
47 |
Correct |
81 ms |
46280 KB |
ok |
48 |
Correct |
110 ms |
48708 KB |
ok |
49 |
Correct |
100 ms |
47208 KB |
ok |
50 |
Correct |
143 ms |
57244 KB |
ok |
51 |
Correct |
156 ms |
54976 KB |
ok |
52 |
Correct |
66 ms |
44492 KB |
ok |
53 |
Correct |
57 ms |
43316 KB |
ok |
54 |
Correct |
82 ms |
44388 KB |
ok |
55 |
Correct |
82 ms |
45772 KB |
ok |
56 |
Correct |
67 ms |
44492 KB |
ok |
57 |
Correct |
136 ms |
55232 KB |
ok |
58 |
Correct |
144 ms |
55244 KB |
ok |
59 |
Correct |
144 ms |
55236 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
0 ms |
348 KB |
ok |
4 |
Correct |
1 ms |
348 KB |
ok |
5 |
Correct |
1 ms |
348 KB |
ok |
6 |
Correct |
1 ms |
348 KB |
ok |
7 |
Correct |
1 ms |
348 KB |
ok |
8 |
Correct |
3 ms |
2396 KB |
ok |
9 |
Correct |
77 ms |
46708 KB |
ok |
10 |
Correct |
1387 ms |
850572 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
0 ms |
348 KB |
ok |
13 |
Correct |
1 ms |
348 KB |
ok |
14 |
Correct |
1 ms |
348 KB |
ok |
15 |
Correct |
0 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
348 KB |
ok |
17 |
Correct |
1 ms |
348 KB |
ok |
18 |
Correct |
0 ms |
348 KB |
ok |
19 |
Correct |
1 ms |
344 KB |
ok |
20 |
Correct |
1 ms |
348 KB |
ok |
21 |
Correct |
0 ms |
348 KB |
ok |
22 |
Correct |
1 ms |
348 KB |
ok |
23 |
Correct |
0 ms |
348 KB |
ok |
24 |
Correct |
0 ms |
348 KB |
ok |
25 |
Correct |
1 ms |
348 KB |
ok |
26 |
Correct |
0 ms |
348 KB |
ok |
27 |
Correct |
0 ms |
348 KB |
ok |
28 |
Correct |
0 ms |
348 KB |
ok |
29 |
Correct |
1 ms |
348 KB |
ok |
30 |
Correct |
0 ms |
348 KB |
ok |
31 |
Correct |
0 ms |
348 KB |
ok |
32 |
Correct |
0 ms |
348 KB |
ok |
33 |
Correct |
0 ms |
348 KB |
ok |
34 |
Correct |
1 ms |
348 KB |
ok |
35 |
Correct |
1 ms |
604 KB |
ok |
36 |
Correct |
1 ms |
468 KB |
ok |
37 |
Correct |
1 ms |
348 KB |
ok |
38 |
Correct |
0 ms |
348 KB |
ok |
39 |
Correct |
0 ms |
348 KB |
ok |
40 |
Correct |
1 ms |
348 KB |
ok |
41 |
Correct |
1 ms |
348 KB |
ok |
42 |
Correct |
0 ms |
348 KB |
ok |
43 |
Correct |
1 ms |
344 KB |
ok |
44 |
Correct |
1 ms |
348 KB |
ok |
45 |
Correct |
1 ms |
344 KB |
ok |
46 |
Correct |
1 ms |
604 KB |
ok |
47 |
Correct |
298 ms |
58900 KB |
ok |
48 |
Correct |
268 ms |
56396 KB |
ok |
49 |
Correct |
366 ms |
64960 KB |
ok |
50 |
Correct |
308 ms |
64796 KB |
ok |
51 |
Correct |
320 ms |
64700 KB |
ok |
52 |
Correct |
81 ms |
46280 KB |
ok |
53 |
Correct |
110 ms |
48708 KB |
ok |
54 |
Correct |
100 ms |
47208 KB |
ok |
55 |
Correct |
143 ms |
57244 KB |
ok |
56 |
Correct |
156 ms |
54976 KB |
ok |
57 |
Correct |
66 ms |
44492 KB |
ok |
58 |
Correct |
57 ms |
43316 KB |
ok |
59 |
Correct |
82 ms |
44388 KB |
ok |
60 |
Correct |
82 ms |
45772 KB |
ok |
61 |
Correct |
67 ms |
44492 KB |
ok |
62 |
Correct |
136 ms |
55232 KB |
ok |
63 |
Correct |
144 ms |
55244 KB |
ok |
64 |
Correct |
144 ms |
55236 KB |
ok |
65 |
Correct |
4242 ms |
947868 KB |
ok |
66 |
Correct |
2123 ms |
870208 KB |
ok |
67 |
Correct |
1353 ms |
828628 KB |
ok |
68 |
Correct |
3746 ms |
950368 KB |
ok |
69 |
Correct |
3894 ms |
954528 KB |
ok |
70 |
Correct |
3978 ms |
952824 KB |
ok |
71 |
Partially correct |
3683 ms |
946236 KB |
partial |
72 |
Correct |
1677 ms |
861056 KB |
ok |
73 |
Correct |
2695 ms |
883880 KB |
ok |
74 |
Correct |
2959 ms |
907848 KB |
ok |
75 |
Correct |
2662 ms |
883580 KB |
ok |
76 |
Correct |
1782 ms |
926740 KB |
ok |
77 |
Correct |
1929 ms |
926180 KB |
ok |
78 |
Correct |
2686 ms |
927840 KB |
ok |
79 |
Correct |
1268 ms |
816572 KB |
ok |
80 |
Correct |
1354 ms |
816836 KB |
ok |
81 |
Correct |
1589 ms |
830140 KB |
ok |
82 |
Correct |
2533 ms |
858540 KB |
ok |
83 |
Correct |
2283 ms |
848572 KB |
ok |
84 |
Correct |
1288 ms |
814524 KB |
ok |
85 |
Correct |
1299 ms |
807300 KB |
ok |
86 |
Correct |
1486 ms |
825732 KB |
ok |
87 |
Correct |
1381 ms |
828280 KB |
ok |
88 |
Correct |
1477 ms |
832468 KB |
ok |
89 |
Correct |
1362 ms |
849592 KB |
ok |
90 |
Correct |
1244 ms |
851000 KB |
ok |
91 |
Correct |
1396 ms |
850452 KB |
ok |
92 |
Correct |
1190 ms |
817560 KB |
ok |
93 |
Correct |
1394 ms |
825012 KB |
ok |
94 |
Correct |
1071 ms |
809348 KB |
ok |
95 |
Correct |
1002 ms |
807624 KB |
ok |
96 |
Correct |
956 ms |
807048 KB |
ok |
97 |
Correct |
974 ms |
807564 KB |
ok |
98 |
Correct |
914 ms |
803020 KB |
ok |
99 |
Correct |
2595 ms |
942748 KB |
ok |
100 |
Correct |
1920 ms |
932664 KB |
ok |
101 |
Correct |
1917 ms |
931484 KB |
ok |
102 |
Correct |
1927 ms |
931940 KB |
ok |
103 |
Correct |
1947 ms |
929704 KB |
ok |
104 |
Correct |
1949 ms |
928924 KB |
ok |
105 |
Correct |
1878 ms |
927404 KB |
ok |
106 |
Correct |
1927 ms |
932268 KB |
ok |
107 |
Correct |
2295 ms |
939764 KB |
ok |
108 |
Correct |
2607 ms |
920592 KB |
ok |
109 |
Correct |
2758 ms |
925092 KB |
ok |