跳转至

线性基

回想高中数学立体几何中基向量的概念,我们可以在三维欧氏空间中找到一组基向量 ,之后空间中任意一个向量都可以由这组基向量表示。换句话说,我们可以 通过有限的基向量来描述无限的三维空间,这足以体现基向量的重要性。

三维欧氏空间是特殊的 线性空间,三维欧氏空间的基向量在线性空间中就被推广为了线性基。

OI 中有关线性基的应用一般只涉及两类线性空间: 维实线性空间 布尔域 线性空间 ,我们会在 应用 一节中详细介绍。若您不熟悉线性代数,则推荐从应用部分开始阅读。

以下会从一般的线性空间出发来介绍线性基,并给出线性基的常见性质。

前置知识:线性空间

线性基是线性空间的一组基,是研究线性空间的重要工具。

定义

称线性空间 的一个极大线性无关组为 的一组 Hamel 基线性基,简称

规定线性空间 的基为空集。

可以证明任意线性空间均存在线性基1,我们定义线性空间 维数 为线性基的元素个数(或势),记作

性质

  1. 对于有限维线性空间 , 设其维数为 , 则:

    1. 中的任意 个向量线性相关。

    2. 中的任意 个线性无关的向量均为 的基。

    3. 中的任意向量均可被向量组 线性表出,则其是 的一个基。

      证明

      任取 中的一组基 , 由已知条件,向量组 可被 线性表出,故

      因此

    4. 中任意线性无关向量组 均可通过插入一些向量使得其变为 的一个基。

  2. (子空间维数公式)令 是关于 的有限维线性空间,且 也是有限维的,则

    证明

    ,,.

    的一组基 , 将其分别扩充为 中的基:.

    接下来只需证明向量组 线性无关即可。

    .

    .

    注意到上式左边在 中,右边在 中,故两边均在 中,因此

    , 进而

  3. 是关于 的有限维线性空间,且 也是有限维的,则下列诸款等价:

    1. .
    2. .
    3. 的一组基, 的一组基,则 的一组基。
    Note

    1,3 两条可推广到无限维线性空间中

例子

考虑 的基。

  1. 如图

    是一组基。

  2. 如图

    是一组基。

  3. 如图

    不是一组基,因为 .

  4. 如图

    不是一组基,因为 .

正交基与单位正交基

若线性空间 的一组基 满足 (即两两正交),则称这组基是 正交基

若线性空间 的一组正交基 还满足 ,则称这组基是 单位正交基

任意有限维线性空间 的基都可以通过 Schmidt 正交化 变换为正交基。

应用

根据前文内容,我们可以利用线性基实现:

  1. 求给定向量组的秩;
  2. 对给定的向量组,找到一组极大线性无关组(或其张成的线性空间的一组基);
  3. 向给定的向量组插入某些向量,在插入操作后的向量组中找到一组极大线性无关组(或其张成的线性空间的一组基);
  4. 对找到的一组极大线性无关组(或基),判断某向量能否被其线性表出;
  5. 对找到的一组极大线性无关组(或基),求其张成的线性空间中的特殊元素(如最大元、最小元等)。

在 OI 中,我们一般将 维实线性空间 下的线性基称为 实数线性基 维布尔域线性空间 下的线性基称为 异或线性基

Tip

中的加法为异或,乘法为与,可以证明 是域。

可以证明代数系统 是线性空间,其中:

即加法是异或,数乘是与。

以异或线性基为例,我们可以根据给定的一组布尔序列 构造出一组异或线性基 ,这组基有如下性质:

  1. 中任意非空子集的异或和不为
  2. 中的任意元素 ,都可在 中取出若干元素使其异或和为
  3. 对任意满足上两条的集合 ,其元素个数不会小于 的元素个数。

我们可以利用异或线性基实现:

  1. 判断一个数能否表示成某数集子集的异或和;
  2. 求一个数表示成某数集子集异或和的方案数;
  3. 求某数集子集的最大/最小/第 大/第 小异或和;
  4. 求一个数在某数集子集异或和中的排名。

构造方法

因为异或线性基与实数线性基没有本质差别,所以接下来以异或线性基为例,实数线性基版本的代码只需做一点简单修改即可。

贪心法

