【HDU3364】Lanterns-高斯消元求异或方程组的解数

测试地址：Lanterns

题目大意：有N个灯和M个开关，每个开关控制若干个灯，开启开关会改变所有其控制的灯的状态，一开始所有灯都不亮，每个开关最多开启一次，问有多少种可能方案使所有灯达到某种状态？

做法：经典的求解异或方程组的问题。列N个方程，有M个未知数，组成一个系数矩阵A，如果开关i控制灯j，则Aj,i为1，否则为0，而第i个方程的常数项就是该灯最后要达到的状态。然而题目并不是要求我们找出一个解，而是求解数，我们就高斯消元将增广矩阵化为上三角矩阵，求出系数矩阵和增广矩阵的秩，如果相等则方程组有解，否则无解。设这个秩为k，我们可以看出有m-k个未知数是自由变元，它们可以任意取值，每一组自由变元的取值都确定其他未知数的值，所以答案就是2^(m-k)。还要注意每次我们求解方程组都会改变原来增广矩阵的值，所以要先把原先的增广矩阵存下来，每次询问再赋值给计算的矩阵。由于这一题m较小，在高斯消元的时候也可以把矩阵的每一行存成一个长整数，这样代码就可以变得异常简洁，但也要注意诸如1<<i这样的求值，实时计算的话当i太大时就会出错，所以应先算出这些2的次幂的值，存进一个数组bit，需要时直接调用即可。

以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define ll __int64
using namespace std;
int T,n,m,q;
ll A[110],B[110],sA[110],bit[110];

ll gauss()
{
  int k=1;
  for(int i=1;i<=m;i++)
  {
    for(int j=k;j<=n;j++)
	  if (A[j]&bit[i-1])
	  {
	    swap(A[j],A[k]),swap(B[j],B[k]);
		for(int u=k+1;u<=n;u++)
		  if (A[u]&bit[i-1]) A[u]=A[u]^A[k],B[u]=B[u]^B[k];
		k++;break;
	  }
  }
  for(int i=k;i<=n;i++)
    if (!A[i]&&B[i]) return 0;
  return bit[m-k+1];
}

int main()
{
  scanf("%d",&T);
  bit[0]=1;
  for(int i=1;i<=55;i++) bit[i]=bit[i-1]*2;
  for(int t=1;t<=T;t++)
  {
    scanf("%d%d",&n,&m);
	printf("Case %d:\n",t);
	for(int i=1;i<=n;i++) sA[i]=0;
	for(int i=1,k;i<=m;i++)
	{
	  scanf("%d",&k);
	  for(int j=1;j<=k;j++)
	  {
	    int a;
		scanf("%d",&a);
		sA[a]+=bit[i-1];
	  }
	}
	scanf("%d",&q);
	for(int i=1;i<=q;i++)
	{
	  for(int j=1;j<=n;j++)
	    scanf("%I64d",&B[j]),A[j]=sA[j];
	  printf("%I64d\n",gauss());
	}
  }
  
  return 0;
}