【题目分析】
高斯消元求线性基。
题目本身不难,但是两种维护线性基的方法引起了我的思考。
void gauss(){
k=n;
F(i,1,n){
F(j,i+1,n) if (a[j]>a[i]) swap(a[i],a[j]);
if (!a[i]) {k=i-1; break;}
D(j,30,0) if (a[i]>>j & 1){
b[i]=j;
F(x,1,n) if (x!=i && a[x]>>j&1) a[x]^=a[i];
break;
}
}
}
——高斯消元求线性基
for (int i=1;i<=n;++i)
for (int j=31;j>=0;--j)
if ((a[i]>>j)&1){
if (!lb[j]) {lb[j]=a[i]; cnt++; break;}
else a[i]^=lb[j];
}
——动态维护线性基
不会高斯消元解Xor方程组的我,直接使用了第二种方式求解,发现直接WA飞了。
(后来一想,居然过了样例)。
那么他们有什么差别呢。
我对拍了许多组,发现他们求出的线性基的大小是相同的。
但是高斯消元的线性基有一个神奇的特征,是使得该位为1的最小的数。(最小的)
那么有必要去写高斯消元吗?
显然不必要,做一个小操作就好了。
于是改了改动态维护线性基的代码,成了这个样子 ↓
for (int i=1;i<=n;++i)
for (int j=31;j>=0;--j)
if ((a[i]>>j)&1){
if (!lb[j]) {lb[j]=a[i]; cnt++; break;}
else a[i]^=lb[j];
}
for (int i=31;i>=0;--i)
if (lb[i])
for (int j=i-1;j>=0;--j)
if ((lb[i]>>j)&1) lb[i]^=lb[j];
——改版
神奇的AC了。线性基与高斯消元的结果相同。
考虑时间复杂度,都是log*n的,自然没什么差别,但是用哪种就是仁者见仁智者见智了。
实际上高斯消元会快一些(达不到复杂度上限),而动态维护线性基是标准的上限(雾)
【代码】
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <set>
#include <map>
#include <string>
#include <algorithm>
#include <vector>
#include <iostream>
#include <queue>
using namespace std;
#define maxn 100005
#define ll long long
int read()
{
int x=0,f=1; char ch=getchar();
while (ch<'0'||ch>'9') {if (ch=='-') f=-1; ch=getchar();}
while (ch>='0'&&ch<='9') {x=x*10+ch-'0'; ch=getchar();}
return x*f;
}
const int mod=10086;
int n,cnt=0,q;
int lb[34],a[maxn];
int power(int a,int b)
{
// printf ("Pow %d ^ %d is ",a,b);
int ret=1;
while (b)
{
if (b&1) (ret*=a)%=mod;
(a*=a)%=mod;
b>>=1;
}
// printf("%d\n",ret);
return ret;
}
bool cmp(int a,int b)
{
return a>b;
}
int main()
{
n=read();
for (int i=1;i<=n;++i) a[i]=read();
// sort(a+1,a+n+1);
// for (int i=1;i<=n;++i) cout<<a[i]<<" "; cout<<endl;
q=read();
// cout<<"query : "<<q<<endl;
for (int i=1;i<=n;++i)
for (int j=31;j>=0;--j)
if ((a[i]>>j)&1){
if (!lb[j]) {lb[j]=a[i]; cnt++; break;}
else a[i]^=lb[j];
}
for (int i=31;i>=0;--i)
if (lb[i])
for (int j=i-1;j>=0;--j)
if ((lb[i]>>j)&1) lb[i]^=lb[j];
// printf("The Xor Base is %d\n",cnt);
int rk=0,x=0,tmp=0;
for (int j=31;j>=0;--j)
if (lb[j]){
tmp++;
// cout<<tmp<<":"<<lb[j]<<endl;
if ((x^lb[j])>q) continue;
x^=lb[j];
// printf("now add %d\n",cnt-tmp);
rk=(rk+power(2,cnt-tmp))%mod;
}
for (int i=1;i<=n-cnt;++i)
rk=(rk*2)%mod;
rk++;
cout<<rk<<endl;
}