HDU---3364：Lanterns【高斯消元】

题意：

开关问题的变形，此题求解集的数量

分析:

构造N * M 的矩阵高斯消元，解集和自由变元的数量有关，因为每个自由变元可以任意取都能得到一组解集，由乘法原理知道：解集的数量为a^x（a为每个自由变元的取值个数，x为自由变元的数目）；其实对于一般的方程，自由变元的个数为 x =（变量数var - 非0系数向量的个数）；所以此题的答案为 2^x

代码：

#include <bits/stdc++.h>

using namespace std;

const int MAXN = 55;
int t,equ,var,k,q,lan,n,m;
int b[MAXN][MAXN],a[MAXN][MAXN];
void init()
{
    for(int i = 0;i < n; ++i){
        for(int j = 0;j < m; ++j)
            a[i][j] = b[i][j];
        a[i][m] = 0;
    }
}
int Gauss()
{
    int r,c;
    for(r = 0,c = 0;r < equ && c < var; ++r,++c)
    {
        int maxr = r;
        for(int i = r + 1;i < equ; ++i){
            if(a[i][c] > a[maxr][c]) maxr = i;
        }
        if(a[maxr][c] == 0){
            r--;
            continue;
        }
        if(maxr != r){
            for(int i = c;i <= var; ++i) swap(a[maxr][i],a[r][i]);
        }
        for(int i = r + 1;i < equ; ++i){
            if(!a[i][c]) continue;
            for(int j = c;j <= var; ++j)
                a[i][j] ^= a[r][j]; 
        }
    }
    for(int i = r;i < equ; ++i){
        if(a[i][c]) return -1;
    }
    return var - r;
}
int main()
{
    cin>>t;
    for(int Case = 1;Case <= t; ++Case)
    {
        cin>>n>>m;
        memset(b,0,sizeof(b));
        equ = n;var = m;
        for(int i = 0; i < m; ++i){
            cin>>k;
            for(int j = 1;j <= k; ++j){
                cin>>lan;
                b[lan-1][i] = 1;
            }
        }
        cin>>q;
        printf("Case %d:\n", Case);
        while(q--){
            init();
            for(int i = 0;i < n; ++i) cin>>a[i][m];
            int num = Gauss();
            if(num == -1) cout<<0<<endl;
            else cout<<(1LL<<num)<<endl;
        }
    }
    return 0;
}