Guass消元 poj 1830

代码是拷别人的。。。。
举个例子
当n=3时
灯1  (a11*x1)^(a12*x2)^(a13*x3)=sta[0]^end[0]
灯2  (a21*x1)^(a22*x2)^(a23*x3)=sta[1]^end[1]
灯3  (a31*x1)^(a32*x2)^(a33*x3)=sta[2]^end[2]
用sta和end数组表示初始状态和最终状态 sta[i]^end[i]的结果就是中间操作数异或的结果，用b[i]=sta[i]^end[i]
xi=0时表示开关i没有操作，xi=1是表示开关i操作一次
aij 表示 灯j会影响灯i  这里需要仔细理解，很重要。。。


用高斯消元法求出求出增广矩阵，
对于增广矩阵
当某一行的前n个数都为0，第n+1个数不为0时，方程无解
r()表示矩阵的秩

若r(A)=r(A,b)=n则有唯一解。若r(a,b)=r(a)<n时，有多组解，个数为2^(n-r)，为什么？我也不知道啊，线性代数没学好啊

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;
const int Max = 31;
int n,a[Max][Max],sta[Max],end[Max];
inline int fab(int x) { return( x > 0 ? x : -x); }
int gauss()   // 高斯消元法
{
   int i,j,k;
   for(i = 0, j = 0; i < n,j < n;)
   {
      int id = i;
      for(k = i+1; k < n; k++)
         if(fab(a[k][j]) > fab(a[id][j]))
            id = k;
      if(id != i)
      {
         for(k = j; k <= n; k++)
            swap(a[i][k],a[id][k]);
      }
      if(a[i][j] == 0) { j++; continue; }
      for(k = i+1; k < n; k++)
      {
         if(a[k][j] == 0) continue;
         for(int l = j; l <= n; l++)
            a[k][l] = a[k][l] ^ a[i][l];  //此处做了特殊处理，不跟传统的高斯一样。
      }
      i++, j++;
   }
   //printf("i=%d\n",i);
   for(int k = i; k < n; k++) if(a[k][n] != 0) return(-1);

   return(1<<(n-i));   // i表示矩阵的秩
}
int main()
{
   int t;
   scanf("%d",&t);
   while(t--)
   {
      scanf("%d",&n);
      for(int i = 0; i < n; i++) scanf("%d",&sta[i]);
      for(int i = 0; i < n; i++) scanf("%d",&end[i]);
      int k,j;
      memset(a,0,sizeof(a));
      while(scanf("%d%d",&k,&j) == 2 && k && j) a[j-1][k-1] = 1;
      for(int i = 0; i < n; i++)
      {
         a[i][n] = sta[i] ^ end[i];
         a[i][i] = 1;     // 灯i对灯i有影响
      }
      int ans = gauss();
      if(ans == -1) printf("Oh,it's impossible~!!\n");
      else printf("%d\n",ans);
   }
   return(0);
}