代码

#include <stdint.h>
#include <sys/auxv.h>
#include <unistd.h>

struct DuckInfo {
  uint64_t abi_version;
  const char *stdin_ptr; uint64_t stdin_size;
  char *stdout_ptr; uint64_t stdout_limit; uint64_t stdout_size;
  char *stderr_ptr; uint64_t stderr_limit; uint64_t stderr_size;
  const char *IB_ptr; uint64_t IB_limit;
  char *OB_ptr; uint64_t OB_limit;
  uint64_t tsc_frequency;
} __attribute__((packed));

#define NL (1<<20)
#define N1 (NL>>6)   /* 16384 */
#define N2 (N1>>6)   /* 256 */
#define N3 (N2>>6)   /* 4 */
#define MASK 0x3FFFFFFFu
#define FULL MASK

static uint32_t limb[NL];
static uint64_t fm1[N1], nm1[N1];
static uint64_t fm2[N2], nm2[N2];
static uint64_t fm3[N3], nm3[N3];
static uint64_t fm4, nm4;

static const char *P;
static uint64_t O;

static inline int ctz(uint64_t x){ return __builtin_ctzll(x); }

/* recompute ancestors of L1 word l1 after fm1[l1]/nm1[l1] changed */
static inline void bubble_l1(int l1){
    int l2 = l1 >> 6;
    uint64_t b = 1ull << (l1 & 63);
    if (fm1[l1] == ~0ull) fm2[l2] |= b; else fm2[l2] &= ~b;
    if (nm1[l1] != 0) nm2[l2] |= b; else nm2[l2] &= ~b;
    int l3 = l2 >> 6;
    b = 1ull << (l2 & 63);
    if (fm2[l2] == ~0ull) fm3[l3] |= b; else fm3[l3] &= ~b;
    if (nm2[l2] != 0) nm3[l3] |= b; else nm3[l3] &= ~b;
    b = 1ull << l3;
    if (fm3[l3] == ~0ull) fm4 |= b; else fm4 &= ~b;
    if (nm3[l3] != 0) nm4 |= b; else nm4 &= ~b;
}

/* materialize L1 word l1 (make fm1[l1]/nm1[l1] authoritative) */
static inline void materialize_word(int l1){
    int l2 = l1 >> 6, l3 = l2 >> 6;
    uint64_t b3 = 1ull << l3;
    if (fm4 & b3){ fm3[l3] = ~0ull; nm3[l3] = ~0ull; }
    else if (!(nm4 & b3)){ fm3[l3] = 0; nm3[l3] = 0; }
    uint64_t b2 = 1ull << (l2 & 63);
    if (fm3[l3] & b2){ fm2[l2] = ~0ull; nm2[l2] = ~0ull; }
    else if (!(nm3[l3] & b2)){ fm2[l2] = 0; nm2[l2] = 0; }
    uint64_t b1 = 1ull << (l1 & 63);
    if (fm2[l2] & b1){ fm1[l1] = ~0ull; nm1[l1] = ~0ull; }
    else if (!(nm2[l2] & b1)){ fm1[l1] = 0; nm1[l1] = 0; }
}

/* add v (< 2^30) to limb pos; returns carry (0/1) */
static inline int point_add(int pos, uint32_t v){
    int l1 = pos >> 6;
    materialize_word(l1);
    uint64_t b = 1ull << (pos & 63);
    uint32_t old;
    if (fm1[l1] & b) old = FULL;
    else if (!(nm1[l1] & b)) old = 0;
    else old = limb[pos];
    uint32_t nv = old + v;
    int ovf = nv >> 30;
    uint32_t n30 = nv & MASK;
    limb[pos] = n30;
    if (n30 == FULL) fm1[l1] |= b; else fm1[l1] &= ~b;
    if (n30 != 0) nm1[l1] |= b; else nm1[l1] &= ~b;
    bubble_l1(l1);
    return ovf;
}

/* subtract v (< 2^30) from limb pos; returns borrow (0/1) */
static inline int point_sub(int pos, uint32_t v){
    int l1 = pos >> 6;
    materialize_word(l1);
    uint64_t b = 1ull << (pos & 63);
    uint32_t old;
    if (fm1[l1] & b) old = FULL;
    else if (!(nm1[l1] & b)) old = 0;
    else old = limb[pos];
    int32_t nv = (int32_t)old - (int32_t)v;
    int borrow = nv < 0;
    uint32_t res = borrow ? (uint32_t)(nv + (1<<30)) : (uint32_t)nv;
    limb[pos] = res;
    if (res == FULL) fm1[l1] |= b; else fm1[l1] &= ~b;
    if (res != 0) nm1[l1] |= b; else nm1[l1] &= ~b;
    bubble_l1(l1);
    return borrow;
}

