[线性代数] hdu 3949 XOR

题意：

给你N个数，要你从中取出若干个进行异或运算，求最后所有可以得到的异或结果中的

第k小的异或值。N<=10^18

这个题是用来研究xor性质的一道好题。

首先我们可以暴力地找出些规律，我们发现不管拿出多少个数，他们能xor到的数，出现的次数都是一样的，并且都是2的倍数。事实上，我们不论用原数ai还是ai xor aj得到的数，去xor起来得到的数都是这些。

这样我们可以想方设法地化简我们用到的基数，也就是找到一组基数使得它们能够xor出所有原数能xor出的数，并且要尽量简单。我们可以模仿高斯消元的过程，就是把每个数当成一行二进制数去消元，然后保留一个类似上三角的东西，从高位向低位每一行保留连续的几个1，这些1行与行之间是不重复的，但每一行之中还会有额外的一些1，这些1是可重复的，就造成了一些数是xor不出来的。

假设我们选出了m行，也就是m个基数，那么我们可以xor出2^m个数，如果总共给你了n个数，那么就有n-m个数是可以被那些基数完全替代的，这些行在高斯消元的过程中就是自由变量。这些行的值为0，也就是它去xor任意一个数，原数都不会改变，这样就产生了2^n-m个重复的数，这样便验证了我们一开始的结论。

其实这组基数还有一个很好的性质，就是xor出的第k大(去重)的数，恰好是k二进制分解后对应位上的基数xor起来的到的数。因为从上到下每个基数代表了从高到低的一些二进制位，那么我们就相当于用这些连续的二进制位去组合形成所有的数。

那对于这个题来说，我们要求第k大的数，我们把m行基数求出来了以后，选出k的各个二进制位对应的基数，xor起来就是答案。

注意long long，注意0的情况。

代码其中gauss部分，可以作为解模2的方程组的优化。

#include"cstdlib"
#include"cstdio"
#include"cstring"
#include"cmath"
#include"queue"
#include"algorithm"
#include"iostream"
using namespace std;
unsigned __int64 a[123450];
unsigned __int64 bit[66];
int n;
void debug()
{
    for(int i=0; i<n; i++) printf("%I64u ",a[i]);
    puts("");
}
int gauss()         //可以用来解模二的方程，加快速度
{
    int row=0;
    int i,j;
    for(i=63; i>=0; i--)
    {
        for(j=row; j<n; j++) if(a[j]&bit[i]) break;
        if(j!=n)
        {
            swap(a[row],a[j]);
            for(j=0; j<n; j++)
            {
                if(j==row) continue;
                if(a[j]&bit[i])
                    a[j]^=a[row];
            }
            row++;
        }
    }
    return row;
}
int main()
{
    int t,cas=1;
    unsigned __int64 tt=1;
    cin>>t;
    for(int i=0; i<64; i++)
    {
        bit[i]=tt;
        tt*=2;
    }

    while(t--)
    {
        int i;
        cin>>n;
        int nofree;
        for(i=0; i<n; i++) scanf("%I64u",&a[i]);
        //debug();
        nofree=gauss();
        //debug();
        int q;
        cin>>q;
        printf("Case #%d:\n",cas++);
        while(q--)
        {
            unsigned __int64 k;
            scanf("%I64u",&k);
            if(nofree<n)
            {
                if(k==1)
                {
                    puts("0");
                    continue;
                }
                else k--;
            }
            if(k>=((unsigned __int64)1<<nofree)) printf("-1\n");
            else
            {
                unsigned __int64 ans=0;
                for(i=0; i<64; i++)
                {
                    if(k&bit[i]) ans^=a[nofree-i-1];
                }
                printf("%I64u\n",ans);
            }
        }


    }
}