[AcWing] 884. 高斯消元解异或线性方程组 （C++实现）高斯消元解异或线性方程组模板题

与高斯消元求解线性方程组差不多，详见 4. 可能有帮助的前置习题

1. 题目

2. 读题（需要重点注意的东西）

思路：

类似线性代数的求通解过程

本题的思路

3. 解法

---------------------------------------------------解法---------------------------------------------------

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 110;


int n;
int a[N][N];


int gauss()
{
    int c, r; // c 列 r 行
    // 枚举列，每次迭代进行如下操作
    for (c = 0, r = 0; c < n; c ++ )
    {
        int t = r;
        // 找到第一个非0行
        for (int i = r; i < n; i ++ )
            if (a[i][c])
                t = i;

        if (!a[t][c]) continue;
        
        // 将这行放到最高一行（交换）
        for (int i = c; i <= n; i ++ ) swap(a[r][i], a[t][i]);
        
        // 用枚举的这一行将下面所有行的首系数消为0
        for (int i = r + 1; i < n; i ++ )
            if (a[i][c])
                for (int j = n; j >= c; j -- )
                    a[i][j] ^= a[r][j];

        r ++ ;
    }
    
    // 如果消完系数后上三角的行数 < n，判断是无穷多解还是无解
    if (r < n)
    {
        for (int i = r; i < n; i ++ )
            if (a[i][n])
                // 无解
                return 2;
        // 无穷多解
        return 1;
    }

    // 否则 r = n，有唯一解，则倒着推一遍，算出所有的解
    for (int i = n - 1; i >= 0; i -- )
        for (int j = i + 1; j < n; j ++ )
            a[i][n] ^= a[i][j] * a[j][n];

    return 0;
}


int main()
{
    cin >> n;

    for (int i = 0; i < n; i ++ )
        for (int j = 0; j < n + 1; j ++ )
            cin >> a[i][j];

    int t = gauss();

    if (t == 0)
    {
        for (int i = 0; i < n; i ++ ) cout << a[i][n] << endl;
    }
    // 有无穷多解的情况
    else if (t == 1) puts("Multiple sets of solutions");
    // 无解的情况
    else puts("No solution");

    return 0;
}

---

4. 可能有帮助的前置习题

• [AcWing]883. 高斯消元解线性方程组（C++实现）高斯消元解线性方程组模板题

5. 所用到的数据结构与算法思想

• 高斯消元

6. 总结

高斯消元解异或线性方程组模板题，与高斯消元求解线性方程组差不多，理解思路并背下代码。