#include "rect.h"
#include <bits/stdc++.h>
using namespace std;
#pragma GCC optimize("O3")
const int N = 2510;
int n,m;
vector<pair<int,int>> gdr[N];
vector<pair<int,int>> gdc[N];
struct Side {
int r1, r2, c1, c2;
};
vector<pair<int,int>> hsides[N][N];
vector<pair<int,int>> vsides[N][N];
vector<pair<int,int>> swp[N];
int fen[N];
void upd(int at, int add) {
at++;
while(at) {
fen[at]+=add;
at-=at&(-at);
}
}
int qry(int at) {
at++;
int sum=0;
while(at<N) {
sum+=fen[at];
at+=at&(-at);
}
return sum;
}
long long count_rectangles(std::vector<std::vector<int> > a) {
n = a.size();
m = a[0].size();
for(int i = 0;i<n;i++) {
vector<int> stk;
for(int j = 0;j<m;j++) {
while(stk.size() and a[i][stk.back()] < a[i][j]){
stk.pop_back();
}
if(stk.size()) {
gdr[i].push_back({stk.back(), j});
}
stk.push_back(j);
}
stk.clear();
for(int j = m-1;j>=0;j--) {
while(stk.size() and a[i][stk.back()] < a[i][j])stk.pop_back();
if(stk.size()) {
if(a[i][j] != a[i][stk.back()])
gdr[i].push_back({j, stk.back()});
}
stk.push_back(j);
}
sort(gdr[i].begin(), gdr[i].end());
}
for(int i = 0;i<m;i++) {
vector<int> stk;
for(int j = 0;j<n;j++) {
while(stk.size() and a[stk.back()][i] < a[j][i]){
stk.pop_back();
}
if(stk.size()) {
gdc[i].push_back({stk.back(), j});
}
stk.push_back(j);
}
stk.clear();
for(int j = n-1;j>=0;j--) {
while(stk.size() and a[stk.back()][i] < a[j][i])stk.pop_back();
if(stk.size()) {
gdc[i].push_back({j, stk.back()});
}
stk.push_back(j);
}
sort(gdc[i].begin(), gdc[i].end());
gdc[i].erase(unique(gdc[i].begin(), gdc[i].end()), gdc[i].end());
}
Side seg;
for(int i = 0;i<n;i++) {
for(auto pr:gdr[i]) {
if(pr.first +1 == pr.second)continue;
if(i!=n-1 and (*lower_bound(gdr[i+1].begin(), gdr[i+1].end(), pr)) == pr)continue;
seg.r1 = pr.first;
seg.r2 = pr.second;
seg.c2 = i;
for(int j = i;j>=0;j--){
auto it = (lower_bound(gdr[j].begin(), gdr[j].end(), pr));
if(it == gdr[j].end())break;
if((*it)!=pr)break;
seg.c1=j;
}
seg.c1--;
seg.c2++;
seg.c1=max(seg.c1, 0);
seg.c2=min(seg.c2, n-1);
for(int j = seg.c1;j<=seg.c2;j++) {
hsides[j][seg.r2].push_back(pair<int,int>(seg.c1, seg.r1));
}
}
}
for(int i = 0;i<m;i++) {
for(auto pr:gdc[i]) {
if(pr.first +1 == pr.second)continue;
if(i!=m-1 and (*lower_bound(gdc[i+1].begin(), gdc[i+1].end(), pr)) == pr)continue;
seg.r1 = pr.first;
seg.r2 = pr.second;
seg.c2 = i;
for(int j = i;j>=0;j--){
auto it = (lower_bound(gdc[j].begin(), gdc[j].end(), pr));
if(it == gdc[j].end())break;
if((*it)!=pr)break;
seg.c1=j;
}
seg.c1--;
seg.c2++;
seg.c1=max(seg.c1, 0);
seg.c2=min(seg.c2, m-1);
for(int j = seg.c1;j<=seg.c2;j++) {
vsides[seg.r2][j].push_back(pair<int,int>(seg.r1, seg.c1));
}
}
}
long long ans=0;
for(int i = 0;i<n;i++) {
for(int j = 0;j<m;j++) {
int tmp=0;
// for(auto k:vsides[i][j]) {
// for(auto k2:hsides[i][j]) {
// if(k2.second >= k.second and k2.first <= k.first) {
// ans++;
// tmp++;
// }
// }
// }
sort(vsides[i][j].begin(), vsides[i][j].end());
sort(hsides[i][j].begin(), hsides[i][j].end());
int p1 = 0;
for(auto k:vsides[i][j]) {
while(p1 < (int)hsides[i][j].size() and hsides[i][j][p1].first <= k.first) {
upd(hsides[i][j][p1++].second, 1);
}
ans+=qry(k.second);
}
while(p1 < (int)hsides[i][j].size()) {
upd(hsides[i][j][p1++].second, 1);
}
for(auto k:hsides[i][j]) {
upd(k.second, -1);
}
}
}
return ans;
}
Compilation message
rect.cpp: In function 'long long int count_rectangles(std::vector<std::vector<int> >)':
rect.cpp:135:8: warning: unused variable 'tmp' [-Wunused-variable]
135 | int tmp=0;
| ^~~
rect.cpp:99:10: warning: 'seg.Side::c1' may be used uninitialized in this function [-Wmaybe-uninitialized]
99 | seg.c1--;
| ~~~~~~^~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
40 ms |
296288 KB |
Output is correct |
2 |
Correct |
54 ms |
296304 KB |
Output is correct |
3 |
Correct |
40 ms |
296528 KB |
Output is correct |
4 |
Correct |
50 ms |
296540 KB |
Output is correct |
5 |
Correct |
44 ms |
296408 KB |
Output is correct |
6 |
Correct |
47 ms |
296532 KB |
Output is correct |
7 |
Correct |
48 ms |
296540 KB |
Output is correct |
8 |
Correct |
41 ms |
296536 KB |
Output is correct |
9 |
Correct |
40 ms |
296540 KB |
Output is correct |
10 |
Correct |
46 ms |
296528 KB |
Output is correct |
11 |
Correct |
59 ms |
296532 KB |
Output is correct |
12 |
Correct |
41 ms |
296532 KB |
Output is correct |
13 |
Correct |
42 ms |
296272 KB |
Output is correct |
14 |
Correct |
74 ms |
296468 KB |
Output is correct |
15 |
Correct |
70 ms |
296272 KB |
Output is correct |
16 |
Correct |
52 ms |
296452 KB |
Output is correct |
17 |
Correct |
60 ms |
296352 KB |
Output is correct |
18 |
Correct |
58 ms |
296404 KB |
Output is correct |
19 |
Correct |
51 ms |
296532 KB |
Output is correct |
20 |
Correct |
65 ms |
296276 KB |
Output is correct |
21 |
Correct |
42 ms |
296280 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
40 ms |
296288 KB |
Output is correct |
2 |
Correct |
54 ms |
296304 KB |
Output is correct |
3 |
Correct |
40 ms |
296528 KB |
Output is correct |
4 |
Correct |
50 ms |
296540 KB |
Output is correct |
5 |
Correct |
44 ms |
296408 KB |
Output is correct |
6 |
Correct |
47 ms |
296532 KB |
Output is correct |
7 |
Correct |
48 ms |
296540 KB |
Output is correct |
8 |
Correct |
41 ms |
296536 KB |
Output is correct |
9 |
Correct |
40 ms |
296540 KB |
Output is correct |
10 |
Correct |
46 ms |
296528 KB |
Output is correct |
11 |
Correct |
59 ms |
296532 KB |
Output is correct |
12 |
Correct |
41 ms |
296532 KB |
Output is correct |
13 |
Correct |
42 ms |
296272 KB |
Output is correct |
14 |
Correct |
74 ms |
296468 KB |
Output is correct |
15 |
Correct |
70 ms |
296272 KB |
Output is correct |
16 |
Correct |
52 ms |
296452 KB |
Output is correct |
17 |
Correct |
60 ms |
296352 KB |
Output is correct |
18 |
Correct |
58 ms |
296404 KB |
Output is correct |
19 |
Correct |
51 ms |
296532 KB |
Output is correct |
20 |
Correct |
65 ms |
296276 KB |
Output is correct |
21 |
Correct |
42 ms |
296280 KB |
Output is correct |
22 |
Correct |
57 ms |
296928 KB |
Output is correct |
23 |
Correct |
49 ms |
296932 KB |
Output is correct |
24 |
Correct |
46 ms |
297012 KB |
Output is correct |
25 |
Correct |
49 ms |
296980 KB |
Output is correct |
26 |
Correct |
56 ms |
297080 KB |
Output is correct |
27 |
Correct |
76 ms |
297300 KB |
Output is correct |
28 |
Correct |
79 ms |
297300 KB |
Output is correct |
29 |
Correct |
51 ms |
296788 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
40 ms |
296288 KB |
Output is correct |
2 |
Correct |
54 ms |
296304 KB |
Output is correct |
3 |
Correct |
40 ms |
296528 KB |
Output is correct |
4 |
Correct |
50 ms |
296540 KB |
Output is correct |
5 |
Correct |
44 ms |
296408 KB |
Output is correct |
6 |
Correct |
47 ms |
296532 KB |
Output is correct |
7 |
Correct |
48 ms |
296540 KB |
Output is correct |
8 |
Correct |
41 ms |
296536 KB |
Output is correct |
9 |
Correct |
40 ms |
296540 KB |
Output is correct |
10 |
Correct |
46 ms |
296528 KB |
Output is correct |
11 |
Correct |
59 ms |
296532 KB |
Output is correct |
12 |
Correct |
41 ms |
296532 KB |
Output is correct |
13 |
Correct |
42 ms |
296272 KB |
Output is correct |
14 |
Correct |
74 ms |
296468 KB |
Output is correct |
15 |
Correct |
70 ms |
296272 KB |
Output is correct |
16 |
Correct |
52 ms |
296452 KB |
Output is correct |
17 |
Correct |
57 ms |
296928 KB |
Output is correct |
18 |
Correct |
49 ms |
296932 KB |
Output is correct |
19 |
Correct |
46 ms |
297012 KB |
Output is correct |
20 |
Correct |
49 ms |
296980 KB |
Output is correct |
21 |
Correct |
56 ms |
297080 KB |
Output is correct |
22 |
Correct |
76 ms |
297300 KB |
Output is correct |
23 |
Correct |
79 ms |
297300 KB |
Output is correct |
24 |
Correct |
51 ms |
296788 KB |
Output is correct |
25 |
Correct |
60 ms |
296352 KB |
Output is correct |
26 |
Correct |
58 ms |
296404 KB |
Output is correct |
27 |
Correct |
51 ms |
296532 KB |
Output is correct |
28 |
Correct |
65 ms |
296276 KB |
Output is correct |
29 |
Correct |
42 ms |
296280 KB |
Output is correct |
30 |
Correct |
87 ms |
299064 KB |
Output is correct |
31 |
Correct |
64 ms |
299092 KB |
Output is correct |
32 |
Correct |
65 ms |
299112 KB |
Output is correct |
33 |
Correct |
72 ms |
300116 KB |
Output is correct |
34 |
Correct |
98 ms |
301520 KB |
Output is correct |
35 |
Correct |
109 ms |
301396 KB |
Output is correct |
36 |
Correct |
111 ms |
301104 KB |
Output is correct |
37 |
Correct |
106 ms |
300916 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
40 ms |
296288 KB |
Output is correct |
2 |
Correct |
54 ms |
296304 KB |
Output is correct |
3 |
Correct |
40 ms |
296528 KB |
Output is correct |
4 |
Correct |
50 ms |
296540 KB |
Output is correct |
5 |
Correct |
44 ms |
296408 KB |
Output is correct |
6 |
Correct |
47 ms |
296532 KB |
Output is correct |
7 |
Correct |
48 ms |
296540 KB |
Output is correct |
8 |
Correct |
41 ms |
296536 KB |
Output is correct |
9 |
Correct |
40 ms |
296540 KB |
Output is correct |
10 |
Correct |
46 ms |
296528 KB |
Output is correct |
11 |
Correct |
59 ms |
296532 KB |
Output is correct |
12 |
Correct |
41 ms |
296532 KB |
Output is correct |
13 |
Correct |
42 ms |
296272 KB |
Output is correct |
14 |
Correct |
74 ms |
296468 KB |
Output is correct |
15 |
Correct |
70 ms |
296272 KB |
Output is correct |
16 |
Correct |
52 ms |
296452 KB |
Output is correct |
17 |
Correct |
57 ms |
296928 KB |
Output is correct |
18 |
Correct |
49 ms |
296932 KB |
Output is correct |
19 |
Correct |
46 ms |
297012 KB |
Output is correct |
20 |
Correct |
49 ms |
296980 KB |
Output is correct |
21 |
Correct |
56 ms |
297080 KB |
Output is correct |
22 |
Correct |
76 ms |
297300 KB |
Output is correct |
23 |
Correct |
79 ms |
297300 KB |
Output is correct |
24 |
Correct |
51 ms |
296788 KB |
Output is correct |
25 |
Correct |
87 ms |
299064 KB |
Output is correct |
26 |
Correct |
64 ms |
299092 KB |
Output is correct |
27 |
Correct |
65 ms |
299112 KB |
Output is correct |
28 |
Correct |
72 ms |
300116 KB |
Output is correct |
29 |
Correct |
98 ms |
301520 KB |
Output is correct |
30 |
Correct |
109 ms |
301396 KB |
Output is correct |
31 |
Correct |
111 ms |
301104 KB |
Output is correct |
32 |
Correct |
106 ms |
300916 KB |
Output is correct |
33 |
Correct |
60 ms |
296352 KB |
Output is correct |
34 |
Correct |
58 ms |
296404 KB |
Output is correct |
35 |
Correct |
51 ms |
296532 KB |
Output is correct |
36 |
Correct |
65 ms |
296276 KB |
Output is correct |
37 |
Correct |
42 ms |
296280 KB |
Output is correct |
38 |
Correct |
205 ms |
332992 KB |
Output is correct |
39 |
Correct |
184 ms |
335784 KB |
Output is correct |
40 |
Correct |
209 ms |
335820 KB |
Output is correct |
41 |
Correct |
235 ms |
338792 KB |
Output is correct |
42 |
Correct |
185 ms |
332624 KB |
Output is correct |
43 |
Correct |
300 ms |
332612 KB |
Output is correct |
44 |
Correct |
202 ms |
332628 KB |
Output is correct |
45 |
Correct |
262 ms |
330840 KB |
Output is correct |
46 |
Correct |
280 ms |
334168 KB |
Output is correct |
47 |
Correct |
291 ms |
338900 KB |
Output is correct |
48 |
Correct |
540 ms |
360492 KB |
Output is correct |
49 |
Correct |
641 ms |
360948 KB |
Output is correct |
50 |
Correct |
304 ms |
328824 KB |
Output is correct |
51 |
Correct |
340 ms |
328732 KB |
Output is correct |
52 |
Correct |
486 ms |
354128 KB |
Output is correct |
53 |
Correct |
440 ms |
354516 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
84 ms |
296784 KB |
Output is correct |
2 |
Correct |
80 ms |
296688 KB |
Output is correct |
3 |
Correct |
89 ms |
296720 KB |
Output is correct |
4 |
Correct |
54 ms |
296276 KB |
Output is correct |
5 |
Correct |
66 ms |
297044 KB |
Output is correct |
6 |
Correct |
47 ms |
297044 KB |
Output is correct |
7 |
Correct |
99 ms |
297040 KB |
Output is correct |
8 |
Correct |
73 ms |
296844 KB |
Output is correct |
9 |
Correct |
47 ms |
297048 KB |
Output is correct |
10 |
Correct |
98 ms |
296528 KB |
Output is correct |
11 |
Correct |
81 ms |
296640 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
60 ms |
296352 KB |
Output is correct |
2 |
Correct |
58 ms |
296404 KB |
Output is correct |
3 |
Correct |
51 ms |
296532 KB |
Output is correct |
4 |
Correct |
65 ms |
296276 KB |
Output is correct |
5 |
Correct |
42 ms |
296280 KB |
Output is correct |
6 |
Correct |
45 ms |
296396 KB |
Output is correct |
7 |
Correct |
1374 ms |
510220 KB |
Output is correct |
8 |
Correct |
3023 ms |
749632 KB |
Output is correct |
9 |
Correct |
3111 ms |
751908 KB |
Output is correct |
10 |
Correct |
2945 ms |
758504 KB |
Output is correct |
11 |
Correct |
444 ms |
411332 KB |
Output is correct |
12 |
Correct |
747 ms |
505584 KB |
Output is correct |
13 |
Correct |
785 ms |
518904 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
40 ms |
296288 KB |
Output is correct |
2 |
Correct |
54 ms |
296304 KB |
Output is correct |
3 |
Correct |
40 ms |
296528 KB |
Output is correct |
4 |
Correct |
50 ms |
296540 KB |
Output is correct |
5 |
Correct |
44 ms |
296408 KB |
Output is correct |
6 |
Correct |
47 ms |
296532 KB |
Output is correct |
7 |
Correct |
48 ms |
296540 KB |
Output is correct |
8 |
Correct |
41 ms |
296536 KB |
Output is correct |
9 |
Correct |
40 ms |
296540 KB |
Output is correct |
10 |
Correct |
46 ms |
296528 KB |
Output is correct |
11 |
Correct |
59 ms |
296532 KB |
Output is correct |
12 |
Correct |
41 ms |
296532 KB |
Output is correct |
13 |
Correct |
42 ms |
296272 KB |
Output is correct |
14 |
Correct |
74 ms |
296468 KB |
Output is correct |
15 |
Correct |
70 ms |
296272 KB |
Output is correct |
16 |
Correct |
52 ms |
296452 KB |
Output is correct |
17 |
Correct |
57 ms |
296928 KB |
Output is correct |
18 |
Correct |
49 ms |
296932 KB |
Output is correct |
19 |
Correct |
46 ms |
297012 KB |
Output is correct |
20 |
Correct |
49 ms |
296980 KB |
Output is correct |
21 |
Correct |
56 ms |
297080 KB |
Output is correct |
22 |
Correct |
76 ms |
297300 KB |
Output is correct |
23 |
Correct |
79 ms |
297300 KB |
Output is correct |
24 |
Correct |
51 ms |
296788 KB |
Output is correct |
25 |
Correct |
87 ms |
299064 KB |
Output is correct |
26 |
Correct |
64 ms |
299092 KB |
Output is correct |
27 |
Correct |
65 ms |
299112 KB |
Output is correct |
28 |
Correct |
72 ms |
300116 KB |
Output is correct |
29 |
Correct |
98 ms |
301520 KB |
Output is correct |
30 |
Correct |
109 ms |
301396 KB |
Output is correct |
31 |
Correct |
111 ms |
301104 KB |
Output is correct |
32 |
Correct |
106 ms |
300916 KB |
Output is correct |
33 |
Correct |
205 ms |
332992 KB |
Output is correct |
34 |
Correct |
184 ms |
335784 KB |
Output is correct |
35 |
Correct |
209 ms |
335820 KB |
Output is correct |
36 |
Correct |
235 ms |
338792 KB |
Output is correct |
37 |
Correct |
185 ms |
332624 KB |
Output is correct |
38 |
Correct |
300 ms |
332612 KB |
Output is correct |
39 |
Correct |
202 ms |
332628 KB |
Output is correct |
40 |
Correct |
262 ms |
330840 KB |
Output is correct |
41 |
Correct |
280 ms |
334168 KB |
Output is correct |
42 |
Correct |
291 ms |
338900 KB |
Output is correct |
43 |
Correct |
540 ms |
360492 KB |
Output is correct |
44 |
Correct |
641 ms |
360948 KB |
Output is correct |
45 |
Correct |
304 ms |
328824 KB |
Output is correct |
46 |
Correct |
340 ms |
328732 KB |
Output is correct |
47 |
Correct |
486 ms |
354128 KB |
Output is correct |
48 |
Correct |
440 ms |
354516 KB |
Output is correct |
49 |
Correct |
84 ms |
296784 KB |
Output is correct |
50 |
Correct |
80 ms |
296688 KB |
Output is correct |
51 |
Correct |
89 ms |
296720 KB |
Output is correct |
52 |
Correct |
54 ms |
296276 KB |
Output is correct |
53 |
Correct |
66 ms |
297044 KB |
Output is correct |
54 |
Correct |
47 ms |
297044 KB |
Output is correct |
55 |
Correct |
99 ms |
297040 KB |
Output is correct |
56 |
Correct |
73 ms |
296844 KB |
Output is correct |
57 |
Correct |
47 ms |
297048 KB |
Output is correct |
58 |
Correct |
98 ms |
296528 KB |
Output is correct |
59 |
Correct |
81 ms |
296640 KB |
Output is correct |
60 |
Correct |
45 ms |
296396 KB |
Output is correct |
61 |
Correct |
1374 ms |
510220 KB |
Output is correct |
62 |
Correct |
3023 ms |
749632 KB |
Output is correct |
63 |
Correct |
3111 ms |
751908 KB |
Output is correct |
64 |
Correct |
2945 ms |
758504 KB |
Output is correct |
65 |
Correct |
444 ms |
411332 KB |
Output is correct |
66 |
Correct |
747 ms |
505584 KB |
Output is correct |
67 |
Correct |
785 ms |
518904 KB |
Output is correct |
68 |
Correct |
60 ms |
296352 KB |
Output is correct |
69 |
Correct |
58 ms |
296404 KB |
Output is correct |
70 |
Correct |
51 ms |
296532 KB |
Output is correct |
71 |
Correct |
65 ms |
296276 KB |
Output is correct |
72 |
Correct |
42 ms |
296280 KB |
Output is correct |
73 |
Correct |
2120 ms |
729540 KB |
Output is correct |
74 |
Correct |
2102 ms |
801304 KB |
Output is correct |
75 |
Correct |
2467 ms |
797008 KB |
Output is correct |
76 |
Correct |
2375 ms |
850448 KB |
Output is correct |
77 |
Correct |
2747 ms |
748804 KB |
Output is correct |
78 |
Correct |
3748 ms |
775132 KB |
Output is correct |
79 |
Correct |
4032 ms |
811380 KB |
Output is correct |
80 |
Runtime error |
4011 ms |
1048576 KB |
Execution killed with signal 9 |
81 |
Halted |
0 ms |
0 KB |
- |