Submission #650398

# Submission time Handle Problem Language Result Execution time Memory
650398 2022-10-13T16:41:18 Z ymm Bubble Sort 2 (JOI18_bubblesort2) C++17
17 / 100
110 ms 9292 KB
#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];
	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 == 0) break;
		if (l > 1) UP(l-1);
		if (r < S) 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);
}

__attribute__((optimize("O3,unroll-loops"),target("avx")))
void mov(int i, int x)
{
	int j = i - i%S;
	for (; srt_ind[j] != i; ++j);
	int dst =   lower_bound(srt_ind + i-i%S, srt_ind + i-i%S + S, x)
	          - srt_ind;
	if (j < dst) {
		--dst;
		Loop (k, j, dst) {
			val_srt[k] = val_srt[k+1];
			srt_ind[k] = srt_ind[k+1];
		}
	} else {
		LoopR (k, dst, j) {
			val_srt[k+1] = val_srt[k];
			srt_ind[k+1] = srt_ind[k];
		}
	}
	val_srt[dst] = x;
	srt_ind[dst] = i;
}

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;
}
# Verdict Execution time Memory Grader output
1 Correct 17 ms 8276 KB Output is correct
2 Correct 29 ms 8276 KB Output is correct
3 Correct 39 ms 8292 KB Output is correct
4 Correct 40 ms 8236 KB Output is correct
5 Correct 39 ms 8276 KB Output is correct
6 Correct 38 ms 8276 KB Output is correct
7 Correct 37 ms 8276 KB Output is correct
8 Correct 39 ms 8284 KB Output is correct
9 Correct 39 ms 8276 KB Output is correct
10 Correct 37 ms 8272 KB Output is correct
11 Correct 38 ms 8280 KB Output is correct
12 Correct 38 ms 8276 KB Output is correct
13 Correct 39 ms 8356 KB Output is correct
14 Correct 38 ms 8276 KB Output is correct
15 Correct 38 ms 8272 KB Output is correct
16 Correct 42 ms 8272 KB Output is correct
17 Correct 38 ms 8276 KB Output is correct
18 Correct 39 ms 8276 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 17 ms 8276 KB Output is correct
2 Correct 29 ms 8276 KB Output is correct
3 Correct 39 ms 8292 KB Output is correct
4 Correct 40 ms 8236 KB Output is correct
5 Correct 39 ms 8276 KB Output is correct
6 Correct 38 ms 8276 KB Output is correct
7 Correct 37 ms 8276 KB Output is correct
8 Correct 39 ms 8284 KB Output is correct
9 Correct 39 ms 8276 KB Output is correct
10 Correct 37 ms 8272 KB Output is correct
11 Correct 38 ms 8280 KB Output is correct
12 Correct 38 ms 8276 KB Output is correct
13 Correct 39 ms 8356 KB Output is correct
14 Correct 38 ms 8276 KB Output is correct
15 Correct 38 ms 8272 KB Output is correct
16 Correct 42 ms 8272 KB Output is correct
17 Correct 38 ms 8276 KB Output is correct
18 Correct 39 ms 8276 KB Output is correct
19 Incorrect 110 ms 8588 KB Output isn't correct
20 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 51 ms 9292 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 17 ms 8276 KB Output is correct
2 Correct 29 ms 8276 KB Output is correct
3 Correct 39 ms 8292 KB Output is correct
4 Correct 40 ms 8236 KB Output is correct
5 Correct 39 ms 8276 KB Output is correct
6 Correct 38 ms 8276 KB Output is correct
7 Correct 37 ms 8276 KB Output is correct
8 Correct 39 ms 8284 KB Output is correct
9 Correct 39 ms 8276 KB Output is correct
10 Correct 37 ms 8272 KB Output is correct
11 Correct 38 ms 8280 KB Output is correct
12 Correct 38 ms 8276 KB Output is correct
13 Correct 39 ms 8356 KB Output is correct
14 Correct 38 ms 8276 KB Output is correct
15 Correct 38 ms 8272 KB Output is correct
16 Correct 42 ms 8272 KB Output is correct
17 Correct 38 ms 8276 KB Output is correct
18 Correct 39 ms 8276 KB Output is correct
19 Incorrect 110 ms 8588 KB Output isn't correct
20 Halted 0 ms 0 KB -