对原集合的每个数 转为二进制,从高位向低位扫,对于第 位是 的,如果 不存在,那么令 并结束扫描,如果存在,令

查询原集合内任意几个元素 的最大值,只需将线性基从高位向低位扫,若 上当前扫到的 答案变大,就把答案异或上

为什么能行呢?因为从高往低位扫,若当前扫到第 位,意味着可以保证答案的第 位为 ,且后面没有机会改变第 位。

查询原集合内任意几个元素 的最小值,就是线性基集合所有元素中最小的那个。

查询某个数是否能被异或出来,类似于插入,如果最后插入的数 被异或成了 ,则能被异或出来。

代码(洛谷 P3812 【模板】线性基
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
#include <bits/stdc++.h>
using ull = unsigned long long;

ull p[64];

void insert(ull x) {
  for (int i = 63; ~i; --i) {
    if (!(x >> i))  // x 的第 i 位是 0
      continue;
    if (!p[i]) {
      p[i] = x;
      break;
    }
    x ^= p[i];
  }
}

int main() {
  int n;
  scanf("%d", &n);
  ull a;
  for (int i = 1; i <= n; ++i) {
    scanf("%llu", &a);
    insert(a);
  }
  ull ans = 0;
  for (int i = 63; ~i; --i) {
    ans = std::max(ans, ans ^ p[i]);
  }
  printf("%llu\n", ans);
  return 0;
}

高斯消元法

高斯消元法相当于从线性方程组的角度去构造线性基,正确性显然。

代码(洛谷 P3812 【模板】线性基
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
#include <bits/stdc++.h>
using ull = unsigned long long;
const int MAXN = 1e5 + 5;

ull deg(ull num, int deg) { return num & (1ull << deg); }

ull a[MAXN];

int main() {
  int n;
  scanf("%d", &n);
  for (int i = 1; i <= n; ++i) scanf("%llu", &a[i]);
  int row = 1;
  for (int col = 63; ~col && row <= n; --col) {
    for (int i = row; i <= n; ++i) {
      if (deg(a[i], col)) {
        std::swap(a[row], a[i]);
        break;
      }
    }
    if (!deg(a[row], col)) continue;
    for (int i = 1; i <= n; ++i) {
      if (i == row) continue;
      if (deg(a[i], col)) {
        a[i] ^= a[row];
      }
    }
    ++row;
  }
  ull ans = 0;
  for (int i = 1; i < row; ++i) {
    ans ^= a[i];
  }
  printf("%llu\n", ans);
  return 0;
}

性质

贪心法构造的线性基具有如下性质:

  • 线性基没有异或和为 的子集。
  • 线性基中各数二进制最高位不同。

高斯消元法构造出的线性基满足如下性质:

  • 高斯消元后的矩阵是一个行简化阶梯形矩阵。

    该性质包含了贪心法构造的线性基满足的两条性质

    如果不理解这条性质的正确性,可以跳转 高斯消元

提供一组样例:

1
2
5
633 211 169 841 1008

二进制表示:

1
2
3
4
5
1001111001
0011010011
0010101001
1101001001
1111110000

贪心法生成的线性基:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
1001111001
0100110000
0011010011
0001111010
0000000000
0000010000
0000000000
0000000000
0000000000
0000000000

高斯消元法生成的线性基:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
1000000011
0100100000
0010101001
0001101010
0000010000
0000000000
0000000000
0000000000
0000000000
0000000000

这是一条非常好的性质,能帮我们更方便的解决很多问题。比如:给定一些数,选其中一些异或起来,求异或最大值,如果用贪心法构造线性基,需要再做一遍贪心,如果 ans 的当前位是 0,那么异或一定会更优,否则当前位如果为 1,则一定不会更优;而使用高斯消元法构造线性基后直接将线性基中所有元素都异或起来输出即可。

对于其他比较经典的问题(查询一个数能否被异或得到,查询能被异或得到的第 大数等),高斯消元法得到的线性基也能更加方便地解决。

时间复杂度

设向量长度为 , 总数为 , 则时间复杂度为 . 其中高斯消元法的常数略大。

若是实数线性基,则时间复杂度为 .

练习题

参考资料与注释

  1. 丘维声,高等代数(下)。清华大学出版社。
  2. Basis (linear algebra) - Wikipedia
  3. Vector Basis -- from Wolfram MathWorld