【HDU3949】XOR（线性基）

Description

给你$n$个不大于$10^{18}$的数，有$q$次询问，每次给定$k$，求这$n$个数中所有可以通过异或运算得到的数中第$k$小的数。

Solution

算是线性基模板题吧。
把$n$个数插入线性基，进行如下处理：
从高位向低位枚举，若一个数除了最高位外还有其他位上有$1$且该位对应的低位的数非$0$，那么将高位的数异或低位上的数。感觉说不清楚，具体看代码中的函数LinearBasis.init()吧。。
然后就可以将线性基中的非$0$的数存下来，根据$k$按位取。

注意处理可以异或成$0$的情况：如果插入一个数时发现该数和当前线性基线性相关，那么就可以通过异或得到 0 0 <script type="math/tex" id="MathJax-Element-579">0</script>。

Source

/************************************************
 * Au: Hany01
 * Date: Mar 11th, 2018
 * Prob: [HDU3949] XOR
 * Email: hany01@foxmail.com
************************************************/

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define fir first
#define sec second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000007)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia

template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

template <typename T> inline T read()
{
    register T _, __; register char c_;
    for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
    for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

inline void File()
{
#ifdef hany01
    freopen("hdu3949.in", "r", stdin);
    freopen("hdu3949.out", "w", stdout);
#endif
}

const int sz = 61;

struct LinearBasis
{
    LL a[sz + 2], b[sz + 2];
    int cnt, mark;

    inline void Clear() { Set(a, 0), mark = 0, cnt = 0; }

    inline void Insert(LL x)
    {
        Fordown(i, sz, 0) {
            if (x >> i & 1) {
                if (!a[i]) { a[i] = x; break; }
                x ^= a[i];
            }
            if (!x) { mark = 1; break; }
        }
    }

    inline void Init()
    {
        Fordown(i, 60, 0) if (a[i]) {
            Fordown(j, i - 1, 0) if (a[i] >> j & 1) a[i] ^= a[j];
            b[cnt ++] = a[i];
        }
        reverse(b, b + cnt);
    }

    inline LL query(LL k)
    {
        if (mark) -- k;
        if (k > (1ll << cnt) - 1) return -1;
        LL Ans = 0;
        rep(i, cnt) if (k >> (LL)i & 1ll) Ans ^= b[i];
        return Ans;
    }

}Base;

int main()
{
    File();

    for (register int T = read<int>(), Case = 0; T --; ){
        printf("Case #%d:\n", ++ Case);
        Base.Clear();
        for (register int n = read<int>(); n --; ) Base.Insert(read<LL>());
        Base.Init();
        for (register int n = read<int>(); n --; ) cout << Base.query(read<LL>()) << endl;
    }

    return 0;
}
//采得百花成蜜后，为谁辛苦为谁甜。
//    -- 罗隐《蜂》