欢迎大家访问我的老师的OJ———caioj.cn
题面描述
思路
用高斯消元求出 a 1 , a 2 , . . . . . . a n a_1,a_2,......a_n a1,a2,......an构成的异或空间的基。不妨设这个基由整数 b 1 , b 2 , . . . b t b_1,b_2,...b_t b1,b2,...bt构成,其中 b 1 > b 2 > . . . > b t b_1>b_2>...>b_t b1>b2>...>bt。从基底中选取若干个整数,显然有 2 t 2^t 2t种取法,因此t维异或空间中共有 2 t 2^t 2t个整数,与这 2 t 2^t 2t中取法一一对应。
因此,高斯消元时,应该有t个主元,分别是 b 1 , b 2 , . . . . . . b t b_1,b_2,......b_t b1,b2,......bt的最高位
显然最高位 c 1 > c 2 > , . . . . . . > c t c_1>c_2>,......>c_t c1>c2>,......>ct.
在简化阶梯形矩阵中,对于每个主元,有且仅有一个位置 ( i , j ) (i,j) (i,j),满足该主元的系数非零。第 j j j列的其他位置都是 0 0 0,第i行的第 1 1 1~ j − 1 j-1 j−1列都是 0 0 0.
我们可以把询问中给出的 k − 1 k-1 k−1进行二进制分解(减掉 1 1 1是为了把二进制下最小的数定义为第 0 0 0小,题目中是第 1 1 1小),若 k − 1 k-1 k−1的第 j ( 0 ≤ j < t ) j(0\le j<t) j(0≤j<t)位等于 1 1 1,就选取 b t − j b_{t-j} bt−j,最后把所有选出的 b b b异或起来,就得到了答案。特别地,若 k ≥ 2 t k\ge2^t k≥2t,则无解。
另外,异或空间中可能存在 a i xor a i a_i \operatorname{xor} a_i aixorai的情况(而题目不允许出现这种情况)
所以异或空间的第k小(一定含0)的整数与题目中所求的第k小的整数稍有不同。
因此我们要加以判断(在简化阶梯型矩阵中,若出现”零行“,出现”零行”,意味着多个数可以异或为 0 0 0,但有时候我们并不能从多个数中异或出 0 0 0,这时候我们就要稍作调整让 k + + k++ k++,使得k变回题目原意的 k k k,因为我们并不能以 0 0 0为第一 小的数,符合题意).
如果出现多个数异或为 0 0 0,则要将 k − 1 k-1 k−1。
否则,则不用。
AC code
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define ull unsigned long long
using namespace std;
const int N=10010;
ull a[N];
int n,m,t,T;
int main()
{
scanf("%d",&T);
for(int C=1;C<=T;C++)
{
scanf("%d",&n);
for(int i=1;i<=n;i++)scanf("%llu",&a[i]);
bool zero=0;t=n;
for(int i=1;i<=n;i++)
{
int x=i;
for(int j=i+1;j<=n;j++)if(a[j]>a[x])x=j;
swap(a[x],a[i]);
if(a[i]==0)//剩下都是自由元
{
t=i-1;zero=1;break;
}
for(int j=63;j>=0;j--)
if(a[i]>>j&1)
{
for(int k=n;k>=1;k--)
if(a[k]>>j&1&&k!=i)
a[k]^=a[i];
break;
}
}
printf("Case #%d:\n",C);
scanf("%d",&m);
while(m--)
{
ull k,ans=0;
scanf("%llu",&k);
if(zero)k--;
if(k>=1llu<<t)puts("-1");
else
{
for(int i=t-1;i>=0;i--)
if(k>>i&1)ans^=a[t-i];
printf("%llu\n",ans);
}
}
}
return 0;
}