【题目分析】

高斯消元求线性基。

题目本身不难,但是两种维护线性基的方法引起了我的思考。

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;
}

  

BZOJ 2844 albus就是要第一个出场 ——高斯消元 线性基的更多相关文章

  1. bzoj 2844: albus就是要第一个出场 高斯消元

    LINK 题意:看题目不如看样例解释.给出有n个数的集合,对这些子集中的数求异或,升序统计所有子集得到的数(重复会被计入),询问一个数x,问这个数出现的第一个位置 思路:在这里要求一个所有可能出现的异 ...

  2. BZOJ 2844: albus就是要第一个出场 [高斯消元XOR 线性基]

    2844: albus就是要第一个出场 题意:给定一个n个数的集合S和一个数x,求x在S的$2^n$个子集从小到大的异或和序列中最早出现的位置 一开始看错题了...人家要求的是x第一次出现位置不是第x ...

  3. 【BZOJ2844】albus就是要第一个出场 高斯消元求线性基

    [BZOJ2844]albus就是要第一个出场 Description 已知一个长度为n的正整数序列A(下标从1开始), 令 S = { x | 1 <= x <= n }, S 的幂集2 ...

  4. BZOJ 2844: albus就是要第一个出场

    2844: albus就是要第一个出场 Time Limit: 6 Sec  Memory Limit: 128 MBSubmit: 1134  Solved: 481[Submit][Status] ...

  5. BZOJ 2844 albus就是要第一个出场(高斯消元)

    题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=2844 题意: 给出一个长度为n的正整数数列A.每次选出A的一个子集进行抑或(空集抑或值为 ...

  6. BZOJ 4004 JLOI2015 装备购买 高斯消元+线性基

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=4004 Description 脸哥最近在玩一款神奇的游戏,这个游戏里有 n 件装备,每件装 ...

  7. bzoj 2844 albus就是要第一个出场 - 线性基

    题目传送门 这是个通往vjudge的虫洞 这是个通往bzoj的虫洞 题目大意 给定集合$S$,现在将任意$A\subseteq S$中的元素求异或和,然后存入一个数组中(下标从1开始),然后从小到大排 ...

  8. bzoj 2844 albus就是要第一个出场 异或和出现次数 线性基

    题目链接 题意 给定\(n\)个数,将其所有的子集(\(2^n\)个)的异或和按升序排列.给出一个询问\(q\),问\(q\)在该序列中第一次出现位置的下标(下标从\(1\)开始). 题解 结论 记其 ...

  9. BZOJ 2844 高斯消元 线性基

    思路: //By SiriusRen #include <cstdio> #include <cstring> #include <algorithm> using ...

随机推荐

  1. ANTLR3完全参考指南读书笔记[03]

    前言 文中第4章内容有点多,有点枯燥,但不坚持一下,之前所做的工作就白做了. 再次确认一下总体目标: protege4编辑器中Class Definition中语法解析和错误提示: Java虚拟机规范 ...

  2. Oracle 基础 游标

    一:游标的基本原理 游标用来处理从数据库中检索的多行记录(使用SELECT语句).利用游标,程序可以逐个地处理和遍历一次检索返回的整个记录集. 为了处理SQL语句,Oracle将在内存中分配一个区域, ...

  3. SQLBindCol---数组输出

    SQLSetStmtAttr(hstmt,SQL_ATTR_ROW_BIND_TYPE,SQL_BIND_BY_COLUMN,0); SQLSetStmtAttr(hstmt,SQL_ATTR_ROW ...

  4. html - table 表格不被撑开,td某些列宽度固定某些列自适应

    table-layout 属性的解释:http://www.w3school.com.cn/cssref/pr_tab_table-layout.asp 1,实现table细边框,设置如下css: t ...

  5. Chrome 插件——二维码生成笔记

    原来用了几个生成二维码的插件,总是遇见各种问题……最后索性自己弄一个,这里顺便记录一下. Chrome 插件很开放!!!你只要拿到了crx文件,然后把文件的后缀名改为zip,就可以解压了,最后一切的一 ...

  6. MapReduce论文学习

    MapReduce和区块链有什么相同的地方? 我的天哪,他俩还有相同的地方呢.我书读的少,你别骗我. 他俩还真有相同点,绝不忽悠. 他俩都有一个高大上的名字. 区块链就是一个分布式数据库,并不是什么神 ...

  7. Python-random模块-59

    random模块: 随机数模块 >>> import random #随机小数 >>> random.random() # 大于0且小于1之间的小数 0.76643 ...

  8. jieba安装

    执行“pip install jieba”后忽略此条提示" You are using pip version 9.0.3, however version 10.0.1 is availa ...

  9. jquery获取父元素或父节点的方法

    jquery获取父元素方法比较多,比如parent(),parents(),closest()这些都能帮你实现查找父元素或节点,下面我们来一一讲解: 先举个例子: <ul class=" ...

  10. CSS的下拉菜单被挡住,修改Z-INDEX也不成功

    CSS的下拉菜单被挡住,修改Z-INDEX也不成功 做了一个鼠标放上去就出现的下拉菜单,但是这个下拉的内容被挡住了. Z-INDEX 是设置不同块的层次的,我修改了问题还是有. 后来发现是必须要把该便 ...