QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#326579 | #8229. 栈 | Xiaohuba | 18 | 24ms | 108572kb | C++23 | 9.6kb | 2024-02-13 14:36:28 | 2024-02-13 14:36:28 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
// #define LOCK_GETCHAR
// #define USE_INT_128
#if __cplusplus < 201400
#warning "Please use c++14 or higher."
#define CONSTEXPR_FUNC
#define ENABLE_IF_INT
#else
#define CONSTEXPR_FUNC constexpr
#define ENABLE_IF_INT , enable_if_t<_is_integer<T>, int> = 0
template <class T> constexpr bool _is_integer = numeric_limits<T>::is_integer;
template <> constexpr bool _is_integer<bool> = false;
template <> constexpr bool _is_integer<char> = false;
#ifdef USE_INT_128
template <> constexpr bool _is_integer<__int128> = true;
template <> constexpr bool _is_integer<__uint128_t> = true;
#endif
template <class T ENABLE_IF_INT>
constexpr T INF = numeric_limits<T>::max() >> 1;
#endif
#if !defined(_WIN32) && !defined(LOCK_GETCHAR)
#define getchar getchar_unlocked
#endif
#define il inline
#define mkp make_pair
#define fi first
#define se second
#define For(i, j, k) for (decltype(j - k) i = (j); i <= (k); ++i) // NOLINT
#define ForDown(i, j, k) for (decltype(j - k) i = (j); i >= (k); --i) // NOLINT
#define pb push_back
#define eb emplace_back
#ifndef ONLINE_JUDGE
#define FileIO(filename) \
freopen(filename ".in", "r", stdin); \
freopen(filename ".out", "w", stdout)
#else
#define FileIO(filename) void(0)
#endif
using ll = long long;
using uint = unsigned int;
using ull = unsigned long long;
using db = double;
using ldb = long double;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
#ifdef USE_INT_128
using lll = __int128_t;
using ulll = __uint128_t;
#endif
// clang-format off
template<typename T> constexpr il T sq(const T & x){ return x * x; }
template<typename T> CONSTEXPR_FUNC il void cmin(T & x, const T &y){ x = min(x, y); }
template<typename T> CONSTEXPR_FUNC il void cmax(T & x, const T &y){ x = max(x, y);}
template<typename T> CONSTEXPR_FUNC il T qpow(T x, ull y, T mod){T ans = 1; x %= mod; while (y) { if(y & 1)(ans *= x) %= mod;(x *= x) %= mod; y >>= 1;} return ans;}
template<typename T> CONSTEXPR_FUNC il T qpow(T x, ull y){T ans = 1; while (y) {if(y & 1) ans *= x;x *= x;y >>= 1;} return ans;}
template<typename T ENABLE_IF_INT> il void read(T &x){ x = 0; int f = 1; int c = getchar(); while(!isdigit(c)) {if (c == '-') f = -1;c = getchar();} while(isdigit(c)) {x = x * 10 + c - '0';c = getchar();} x *= f;}
template<typename T, typename ... Args> il void read(T &x, Args &... y){ read(x); read(y...); }
// clang-format on
// File head end
namespace {
class WBLT {
static il constexpr db alpha = 1.0 - 1.4142136 / 2;
static il constexpr ll INITIAL_SIZE = 3e6 + 5;
struct Node {
int lc, rc, sz;
ll sum, weight;
Node() : lc(0), rc(0), sz(0), sum(0), weight(0) {}
};
static il array<Node, INITIAL_SIZE> T;
static il int cnt = 0;
static il int new_nd() {
// if (cnt + 1 == T.size())
// T.resize(T.size() * 1.5);
return ++cnt;
}
#define lc(p) (T[p].lc)
#define rc(p) (T[p].rc)
#define sz(p) (T[p].sz)
#define weight(p) (T[p].weight)
static il void pu(int p) {
assert(lc(p) && rc(p));
sz(p) = sz(lc(p)) + sz(rc(p));
T[p].sum = T[lc(p)].sum + T[rc(p)].sum;
T[p].weight = T[lc(p)].weight + T[rc(p)].weight;
}
static il int __merge(int lc, int rc) {
int p = new_nd();
return lc(p) = lc, rc(p) = rc, pu(p), p;
}
static int merge(int p, int q) {
if (!p || !q)
return p | q;
auto val = alpha * (sz(p) + sz(q));
if (min(sz(p), sz(q)) >= val)
return __merge(p, q);
else if (sz(p) >= sz(q)) {
if (sz(lc(p)) >= val) {
int pre = p;
p = new_nd(), T[p] = T[pre];
return rc(p) = merge(rc(p), q), pu(p), p;
} else {
int u = merge(lc(p), lc(rc(p))), v = merge(rc(rc(p)), q);
return merge(u, v);
}
} else {
if (sz(rc(q)) >= val) {
int pre = q;
q = new_nd(), T[q] = T[pre];
return lc(q) = merge(p, lc(q)), pu(q), q;
} else {
int u = merge(p, lc(lc(q))), v = merge(rc(lc(q)), rc(q));
return merge(u, v);
}
}
}
static pii split(int p, ll rk) {
assert(p);
if (!rk)
return {0, p};
else if (sz(p) == 1) {
if (rk == weight(p))
return {p, 0};
else {
assert(rk < weight(p));
int u = new_nd(), v = new_nd();
ll val = T[p].sum / weight(p);
sz(u) = sz(v) = 1;
weight(u) = rk, T[u].sum = val * weight(u);
weight(v) = weight(p) - rk, T[v].sum = val * weight(v);
return {u, v};
}
}
if (rk <= weight(lc(p))) {
auto [u, v] = split(lc(p), rk);
return {u, merge(v, rc(p))};
} else {
auto [u, v] = split(rc(p), rk - weight(lc(p)));
return {merge(lc(p), u), v};
}
}
static ll query(int p, ll ql, ll qr) {
// if (!(ql >= 1 && ql <= qr))
// cerr << ql << ' ' << qr << '\n';
assert(ql >= 1 && ql <= qr);
if (!p)
return 0;
else if (sz(p) == 1) {
ll cnt = min(qr, weight(p)) - ql + 1, val = T[p].sum / T[p].weight;
return cnt * val;
} else if (ql == 1 && qr >= weight(p))
return T[p].sum;
ll ans = 0, mid = weight(lc(p));
if (ql <= mid)
ans += query(lc(p), ql, qr);
if (qr > mid)
ans += query(rc(p), max(1ll, ql - mid), qr - mid);
return ans;
}
int rt = 0;
WBLT(int _rt) : rt(_rt) {}
public:
WBLT() : rt(0) {}
il void join(WBLT rhs) { rt = merge(rt, rhs.rt); }
static il WBLT join2(WBLT x, WBLT y) { return {merge(x.rt, y.rt)}; }
il void push(ll val, ll cnt) {
int p = new_nd();
sz(p) = 1, weight(p) = cnt, T[p].sum = val * weight(p);
rt = merge(rt, p);
}
il void pop_back(ll cnt) {
if (weight(rt) <= cnt)
return rt = 0, void();
auto [u, v] = split(rt, weight(rt) - cnt);
assert(weight(u) == weight(rt) - cnt);
rt = u;
}
il void pop_front(ll cnt) {
if (weight(rt) <= cnt)
return rt = 0, void();
auto [u, v] = split(rt, cnt);
rt = v;
}
il ll qry_sum(ll st, ll ed) const {
if (st > this->size())
return 0;
else {
// cerr << st << ' ' << min(ed, weight(rt)) << '\n';
return query(rt, st, min(ed, weight(rt)));
}
}
il ll size() const { return weight(rt); }
il void clear() { rt = 0; }
#undef lc
#undef rc
#undef sz
#undef weight
};
namespace SGT {
constexpr ll MAXN = 1e5 + 5;
using tag_t = tuple<short, ll, WBLT>;
il short &_type(tag_t &x) { return get<0>(x); }
struct Node {
int l, r;
array<tag_t, 2> tags;
Node() : l(0), r(0), tags() {}
} static T[MAXN << 1];
#define mid(p) ((T[p].l + T[p].r) >> 1)
#define lc(p) (mid(p) << 1)
#define rc(p) (mid(p) << 1 | 1)
il void f(int p, tag_t tg) {
int id = !!get<0>(T[p].tags[1]), tp1 = _type(T[p].tags[id]), tp2 = _type(tg);
tag_t &cur = T[p].tags[id];
assert(!id || tp1 == 2);
if (tp2 == 1) { // pop
ll cnt = get<1>(tg), val = get<2>(cur).size();
if (tp1 == 1)
get<1>(cur) += cnt;
else if (val >= cnt)
get<2>(cur).pop_back(cnt);
else {
get<2>(cur).clear();
if (id) {
assert(get<0>(T[p].tags[0]) == 1);
cur = {};
}
_type(T[p].tags[0]) = 1;
get<1>(T[p].tags[0]) += cnt - val;
}
} else if (tp2 == 2) { // push
if (tp1 != 1)
_type(cur) = 2, get<2>(cur).join(get<2>(tg));
else {
assert(id == 0);
T[p].tags[1] = tg;
assert(get<0>(T[p].tags[1]) == 2);
}
}
}
il void pd(int p) {
f(lc(p), T[p].tags[0]), f(rc(p), T[p].tags[0]);
f(lc(p), T[p].tags[1]), f(rc(p), T[p].tags[1]);
T[p].tags[0] = T[p].tags[1] = {};
}
void build(int p, int l, int r) {
T[p].l = l, T[p].r = r;
if (l == r)
return;
int mid = (l + r) >> 1;
build(lc(p), l, mid), build(rc(p), mid + 1, r);
}
void push(int p, int ql, int qr, int val, int cnt) {
int l = T[p].l, r = T[p].r;
if (ql <= l && qr >= r) {
WBLT tag;
tag.push(val, cnt);
f(p, make_tuple(2, 0, tag));
return;
}
pd(p);
if (ql <= mid(p))
push(lc(p), ql, qr, val, cnt);
if (qr > mid(p))
push(rc(p), ql, qr, val, cnt);
}
void pop(int p, int ql, int qr, ll cnt) {
int l = T[p].l, r = T[p].r;
if (ql <= l && qr >= r)
return f(p, make_tuple(1, cnt, WBLT{}));
pd(p);
if (ql <= mid(p))
pop(lc(p), ql, qr, cnt);
if (qr > mid(p))
pop(rc(p), ql, qr, cnt);
}
ll qry(int p, int pos, ll L, ll R) {
if (T[p].l == T[p].r) {
int id = (get<0>(T[p].tags[0]) == 1);
// cerr << "> " << id << ' ' << get<1>(T[p].tags[0]) << ' '
// << get<2>(T[p].tags[id]).size() << '\n';
return get<2>(T[p].tags[id]).qry_sum(L, R);
}
pd(p);
if (pos <= mid(p))
return qry(lc(p), pos, L, R);
else
return qry(rc(p), pos, L, R);
}
#undef mid
#undef lc
#undef rc
}; // namespace SGT
int n, m;
il void Main() {
read(n, m);
SGT::build(1, 1, n);
For(i, 1, m) {
int op, x = 0;
ll y = 0, z = 0, w = 0;
read(op, x, y, z);
if (op == 1) {
read(w);
SGT::push(1, x, y, w, z);
} else if (op == 2)
SGT::pop(1, x, y, z);
else if (op == 3)
printf("%lld\n", SGT::qry(1, x, y, z));
// if (i > 16)
// cout << op << ' ' << x << ' ' << y << ' ' << z << '\n';
}
// WBLT tr;
// tr.push(2, 3);
// cout << tr.size() << ' ' << tr.qry_sum(2, 4) << '\n';
}
} // namespace
signed main() { return Main(), 0; }
/*
[1 10]
[1 5][6 10]
[1 3][4 5][6 8][9 10]
[1 2][3][4][5][6 7][8][9][10]
[1][2] [6][7]
*/
详细
Subtask #1:
score: 18
Accepted
Test #1:
score: 18
Accepted
time: 3ms
memory: 108444kb
input:
4907 4910 2 763 3330 1 3 307 1 1 1 2262 3430 22699 89397 1 1915 4000 51541 67587 2 212 2990 9763 2 1086 2162 1 2 1813 4496 16760 1 51 2796 68005 99390 1 1267 1519 74236 66178 3 1768 23808 54314 2 900 4122 27758 3 3287 17350 28989 2 3277 4024 3633 2 444 4866 1 2 353 4219 1061 1 987 3141 99906 17320 2...
output:
0 3032090730 903396180 471569175 200648623 98486697 647114751 123945 50793012 61782451 0 0 0 762429740 321140700 871619914 536311874 5361094892 0 1792521566 6640518748 2415375780 249435711 225987900 5250788038 1145132507 140071334 0 118545795 3086405469 5646099271 84280112 1232466642 4992966775 7968...
result:
ok 1622 numbers
Test #2:
score: 0
Accepted
time: 8ms
memory: 108512kb
input:
4992 4958 2 2965 3892 1 3 2141 1 1 3 4963 1 1 3 2298 1 1 3 2236 1 1 1 3393 4668 65533 8224 1 884 2343 86158 41289 3 4617 12174 63763 2 898 4236 44143 2 2860 4246 1 2 2696 3014 1 2 496 1635 15779 3 2230 8262 39805 2 3287 3627 5350 2 3443 4900 19874 1 535 1622 26926 88406 1 3272 3747 94290 19081 2 812...
output:
0 0 0 0 424276160 1302420216 0 393525459 0 188112684 0 38587680 696225296 717180100 2193537294 297696966 0 0 0 124461621 26876032 1609925499 0 3681040200 51602516 1502016 0 8636592 1138256753 196684293 0 16126264 959145423 58640451 1945754097 2949696960 0 3577791360 2029416288 2361290004 5882833609 ...
result:
ok 1597 numbers
Test #3:
score: 0
Accepted
time: 16ms
memory: 108572kb
input:
4980 4984 1 183 4891 75896 45281 2 767 3528 1367 3 2313 45535 49112 2 529 4006 1568 2 2114 2971 3819 3 3237 30655 31381 1 2074 2176 51631 35339 3 1602 95 16082 2 1340 3727 9249 2 105 1907 11928 3 2461 15189 33999 2 1450 1956 4721 1 700 4760 3043 92859 2 329 2992 6515 3 1295 10832 40486 2 3477 4447 8...
output:
162015418 32919287 723952628 851780891 1342808055 38307726 4701651115 903944603 240532672 652952020 1168430924 2253203546 3542990917 5872603595 305017015 657095398 25321688 1834305604 0 256832266 2732654054 1757936801 1158797383 656866283 3470700279 694675745 1042863834 76284096 6705704850 475899629...
result:
ok 1645 numbers
Test #4:
score: 0
Accepted
time: 16ms
memory: 108564kb
input:
4976 4948 2 858 1218 1 1 780 1528 70910 12344 1 681 4398 25997 59182 1 4564 4692 72420 96925 1 1124 2604 98159 98651 3 4734 1 1 2 1921 3430 3373 1 3805 3909 56118 23383 2 1471 2065 23679 2 1052 1154 30740 1 1098 2180 13716 29728 1 1094 3585 2073 93894 1 2024 4201 39023 1713 3 1571 21453 96893 3 1297...
output:
0 7262943486 185110624 53327400 957813600 383014415 1539405803 896522316 1454164560 7158196459 479198625 1943839360 1189657450 23355822139 2684778350 183742084 6400082784 2310401225 2082631008 5644811789 1875949890 3185562597 7185156304 3147144197 1588457333 676240200 1122598010 8758314557 100699296...
result:
ok 1663 numbers
Test #5:
score: 0
Accepted
time: 24ms
memory: 108512kb
input:
4944 4934 1 468 4845 87517 63656 3 4756 22899 79177 1 761 1054 45331 86403 1 24 2806 46189 11490 1 2602 4446 12528 14317 3 2601 51537 65051 1 1502 3573 79699 84830 3 1998 35405 151264 1 2400 4041 95071 83748 1 2050 3772 23643 53614 3 2261 51072 236192 2 1317 1908 6197 2 949 2543 30190 1 1457 4573 33...
output:
3582496024 860310840 5337461878 10833286574 1397502876 3735482073 4207877274 17671620 10854427218 1484917319 5462491276 1497165465 1453546510 1672688712 1158344316 1014734250 3797802047 15668090927 14634073116 32337553147 2159971110 12088416736 90924880 1410366456 13829776128 12126485158 18393654569...
result:
ok 829 numbers
Subtask #2:
score: 0
Runtime Error
Test #6:
score: 0
Runtime Error
input:
99999 99998 1 5026 18575 27178 90423 3 30623 1 1 3 76936 1 1 1 77021 95683 84664 24734 1 46085 74886 40512 11266 3 5048 8594 22468 1 53318 77721 97151 70784 1 70645 91192 37556 13013 1 56752 56940 91812 62887 1 7928 34576 87339 69404 3 74875 32807 100970 3 22338 17221 25771 3 21421 20602 57957 3 717...
output:
0 0 1254619125 4366274868 593473604 2592655824 3657975552 5652513833 110091352 1226646296 1989326852 763582808 8205318671 1659086055 3012598941 20085582585 3242801176 17381308704 24555397019 4722824224 20308857160 899316516 38935050954 988382364 13341823621 11397759491 2449683584 5875277101 80572355...
result:
Subtask #3:
score: 0
Runtime Error
Test #12:
score: 0
Runtime Error
input:
100000 99993 1 47773 70467 16065 1 2 52349 78446 2304 3 40821 1 1 1 40216 93069 78144 1 1 41089 43671 76025 1 2 35263 68629 31066 3 79881 13534 57327 3 5556 1 1 2 21962 38192 1 1 664 58116 9417 1 3 28089 6039 7989 2 88500 90302 9946 3 63215 49410 60770 2 11069 89527 57581 2 70303 97603 12363 1 3420 ...
output:
0 43794 0 1951 11361 129 898 29245 7969 1947 34972 16405 59952 123666 24537 68209 34537 0 32225 37527 0 31810 16824 96178 14285 300941 57614 1602 129470 61935 4068 114182 17609 152949 26099 179359 250368 4796 183349 125791 17414 61871 42058 0 2698 183297 23029 54464 0 26259 204595 35604 0 0 18437 29...
result:
Subtask #4:
score: 0
Runtime Error
Test #17:
score: 0
Runtime Error
input:
99999 99996 3 77889 1 10000000000 1 6316 86327 89644 386 3 9260 1 10000000000 2 2603 47234 69717 2 20260 73011 19290 2 62477 81233 26127 1 50140 68508 37004 98794 2 14449 22788 16063 1 43860 84932 50375 21777 1 67345 94584 28202 66610 2 661 68654 1 1 14411 94422 82738 61196 1 16563 94416 4920 38408 ...
output:
0 34602584 0 0 27739639583 1363823412 0 1902514434 1902514434 2147553884 1902514434 15794375094 0 4192446657 15797478185 13141921145 0 6351944090 5775183021 363222594 1995572111 2193350882 0 6843261316 5508935691 250667843 0 14181223499 7734049978 21958753162 12852564544 4496343819 15011219087 11331...
result:
Subtask #5:
score: 0
Skipped
Dependency #1:
100%
Accepted
Dependency #2:
0%