poj 1830(高斯消元求解。。。。)

理解了方程的一个细节。。。程序中给出。。。。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int map1[35][35], ans, m, n;
int a[35],b[35];
int gauss() {
    int i,j,k,t;
    for(i=0,j=0; i<m && j<m; i++,j++) {
        k=i;
        while(!map1[k][j] && k<m)
            k++;
        if(k==m){
            i--;continue;
        }
        if(k!=i)
            for(t = j; t <= m; t++)
                swap(map1[i][t],map1[k][t]);
        for(k=i+1; k < m; k++)
            if(map1[k][j])
                for(t=j;t<=m;t++)
                    map1[k][t]^=map1[i][t];
    }
    for(k=i;k<m;k++)
        if(map1[k][m])
            return -1;
    if(i==m) return 1;
    else
        return 1<<(m-i);
}

int main() {
    int cas, i, j, k, ans;
    scanf("%d", &cas);
    while (cas--) {
        scanf("%d", &m);
        memset(map1, 0, sizeof (map1));
        for(i=0;i<m;i++)
            scanf("%d",&a[i]);
        for(i=0;i<m;i++)
            scanf("%d",&b[i]);
        for(i=0;i<m;i++){
            map1[i][m]=(a[i]^b[i]);
            map1[i][i]=1;
        }
        while(scanf("%d%d",&j,&k) && (j||k))
            map1[k-1][j-1]=1; // 被动的那个放前面，不然会出错。。
        ans=gauss();
        if(ans==-1)
            printf("Oh,it's impossible~!!\n");
        else
            printf("%d\n",ans);
    }
    return 0;
}