#include "bubblesort2.h"
#include <bits/stdc++.h>
#define Loop(x,l,r) for (ll x = (l); x < (ll)(r); ++x)
#define LoopR(x,l,r) for (ll x = (r)-1; x >= (ll)(l); --x)
typedef long long ll;
typedef std::pair<int, int> pii;
typedef std::pair<ll , ll > pll;
using namespace std;
const int S = 2048;
const int N = 500'032 + S;
int seg_mx[N/S][S*2];
int seg_cnt[N/S][S*2];
int val[N], cnt[N], val_srt[N], srt_ind[N];
int n;
void add_seg(int l, int r, int x, int t)
{
#define UP(i) seg_mx[t][i] = max(seg_mx[t][2*(i)],\
seg_mx[t][2*(i)+1])\
+ seg_cnt[t][i];
vector<tuple<int,int,int>> bfs = {{1, 0, S}};
for (int j = 0; j < bfs.size(); ++j) {
auto [i, b, e] = bfs[j];
if (l <= b && e <= r) {
seg_cnt[t][i] += x;
seg_mx[t][i] += x;
} else {
if (l < (b+e)/2)
bfs.emplace_back(2*i+0, b, (b+e)/2);
if ((b+e)/2 < r)
bfs.emplace_back(2*i+1, (b+e)/2, e);
}
}
while (get<0>(bfs.back()) >= S)
bfs.pop_back();
LoopR (j,0,bfs.size()) {
int i = get<0>(bfs[j]);
UP(i);
}
return;
l += S; r += S;
for (;;) {
if (l != r && (l&1)) {
seg_mx[t][l] += x;
seg_cnt[t][l] += x;
++l;
}
if (l != r && (r&1)) {
--r;
seg_mx[t][r] += x;
seg_cnt[t][r] += x;
}
l /= 2; r /= 2;
if (l == 1) {
UP(r);
break;
}
UP(l-1); UP(r);
}
}
void build_seg(int t)
{
Loop (i,S,2*S) {
seg_cnt[t][i] = cnt[srt_ind[S*t + i-S]];
seg_mx[t][i] = seg_cnt[t][i];
}
LoopR (i,1,S) {
seg_mx[t][i] = max(seg_mx[t][2*i], seg_mx[t][2*i+1]);
seg_cnt[t][i] = 0;
}
}
void flush_seg(int t)
{
Loop (i,1,S) {
seg_cnt[t][2*i] += seg_cnt[t][i];
seg_cnt[t][2*i+1] += seg_cnt[t][i];
}
Loop (i,S,2*S)
cnt[srt_ind[t*S + i-S]] = seg_cnt[t][i];
}
int get_gt(int l, int x)
{
return S - ( upper_bound(val_srt + l, val_srt + l + S, x)
- (val_srt + l));
}
void up_range(int l, int xl, int xr, int val)
{
xl = lower_bound(val_srt + l, val_srt + l + S, xl) - val_srt - l;
xr = lower_bound(val_srt + l, val_srt + l + S, xr) - val_srt - l;
if (xl < xr)
add_seg(xl, xr, val, l/S);
}
void mov(int i, int x)
{
int j = i - i%S;
for (; srt_ind[j] != i; ++j);
val_srt[j] = x;
while ((j+1)%S && val_srt[j] > val_srt[j+1]) {
swap(val_srt[j], val_srt[j+1]);
swap(srt_ind[j], srt_ind[j+1]);
++j;
}
while (j%S && val_srt[j] < val_srt[j-1]) {
swap(val_srt[j], val_srt[j-1]);
swap(srt_ind[j], srt_ind[j-1]);
--j;
}
}
void up(int i, int x)
{
int pre = val[i];
int cnt_gt = 0;
for (int l = 0; l < i - i%S; l += S)
cnt_gt += get_gt(l, x);
flush_seg(i/S);
val[i] = x;
for (int j = i - i%S; j < i; ++j)
cnt_gt += val[j] > x;
cnt[i] = cnt_gt;
for (int j = i + 1; j < i - i%S + S; ++j)
cnt[j] += (x > val[j]) - (pre > val[j]);
mov(i, x);
build_seg(i/S);
int val = pre < x? 1: -1;
int xl = pre < x? pre: x;
int xr = pre < x? x: pre;
for (int l = i - i%S + S; l < n; l += S)
up_range(l, xl, xr, val);
}
int get()
{
int mx = 0;
for (int l = 0; l < n; l += S)
mx = max(seg_mx[l/S][1], mx);
return mx;
}
namespace init {
int fen[N];
void fen_add(int i, int x)
{
++i;
while (i < N) {
fen[i] += x;
i += i & -i;
}
}
int fen_get(int r)
{
int ans = 0;
while (r > 0) {
ans += fen[r];
r -= r & -r;
}
return ans;
}
void init(vector<int> A)
{
static pii pii_srt[N];
n = A.size();
Loop (i,0,n) {
val[i] = A[i];
pii_srt[i] = {val[i], i};
}
Loop (i,n,N) {
val[i] = 1e9+10;
pii_srt[i] = {val[i], i};
}
for (int l = 0; l < n; l += S) {
sort(pii_srt+l, pii_srt+l+S);
Loop (i,l,l+S) {
val_srt[i] = pii_srt[i].first;
srt_ind[i] = pii_srt[i].second;
}
}
vector<int> cmp(val, val+N);
sort(cmp.begin(), cmp.end());
cmp.resize(unique(cmp.begin(), cmp.end()) - cmp.begin());
int mx = cmp.size() - 1;
Loop (i,0,n) {
int j = lower_bound(cmp.begin(), cmp.end(), val[i])
- cmp.begin();
cnt[i] = fen_get(mx-j);
fen_add(mx-j, 1);
}
for (int l = 0; l < n; l += S)
build_seg(l/S);
}
}
vector<int> countScans(vector<int> A, vector<int> X, vector<int> V)
{
init::init(A);
vector<int> ans;
Loop (i,0,X.size()) {
up(X[i], V[i]);
ans.push_back(get());
}
return ans;
}
Compilation message
bubblesort2.cpp: In function 'void add_seg(int, int, int, int)':
bubblesort2.cpp:23:20: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::tuple<int, int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
23 | for (int j = 0; j < bfs.size(); ++j) {
| ~~^~~~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
19 ms |
8276 KB |
Output is correct |
2 |
Correct |
29 ms |
8276 KB |
Output is correct |
3 |
Correct |
44 ms |
8276 KB |
Output is correct |
4 |
Correct |
44 ms |
8260 KB |
Output is correct |
5 |
Correct |
40 ms |
8276 KB |
Output is correct |
6 |
Correct |
40 ms |
8276 KB |
Output is correct |
7 |
Correct |
39 ms |
8276 KB |
Output is correct |
8 |
Correct |
45 ms |
8276 KB |
Output is correct |
9 |
Correct |
44 ms |
8276 KB |
Output is correct |
10 |
Correct |
45 ms |
8280 KB |
Output is correct |
11 |
Correct |
45 ms |
8352 KB |
Output is correct |
12 |
Correct |
40 ms |
8280 KB |
Output is correct |
13 |
Correct |
44 ms |
8276 KB |
Output is correct |
14 |
Correct |
42 ms |
8276 KB |
Output is correct |
15 |
Correct |
45 ms |
8276 KB |
Output is correct |
16 |
Correct |
44 ms |
8276 KB |
Output is correct |
17 |
Correct |
43 ms |
8276 KB |
Output is correct |
18 |
Correct |
43 ms |
8272 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
19 ms |
8276 KB |
Output is correct |
2 |
Correct |
29 ms |
8276 KB |
Output is correct |
3 |
Correct |
44 ms |
8276 KB |
Output is correct |
4 |
Correct |
44 ms |
8260 KB |
Output is correct |
5 |
Correct |
40 ms |
8276 KB |
Output is correct |
6 |
Correct |
40 ms |
8276 KB |
Output is correct |
7 |
Correct |
39 ms |
8276 KB |
Output is correct |
8 |
Correct |
45 ms |
8276 KB |
Output is correct |
9 |
Correct |
44 ms |
8276 KB |
Output is correct |
10 |
Correct |
45 ms |
8280 KB |
Output is correct |
11 |
Correct |
45 ms |
8352 KB |
Output is correct |
12 |
Correct |
40 ms |
8280 KB |
Output is correct |
13 |
Correct |
44 ms |
8276 KB |
Output is correct |
14 |
Correct |
42 ms |
8276 KB |
Output is correct |
15 |
Correct |
45 ms |
8276 KB |
Output is correct |
16 |
Correct |
44 ms |
8276 KB |
Output is correct |
17 |
Correct |
43 ms |
8276 KB |
Output is correct |
18 |
Correct |
43 ms |
8272 KB |
Output is correct |
19 |
Correct |
133 ms |
8588 KB |
Output is correct |
20 |
Correct |
161 ms |
8620 KB |
Output is correct |
21 |
Correct |
140 ms |
8660 KB |
Output is correct |
22 |
Correct |
142 ms |
8660 KB |
Output is correct |
23 |
Correct |
156 ms |
8660 KB |
Output is correct |
24 |
Correct |
149 ms |
8612 KB |
Output is correct |
25 |
Correct |
147 ms |
8660 KB |
Output is correct |
26 |
Correct |
150 ms |
8660 KB |
Output is correct |
27 |
Correct |
153 ms |
8620 KB |
Output is correct |
28 |
Correct |
147 ms |
8612 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
66 ms |
9300 KB |
Output is correct |
2 |
Correct |
603 ms |
10136 KB |
Output is correct |
3 |
Correct |
1245 ms |
10940 KB |
Output is correct |
4 |
Correct |
1263 ms |
10956 KB |
Output is correct |
5 |
Correct |
1174 ms |
10996 KB |
Output is correct |
6 |
Correct |
1193 ms |
10956 KB |
Output is correct |
7 |
Correct |
1073 ms |
10956 KB |
Output is correct |
8 |
Correct |
1099 ms |
10956 KB |
Output is correct |
9 |
Correct |
1097 ms |
10956 KB |
Output is correct |
10 |
Correct |
745 ms |
10988 KB |
Output is correct |
11 |
Correct |
743 ms |
10996 KB |
Output is correct |
12 |
Correct |
739 ms |
10996 KB |
Output is correct |
13 |
Correct |
737 ms |
10996 KB |
Output is correct |
14 |
Correct |
707 ms |
10912 KB |
Output is correct |
15 |
Correct |
725 ms |
10880 KB |
Output is correct |
16 |
Correct |
723 ms |
10996 KB |
Output is correct |
17 |
Correct |
727 ms |
10956 KB |
Output is correct |
18 |
Correct |
712 ms |
10956 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
19 ms |
8276 KB |
Output is correct |
2 |
Correct |
29 ms |
8276 KB |
Output is correct |
3 |
Correct |
44 ms |
8276 KB |
Output is correct |
4 |
Correct |
44 ms |
8260 KB |
Output is correct |
5 |
Correct |
40 ms |
8276 KB |
Output is correct |
6 |
Correct |
40 ms |
8276 KB |
Output is correct |
7 |
Correct |
39 ms |
8276 KB |
Output is correct |
8 |
Correct |
45 ms |
8276 KB |
Output is correct |
9 |
Correct |
44 ms |
8276 KB |
Output is correct |
10 |
Correct |
45 ms |
8280 KB |
Output is correct |
11 |
Correct |
45 ms |
8352 KB |
Output is correct |
12 |
Correct |
40 ms |
8280 KB |
Output is correct |
13 |
Correct |
44 ms |
8276 KB |
Output is correct |
14 |
Correct |
42 ms |
8276 KB |
Output is correct |
15 |
Correct |
45 ms |
8276 KB |
Output is correct |
16 |
Correct |
44 ms |
8276 KB |
Output is correct |
17 |
Correct |
43 ms |
8276 KB |
Output is correct |
18 |
Correct |
43 ms |
8272 KB |
Output is correct |
19 |
Correct |
133 ms |
8588 KB |
Output is correct |
20 |
Correct |
161 ms |
8620 KB |
Output is correct |
21 |
Correct |
140 ms |
8660 KB |
Output is correct |
22 |
Correct |
142 ms |
8660 KB |
Output is correct |
23 |
Correct |
156 ms |
8660 KB |
Output is correct |
24 |
Correct |
149 ms |
8612 KB |
Output is correct |
25 |
Correct |
147 ms |
8660 KB |
Output is correct |
26 |
Correct |
150 ms |
8660 KB |
Output is correct |
27 |
Correct |
153 ms |
8620 KB |
Output is correct |
28 |
Correct |
147 ms |
8612 KB |
Output is correct |
29 |
Correct |
66 ms |
9300 KB |
Output is correct |
30 |
Correct |
603 ms |
10136 KB |
Output is correct |
31 |
Correct |
1245 ms |
10940 KB |
Output is correct |
32 |
Correct |
1263 ms |
10956 KB |
Output is correct |
33 |
Correct |
1174 ms |
10996 KB |
Output is correct |
34 |
Correct |
1193 ms |
10956 KB |
Output is correct |
35 |
Correct |
1073 ms |
10956 KB |
Output is correct |
36 |
Correct |
1099 ms |
10956 KB |
Output is correct |
37 |
Correct |
1097 ms |
10956 KB |
Output is correct |
38 |
Correct |
745 ms |
10988 KB |
Output is correct |
39 |
Correct |
743 ms |
10996 KB |
Output is correct |
40 |
Correct |
739 ms |
10996 KB |
Output is correct |
41 |
Correct |
737 ms |
10996 KB |
Output is correct |
42 |
Correct |
707 ms |
10912 KB |
Output is correct |
43 |
Correct |
725 ms |
10880 KB |
Output is correct |
44 |
Correct |
723 ms |
10996 KB |
Output is correct |
45 |
Correct |
727 ms |
10956 KB |
Output is correct |
46 |
Correct |
712 ms |
10956 KB |
Output is correct |
47 |
Correct |
6164 ms |
17216 KB |
Output is correct |
48 |
Execution timed out |
9083 ms |
34436 KB |
Time limit exceeded |
49 |
Halted |
0 ms |
0 KB |
- |