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)

下面是代码：

#include 
#include 
#include 
#include 
#include 

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;ia[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>T;
    while (T>0)
    {
          cin>>n;
          memset(a,0,sizeof(a));
          for (i=0;i>s[i];
          for (i=0;i>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~!!"<