/* first limb >= p that is not full */
static inline int find_nonfull(int p){
    int l1 = p >> 6, l2 = p >> 12, l3 = p >> 18, off = p & 63;
    materialize_word(l1);
    uint64_t m1 = ~fm1[l1] >> off;
    if (m1) return (l1 << 6) + off + ctz(m1);
    int sh2 = (l1 & 63) + 1;
    uint64_t m2 = (sh2 == 64) ? 0 : (~fm2[l2] >> sh2);
    if (m2){
        int l1n = (l2 << 6) + sh2 + ctz(m2);
        uint64_t b = 1ull << (l1n & 63);
        if (fm2[l2] & b){ fm1[l1n] = ~0ull; nm1[l1n] = ~0ull; }
        else if (!(nm2[l2] & b)){ fm1[l1n] = 0; nm1[l1n] = 0; }
        return (l1n << 6) + ctz(~fm1[l1n]);
    }
    int sh3 = (l2 & 63) + 1;
    uint64_t m3 = (sh3 == 64) ? 0 : (~fm3[l3] >> sh3);
    if (m3){
        int l2n = (l3 << 6) + sh3 + ctz(m3);
        uint64_t b2 = 1ull << (l2n & 63);
        if (fm3[l3] & b2){ fm2[l2n] = ~0ull; nm2[l2n] = ~0ull; }
        else if (!(nm3[l3] & b2)){ fm2[l2n] = 0; nm2[l2n] = 0; }
        int r2 = ctz(~fm2[l2n]);
        int l1n = (l2n << 6) + r2;
        uint64_t b1 = 1ull << (l1n & 63);
        if (fm2[l2n] & b1){ fm1[l1n] = ~0ull; nm1[l1n] = ~0ull; }
        else if (!(nm2[l2n] & b1)){ fm1[l1n] = 0; nm1[l1n] = 0; }
        return (l1n << 6) + ctz(~fm1[l1n]);
    }
    uint64_t m4 = ~fm4 >> (l3 + 1);
    if (m4){
        int l3n = l3 + 1 + ctz(m4);
        uint64_t b3 = 1ull << l3n;
        if (fm4 & b3){ fm3[l3n] = ~0ull; nm3[l3n] = ~0ull; }
        else if (!(nm4 & b3)){ fm3[l3n] = 0; nm3[l3n] = 0; }
        int r3 = ctz(~fm3[l3n]);
        int l2n = (l3n << 6) + r3;
        uint64_t b2 = 1ull << (l2n & 63);
        if (fm3[l3n] & b2){ fm2[l2n] = ~0ull; nm2[l2n] = ~0ull; }
        else if (!(nm3[l3n] & b2)){ fm2[l2n] = 0; nm2[l2n] = 0; }
        int r2 = ctz(~fm2[l2n]);
        int l1n = (l2n << 6) + r2;
        uint64_t b1 = 1ull << (l1n & 63);
        if (fm2[l2n] & b1){ fm1[l1n] = ~0ull; nm1[l1n] = ~0ull; }
        else if (!(nm2[l2n] & b1)){ fm1[l1n] = 0; nm1[l1n] = 0; }
        return (l1n << 6) + ctz(~fm1[l1n]);
    }
    return NL;
}

