线性基计数(k大异或和)
链接:#114. k 大异或和 - 题目 - LibreOJ (loj.ac)
题意:给定可重集 S,多次询问求第 k 小的真子集异或和。
题解:重新构造线性基,使 p [ i ] p[i] p[i] 的第 j 位在 p [ j ] p[j] p[j] 存在的情况下不为 1,这样异或 p [ j ] p[j] p[j] 一定使数值变大,抽象化为二进制便是取 p [ j ] p[j] p[j] 相当于 j 位取 1,则可以将 k 进行按位查看是否异或得结果。
tips:当线性基个数 c n t cnt cnt 小于 n 时,说明非满秩,所有结果中可以异或出 0,这样看来第 1 小得为 0 而非 p [ 0 ] p[0] p[0],所以对k 减一。总的异或结果是不超过 2 c n t − 1 2^{cnt}-1 2cnt−1 的。
//#pragma GCC optimize("O3")
//#pragma GCC optimize("unroll-loops")
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=55;
ll p[maxn],f[maxn],cnt,n,m,ans,u,v,w;
inline void insert(ll x)
{
for(int i=50;i>=0;i--)
{
if(x>>i&1)
{
if(!p[i]){p[i]=x; break;}
x^=p[i];
}
}
return;
}
inline void rebuild()
{
for(int i=0;i<=50;i++)
{
for(int j=i-1;j>=0;j--)
{
if(p[i]>>j&1)p[i]^=p[j];
}
if(p[i])f[cnt++]=p[i];
}
return;
}
inline ll query(ll k)
{
if(cnt!=n)k--;
if(k>=(1ll<<cnt))return -1;
ll ans=0;
for(int i=0;i<cnt;i++)ans^=f[i]*(k>>i&1);
return ans;
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0); cout.tie(0);
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>u,insert(u);
}
rebuild(),cin>>m;
for(int i=1;i<=m;i++)
{
ll k; cin>>k,cout<<query(k)<<"\n";
}
return 0;
}