QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#326548 | #8229. 栈 | Xiaohuba | 0 | 110ms | 107696kb | C++23 | 9.3kb | 2024-02-13 13:34:25 | 2024-02-13 13:34:25 |
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 {
// constexpr ll MAXN = ...;
class WBLT {
static il constexpr db alpha = 1.0 - 1.4142136 / 2;
static il constexpr ll INITIAL_SIZE = 3e6;
struct Node {
int lc, rc, sz;
ll sum, weight;
Node() : lc(0), rc(0), sz(0), sum(0), weight(0) {}
};
static il vector<Node> T{INITIAL_SIZE};
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<1>(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;
}
}
}
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 << "> " << 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;
ll y, z, w;
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));
}
// WBLT tr;
// tr.push(2, 3);
// cout << tr.size() << ' ' << tr.qry_sum(2, 4) << '\n';
}
} // namespace
signed main() { return Main(), 0; }
/*
[1 , 4]
[1, 2][3, 4]
[1][2][3][4]
*/
Details
Tip: Click on the bar to expand more detailed information
Subtask #1:
score: 0
Wrong Answer
Test #1:
score: 0
Wrong Answer
time: 12ms
memory: 107572kb
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 203015110 200648623 98486697 214102593 123945 0 61782451 0 58379 58379 762429740 321140700 638060258 0 4594954595 58379 0 2678303420 209997360 838785675 0 3948193587 1254619275 0 47775 118545795 2323313145 515151640 84280112 0 3057798043 1276548448 926862294 47775 31720 135678...
result:
wrong answer 4th numbers differ - expected: '471569175', found: '203015110'
Subtask #2:
score: 0
Memory Limit Exceeded
Test #6:
score: 0
Memory Limit Exceeded
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
Wrong Answer
Test #12:
score: 0
Wrong Answer
time: 110ms
memory: 107696kb
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 0 0 75705 24537 34599 0 1 32225 37527 1 0 0 20502 14285 130253 57614 1602 0 8867 4068 65762 0 152949 26099 47500 135190 4796 20259 11098 17414 61871 42058 1 2698 183297 23029 52052 0 26259 0 0 1 1 0 0 108 1 0 6278 0 30138 17936 0 12422 36659 34954 5...
result:
wrong answer 12th numbers differ - expected: '16405', found: '0'
Subtask #4:
score: 0
Wrong Answer
Test #17:
score: 0
Wrong Answer
time: 105ms
memory: 107560kb
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 358059420 1608884520 1608884520 1015595945 1608884520 489109468 293629914 4192446657 1365740406 3353506353 4192446657 5216355443 1308249668 77203035 77203035 77203035 0 6843261316 1928348556 2725517909 77203035 411945666 7138242390 6542442770 12852564544 4496343...
result:
wrong answer 7th numbers differ - expected: '0', found: '358059420'
Subtask #5:
score: 0
Skipped
Dependency #1:
0%