如何实现二元方程的高斯消元答案

【问题标题】：how to implementing gaussian elimination for binary equations如何实现二元方程的高斯消元
【发布时间】：2012-07-14 12:54:38
【问题描述】：

我有这个方程组1=x⊕y⊕z
1=x⊕y⊕w
0=x⊕w⊕z
1=w⊕y⊕z

我正在尝试实施高斯消除来解决该系统，如 here 所述，用 XOR 替换除法、减法和乘法，但它给出了我的错误答案..正确答案是 (x,y,z,w)=(0,1,0,0)
我做错了什么？

public static void ComputeCoefficents(byte[,] X, byte[] Y)
    {
        int I, J, K, K1, N;
        N = Y.Length;
        for (K = 0; K < N; K++)
        {
            K1 = K + 1;
            for (I = K; I < N; I++)
            {
                if (X[I, K] != 0)
                {
                    for (J = K1; J < N; J++)
                    {
                        X[I, J] /= X[I, K];
                    }
                    //Y[I] /= X[I, K];
                    Y[I] ^= X[I, K];

                }
            }
            for (I = K1; I < N; I++)
            {
                if (X[I, K] != 0)
                {
                    for (J = K1; J < N; J++)
                    {
                        X[I, J] ^= X[K, J];
                    }
                    Y[I] ^= Y[K];
                }
            }
        }
        for (I = N - 2; I >= 0; I--)
        {
            for (J = N - 1; J >= I + 1; J--)
            {
                //Y[I] -= AndOperation(X[I, J], Y[J]);
                Y[I] ^= (byte)(X[I, J]* Y[J]);

            }
        }
    }

【问题讨论】：

似乎有些可疑。乘以 1 应该使值保持不变 - 如果您将其替换为 xor，乘以 1 将反转所有内容。通常，如果用逻辑运算符替换算术，则加法将替换为or 或xor。 or 没有简单的逆，但xor 是它自己的逆。乘法替换为and。 and 没有简单的逆运算，但在这种情况下您不需要除法。你甚至不需要乘法。
此算法仅在特定条件下有效。它无论如何都不稳定。您可以离开第一个 for 循环，因为除法只会执行 1，这不会导致任何变化。那么，在我看来，剩下的for (J = K1; J < N; J++)一定是for (J = 0; J < N; J++)。但同样，该算法并不稳定，您应该真正使用库来完成该任务。或者，如果您的问题不那么普遍，也许我们可以编写一个算法。

标签： c# math binary

【解决方案1】：

我认为您正在尝试为此应用高斯消除模型 2。

一般来说，如果你的方程是这样的，你可以做高斯消除 mod k

a_1 * x + b_1 * y + c_1 * z = d_1
a_2 * x + b_2 * y + c_2 * z = d_2
a_3 * x + b_3 * y + c_3 * z = d_3
a_4 * x + b_4 * y + c_4 * z = d_4

在 Z2 中 * 是 and 和 + 是 xor，所以你可以使用高斯消元法来求解形式的方程

x (xor) y (xor) z   = 1
x (xor) y (xor) w   = 1 
x (xor) z (xor) w   = 0
y (xor) z (xor) w   = 1

让我们手动使用高斯消元来完成这个方程。

对应的增广矩阵为：

 1 1 1 0 | 1
 1 1 0 1 | 1
 1 0 1 1 | 0
 0 1 1 1 | 1

 1 1 1 0 | 1
 0 0 1 1 | 0   (R2 = R2 + R1)
 0 1 0 1 | 1   (R3 = R3 + R1)
 0 1 1 1 | 1

 1 1 1 0 | 1
 0 1 1 1 | 1   (R2 = R4)
 0 1 0 1 | 1   
 0 0 1 1 | 0   (R4 = R2)

 1 0 0 1 | 0   (R1 = R1 + R2)
 0 1 1 1 | 1   
 0 0 1 0 | 0   (R3 = R3 + R2)   
 0 0 1 1 | 0   

 1 0 0 1 | 0
 0 1 0 1 | 1   (R2 = R2 + R3)  
 0 0 1 0 | 0      
 0 0 0 1 | 0   (R4 = R4 + R3)

 1 0 0 0 | 0   (R1 = R1 + R4)
 0 1 0 0 | 1   (R2 = R2 + R4)  
 0 0 1 0 | 0      
 0 0 0 1 | 0

给出 (x,y,z,w) = (0,1,0,0) 的解。

但这需要行旋转 - 我在您的代码中看不到。

您的代码中还有一些乘法和除法可能不需要存在。我希望代码如下所示：（您需要修复 TODO）。

public static void ComputeCoefficents(byte[,] X, byte[] Y) {
  int I, J, K, K1, N;
  N = Y.Length;

  for (K = 0; K < N; K++) {
    //First ensure that we have a non-zero entry in X[K,K]
    if( X[K,K] == 0 ) {
      for(int i = 0; i<N ; ++i ) { 
        if(X[i,K] != 0 ) {
             for( ... ) //TODO: A loop to swap the entries
             //TODO swap entries in Y too
           }
      }
    if( X[K,K] == 0 ) {
       // TODO: Handle the case where we have a zero column 
       //      - for now we just move on to the next column
       //      - This means we have no solutions or multiple 
       //        solutions
       continue
    }

    // Do full row elimination.
    for( int I = 0; I<N; ++I)
    {
       if( I!=K ){ //Don't self eliminate
         if( X[I,K] ) { 
           for( int J=K; J<N; ++J ) { X[I,J] = X[I,J] ^ X[K,J]; }
           Y[J] = Y[J] ^ Y[K];
         }
       }
    }
  }

  //Now assuming we didnt hit any zero columns Y should be our solution.

}

【讨论】：