高斯消元解异或线性方程组

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 110;

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

int gauss()
{
    int c, r; 
    for (c = 0, r = 0; c < n; c ++ )
    {
        int t = r;
        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]);
        
        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 ++ ;
    }
   
    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;
}