题意
传送门 HDU 3949
题解
处理集合异或的值域,考虑线性基。在不允许出现 a i ⊕ a i a_i\oplus a_i ai⊕ai 运算的限制下,线性基的非零子集的异或即集合的异或非零值域,若集合中出现可以由线性基表达的值,或说线性基的规模小于集合规模,那么集合异或的值域包含零,需要特殊处理。
考虑值域第 k k k 大,从高位向低位用类似二分的方法处理。假设已处理好大于 i i i 位的线性基,若第 i i i 位存在线性基, [ 0 , i − 1 ] [0,i-1] [0,i−1] 位存在 n n n 个线性基,那么可以构成 2 n 2^n 2n 个包含空集且小于包含第 i i i 位线性基的值,根据 k k k 值判断是否使用第 i i i 位线性基。
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
#define maxl 60
typedef long long ll;
int T, N, Q, cnt[maxl + 1];
ll B[maxl + 1];
void insert(ll n)
{
for (int i = maxl; i >= 0; --i)
{
if (n >> i & 1)
{
if (B[i])
n ^= B[i];
else
{
for (int j = 0; j < i; ++j)
if (n >> j & 1)
n ^= B[j];
for (int j = i + 1; j <= maxl; ++j)
if (B[j] >> i & 1)
B[j] ^= n;
B[i] = n;
break;
}
}
}
}
int main()
{
scanf("%d", &T);
for (int c = 1; c <= T; ++c)
{
printf("Case #%d:\n", c);
memset(B, 0, sizeof(B));
scanf("%d", &N);
for (int i = 0; i < N; ++i)
{
ll a;
scanf("%lld", &a);
insert(a);
}
for (int i = 1; i <= maxl; ++i)
cnt[i] = cnt[i - 1] + (B[i - 1] ? 1 : 0);
int noZero = cnt[maxl] < N ? 0 : 1;
scanf("%d", &Q);
while (Q--)
{
ll k, res = 0;
scanf("%lld", &k);
if (k > (1LL << cnt[maxl]) - noZero)
{
puts("-1");
continue;
}
for (int i = maxl; i >= 0; --i)
{
if (!B[i])
continue;
ll n = 1LL << cnt[i];
if (res)
{
if (k > n)
res ^= B[i], k -= n;
}
else if (k > n - noZero)
res ^= B[i], k -= n - noZero;
}
printf("%lld\n", res);
}
}
return 0;
}
更简单地,将线性基的元素看做二进制的各位,第 i i i 位的 1 , 0 1,0 1,0 分别代表是否选取按位数排序的第 i i i 个线性基进行异或。那么当 k k k 二进制位取 1 1 1 时答案异或这一位对应的线性基即可。由 n n n 个元素组成的非空子集数为 2 n − 1 2^n-1 2n−1,特殊处理值域出现零的情况。
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
#define maxl 60
typedef long long ll;
int T, N, Q, cnt[maxl + 1];
ll B[maxl + 1];
void insert(ll n)
{
for (int i = maxl; i >= 0; --i)
{
if (n >> i & 1)
{
if (B[i])
n ^= B[i];
else
{
for (int j = 0; j < i; ++j)
if (n >> j & 1)
n ^= B[j];
for (int j = i + 1; j <= maxl; ++j)
if (B[j] >> i & 1)
B[j] ^= n;
B[i] = n;
break;
}
}
}
}
int main()
{
scanf("%d", &T);
for (int c = 1; c <= T; ++c)
{
printf("Case #%d:\n", c);
memset(B, 0, sizeof(B));
scanf("%d", &N);
for (int i = 0; i < N; ++i)
{
ll a;
scanf("%lld", &a);
insert(a);
}
for (int i = 1; i <= maxl; ++i)
cnt[i] = cnt[i - 1] + (B[i - 1] ? 1 : 0);
int hasZero = cnt[maxl] < N ? 1 : 0;
scanf("%d", &Q);
while (Q--)
{
ll k, res = 0;
scanf("%lld", &k);
k -= hasZero;
if (k > (1LL << cnt[maxl]) - 1)
{
puts("-1");
continue;
}
for (int i = maxl; i >= 0; --i)
{
if (!B[i])
continue;
if (k >> cnt[i] & 1)
res ^= B[i];
}
printf("%lld\n", res);
}
}
return 0;
}