/tmp/solutions/build/predecessor_problem-main.cpp:
1| |#include <common.h>
2| |prelude;
3| |
4| |namespace {
5| |
6| |u64 root;
7| |u64 lv18[39];
8| |u64 lv12[2442];
9| |u64 leaf[156250];
10| |
11| 0|void set(u32 x) {
12| 0| leaf[x >> 6] |= 1ull << (x & 63);
13| 0| lv12[x >> 12] |= 1ull << ((x >> 6) & 63);
14| 0| lv18[x >> 18] |= 1ull << ((x >> 12) & 63);
15| 0| root |= 1ull << (x >> 18);
16| 0|}
17| |
18| 0|void unset(u32 x) {
19| 0| if (!(leaf[x >> 6] &= ~(1ull << (x & 63))))
------------------
| Branch (19:7): [True: 0.00%, False: 0.00%]
------------------
20| 0| if (!(lv12[x >> 12] &= ~(1ull << ((x >> 6) & 63))))
------------------
| Branch (20:9): [True: 0.00%, False: 0.00%]
------------------
21| 0| if (!(lv18[x >> 18] &= ~(1ull << ((x >> 12) & 63))))
------------------
| Branch (21:11): [True: 0.00%, False: 0.00%]
------------------
22| 0| root &= ~(1ull << (x >> 18));
23| 0|}
24| |
25| 0|bool get(u32 x) { return leaf[x >> 6] & (1ull << (x & 63)); }
26| |
27| 499k|u32 nxt(u32 x) {
28| 499k| if (u64 tmp = leaf[x >> 6] & (-1ull << (x & 63)); !tmp) {
------------------
| Branch (28:53): [True: 100.00%, False: 0.00%]
------------------
29| 499k| if (u64 tmp = lv12[x >> 12] & (-2ull << ((x >> 6) & 63)); !tmp) {
------------------
| Branch (29:63): [True: 99.80%, False: 0.20%]
------------------
30| 498k| if (u64 tmp = lv18[x >> 18] & (-2ull << ((x >> 12) & 63)); !tmp) {
------------------
| Branch (30:66): [True: 91.29%, False: 8.71%]
------------------
31| 455k| if (u64 tmp = root & (-2ull << (x >> 18)); !tmp) {
------------------
| Branch (31:52): [True: 0.79%, False: 99.21%]
------------------
32| 3.61k| return -1;
33| 451k| } else {
34| 451k| u32 ans = __builtin_ctzll(tmp);
35| 451k| ans = ans << 6 | __builtin_ctzll(lv18[ans]);
36| 451k| ans = ans << 6 | __builtin_ctzll(lv12[ans]);
37| 451k| ans = ans << 6 | __builtin_ctzll(leaf[ans]);
38| 451k| return ans;
39| 451k| }
40| 455k| } else {
41| 43.4k| u32 ans = (x >> 18) << 6 | __builtin_ctzll(tmp);
42| 43.4k| ans = ans << 6 | __builtin_ctzll(lv12[ans]);
43| 43.4k| ans = ans << 6 | __builtin_ctzll(leaf[ans]);
44| 43.4k| return ans;
45| 43.4k| }
46| 498k| } else {
47| 1.02k| u32 ans = (x >> 12) << 6 | __builtin_ctzll(tmp);
48| 1.02k| ans = ans << 6 | __builtin_ctzll(leaf[ans]);
49| 1.02k| return ans;
50| 1.02k| }
51| 499k| } else {
52| 13| u32 ans = (x >> 6) << 6 | __builtin_ctzll(tmp);
53| 13| return ans;
54| 13| }
55| 499k|}
56| |
57| 500k|u32 pre(u32 x) {
58| 500k| if (u64 tmp = leaf[x >> 6] & ~(-2ull << (x & 63)); !tmp) {
------------------
| Branch (58:54): [True: 100.00%, False: 0.00%]
------------------
59| 500k| if (u64 tmp = lv12[x >> 12] & ~(-1ull << ((x >> 6) & 63)); !tmp) {
------------------
| Branch (59:64): [True: 99.93%, False: 0.07%]
------------------
60| 500k| if (u64 tmp = lv18[x >> 18] & ~(-1ull << ((x >> 12) & 63)); !tmp) {
------------------
| Branch (60:67): [True: 90.68%, False: 9.32%]
------------------
61| 453k| if (u64 tmp = root & ~(-1ull << (x >> 18)); !tmp) {
------------------
| Branch (61:53): [True: 20.42%, False: 79.58%]
------------------
62| 92.5k| return -1;
63| 360k| } else {
64| 360k| u32 ans = 63 - __builtin_clzll(tmp);
65| 360k| ans = ans << 6 | 63 - __builtin_clzll(lv18[ans]);
66| 360k| ans = ans << 6 | 63 - __builtin_clzll(lv12[ans]);
67| 360k| ans = ans << 6 | 63 - __builtin_clzll(leaf[ans]);
68| 360k| return ans;
69| 360k| }
70| 453k| } else {
71| 46.5k| u32 ans = (x >> 18) << 6 | 63 - __builtin_clzll(tmp);
72| 46.5k| ans = ans << 6 | 63 - __builtin_clzll(lv12[ans]);
73| 46.5k| ans = ans << 6 | 63 - __builtin_clzll(leaf[ans]);
74| 46.5k| return ans;
75| 46.5k| }
76| 500k| } else {
77| 370| u32 ans = (x >> 12) << 6 | 63 - __builtin_clzll(tmp);
78| 370| ans = ans << 6 | 63 - __builtin_clzll(leaf[ans]);
79| 370| return ans;
80| 370| }
81| 500k| } else {
82| 16| u32 ans = (x >> 6) << 6 | 63 - __builtin_clzll(tmp);
83| 16| return ans;
84| 16| }
85| 500k|}
86| |
87| |} // namespace
88| |
89| 1|int main() {
90| 1| rd rd;
91| 1| wt wt;
92| 1| int n = rd.uh();
93| 1| int q = rd.uh();
94| 1|#ifdef LOCAL
95| 1| root = 0;
96| 1| std::memset(lv18, 0, sizeof(lv18));
97| 1| std::memset(lv12, 0, sizeof(lv12));
98| 1| std::memset(leaf, 0, sizeof(leaf));
99| 1|#endif
100| 1| int i = 0;
101| 156k| for (; i + 64 < n; i += 64) {
^156k
------------------
| Branch (101:10): [True: 100.00%, False: 0.00%]
------------------
102| 156k| u32 a = _mm256_movemask_epi8(_mm256_cmpeq_epi8(
103| 156k| _mm256_loadu_si256(reinterpret_cast<const __m256i_u *>(rd.p + i)),
104| 156k| _mm256_set1_epi8('1')));
105| 156k| u32 b = _mm256_movemask_epi8(_mm256_cmpeq_epi8(
106| 156k| _mm256_loadu_si256(reinterpret_cast<const __m256i_u *>(rd.p + i + 32)),
107| 156k| _mm256_set1_epi8('1')));
108| 156k| leaf[i >> 6] = a | u64(b) << 32;
109| 156k| }
110| 65| for (; i < n; ++i) leaf[i >> 6] |= u64(rd.p[i] - '0') << (i & 63);
^64 ^64
------------------
| Branch (110:10): [True: 98.46%, False: 1.54%]
------------------
111| 156k| for (int i = 0; i <= (n - 1) / 64; ++i)
^156k
------------------
| Branch (111:19): [True: 100.00%, False: 0.00%]
------------------
112| 156k| lv12[i >> 6] |= u64(!!leaf[i]) << (i & 63);
113| 2.44k| for (int i = 0; i <= (n - 1) / 64 / 64; ++i)
^2.44k
------------------
| Branch (113:19): [True: 99.96%, False: 0.04%]
------------------
114| 2.44k| lv18[i >> 6] |= u64(!!lv12[i]) << (i & 63);
115| 40| for (int i = 0; i <= (n - 1) / 64 / 64 / 64; ++i)
^39
------------------
| Branch (115:19): [True: 97.50%, False: 2.50%]
------------------
116| 39| root |= u64(!!lv18[i]) << (i & 63);
117| 1| rd.p += n + 1;
118| 1.00M| while (q--) {
------------------
| Branch (118:10): [True: 100.00%, False: 0.00%]
------------------
119| 1.00M| let t = rd.u1();
120| 1.00M| let k = rd.uh();
121| 1.00M| if (t == 0) {
------------------
| Branch (121:9): [True: 0.00%, False: 100.00%]
------------------
122| 0| set(k);
123| 0| }
124| 1.00M| if (t == 1) {
------------------
| Branch (124:9): [True: 0.00%, False: 100.00%]
------------------
125| 0| unset(k);
126| 0| }
127| 1.00M| if (t == 2) {
------------------
| Branch (127:9): [True: 0.00%, False: 100.00%]
------------------
128| 0| let x = get(k);
129| 0| wt.uw(x);
130| 0| }
131| 1.00M| if (t == 3) {
------------------
| Branch (131:9): [True: 49.95%, False: 50.05%]
------------------
132| 499k| let x = nxt(k);
133| 499k| if (~x) {
------------------
| Branch (133:11): [True: 99.28%, False: 0.72%]
------------------
134| 495k| wt.uw(x);
135| 495k| } else {
136| 3.61k| wt.puts(" -1", 3);
137| 3.61k| }
138| 499k| }
139| 1.00M| if (t == 4) {
------------------
| Branch (139:9): [True: 50.05%, False: 49.95%]
------------------
140| 500k| let x = pre(k);
141| 500k| if (~x) {
------------------
| Branch (141:11): [True: 81.50%, False: 18.50%]
------------------
142| 407k| wt.uw(x);
143| 407k| } else {
144| 92.5k| wt.puts(" -1", 3);
145| 92.5k| }
146| 500k| }
147| 1.00M| }
148| 1| return 0;
149| 1|}