QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#669776 | #9242. An Easy Geometry Problem | Andycraft | WA | 316ms | 37412kb | C++20 | 4.7kb | 2024-10-23 19:38:07 | 2024-10-23 19:38:08 |
Judging History
answer
#define DEBUG 0
#include <iostream>
#include <algorithm>
#include <cassert>
#include <vector>
typedef long long LL;
typedef unsigned long long ULL;
template <class T> using Arr = std::vector<T>;
const int base = 998244353;
typedef long long LL;
#define MID ((l + r) >> 1)
#define LS (p << 1)
#define RS (LS | 1)
const long long MOD = 100055128505716009l;
// struct ULL {
// long long _;
// ULL(long long __) : _(__) {}
// ULL() : _(0) {}
// friend ULL operator+(ULL u, ULL v) { return (u._ + v._) % MOD; }
// friend ULL &operator+=(ULL &u, ULL v) { return u = ((u._ + v._) % MOD); }
// friend ULL &operator*=(ULL &u, ULL v) { return u = u * v; }
// friend ULL operator-(ULL u, ULL v) { return (u._ - v._) % MOD; }
// friend ULL operator*(ULL u, ULL v) { return (ULL)((__int128(u._) * v._) % MOD); }
// friend bool operator!=(ULL u, ULL v) { return (u._ % MOD + MOD) % MOD != (v._ % MOD + MOD) % MOD; }
// friend bool operator==(ULL u, ULL v) { return (u._ % MOD + MOD) % MOD == (v._ % MOD + MOD) % MOD; }
// };
const int MAXN = 200002;
const int B = 5e8;
int n, q, K, b;
int a[MAXN];
ULL pw[MAXN], ppw[MAXN];
ULL tr1[MAXN << 2], tag1[MAXN << 2], add1[MAXN << 2], hash1[MAXN];
ULL tr2[MAXN << 2], tag2[MAXN << 2], add2[MAXN << 2], hash2[MAXN];
void build(ULL tr[], ULL hash[], int p, int l, int r) {
if (l == r) {
tr[p] = hash[l];
return;
}
build(tr, hash, LS, l, MID);
build(tr, hash, RS, MID + 1, r);
}
ULL cons(int l, int r) {
return ppw[r - l] * pw[l];
}
ULL ask(ULL tr[], ULL tag[], ULL add[], int p, int l, int r, int ps) {
if (ps < l || ps > r)
return 0;
ULL delta = tag[p] * cons(l, ps) + add[p];
if (l == r)
return tr[p] + delta;
return delta + (ask(tr, tag, add, LS, l, MID, ps) + ask(tr, tag, add, RS, MID + 1, r, ps));
}
ULL ask_int(int l, int r, ULL tr[], ULL tag[], ULL add[]) {
ULL ls = ask(tr, tag, add, 1, 1, n, l - 1);
ULL rs = ask(tr, tag, add, 1, 1, n, r);
// printf("ask_int %d %d %llu %llu\n", l, r, ls, rs);
return rs - ls;
}
void add_add(ULL add[], int p, int l, int r, int lp, int rp, ULL v) {
if (rp < l || lp > r)
return;
if (lp <= l && r <= rp) {
add[p] += v;
return;
}
add_add(add, LS, l, MID, lp, rp, v);
add_add(add, RS, MID + 1, r, lp, rp, v);
}
void add_tag(ULL add[], ULL tag[], int p, int l, int r, int lp, int rp, ULL v) {
if (rp < l || lp > r)
return;
if (lp <= l && r <= rp) {
tag[p] += v;
return;
}
add_tag(add, tag, LS, l, MID, lp, rp, v);
add_tag(add, tag, RS, MID + 1, r, lp, rp, v);
if (rp >= MID + 1 && lp <= MID)
add_add(add, RS, MID + 1, r, MID + 1, rp, v * cons(lp, MID));
}
int main() {
std::ios::sync_with_stdio(false);
std::cin >> n >> q >> K >> b;
for (int i = 1; i <= n; ++i) {
std::cin >> a[i];
a[i] = a[i] * 2 - K * i;
a[i] += B;
assert(a[i] > 0);
}
pw[0] = 1;
ppw[0] = 1;
for (int i = 1; i < MAXN; ++i) {
pw[i] = pw[i - 1] * base;
ppw[i] = ppw[i - 1] + pw[i];
}
for (int i = 1; i <= n; ++i)
hash1[i] = hash1[i - 1] + pw[i] * (a[i] + 2 * b);
for (int i = 1; i <= n; ++i)
hash2[i] = hash2[i - 1] + pw[i] * a[n - i + 1];
build(tr1, hash1, 1, 1, n);
build(tr2, hash2, 1, 1, n);
// puts("HELLO");
#if DEBUG
for (int i = 1; i <= n; ++i)
printf("%llu ", hash1[i]);
puts("");
printf("! %llu %llu\n", hash1[2] * pw[5], (hash2[6] - hash2[4]) * pw[1]);
for (int i = 1; i <= n; ++i)
printf("%llu ", hash2[i]);
puts("");
for (int i = 1; i <= n; ++i)
printf("%llu ", ask(tr1, tag1, add1, 1, 1, n, i));
puts("");
#endif
for (int i = 1; i <= q; ++i) {
int op;
std::cin >> op;
if (op == 1) {
int l, r, v;
std::cin >> l >> r >> v;
v *= 2;
add_tag(add1, tag1, 1, 1, n, l, r, v);
add_tag(add2, tag2, 1, 1, n, n - r + 1, n - l + 1, v);
add_add(add1, 1, 1, n, r + 1, n, v * cons(l, r));
add_add(add2, 1, 1, n, n - l + 2, n, v * cons(n - r + 1, n - l + 1));
} else {
int x;
std::cin >> x;
int h = 0;
for (; x - (1 << h) >= 1 && x + (1 << h) <= n; ++h) {
int lp = x - (1 << h), rp = x + (1 << h);
ULL l = ask_int(lp, x - 1, tr1, tag1, add1);
ULL r = ask_int(n - rp + 1, n - (x + 1) + 1, tr2, tag2, add2);
// printf("== %llu %llu\n", l, r);
// l += 2 * b * cons(lp, x - 1);
l *= pw[n - rp + 1];
r *= pw[lp];
// printf("= %llu %llu\n", l, r);
if (l != r)
break;
}
if (h == 0) {
puts("0");
continue;
}
int ans = 1 << (h - 1);
for (--h; h >= 0; --h) {
int lp = x - (ans + (1 << h));
int rp = x + (ans + (1 << h));
if (lp < 1 || rp > n)
continue;
ULL l = ask_int(lp, x - 1, tr1, tag1, add1);
ULL r = ask_int(n - rp + 1, n - (x + 1) + 1, tr2, tag2, add2);
// l += 2 * b * cons(lp, x - 1);
l *= pw[n - rp + 1];
r *= pw[lp];
// printf("# %llu %llu\n", l, r);
if (l == r)
ans += (1 << h);
}
printf("%d\n", ans);
}
}
return 0;
}
詳細信息
Test #1:
score: 100
Accepted
time: 2ms
memory: 12960kb
input:
6 6 6 2 1 5 9 10 15 18 2 2 1 3 3 -3 2 2 1 3 4 3 2 3 2 4
output:
1 0 2 0
result:
ok 4 number(s): "1 0 2 0"
Test #2:
score: 0
Accepted
time: 6ms
memory: 13000kb
input:
5000 5000 2 0 -329 -328 -327 -326 -325 -324 -323 -322 -321 -320 -319 -318 -317 -316 -315 -314 -313 -312 -311 -310 -309 -308 -307 -306 -305 -304 -303 -302 -301 -300 -299 -298 -297 -296 -295 -294 -293 -292 -291 -290 -289 -288 -287 -286 -285 -284 -283 -282 -281 -280 -279 -278 -277 -276 -275 -274 -273 -...
output:
2 304 73 29 61 292 139 48 17 99 6 5 53 93 3 91 65 29 33 306 21 24 17 21 281 12 16 1 33 7 18 96 7 40 39 13 7 46 43 16 1 72 33 16 22 5 6 189 27 1 35 107 43 34 3 27 20 21 44 56 96 36 2 27 22 30 32 6 5 105 27 37 12 58 2 21 154 17 110 57 3 7 33 15 24 94 68 25 1 14 10 4 10 2 25 39 36 33 164 11 19 181 11 3...
result:
ok 3337 numbers
Test #3:
score: 0
Accepted
time: 11ms
memory: 20304kb
input:
5000 5000 2 0 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 86...
output:
362 82 14 234 140 5 44 136 22 43 29 96 59 23 25 61 193 22 39 39 23 53 48 76 100 58 120 24 12 106 32 48 73 63 116 16 136 10 28 15 84 30 65 1 54 15 16 70 1 95 74 14 17 20 36 254 22 29 70 172 106 2 25 8 98 35 169 16 2 2 99 10 36 40 3 69 272 170 219 12 79 26 78 100 10 167 140 70 34 17 23 21 55 10 6 17 6...
result:
ok 3313 numbers
Test #4:
score: 0
Accepted
time: 10ms
memory: 12564kb
input:
5000 5000 2 0 -456 -455 -454 -453 -452 -451 -450 -449 -448 -447 -446 -445 -444 -443 -442 -441 -440 -439 -438 -437 -436 -435 -434 -433 -432 -431 -430 -429 -428 -427 -426 -425 -424 -423 -422 -421 -420 -419 -418 -417 -416 -415 -414 -413 -412 -411 -410 -409 -408 -407 -406 -405 -404 -403 -402 -401 -400 -...
output:
8 75 80 408 385 73 37 402 338 43 11 163 3 7 80 0 339 47 384 8 10 47 162 307 30 28 36 14 27 126 271 151 4 11 11 9 92 154 2 15 28 160 205 12 59 79 114 23 22 141 7 12 31 42 120 0 34 2 167 157 76 32 20 298 47 104 76 33 49 34 1 40 16 1 28 7 4 55 14 8 68 17 7 117 1 14 14 80 44 8 45 49 65 15 49 56 50 40 14...
result:
ok 3296 numbers
Test #5:
score: -100
Wrong Answer
time: 316ms
memory: 37412kb
input:
200000 199999 -195 -119 -267 -146 191 -456 835 265 -226 -264 160 -101 739 -988 -967 890 -753 -854 514 491 -733 662 681 -362 804 -714 -1000 -790 931 -450 212 94 239 638 400 -167 -360 18 606 256 445 695 -509 643 -892 213 -32 42 400 733 -667 -986 225 493 -699 547 409 -35 394 920 -163 -908 -576 921 -997...
output:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...
result:
wrong answer 47324th numbers differ - expected: '1', found: '0'