#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 = 7e6;
#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(l, r);
int nj = sp1[i].query(l, r);
int d = r-l+1, h = (nj - ni + 1);
fx.push_back({ d * h, { i, j }});
}
}
sort(fx.rbegin(), fx.rend());
for(auto& u : fx)
{
int i = u.second.first, j = u.second.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) {
return ans;
}
}
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 |
0 ms |
348 KB |
ok |
4 |
Correct |
0 ms |
348 KB |
ok |
5 |
Correct |
0 ms |
348 KB |
ok |
6 |
Correct |
0 ms |
348 KB |
ok |
7 |
Correct |
3 ms |
2396 KB |
ok |
8 |
Correct |
94 ms |
46796 KB |
ok |
9 |
Correct |
1607 ms |
850552 KB |
ok |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
348 KB |
ok |
2 |
Correct |
0 ms |
348 KB |
ok |
3 |
Correct |
1 ms |
344 KB |
ok |
4 |
Correct |
0 ms |
344 KB |
ok |
5 |
Correct |
0 ms |
348 KB |
ok |
6 |
Correct |
0 ms |
348 KB |
ok |
7 |
Correct |
0 ms |
348 KB |
ok |
8 |
Correct |
0 ms |
348 KB |
ok |
9 |
Correct |
0 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
344 KB |
ok |
12 |
Correct |
0 ms |
348 KB |
ok |
13 |
Correct |
0 ms |
344 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 |
344 KB |
ok |
5 |
Correct |
0 ms |
344 KB |
ok |
6 |
Correct |
0 ms |
348 KB |
ok |
7 |
Correct |
0 ms |
348 KB |
ok |
8 |
Correct |
0 ms |
348 KB |
ok |
9 |
Correct |
0 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
0 ms |
344 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
14 |
Correct |
0 ms |
344 KB |
ok |
15 |
Correct |
0 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
348 KB |
ok |
17 |
Correct |
0 ms |
348 KB |
ok |
18 |
Correct |
1 ms |
344 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 |
0 ms |
348 KB |
ok |
23 |
Correct |
1 ms |
600 KB |
ok |
24 |
Correct |
0 ms |
348 KB |
ok |
25 |
Correct |
0 ms |
348 KB |
ok |
26 |
Correct |
1 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 |
344 KB |
ok |
7 |
Correct |
0 ms |
344 KB |
ok |
8 |
Correct |
0 ms |
348 KB |
ok |
9 |
Correct |
0 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
0 ms |
348 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
14 |
Correct |
0 ms |
344 KB |
ok |
15 |
Correct |
0 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
344 KB |
ok |
17 |
Correct |
0 ms |
348 KB |
ok |
18 |
Correct |
0 ms |
348 KB |
ok |
19 |
Correct |
0 ms |
348 KB |
ok |
20 |
Correct |
1 ms |
344 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 |
0 ms |
348 KB |
ok |
25 |
Correct |
1 ms |
600 KB |
ok |
26 |
Correct |
0 ms |
348 KB |
ok |
27 |
Correct |
0 ms |
348 KB |
ok |
28 |
Correct |
1 ms |
348 KB |
ok |
29 |
Correct |
0 ms |
348 KB |
ok |
30 |
Correct |
1 ms |
604 KB |
ok |
31 |
Correct |
1 ms |
604 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 |
0 ms |
348 KB |
ok |
37 |
Correct |
0 ms |
348 KB |
ok |
38 |
Correct |
0 ms |
348 KB |
ok |
39 |
Correct |
1 ms |
348 KB |
ok |
40 |
Correct |
0 ms |
348 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 |
0 ms |
348 KB |
ok |
5 |
Correct |
0 ms |
348 KB |
ok |
6 |
Correct |
1 ms |
344 KB |
ok |
7 |
Correct |
0 ms |
344 KB |
ok |
8 |
Correct |
0 ms |
348 KB |
ok |
9 |
Correct |
0 ms |
348 KB |
ok |
10 |
Correct |
0 ms |
348 KB |
ok |
11 |
Correct |
0 ms |
348 KB |
ok |
12 |
Correct |
0 ms |
348 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
14 |
Correct |
0 ms |
344 KB |
ok |
15 |
Correct |
0 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
344 KB |
ok |
17 |
Correct |
0 ms |
348 KB |
ok |
18 |
Correct |
0 ms |
348 KB |
ok |
19 |
Correct |
0 ms |
348 KB |
ok |
20 |
Correct |
1 ms |
344 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 |
0 ms |
348 KB |
ok |
25 |
Correct |
1 ms |
600 KB |
ok |
26 |
Correct |
0 ms |
348 KB |
ok |
27 |
Correct |
0 ms |
348 KB |
ok |
28 |
Correct |
1 ms |
348 KB |
ok |
29 |
Correct |
0 ms |
348 KB |
ok |
30 |
Correct |
1 ms |
604 KB |
ok |
31 |
Correct |
1 ms |
604 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 |
0 ms |
348 KB |
ok |
37 |
Correct |
0 ms |
348 KB |
ok |
38 |
Correct |
0 ms |
348 KB |
ok |
39 |
Correct |
1 ms |
348 KB |
ok |
40 |
Correct |
0 ms |
348 KB |
ok |
41 |
Correct |
1 ms |
604 KB |
ok |
42 |
Correct |
243 ms |
58644 KB |
ok |
43 |
Correct |
196 ms |
56104 KB |
ok |
44 |
Correct |
249 ms |
64780 KB |
ok |
45 |
Correct |
236 ms |
64964 KB |
ok |
46 |
Correct |
254 ms |
64704 KB |
ok |
47 |
Correct |
76 ms |
46268 KB |
ok |
48 |
Correct |
90 ms |
48844 KB |
ok |
49 |
Correct |
84 ms |
47304 KB |
ok |
50 |
Correct |
119 ms |
57440 KB |
ok |
51 |
Correct |
127 ms |
55000 KB |
ok |
52 |
Correct |
62 ms |
44496 KB |
ok |
53 |
Correct |
56 ms |
43156 KB |
ok |
54 |
Correct |
60 ms |
44368 KB |
ok |
55 |
Correct |
74 ms |
45772 KB |
ok |
56 |
Correct |
65 ms |
44556 KB |
ok |
57 |
Correct |
118 ms |
55204 KB |
ok |
58 |
Correct |
124 ms |
55236 KB |
ok |
59 |
Correct |
130 ms |
55112 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 |
0 ms |
348 KB |
ok |
7 |
Correct |
0 ms |
348 KB |
ok |
8 |
Correct |
3 ms |
2396 KB |
ok |
9 |
Correct |
94 ms |
46796 KB |
ok |
10 |
Correct |
1607 ms |
850552 KB |
ok |
11 |
Correct |
1 ms |
344 KB |
ok |
12 |
Correct |
0 ms |
344 KB |
ok |
13 |
Correct |
0 ms |
348 KB |
ok |
14 |
Correct |
0 ms |
348 KB |
ok |
15 |
Correct |
0 ms |
348 KB |
ok |
16 |
Correct |
0 ms |
348 KB |
ok |
17 |
Correct |
0 ms |
348 KB |
ok |
18 |
Correct |
0 ms |
348 KB |
ok |
19 |
Correct |
0 ms |
344 KB |
ok |
20 |
Correct |
0 ms |
348 KB |
ok |
21 |
Correct |
0 ms |
344 KB |
ok |
22 |
Correct |
0 ms |
348 KB |
ok |
23 |
Correct |
0 ms |
348 KB |
ok |
24 |
Correct |
0 ms |
348 KB |
ok |
25 |
Correct |
1 ms |
344 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 |
0 ms |
348 KB |
ok |
30 |
Correct |
1 ms |
600 KB |
ok |
31 |
Correct |
0 ms |
348 KB |
ok |
32 |
Correct |
0 ms |
348 KB |
ok |
33 |
Correct |
1 ms |
348 KB |
ok |
34 |
Correct |
0 ms |
348 KB |
ok |
35 |
Correct |
1 ms |
604 KB |
ok |
36 |
Correct |
1 ms |
604 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 |
0 ms |
348 KB |
ok |
42 |
Correct |
0 ms |
348 KB |
ok |
43 |
Correct |
0 ms |
348 KB |
ok |
44 |
Correct |
1 ms |
348 KB |
ok |
45 |
Correct |
0 ms |
348 KB |
ok |
46 |
Correct |
1 ms |
604 KB |
ok |
47 |
Correct |
243 ms |
58644 KB |
ok |
48 |
Correct |
196 ms |
56104 KB |
ok |
49 |
Correct |
249 ms |
64780 KB |
ok |
50 |
Correct |
236 ms |
64964 KB |
ok |
51 |
Correct |
254 ms |
64704 KB |
ok |
52 |
Correct |
76 ms |
46268 KB |
ok |
53 |
Correct |
90 ms |
48844 KB |
ok |
54 |
Correct |
84 ms |
47304 KB |
ok |
55 |
Correct |
119 ms |
57440 KB |
ok |
56 |
Correct |
127 ms |
55000 KB |
ok |
57 |
Correct |
62 ms |
44496 KB |
ok |
58 |
Correct |
56 ms |
43156 KB |
ok |
59 |
Correct |
60 ms |
44368 KB |
ok |
60 |
Correct |
74 ms |
45772 KB |
ok |
61 |
Correct |
65 ms |
44556 KB |
ok |
62 |
Correct |
118 ms |
55204 KB |
ok |
63 |
Correct |
124 ms |
55236 KB |
ok |
64 |
Correct |
130 ms |
55112 KB |
ok |
65 |
Correct |
4358 ms |
1005180 KB |
ok |
66 |
Correct |
1891 ms |
866984 KB |
ok |
67 |
Correct |
1184 ms |
824760 KB |
ok |
68 |
Correct |
3965 ms |
1010192 KB |
ok |
69 |
Correct |
4273 ms |
1008800 KB |
ok |
70 |
Correct |
4146 ms |
1008800 KB |
ok |
71 |
Correct |
3349 ms |
954548 KB |
ok |
72 |
Correct |
1358 ms |
857244 KB |
ok |
73 |
Correct |
1999 ms |
880260 KB |
ok |
74 |
Correct |
2185 ms |
903596 KB |
ok |
75 |
Correct |
2020 ms |
879528 KB |
ok |
76 |
Correct |
1891 ms |
997024 KB |
ok |
77 |
Correct |
1931 ms |
998560 KB |
ok |
78 |
Partially correct |
2255 ms |
940092 KB |
partial |
79 |
Correct |
1144 ms |
812360 KB |
ok |
80 |
Correct |
1087 ms |
812476 KB |
ok |
81 |
Correct |
1273 ms |
826036 KB |
ok |
82 |
Correct |
1868 ms |
858652 KB |
ok |
83 |
Correct |
1617 ms |
849328 KB |
ok |
84 |
Correct |
971 ms |
814464 KB |
ok |
85 |
Correct |
927 ms |
807320 KB |
ok |
86 |
Correct |
1052 ms |
825780 KB |
ok |
87 |
Correct |
1086 ms |
828336 KB |
ok |
88 |
Correct |
1256 ms |
831660 KB |
ok |
89 |
Correct |
1220 ms |
849836 KB |
ok |
90 |
Correct |
1130 ms |
851076 KB |
ok |
91 |
Correct |
1197 ms |
849324 KB |
ok |
92 |
Correct |
1145 ms |
817532 KB |
ok |
93 |
Correct |
1228 ms |
824192 KB |
ok |
94 |
Correct |
998 ms |
808840 KB |
ok |
95 |
Correct |
919 ms |
807836 KB |
ok |
96 |
Correct |
918 ms |
808136 KB |
ok |
97 |
Correct |
965 ms |
807564 KB |
ok |
98 |
Correct |
917 ms |
803528 KB |
ok |
99 |
Correct |
3128 ms |
946708 KB |
ok |
100 |
Partially correct |
1985 ms |
936276 KB |
partial |
101 |
Partially correct |
1856 ms |
932260 KB |
partial |
102 |
Partially correct |
1931 ms |
935780 KB |
partial |
103 |
Partially correct |
1937 ms |
940792 KB |
partial |
104 |
Partially correct |
1990 ms |
942756 KB |
partial |
105 |
Partially correct |
2068 ms |
943260 KB |
partial |
106 |
Partially correct |
2145 ms |
945764 KB |
partial |
107 |
Partially correct |
2173 ms |
945276 KB |
partial |
108 |
Correct |
2777 ms |
939980 KB |
ok |
109 |
Correct |
2829 ms |
947344 KB |
ok |