/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.85%, False: 0.15%]
------------------
30| 499k| if (u64 tmp = lv18[x >> 18] & (-2ull << ((x >> 12) & 63)); !tmp) {
------------------
| Branch (30:66): [True: 92.74%, False: 7.26%]
------------------
31| 462k| if (u64 tmp = root & (-2ull << (x >> 18)); !tmp) {
------------------
| Branch (31:52): [True: 28.69%, False: 71.31%]
------------------
32| 132k| return -1;
33| 330k| } else {
34| 330k| u32 ans = __builtin_ctzll(tmp);
35| 330k| ans = ans << 6 | __builtin_ctzll(lv18[ans]);
36| 330k| ans = ans << 6 | __builtin_ctzll(lv12[ans]);
37| 330k| ans = ans << 6 | __builtin_ctzll(leaf[ans]);
38| 330k| return ans;
39| 330k| }
40| 462k| } else {
41| 36.2k| u32 ans = (x >> 18) << 6 | __builtin_ctzll(tmp);
42| 36.2k| ans = ans << 6 | __builtin_ctzll(lv12[ans]);
43| 36.2k| ans = ans << 6 | __builtin_ctzll(leaf[ans]);
44| 36.2k| return ans;
45| 36.2k| }
46| 499k| } else {
47| 738| u32 ans = (x >> 12) << 6 | __builtin_ctzll(tmp);
48| 738| ans = ans << 6 | __builtin_ctzll(leaf[ans]);
49| 738| return ans;
50| 738| }
51| 499k| } else {
52| 7| u32 ans = (x >> 6) << 6 | __builtin_ctzll(tmp);
53| 7| return ans;
54| 7| }
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.86%, False: 0.14%]
------------------
60| 499k| if (u64 tmp = lv18[x >> 18] & ~(-1ull << ((x >> 12) & 63)); !tmp) {
------------------
| Branch (60:67): [True: 89.13%, False: 10.87%]
------------------
61| 445k| if (u64 tmp = root & ~(-1ull << (x >> 18)); !tmp) {
------------------
| Branch (61:53): [True: 4.21%, False: 95.79%]
------------------
62| 18.7k| return -1;
63| 426k| } else {
64| 426k| u32 ans = 63 - __builtin_clzll(tmp);
65| 426k| ans = ans << 6 | 63 - __builtin_clzll(lv18[ans]);
66| 426k| ans = ans << 6 | 63 - __builtin_clzll(lv12[ans]);
67| 426k| ans = ans << 6 | 63 - __builtin_clzll(leaf[ans]);
68| 426k| return ans;
69| 426k| }
70| 445k| } else {
71| 54.2k| u32 ans = (x >> 18) << 6 | 63 - __builtin_clzll(tmp);
72| 54.2k| ans = ans << 6 | 63 - __builtin_clzll(lv12[ans]);
73| 54.2k| ans = ans << 6 | 63 - __builtin_clzll(leaf[ans]);
74| 54.2k| return ans;
75| 54.2k| }
76| 499k| } else {
77| 710| u32 ans = (x >> 12) << 6 | 63 - __builtin_clzll(tmp);
78| 710| ans = ans << 6 | 63 - __builtin_clzll(leaf[ans]);
79| 710| return ans;
80| 710| }
81| 500k| } else {
82| 17| u32 ans = (x >> 6) << 6 | 63 - __builtin_clzll(tmp);
83| 17| return ans;
84| 17| }
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.98%, False: 50.02%]
------------------
132| 499k| let x = nxt(k);
133| 499k| if (~x) {
------------------
| Branch (133:11): [True: 73.44%, False: 26.56%]
------------------
134| 367k| wt.uw(x);
135| 367k| } else {
136| 132k| wt.puts(" -1", 3);
137| 132k| }
138| 499k| }
139| 1.00M| if (t == 4) {
------------------
| Branch (139:9): [True: 50.02%, False: 49.98%]
------------------
140| 500k| let x = pre(k);
141| 500k| if (~x) {
------------------
| Branch (141:11): [True: 96.26%, False: 3.74%]
------------------
142| 481k| wt.uw(x);
143| 481k| } else {
144| 18.7k| wt.puts(" -1", 3);
145| 18.7k| }
146| 500k| }
147| 1.00M| }
148| 1| return 0;
149| 1|}