POJ1830 开关问题 高斯消元

题目链接：http://poj.org/problem?id=1830

题目大意：有N个相同的开关，每个开关都与某些开关有着联系，每当你打开或者关闭某个开关的时候，其他的与此开关相关联的开关也会相应地发生变化，即这些相联系的开关的状态如果原来为开就变为关，如果为关就变为开。你的目标是经过若干次开关操作后使得最后N个开关达到一个特定的状态。对于任意一个开关，最多只能进行一次开关操作。你的任务是，计算有多少种可以达到指定状态的方法。（不计开关操作的顺序）。

分析：每个开关都有动和不动两个选择，所以可以用1 来表示动 ，0表示不动。

假如三个开关， 然后用一组向量来表示最终结果Li 表示第i个开关最终是否动了，即向量组 ( L1, L2, L3)

用（x1,x2,x3）来表示最终是动还是不动，那么就有

a11*x1 + a12*x2 +a13*x3 = L1 mod  2

a21*x1 + a22*x2 +a23*x3 = L2 mod  2

a31*x1 + a32*x2 +a33*x3 = L3 mod  2

如果 i 灯能是j 灯亮 ，那么a[j][i ] ==1 

实现代码如下：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=29;
int a[maxn][maxn+1],x[maxn],equ,var;
int pre[maxn]; //纪录开关起始的状态
void Init()
{
    int tmp,t1,t2;
    memset(a,0,sizeof(a));
    memset(x,0,sizeof(x));
    for(int i=0;i<var;i++)
      scanf("%d",&pre[i]);
    for(int i=0;i<var;i++)
    {
        scanf("%d",&tmp);
        a[i][var]=pre[i]^tmp;
        a[i][i]=1;
    }
    while(scanf("%d%d",&t1,&t2))
    {
        if(t1==0&&t2==0) break;
        a[t2-1][t1-1]=1;
    }
}
void Debug()
{
    for(int i=0;i<equ;i++)
    {
        for(int j=0;j<=var;j++)
          cout<<a[i][j]<<" ";
        cout<<endl;
    }
    cout<<endl;
}
int gcd(int a,int b)
{
    if(a<0) return gcd(-a,b);
    if(b<0) return gcd(a,-b);
    return b==0?a:gcd(b,a%b);
}
void Gauss()
{
    int k,col=0;
    for(k=0;k<equ&&col<var;k++,col++)
    {
        int mx=k;
        for(int i=k+1;i<equ;i++)
           if(a[i][col]>a[mx][col]) mx=i;
        if(mx!=k)
        {
            for(int i=k;i<var+1;i++)
               swap(a[k][i],a[mx][i]);
        }
        if(!a[k][col])
        {
            k--; continue;
        }
        for(int i=k+1;i<equ;i++)
        if(a[i][col])
        {
            int lcm=a[k][col]/gcd(a[k][col],a[i][col])*a[i][col];
            int ta=lcm/a[i][col],tb=lcm/a[k][col];
            for(int j=col;j<var+1;j++)
               a[i][j]=((a[i][j]*ta-a[k][j]*tb)%2+2)%2;
        }
    }
    //Debug();
    for(int i=k;i<equ;i++)
      if(a[i][col])
      {
          puts("Oh,it's impossible~!!");
          return ;
      }
    printf("%d\n",1<<(var-k));
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&equ);
        var=equ;
        Init();
        Gauss();
    }
    return 0;
}