线性基计数(K大异或和)

线性基计数(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 2cnt1 的。

//#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;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值