#include <bits/stdc++.h>
using namespace std;
const int di[] = {1, -1, 0, 0};
const int dj[] = {0, 0, 1, -1};
vector<vector<long long>> f[1 << 4];
int logs[50005];
int mx[50005][250];
/* int Query(int x1, int y1, int x2, int y2) {
int kx = logs[x2 - x1 + 1];
int ky = logs[y2 - y1 + 1];
return max({st[kx][ky][x1][y1], st[kx][ky][x2 - (1 << kx) + 1][y1], st[kx][ky][x1][y2 - (1 << ky) + 1],
st[kx][ky][x2 - (1 << kx) + 1][y2 - (1 << ky) + 1]});
} */
long long GetSum(int x1, int y1, int x2, int y2, int mask) {
long long res = f[mask][x2][y2];
if (x1 > 0) {
res -= f[mask][x1 - 1][y2];
}
if (y1 > 0) {
res -= f[mask][x2][y1 - 1];
}
if (x1 > 0 && y1 > 0) {
res += f[mask][x1 - 1][y1 - 1];
}
return res;
}
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int n, m;
cin >> n >> m;
vector<vector<int>> a(n, vector<int>(m));
for (int i = 0; i < n; i++) {
for (int j = 0; j < m; j++) {
cin >> a[i][j];
}
}
if (n < m) {
vector<vector<int>> b(m, vector<int>(n));
for (int i = 0; i < n; i++) {
for (int j = 0; j < m; j++) {
b[j][i] = a[i][j];
}
}
a = b;
swap(n, m);
}
for (int mask = 0; mask < (1 << 4); mask++) {
f[mask] = vector<vector<long long>>(n, vector<long long>(m));
}
for (int i = 0; i < n; i++) {
for (int j = 0; j < m; j++) {
for (int mask = 0; mask < (1 << 4); mask++) {
if (i > 0) {
f[mask][i][j] += f[mask][i - 1][j];
}
if (j > 0) {
f[mask][i][j] += f[mask][i][j - 1];
}
if (i > 0 && j > 0) {
f[mask][i][j] -= f[mask][i - 1][j - 1];
}
f[mask][i][j] -= a[i][j];
int mx = 0;
for (int d = 0; d < 4; d++) {
if (mask >> d & 1) {
int ni = i + di[d];
int nj = j + dj[d];
if (0 <= ni && ni < n && 0 <= nj && nj < m && a[ni][nj] < a[i][j]) {
mx = max(mx, a[ni][nj]);
}
}
}
f[mask][i][j] += mx;
}
}
}
int res = 0;
for (int x1 = 0; x1 < n; ++x1) {
for (int y1 = 0; y1 < m; ++y1) {
for (int x2 = x1; x2 < n; ++x2) {
for (int y2 = y1; y2 < m; ++y2) {
mx[x2][y2] = a[x2][y2];
if (x1 == x2 && y1 == y2) {
res += 1;
continue;
}
if (x1 < x2) {
mx[x2][y2] = max(mx[x2][y2], mx[x2 - 1][y2]);
}
if (y1 < y2) {
mx[x2][y2] = max(mx[x2][y2], mx[x2][y2 - 1]);
}
if (x1 < x2 && y1 < y2) {
mx[x2][y2] = max(mx[x2][y2], mx[x2 - 1][y2 - 1]);
}
long long s = 0;
if (x1 == x2) {
s += GetSum(x1, y1, x1, y1, 1 << 2);
s += GetSum(x1, y2, x1, y2, 1 << 3);
if (y1 + 1 < y2) {
s += GetSum(x1, y1 + 1, x1, y2 - 1, (1 << 2) + (1 << 3));
}
} else if (y1 == y2) {
s += GetSum(x1, y1, x1, y1, 1 << 0);
s += GetSum(x2, y1, x2, y1, 1 << 1);
if (x1 + 1 < x2) {
s += GetSum(x1 + 1, y1, x2 - 1, y1, (1 << 0) + (1 << 1));
}
} else {
s += GetSum(x1, y1, x1, y1, (1 << 0) + (1 << 2));
s += GetSum(x1, y2, x1, y2, (1 << 0) + (1 << 3));
s += GetSum(x2, y1, x2, y1, (1 << 1) + (1 << 2));
s += GetSum(x2, y2, x2, y2, (1 << 1) + (1 << 3));
if (x1 + 1 < x2) {
s += GetSum(x1 + 1, y1, x2 - 1, y1, (1 << 0) + (1 << 1) + (1 << 2));
s += GetSum(x1 + 1, y2, x2 - 1, y2, (1 << 0) + (1 << 1) + (1 << 3));
}
if (y1 + 1 < y2) {
s += GetSum(x1, y1 + 1, x1, y2 - 1, (1 << 2) + (1 << 3) + (1 << 0));
s += GetSum(x2, y1 + 1, x2, y2 - 1, (1 << 2) + (1 << 3) + (1 << 1));
}
if (x1 + 1 < x2 && y1 + 1 < y2) {
s += GetSum(x1 + 1, y1 + 1, x2 - 1, y2 - 1, (1 << 0) + (1 << 1) + (1 << 2) + (1 << 3));
}
}
if (s == -mx[x2][y2]) {
res += 1;
}
}
}
}
}
cout << res << '\n';
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
344 KB |
Output is correct |
2 |
Execution timed out |
5083 ms |
96308 KB |
Time limit exceeded |
3 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
460 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
460 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
13 ms |
3160 KB |
Output is correct |
8 |
Correct |
12 ms |
3164 KB |
Output is correct |
9 |
Correct |
13 ms |
600 KB |
Output is correct |
10 |
Correct |
10 ms |
604 KB |
Output is correct |
11 |
Correct |
9 ms |
1904 KB |
Output is correct |
12 |
Correct |
11 ms |
1884 KB |
Output is correct |
13 |
Correct |
10 ms |
724 KB |
Output is correct |
14 |
Correct |
5 ms |
604 KB |
Output is correct |
15 |
Correct |
11 ms |
604 KB |
Output is correct |
16 |
Correct |
10 ms |
748 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
460 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
13 ms |
3160 KB |
Output is correct |
8 |
Correct |
12 ms |
3164 KB |
Output is correct |
9 |
Correct |
13 ms |
600 KB |
Output is correct |
10 |
Correct |
10 ms |
604 KB |
Output is correct |
11 |
Correct |
9 ms |
1904 KB |
Output is correct |
12 |
Correct |
11 ms |
1884 KB |
Output is correct |
13 |
Correct |
10 ms |
724 KB |
Output is correct |
14 |
Correct |
5 ms |
604 KB |
Output is correct |
15 |
Correct |
11 ms |
604 KB |
Output is correct |
16 |
Correct |
10 ms |
748 KB |
Output is correct |
17 |
Correct |
245 ms |
13916 KB |
Output is correct |
18 |
Correct |
219 ms |
1604 KB |
Output is correct |
19 |
Correct |
221 ms |
2136 KB |
Output is correct |
20 |
Correct |
223 ms |
1616 KB |
Output is correct |
21 |
Correct |
210 ms |
1368 KB |
Output is correct |
22 |
Correct |
211 ms |
1368 KB |
Output is correct |
23 |
Correct |
209 ms |
1372 KB |
Output is correct |
24 |
Correct |
170 ms |
1372 KB |
Output is correct |
25 |
Correct |
219 ms |
1552 KB |
Output is correct |
26 |
Correct |
225 ms |
1372 KB |
Output is correct |
27 |
Correct |
219 ms |
1368 KB |
Output is correct |
28 |
Correct |
221 ms |
1372 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
460 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
13 ms |
3160 KB |
Output is correct |
8 |
Correct |
12 ms |
3164 KB |
Output is correct |
9 |
Correct |
13 ms |
600 KB |
Output is correct |
10 |
Correct |
10 ms |
604 KB |
Output is correct |
11 |
Correct |
9 ms |
1904 KB |
Output is correct |
12 |
Correct |
11 ms |
1884 KB |
Output is correct |
13 |
Correct |
10 ms |
724 KB |
Output is correct |
14 |
Correct |
5 ms |
604 KB |
Output is correct |
15 |
Correct |
11 ms |
604 KB |
Output is correct |
16 |
Correct |
10 ms |
748 KB |
Output is correct |
17 |
Correct |
245 ms |
13916 KB |
Output is correct |
18 |
Correct |
219 ms |
1604 KB |
Output is correct |
19 |
Correct |
221 ms |
2136 KB |
Output is correct |
20 |
Correct |
223 ms |
1616 KB |
Output is correct |
21 |
Correct |
210 ms |
1368 KB |
Output is correct |
22 |
Correct |
211 ms |
1368 KB |
Output is correct |
23 |
Correct |
209 ms |
1372 KB |
Output is correct |
24 |
Correct |
170 ms |
1372 KB |
Output is correct |
25 |
Correct |
219 ms |
1552 KB |
Output is correct |
26 |
Correct |
225 ms |
1372 KB |
Output is correct |
27 |
Correct |
219 ms |
1368 KB |
Output is correct |
28 |
Correct |
221 ms |
1372 KB |
Output is correct |
29 |
Execution timed out |
5101 ms |
96200 KB |
Time limit exceeded |
30 |
Halted |
0 ms |
0 KB |
- |