0x35 Xor (线性基,高斯消元)

首先求出线性基 然后根据类似二进制的原理 观察 1 1 1 的个数 然后可以算出我凑出的异或值是第几大的

注意 当线性基个数小于 n n n 时会存在 0 0 0 所以要将询问第 k k k 小改为询问第 K − 1 K-1 K1

然后板子基本就是高斯消元?

反正就是用 O ( l o g w ) O(log w) O(logw) 级别的数,代替了原来的 n n n 个数。在什么意义上这两个东西相等呢?就是这 l o g w log w logw 个数,通过线性组合可以得到所有的 n n n 个数。一般用于求解异或和最值的问题。

此题是求第 k k k 小的异或和。

根据线性基的性质,我们可以得到第 k k k 小的就是:把 m m m 个线性基从小到大排序,分别叫 v 0 , v 1 , … v m − 1 v_0,v_1,…v_{m−1} v0,v1,vm1 k k k 的二进制表示为( b x … . b 1 b 0 b_x….b_1b_0 bx.b1b0),则答案就是 ∑ x i = 0 v i ∗ b i \sum\limits_{x_i=0} {v_i∗b_i} xi=0vibi

注意特判0

#include <map>
#include <set>
#include <ctime>
#include <queue>
#include <stack>
#include <cmath>
#include <vector>
#include <bitset>
#include <cstdio>
#include <cctype>
#include <string>
#include <numeric>
#include <cstring>
#include <cassert>
#include <climits>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std ;
#define rep(i, a, b) for (int i = (a); i <= (b); i++)
#define per(i, a, b) for (int i = (a); i >= (b); i--)
#define loop(s, v, it) for (s::iterator it = v.begin(); it != v.end(); it++)
#define cont(i, x) for (int i = head[x]; i; i = e[i].nxt)
#define clr(a) memset(a, 0, sizeof(a))
#define ass(a, sum) memset(a, sum, sizeof(a))
#define lowbit(x) (x & -x)
#define all(x) x.begin(), x.end()
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define iv inline void
#define enter cout << endl
#define siz(x) ((int)x.size())
#define file(s) freopen(s".in", "r", stdin), freopen(s."out", "w", stdout)
typedef long long ll ;
typedef unsigned long long ull ;
typedef pair <int, int> pii ;
typedef vector <int> vi ;
typedef vector <pii> vii ;
typedef queue <int> qi ;
typedef set <int> si ;
typedef map <int, int> mii ;
typedef map <string, int> msi ;
const int N = 100010 ;
const int INF = 0x3f3f3f3f ;
const int iinf = 1 << 30 ;
const ll linf = 2e18 ;
const int MOD = 1000000007 ;
const double eps = 1e-7 ;
void print(int x) { cout << x << endl ; exit(0) ; }
void PRINT(string x) { cout << x << endl ; exit(0) ; }
void douout(double x){ printf("%lf\n", x + 0.0000000001) ; }

ull a[N] ;
int n, m, t, T ;

signed main(){
    scanf("%d", &T) ;
    rep(Case, 1, T) {
    	scanf("%d", &n) ;
		rep(i, 1, n) scanf("%llu", &a[i]) ;
		bool zero = 0 ;
		t = n ;
		rep(i, 1, n) {
			rep(j, i + 1, n) if (a[j] > a[i]) swap(a[j], a[i]) ;
			if (!a[i]) {
				zero = 1 ;
				t = i - 1 ;
				break ;
			}
			per(k, 63, 0)
			if (a[i] >> k & 1) {
				rep(j, 1, n) if (i != j && (a[j] >> k & 1)) a[j] ^= a[i] ;
				break ;
			}
		}
		scanf("%d", &m) ;
		printf("Case #%d:\n", Case) ;
		while (m--) {
			ull k, ans = 0 ; scanf("%llu", &k) ;
			if (zero) k-- ;
			if (k >= 1llu << t) puts("-1") ;
			else {
				per(i, t - 1, 0) if (k >> i & 1) ans ^= a[t - i] ;
				printf("%llu\n", ans) ;
			}
		}
	}


	return 0 ;
}

/*
写代码时请注意:
	1.ll?数组大小,边界?数据范围?
	2.精度?
	3.特判?
	4.至少做一些
思考提醒:
	1.最大值最小->二分?
	2.可以贪心么?不行dp可以么
	3.可以优化么
	4.维护区间用什么数据结构?
	5.统计方案是用dp?模了么?
	6.逆向思维?
*/


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值