QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#279173 | #7904. Rainbow Subarray | Wood | TL | 1271ms | 255276kb | C++23 | 3.0kb | 2023-12-08 13:11:15 | 2023-12-08 13:11:16 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
using i64 = long long;
#ifdef LOCAL
#include "algo/debug.h"
#else
#define debug(...) 42
#endif
struct Node {
Node *l = nullptr;
Node *r = nullptr;
int cnt = 0;
i64 sum = 0;
Node() : l{}, r{}, cnt{}, sum{} {}
};
Node *add(Node *t, int l, int r, int x, int v) {
Node *nt = new Node();
if (t) {
*nt = *t;
}
nt->cnt += 1;
nt->sum += v;
if (r - l > 1) {
int m = (l + r) / 2;
if (m > x) {
nt->l = add(nt->l, l, m, x, v);
} else {
nt->r = add(nt->r, m, r, x, v);
}
}
return nt;
}
int kth(Node *t1, Node *t2, int l, int r, int k) {
if (r - l == 1) {
return l;
}
int cnt = (t2 && t2->l ? t2->l->cnt : 0) - (t1 && t1->l ? t1->l->cnt : 0);
int m = (l + r) / 2;
if (cnt >= k) {
return kth(t1 ? t1->l : nullptr, t2 ? t2->l : nullptr, l, m, k);
} else {
return kth(t1 ? t1->r : nullptr, t2 ? t2->r : nullptr, m, r, k - cnt);
}
}
pair<int, i64> operator+(const pair<int, i64> &lhs, const pair<int, i64> &rhs) {
return make_pair(lhs.first + rhs.first, lhs.second + rhs.second);
}
pair<int, i64> query(Node *t1, Node *t2, int l, int r, int x, int y) {
if (l >= x && r <= y) {
return make_pair((t2 ? t2->cnt : 0) - (t1 ? t1->cnt : 0), (t2 ? t2->sum : 0) - (t1 ? t1->sum : 0));
}
if (l >= y || r <= x) {
return make_pair(0, 0);
}
int m = (l + r) / 2;
return query(t1 ? t1->l : nullptr, t2 ? t2->l : nullptr, l, m, x, y) + query(t1 ? t1->r : nullptr, t2 ? t2->r : nullptr, m, r, x, y);
}
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int tt;
cin >> tt;
while (tt--) {
int n;
i64 k;
cin >> n >> k;
vector<int> a(n);
for (int i = 0; i < n; i++) {
cin >> a[i];
a[i] -= i;
}
auto v = a;
sort(v.begin(), v.end());
v.erase(unique(v.begin(), v.end()), v.end());
const int m = (int) v.size();
vector<int> b(n);
for (int i = 0; i < n; i++) {
b[i] = lower_bound(v.begin(), v.end(), a[i]) - v.begin();
}
vector<Node *> node(n + 1);
for (int i = 0; i < n; i++) {
node[i + 1] = add(node[i], 0, m, b[i], a[i]);
}
debug(a);
auto calc = [&](int x) {
i64 ans = 1e18;
const int s = (x + 1) / 2;
for (int i = 0; i + x <= n; i++) {
int p = kth(node[i], node[i + x], 0, m, s);
// debug(v[p]);
auto [c1, s1] = query(node[i], node[i + x], 0, m, 0, p);
auto [c2, s2] = query(node[i], node[i + x], 0, m, p, m);
i64 res = 1LL * v[p] * c1 - s1 + s2 - 1LL * v[p] * c2;
ans = min(ans, res);
}
// debug(x, ans);
return ans;
};
int low = 1, high = n;
int ans = 1;
while (low <= high) {
int mid = (low + high) / 2;
if (calc(mid) <= k) {
ans = mid;
low = mid + 1;
} else {
high = mid - 1;
}
}
cout << ans << '\n';
}
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 1ms
memory: 3492kb
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: 1271ms
memory: 255276kb
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: -100
Time Limit Exceeded
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...