给孩子们一个手写的bitset快多了

struct bitst {
    static constexpr int SZ = SIZ;
    alignas(32) u64 a[SZ];

    void reset() {
        memset(a, 0, sizeof(a));
    }
    void set(int x) {
        a[x >> 6] |= 1ull << (x & 63);
    }
    void flip(int x) {
        a[x >> 6] ^= 1ull << (x & 63);
    }
    int val(int x) {
        return (a[x >> 6] >> (x & 63)) & 1;
    }

    void operator|=(const bitst &b) {
        for (int i=0;i<SZ;i+=4) {
            __m256i x = _mm256_load_si256((__m256i*)(a+i));
            __m256i y = _mm256_load_si256((__m256i*)(b.a+i));
            __m256i z = _mm256_or_si256(x,y);
            _mm256_store_si256((__m256i*)(a+i),z);
        }
    }

    void operator^=(const bitst &b) {
        for (int i=0;i<SZ;i+=4) {
            __m256i x = _mm256_load_si256((__m256i*)(a+i));
            __m256i y = _mm256_load_si256((__m256i*)(b.a+i));
            __m256i z = _mm256_xor_si256(x,y);
            _mm256_store_si256((__m256i*)(a+i),z);
        }
    }

    void operator&=(const bitst &b) {
        for (int i=0;i<SZ;i+=4) {
            __m256i x = _mm256_load_si256((__m256i*)(a+i));
            __m256i y = _mm256_load_si256((__m256i*)(b.a+i));
            __m256i z = _mm256_and_si256(x,y);
            _mm256_store_si256((__m256i*)(a+i),z);
        }
    }
};