高斯消元解異或方程組 模板

Code：

const int RNS = 55;
const int CNS = 55;
struct matrix{
    int row, col;
    int a[RNS][CNS], b[RNS];

    void init(int _row, int _col)
    {
        row = _row, col = _col;
        memset(a, 0, sizeof(a));
        memset(b, 0, sizeof(b));
    }
};

int guass_xor(matrix e)//m*n的矩陣，異或消元，返回自由變量數 無解返回-1
{
    int m = e.row, n = e.col;
    int i,j,k,r,u;
    for(i = j = 0; i<m && j<n; j++)
    {
        for(r=i, k=i;k<m;k++) if(e.a[k][j]) {r=k; break;}

        if(e.a[r][j])
        {
            if(r!=i)//提高數值穩定性
            {
                for(k=0;k<=n;k++) swap(e.a[r][k],e.a[i][k]);
                swap(e.b[r], e.b[i]);
            }
            for(u=i+1; u<m; u++)
            {
                if(e.a[u][j])
                {
                    for(k=i;k<=n;k++) e.a[u][k] ^= e.a[i][k];
                    e.b[u] ^= e.b[i];
                }
            }
            i++;//非全0式子數，即有界變量數
        }
    }

    for(k = i; k < m; k++) if(e.b[k]) return -1;
    return (n - i);
}