异或线性基
我们考虑这样一个问题:给定
我们考虑用一个集合
当我们要加入
所以只要
但是我们能观察到这样最多只会翻倍
但是这样就不能快速判断
很容易想到
我们对于
而判断
由于线性基中的元素不能相互异或出来,所以能异或出
空间复杂度
代码
#include<bits/stdc++.h>
using namespace std;
using ll = long long;
const int MAXN = 65;
int n;
ll d[MAXN], ans;
void add(ll x) {
for(ll i = 61; i >= 0; --i) {
if(x & (1ll << i)) {
if(d[i]) {
x ^= d[i];
}else {
d[i] = x;
return;
}
}
}
}
int main() {
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cin >> n;
for(int i = 1; i <= n; ++i) {
ll x;
cin >> x;
add(x);
}
for(int i = 61; i >= 0; --i) {
ans = max(ans, (ans ^ d[i]));
}
cout << ans;
return 0;
}
Luogu P3812
题目描述
给定
思路
首先对这些数建线性基,接着从高到低枚举每一位,如果
空间复杂度
代码
#include<bits/stdc++.h>
using namespace std;
using ll = long long;
const int MAXN = 65;
int n;
ll d[MAXN], ans;
void add(ll x) {
for(ll i = 61; i >= 0; --i) {
if(x & (1ll << i)) {
if(d[i]) {
x ^= d[i];
}else {
d[i] = x;
return;
}
}
}
}
int main() {
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cin >> n;
for(int i = 1; i <= n; ++i) {
ll x;
cin >> x;
add(x);
}
for(int i = 61; i >= 0; --i) {
ans = max(ans, (ans ^ d[i]));
}
cout << ans;
return 0;
}
AT ABC141 F
题目描述
给定
思路
我们令
空间复杂度
代码
#include<bits/stdc++.h>
using namespace std;
using ll = long long;
int n;
ll res, ans, d[60];
void Insert(ll x) {
for(int i = 59; i >= 0; --i) {
if((x >> i) & 1) {
if(!d[i]) {
d[i] = x;
return;
}
x ^= d[i];
}
}
}
int main() {
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cin >> n;
for(int i = 1; i <= n; ++i) {
ll x;
cin >> x;
res ^= x;
Insert(x);
}
for(int i = 59; i >= 0; --i) {
if(!((res >> i) & 1)) {
ans = max(ans, ans ^ d[i]);
}else {
Insert(d[i] ^ (1ll << i));
}
}
cout << ans + (res ^ ans);
return 0;
}
CF 895 C
题目描述
给定
思路
我们可以用一个数每个质因子的次幂的奇偶性来表示一个数,而一个完全平方数就是
在异或线性基中的元素是不可能异或出
空间复杂度
代码
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 100001, MAXP = 19, MOD = int(1e9) + 7, p[MAXP] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67};
int n, a[MAXN], num[71], d[MAXN], ans;
bool op = 1;
int C(int x) {
int res = 0;
for(int i = 2; i * i <= x; ++i) {
while(x % i == 0) {
res ^= (1 << num[i]);
x /= i;
}
}
if(x > 1) {
res ^= (1 << num[x]);
}
return res;
}
void Insert(int x) {
for(int i = MAXP - 1; i >= 0; --i) {
if((x >> i) & 1) {
if(!d[i]) {
d[i] = x, ans--;
break;
}else {
x ^= d[i];
}
}
}
}
int Pow(int a, int b) {
int res = 1;
for(; b; a = 1ll * a * a % MOD, b >>= 1) {
if(b & 1) {
res = 1ll * res * a % MOD;
}
}
return res;
}
int main() {
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
for(int i = 0; i < MAXP; ++i) {
num[p[i]] = i;
}
cin >> n;
ans = n;
for(int i = 1, x; i <= n; ++i) {
cin >> x;
Insert(C(x));
}
cout << ((Pow(2, ans) - 1) % MOD + MOD) % MOD;
return 0;
}
CF 1902 F
题目描述
给定一颗树,每个结点都有一个权值
思路
我们对每个点都记录一个线性基,这个线性基是从当前结点往根节点插入的。接着考虑怎么把线性基从父亲转移到儿子。
很容易想到父亲没有插入线性基的元素,儿子肯定也不会插入。所以记录插入线性基的元素,依次枚举并插入即可。
查询时枚举插入线性基的元素,如果在路径上就插入。再进行判断。
由于至多只有
代码
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 200001;
struct xor_linear_basis {
int d[20];
xor_linear_basis() {
fill(d, d + 20, 0);
}
bool Insert(int x) {
for(int i = 20; i >= 0; --i) {
if((x >> i) & 1) {
if(!d[i]) {
d[i] = x;
return 1;
}
x ^= d[i];
}
}
return 0;
}
bool In(int x) {
for(int i = 20; i >= 0; --i) {
if((x >> i) & 1) {
x ^= d[i];
}
}
return !x;
}
}d[MAXN];
int n, a[MAXN], q, f[19][MAXN], dep[MAXN];
vector<int> e[MAXN], ve[MAXN];
void dfs(int u, int fa) {
if(d[u].Insert(a[u])) {
ve[u].emplace_back(u);
}
for(int v : ve[fa]) {
if(d[u].Insert(a[v])) {
ve[u].emplace_back(v);
}
}
dep[u] = dep[fa] + 1, f[0][u] = fa;
for(int i = 1; i <= 18; ++i) {
f[i][u] = f[i - 1][f[i - 1][u]];
}
for(int v : e[u]) {
if(v != fa) {
dfs(v, u);
}
}
}
int LCA(int u, int v) {
if(dep[u] < dep[v]) {
swap(u, v);
}
int d = dep[u] - dep[v];
for(int i = 18; i >= 0; --i) {
if((d >> i) & 1) {
u = f[i][u];
}
}
if(u == v) {
return u;
}
for(int i = 18; i >= 0; --i) {
if(f[i][u] != f[i][v]) {
u = f[i][u], v = f[i][v];
}
}
return f[0][u];
}
int main() {
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
cin >> n;
for(int i = 1; i <= n; ++i) {
cin >> a[i];
}
for(int i = 1, u, v; i < n; ++i) {
cin >> u >> v;
e[u].emplace_back(v);
e[v].emplace_back(u);
}
dfs(1, 0);
cin >> q;
for(int i = 1, u, v, k; i <= q; ++i) {
cin >> u >> v >> k;
int l = LCA(u, v);
xor_linear_basis res;
for(int x : ve[u]) {
if(dep[x] >= dep[l]) {
res.Insert(a[x]);
}
}
for(int x : ve[v]) {
if(dep[x] >= dep[l]) {
res.Insert(a[x]);
}
}
cout << (res.In(k) ? "YES\n" : "NO\n");
}
return 0;
}
本文作者:Yaosicheng124
本文链接:https://www.cnblogs.com/yaosicheng124/p/18426024
版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步