POJ1830，01矩阵高斯消元

拿到这道题后，没啥思路，后来看了几份解题报告后，才发现这题用线性方程组来解，终于知道了线性代数在ACM的用武之地了。由于是中文题，题目解释略。直接说解题思路：设A矩阵是灯泡关系矩阵，A[[I][J]表明：对第J个灯实行一次操作后，第I个灯的状态也随之改变。B为结果矩阵，也就是灯的状态是否改变，改变为1，不变为0.这样AX=B求X即为题目中的解。本题求的是解的个数，有线性代数知识所得：当AX=B增广矩阵的秩为n时，若原矩阵的秩不为N，则无解。若r(A)=r(A,b)=n则有唯一解。若r(a,b)=r(a)<n时，有多组解，个数为2^(n-r).特别注意A[X][Y]的含义，写代码时很有可能搞反。。。

下面是代码：

#include <iostream>
#include <cmath>
#include <cstdlib>
#include <cstring>
#include <cstdio>

using namespace std;

int a[100][100],s[100],e[100];
int n;

void swap(int &a,int &b)
{
     int tmp=a;
     a=b;
     b=tmp;
}

int gauss()
{
    int i,j,k,l;
    for (i=0,j=0;i<n,j<n;i++,j++)
    {
        int row=i;
        for (k=i+1;k<n;k++)
            if (a[k][j]>a[row][j]) row=k;
        if (row!=i)
        {
                   for (k=j;k<=n;k++) swap(a[i][k],a[row][k]);
        }
        if (a[i][j]==0)
        {
                       i--;
                       continue;
        }        
        for (k=i+1;k<n;k++)
        {
            if (a[k][j]==0) continue;
            for (l=j;l<=n;l++) a[k][l]=(a[k][l]-a[i][l]+2)%2;
        }
    }
    for (k=i;k<n;k++)
        if (a[k][n]!=0) return -1;
    if (i==n) return 1;
    else return pow(2.0,n-i); 
}

int main()
{
    int T,i,k,j,r;
    cin>>T;
    while (T>0)
    {
          cin>>n;
          memset(a,0,sizeof(a));
          for (i=0;i<n;i++) cin>>s[i];
          for (i=0;i<n;i++) 
          {
              cin>>e[i];
              a[i][n]=(e[i]-s[i]+2)%2;
              a[i][i]=1;
          }
          do
          {
              int x,y;
              cin>>x>>y;
              if (x+y==0) break;
              a[y-1][x-1]=1;
          }
          while (true);
          if ((r=gauss())<0)
             cout<<"Oh,it's impossible~!!"<<endl;
          else cout<<r<<endl;
          T--;          
    }
    return 0;
}