/* first limb >= p that is not zero */
static inline int find_nonzero(int p){
    int l1 = p >> 6, l2 = p >> 12, l3 = p >> 18, off = p & 63;
    materialize_word(l1);
    uint64_t m1 = nm1[l1] >> off;
    if (m1) return (l1 << 6) + off + ctz(m1);
    int sh2 = (l1 & 63) + 1;
    uint64_t m2 = (sh2 == 64) ? 0 : (nm2[l2] >> sh2);
    if (m2){
        int l1n = (l2 << 6) + sh2 + ctz(m2);
        uint64_t b = 1ull << (l1n & 63);
        if (fm2[l2] & b){ fm1[l1n] = ~0ull; nm1[l1n] = ~0ull; }
        else if (!(nm2[l2] & b)){ fm1[l1n] = 0; nm1[l1n] = 0; }
        return (l1n << 6) + ctz(nm1[l1n]);
    }
    int sh3 = (l2 & 63) + 1;
    uint64_t m3 = (sh3 == 64) ? 0 : (nm3[l3] >> sh3);
    if (m3){
        int l2n = (l3 << 6) + sh3 + ctz(m3);
        uint64_t b2 = 1ull << (l2n & 63);
        if (fm3[l3] & b2){ fm2[l2n] = ~0ull; nm2[l2n] = ~0ull; }
        else if (!(nm3[l3] & b2)){ fm2[l2n] = 0; nm2[l2n] = 0; }
        int r2 = ctz(nm2[l2n]);
        int l1n = (l2n << 6) + r2;
        uint64_t b1 = 1ull << (l1n & 63);
        if (fm2[l2n] & b1){ fm1[l1n] = ~0ull; nm1[l1n] = ~0ull; }
        else if (!(nm2[l2n] & b1)){ fm1[l1n] = 0; nm1[l1n] = 0; }
        return (l1n << 6) + ctz(nm1[l1n]);
    }
    uint64_t m4 = nm4 >> (l3 + 1);
    if (m4){
        int l3n = l3 + 1 + ctz(m4);
        uint64_t b3 = 1ull << l3n;
        if (fm4 & b3){ fm3[l3n] = ~0ull; nm3[l3n] = ~0ull; }
        else if (!(nm4 & b3)){ fm3[l3n] = 0; nm3[l3n] = 0; }
        int r3 = ctz(nm3[l3n]);
        int l2n = (l3n << 6) + r3;
        uint64_t b2 = 1ull << (l2n & 63);
        if (fm3[l3n] & b2){ fm2[l2n] = ~0ull; nm2[l2n] = ~0ull; }
        else if (!(nm3[l3n] & b2)){ fm2[l2n] = 0; nm2[l2n] = 0; }
        int r2 = ctz(nm2[l2n]);
        int l1n = (l2n << 6) + r2;
        uint64_t b1 = 1ull << (l1n & 63);
        if (fm2[l2n] & b1){ fm1[l1n] = ~0ull; nm1[l1n] = ~0ull; }
        else if (!(nm2[l2n] & b1)){ fm1[l1n] = 0; nm1[l1n] = 0; }
        return (l1n << 6) + ctz(nm1[l1n]);
    }
    return NL;
}

/* set limbs [l, r) to zero (they are all full) */
static inline void zero_range(int l, int r){
    int l1a = l >> 6, l1b = (r - 1) >> 6;
    if (l1a == l1b){
        materialize_word(l1a);
        int lo = l & 63, hi = (r - 1) & 63;
        uint64_t mask = (lo == 0 && hi == 63) ? ~0ull : ((1ull << (hi - lo + 1)) - 1) << lo;
        fm1[l1a] &= ~mask;
        nm1[l1a] &= ~mask;
        bubble_l1(l1a);
    } else {
        materialize_word(l1a);
        uint64_t mask = ~0ull << (l & 63);
        fm1[l1a] &= ~mask;
        nm1[l1a] &= ~mask;
        bubble_l1(l1a);
        for (int b = l1a + 1; b < l1b; b++){
            fm1[b] = 0; nm1[b] = 0;
            bubble_l1(b);
        }
        materialize_word(l1b);
        int hi = (r - 1) & 63;
        mask = (hi == 63) ? ~0ull : ((1ull << (hi + 1)) - 1);
        fm1[l1b] &= ~mask;
        nm1[l1b] &= ~mask;
        bubble_l1(l1b);
    }
}

/* set limbs [l, r) to full (they are all zero) */
static inline void full_range(int l, int r){
    int l1a = l >> 6, l1b = (r - 1) >> 6;
    if (l1a == l1b){
        materialize_word(l1a);
        int lo = l & 63, hi = (r - 1) & 63;
        uint64_t mask = (lo == 0 && hi == 63) ? ~0ull : ((1ull << (hi - lo + 1)) - 1) << lo;
        fm1[l1a] |= mask;
        nm1[l1a] |= mask;
        bubble_l1(l1a);
    } else {
        materialize_word(l1a);
        uint64_t mask = ~0ull << (l & 63);
        fm1[l1a] |= mask;
        nm1[l1a] |= mask;
        bubble_l1(l1a);
        for (int b = l1a + 1; b < l1b; b++){
            fm1[b] = ~0ull; nm1[b] = ~0ull;
            bubble_l1(b);
        }
        materialize_word(l1b);
        int hi = (r - 1) & 63;
        mask = (hi == 63) ? ~0ull : ((1ull << (hi + 1)) - 1);
        fm1[l1b] |= mask;
        nm1[l1b] |= mask;
        bubble_l1(l1b);
    }
}

