测试地址: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;
}