POJ 1830 开关问题

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;
const int N=51;
int T,n,m,start[N],end[N],a[N][N];//增广矩阵 
int Gauss(int equ,int var){
    int i,j,k,max_r,col=0;
    for(k=0;k<equ&&col<var;k++,col++){
        max_r=k;//当前这列绝对值最大的行  
        for(i=k+1;i<equ;i++){
            if(abs(a[i][col])>abs(a[max_r][col])){
                max_r=i;
            } 
        } 
        if(max_r!=k){//将该行与当前行交换
            for(j=k;j<=var;j++){
                swap(a[k][j],a[max_r][j]);
            } 
        } 
        if(!a[k][col]){//当前col列k行以下全为0(包括k行) 
            k--;continue;
        }
        for(i=k+1;i<equ;i++){
            if(a[i][col]){//如果不是0，则消除
                for(j=col;j<=var;j++){
                    a[i][j]^=a[k][j];
                } 
            }
        }
    }
    for(i=k;i<equ;i++){
        if(a[i][col]){//无解 
            return -1;
        } 
    } 
    return var-k;//var-k个变元,只有0/1两种取值 
}
int main(){
    for(scanf("%d",&T);T--;){
        scanf("%d",&n);memset(a,0,sizeof a);
        for(int i=0;i<n;i++) scanf("%d",&start[i]);
        for(int i=0;i<n;i++) scanf("%d",&end[i]);
        for(int u,v;scanf("%d%d",&u,&v)!=EOF&&u&&v;a[v-1][u-1]=1);//每行两个数u v，表示如果操作第u个开关，第v个开关的状态也会变化”，说明v受到u的影响,a[v-1][u-1]=1 
        for(int i=0;i<n;i++) a[i][i]=1;//自己可以改变自己
        for(int i=0;i<n;i++) a[i][n]=start[i]^end[i];
        int ans=Gauss(n,n);
        if(ans==-1) puts("Oh,it‘s impossible~!!");
        else printf("%d\n",1<<ans);
    }
    return 0;
}