POJ 1830 开关问题

POJ_1830

    可以用高斯消元求矩阵的秩，这样就知道了有x个变元，而这些变元可以取任意值，对于每种情况，其他元会有唯一的解，所以一共就有2^x种方案。

#include<stdio.h>
#include<string.h>
#include<iostream>
#define MAXD 35
using namespace std;
int N, mat[MAXD][MAXD];
void init()
{
    int i, j, k;
    memset(mat, 0, sizeof(mat));
    scanf("%d", &N);
    for(i = 0; i < N; i ++)
        scanf("%d", &mat[i][N]);
    for(i = 0; i < N; i ++)
    {
        scanf("%d", &k);
        mat[i][N] ^= k;
        mat[i][i] = 1;
    }
    for(;;)
    {
        scanf("%d%d", &i, &j);
        if(!i && !j)
            break;
        mat[j - 1][i - 1] = 1;
    }
}
int gauss()
{
    int i, j, k, x, y;
    for(i = j = 0; j < N; i ++, j ++)
    {
        if(mat[i][j] == 0)
        {
            for(k = i + 1; k < N; k ++)
                if(mat[k][j])
                    break;
            if(k == N)
            {
                -- i;
                continue;
            }
            for(y = j; y <= N; y ++)
                swap(mat[i][y], mat[k][y]);
        }
        for(x = i + 1; x < N; x ++)
            if(mat[x][j])
            {
                for(y = j; y <= N; y ++)
                    mat[x][y] ^= mat[i][y];
            }
    }
    for(x = i; x < N; x ++)
        if(mat[x][N])
            return -1;
    return 1 << (N - i);
}
void solve()
{
    int ans = gauss();
    if(ans == -1)
        printf("Oh,it's impossible~!!\n");
    else
        printf("%d\n", ans);
}
int main()
{
    int t;
    scanf("%d", &t);
    while(t --)
    {
        init();
        solve();
    }
    return 0;
}