POJ 1830 开关问题 解题报告（高斯消元法）

    题目大意：中文题……

    解题报告：问我们有多少种可行方法。仔细想一下，如果矩阵有唯一解，那么只有一种可行方法。如果有一个变元，变元取0或者1对其他非变元的值会有影响，但仍然都是可行解。换句话说，变元有多少种，就有2的多少次方个解。简单算一下就知道了。

    另外，给定的是开关的初状态和末状态，我们可以记录哪些状态是改变的，记为1。未改变的记为0。用高斯消元法计算变元数量即可。

    代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cmath>
using namespace std;

const int maxn = 100+10;
typedef int Matrix[maxn][maxn];

int Gauss(Matrix a, int equ, int var)
{
    int k=0, col=0;
    for(;k<equ && col<var;col++)
    {
        int r=k;
        for(int i=k;i<equ;i++) if(a[i][col])
            r=i;

        if(a[r][col]==0)
            continue;

        if(r!=k) for(int j=col;j<=var;j++)
            swap(a[k][j],a[r][j]);

        for(int i=k+1;i<equ;i++) if(a[i][col])
            for(int j=col;j<=var;j++)
                a[i][j]^=a[k][j];
        k++;
    }

    for(int i=k;i<equ;i++) if(a[i][col])
        return -1;
    return var-k;
}

void work()
{
    int n;
    scanf("%d", &n);

    int sw[30], tmp;
    for(int i=0;i<n;i++)
        scanf("%d", sw+i);
    for(int i=0;i<n;i++)
        scanf("%d", &tmp), sw[i]^=tmp;

    Matrix A;
    memset(A, 0, sizeof(A));
    for(int i=0;i<n;i++)
        A[i][i]=1, A[i][n]=sw[i];

    int u, v;
    while(~scanf("%d%d", &u, &v) && (u||v))
        A[v-1][u-1]=1;
    int ans = Gauss(A, n, n);

    if(ans==-1)
        puts("Oh,it's impossible~!!");
    else
        printf("%d\n", 1<<(ans));
}

int main()
{
    int T;
    scanf("%d", &T);
    while(T--)
        work();
}