hdu3949 XOR(線性基)

HDU3949 XOR(線性基)

題目：

給你\(n\)個數，從其中隨便取任意數問你第\(k\)小的異或和是多少。

題解：

這題是線性基的應用之一。咱們知道一個集合的線性基能夠異或出這個集合的全部異或和，而且方法惟一。對於一個數x可否被異或出來，咱們能夠這樣作，假設x的最高位爲r，那麼在線性基裏面找到最高爲也爲r的數，讓x異或r。而後不斷重複這個操做，若是最後x能爲0那麼確定是能的。而後咱們如今想一想怎麼經過線性基求第k小的異或和。
首先，線性基能夠當作一個最大線性無關組，也就是說最高位以上都是0。如今假設最大線性基是這樣的：

00100010 - - - - - - 4
00011000 - - - - - - 3
00000110 - - - - - - 2
00000001 - - - - - - 1

那麼異或和排序顯然是（用編號表示）：
\(1\Rightarrow2\Rightarrow(2\bigoplus1)\Rightarrow3\Rightarrow(3\bigoplus1)\Rightarrow(3\bigoplus1\bigoplus2)\Rightarrow......\Rightarrow(4\bigoplus1\bigoplus2\bigoplus3)\)
這個排序顯然是按照二進制最高位排序加上去的。可是若是想按這樣的規律數顯然是不行的。不過，咱們能夠聯繫下上面說的一個數x可否組成的原理。若是我將線性基用離散化的思想處理，而後再按照上面的排序思想。就能夠把題目轉化乘k可否由離散後的線性基組成。
舉個例子：咱們假設要求第k小的數，先上面的線性基離散化就能夠當作

1000
0100
0010
0001

而後，這個新的線性基能夠組成\(2^5-1\)排名內的全部數。也就是原線性基全部能夠組成的異或的數量。k必定在這裏面，不然就不存在。所以，咱們只須要分解k的全部1就好了。

代碼：

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
long long a[maxn],p[maxn],tot;
void Guass(int n)
{
    memset(p,0,sizeof(p));
    for(int i=1;i<=n;i++)
    {
        for(int j=63;j>=0;j--)
        {
            if((a[i]>>j)&1)
            {
                if(p[j]) a[i]^=p[j];
                else
                {
                    p[j]=a[i];
                    break;
                }
            }
        }
    }
    for(int i=63;i>=0;i--)
    {
        if(!p[i])continue;
        for(int j=i+1;j<=62;j++)
        {
            if((p[j]>>i)&1) p[j]^=p[i];
        }
    }
    tot=0;
    for(int i=0;i<=63;i++) if(p[i]) p[tot++]=p[i];
}
int main()
{
    int T,i,j,n,Q;
    long long k;
    scanf("%d",&T);
    for(int s=1;s<=T;s++)
    {
        printf("Case #%d:\n",s);
        scanf("%d",&n);
        for(i=1;i<=n;i++) scanf("%I64d",&a[i]);
        Guass(n);
        scanf("%d",&Q);
        while(Q--)
        {
            scanf("%I64d",&k);
            if(n!=tot)k--;
            if(k>=(1ll<<tot))printf("-1\n");
            else
            {
                long long ans=0;
                for(i=0;i<=63;i++) if((k>>i)&1) ans^=p[i];
                printf("%I64d\n",ans);
             } 
        }
    }
}