static inline void add_val(int q, int r, uint32_t a){
    uint32_t lo = (uint32_t)((uint64_t)a << r) & MASK;
    uint32_t hi = (r == 0) ? 0u : (uint32_t)(a >> (30 - r));
    int c0 = point_add(q, lo);
    int c1 = point_add(q + 1, hi + c0);
    if (c1){
        int start = q + 2;
        int p = find_nonfull(start);
        if (p > start) zero_range(start, p);
        point_add(p, 1);
    }
}

static inline void sub_val(int q, int r, uint32_t a){
    uint32_t lo = (uint32_t)((uint64_t)a << r) & MASK;
    uint32_t hi = (r == 0) ? 0u : (uint32_t)(a >> (30 - r));
    int b0 = point_sub(q, lo);
    int b1 = point_sub(q + 1, hi + b0);
    if (b1){
        int start = q + 2;
        int p = find_nonzero(start);
        if (p > start) full_range(start, p);
        point_sub(p, 1);
    }
}

static inline uint32_t read_limb(int pos){
    int l1 = pos >> 6, l2 = pos >> 12, l3 = pos >> 18;
    uint64_t b3 = 1ull << l3;
    if (fm4 & b3) return FULL;
    if (!(nm4 & b3)) return 0;
    uint64_t b2 = 1ull << (l2 & 63);
    if (fm3[l3] & b2) return FULL;
    if (!(nm3[l3] & b2)) return 0;
    uint64_t b1 = 1ull << (l1 & 63);
    if (fm2[l2] & b1) return FULL;
    if (!(nm2[l2] & b1)) return 0;
    uint64_t b = 1ull << (pos & 63);
    if (fm1[l1] & b) return FULL;
    if (!(nm1[l1] & b)) return 0;
    return limb[pos];
}

/* ---- parse ---- */

static inline int rd(){
    while (*P <= ' ') P++;
    int v = 0;
    while (*P > ' ') v = v*10 + (*P - '0'), P++;
    return v;
}

__attribute__((noreturn)) static void done(uint64_t olen, struct DuckInfo *di, char *out, int use_di){
    if (use_di) di->stdout_size = olen;
    else (void)!write(1, out, olen);
    asm volatile("mov $60, %%eax; xor %%edi, %%edi; syscall" ::: "rax","rdi","memory");
    __builtin_unreachable();
}

int main(){
    struct DuckInfo *di = (struct DuckInfo*)getauxval(0x6b637564);
    static char lbuf[1<<24];
    static char obuf[1<<22];
    char *out;
    int use_di = 0;
    if (di && di->abi_version >= 1 && di->stdin_ptr && di->stdout_ptr){
        P = di->stdin_ptr;
        out = di->stdout_ptr;
        use_di = 1;
    } else {
        long n2 = 0, t;
        while (n2 < (long)sizeof(lbuf) && (t = read(0, lbuf + n2, sizeof(lbuf) - n2)) > 0) n2 += t;
        lbuf[n2] = 0;
        P = lbuf;
        out = obuf;
    }

    int n = rd();
    rd(); rd(); rd();
    P++;

    for (int i = 0; i < n; i++){
        int op = *P++ - '0';
        P++;
        if (op == 1){
            int neg = 0;
            if (*P == '-'){ neg = 1; P++; }
            int a = 0;
            while (*P > ' '){ a = a*10 + (*P - '0'); P++; }
            P++;
            int b = 0;
            while (*P > ' '){ b = b*10 + (*P - '0'); P++; }
            P++;
            int q = b / 30;
            int r = b % 30;
            if (neg) sub_val(q, r, (uint32_t)a);
            else if (a) add_val(q, r, (uint32_t)a);
        } else {
            int k = 0;
            while (*P > ' '){ k = k*10 + (*P - '0'); P++; }
            P++;
            int q = k / 30;
            int r = k % 30;
            uint32_t v = read_limb(q);
            out[O++] = (char)('0' + ((v >> r) & 1));
            out[O++] = '\n';
        }
    }
    done(O, di, out, use_di);
}

评测结果