QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#626603#7904. Rainbow SubarraywsxcbAC ✓1312ms249804kbC++172.8kb2024-10-10 10:40:332024-10-10 10:40:34

Judging History

你现在查看的是最新测评结果

  • [2024-10-10 10:40:34]
  • 评测
  • 测评结果:AC
  • 用时:1312ms
  • 内存:249804kb
  • [2024-10-10 10:40:33]
  • 提交

answer

#pragma GCC optimize(3,"Ofast","inline")
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define pb push_back
#define fi first
#define se second
//#define int long long
const int N = 5e5 + 5;
int a[N], t[N], cnt = 0, len;

struct node {
	int ls, rs, cnt = 0;
	ll sum = 0;
} tr[N * 20];
int id = 0;

int updata(int p, int pre, int l, int r, int pos, int k) {
	p = ++id;
	//cout << id << '\n';
	tr[p]=tr[pre];
	tr[p].sum += k;
	tr[p].cnt += 1;
	//cout << l << ' ' << r << ' ' << tr[p].sum << ' ' << tr[p].cnt << '\n';
	if (l == r) {
		return p;
	}
	int mid = (l + r) / 2;
	if (pos <= mid)
		tr[p].ls = updata(tr[p].ls, tr[pre].ls, l, mid, pos, k);
	else
		tr[p].rs = updata(tr[p].rs, tr[pre].rs, mid + 1, r, pos, k);
	return p;
}
ll ct,st,ql,qr;
void get(int p, int l, int r) {
	if (ql <= l && r <= qr) {
		ct+=tr[p].cnt;
		st+=tr[p].sum;
		return ;
	}
	int mid = (l + r) / 2;
	if (mid >= ql) {
		get(tr[p].ls, l, mid);
	}
	if (mid < qr) {
		get(tr[p].rs, mid + 1, r);
	}
}
int rt[N];

pair<int, ll> query1(int l, int r, int x) {	
	st=ct=0;
	ql=x+1,qr=len;
	get(rt[r], 1, len);
	ll c1=ct,s1=st;
	st=ct=0;
	get(rt[l - 1], 1, len);
	return {c1 - ct, s1 - st};
}

int query2(int p1, int p2, int l, int r, int k) {
	if (l == r)
		return l;
	int mid = (l + r) / 2;
	int c = tr[tr[p1].ls].cnt - tr[tr[p2].ls].cnt;
	if (c >= k)
		return query2(tr[p1].ls, tr[p2].ls, l, mid, k);
	return query2(tr[p1].rs, tr[p2].rs, mid + 1, r, k - c);
}

int query2(int l, int r) {
	int mid = (r - l + 2) / 2;
	return query2(rt[r], rt[l - 1], 1, len, mid);
}
int pos[N];
ll pre[N];

void solve() {
	ll n, k;
	cin >> n >> k;
	for (int i = 1; i <= n; i++) {
		cin >> a[i];
		a[i] -= i;
		t[++cnt] = a[i];
		pre[i] = pre[i - 1] + a[i];
	}
	sort(t + 1, t + 1 + cnt);
	len = unique(t + 1, t + 1 + cnt) - t - 1;
	for (int i = 1; i <= n; i++) {
		//cout << a[i] << ' ';
		pos[i] = lower_bound(t + 1, t + 1 + len, a[i]) - t;
		//cout << pos[i] << ' ';
		//continue;
		rt[i] = updata(rt[i], rt[i - 1], 1, len, pos[i], a[i]);
	}

	auto check = [&](int x) {
		for (int i = x ; i <= n ; i++) {
			int l = i - x + 1, r = i;
			int zh = query2(i - x + 1, i);
			ll xt = t[zh];
			auto [c1, s1] = query1(l, r, zh);
			ll c2 = r - l + 1 - c1, s2 = pre[r] - pre[l - 1] - s1;
			ll ans = c2 * xt - s2 + (s1 - c1 * xt);

			if (ans <= k)
				return true;
		}
		return false;
	};

	int le = 1, ri = n, ans = 1;
	while (le <= ri) {
		int mid = (le + ri) >> 1;
		//cout << ri << ' ' << mid << ' ' << le << '\n';
		if (check(mid)) {
			le = mid + 1;
			ans = mid;
		} else
			ri = mid - 1;
	}

	cout << ans << '\n';
	cnt = 0;
	id = 0;
}

signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int t = 1;
	cin >> t;
	while (t--)
		solve();
	return 0;
}

这程序好像有点Bug,我给组数据试试?

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 12ms
memory: 239012kb

input:

5
7 5
7 2 5 5 4 11 7
6 0
100 3 4 5 99 100
5 6
1 1 1 1 1
5 50
100 200 300 400 500
1 100
3

output:

4
3
5
1
1

result:

ok 5 lines

Test #2:

score: 0
Accepted
time: 429ms
memory: 241036kb

input:

11102
2 167959139
336470888 134074578
5 642802746
273386884 79721198 396628655 3722503 471207868
6 202647942
268792718 46761498 443917727 16843338 125908043 191952768
2 717268783
150414369 193319712
6 519096230
356168102 262263554 174936674 407246545 274667941 279198849
9 527268921
421436316 3613460...

output:

1
4
3
2
6
5
7
2
4
1
4
1
1
3
2
2
7
8
7
7
1
7
6
2
4
3
1
6
7
7
3
4
3
9
3
8
6
6
3
1
6
3
1
2
4
6
4
6
4
1
4
7
1
6
3
5
6
6
1
7
5
3
1
6
4
5
3
2
2
6
2
3
10
1
4
3
2
4
5
1
7
5
5
5
8
5
3
6
3
5
5
8
5
4
5
2
1
5
2
3
3
4
8
1
3
1
2
2
8
3
1
6
8
1
8
4
5
6
6
8
4
8
3
2
8
4
5
6
2
6
2
4
1
5
4
5
3
2
4
1
2
1
4
5
8
3
7
3
3
3...

result:

ok 11102 lines

Test #3:

score: 0
Accepted
time: 1123ms
memory: 249744kb

input:

1
500000 17244641009859
54748096 75475634 204928248 276927808 84875072 103158867 27937890 322595515 186026685 45468307 69240390 139887597 188586447 373764525 121365644 310156469 185188306 60350786 211308832 384695957 370562147 208427221 35937909 267590963 126478310 275357775 55361535 335993561 36696...

output:

172998

result:

ok single line: '172998'

Test #4:

score: 0
Accepted
time: 1312ms
memory: 249804kb

input:

1
500000 43483524125987
923264237 374288891 535590429 751244358 124321145 232930851 266089174 543529670 773363571 319728747 580543238 582720391 468188689 490702144 598813561 138628383 284660056 733781508 155605777 931759705 245485733 723534730 257812292 794937524 596788519 188451996 981010588 144836...

output:

174260

result:

ok single line: '174260'

Test #5:

score: 0
Accepted
time: 122ms
memory: 249688kb

input:

1
499999 62499749999
496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081 496813081...

output:

499998

result:

ok single line: '499998'

Extra Test:

score: 0
Extra Test Passed