[高斯消元] poj 1830 开关问题

来源：互联网发布：淘宝运营是做什么的编辑：程序博客网时间：2024/05/16 05:09

题意：给定N(N < 29)个开关，每个开关打开和关闭的时候会引起另外一个开关的变化，本来为打开的会变成关闭，本来关闭的会变成打开。给定N个开关的初始状态和终止状态，以及关联的开关关系，求共有多少种方案从初始状态变成终止状态（不计顺序，并且每个开关只能操作至多一次）。

题解：由于开关只有打开和关闭两种状态，所以对于每个开关的打开和关闭，组合一下总共有2^N种情况，枚举所有情况判可行性，对于这个数据量来说是不现实的，需要想办法优化。

我们用X[i]来表示第i个开关的操作状态（1表示操作，0表示不操作）。

第i个开关会被哪些开关影响是可以知道的（这个关系在输入的时候会给出），假设影响第i个开关的开关列表为L[i][0], L[i][1], L[i][2]... 第i个开关的起始状态为S[i]，终止状态为E[i]，则可以列出N个方程：

( X[0] * A[i][0] + X[1] * A[i][1] + ... + X[n-1] * A[i][n-1] ) % 2 = (E[i] - S[i]);

每个方程代表一个开关被它本身以及它的关联开关操作后的最终状态，系数矩阵A[i][j]代表了开关之间的连带关系:

1) 如果第j个开关的操作能够影响第i个开关的状态，那么A[i][j] = 1;

2) 如果第j个开关的操作不影响第i个开关的状态，那么A[i][j] = 0;

3) 特殊的A[i][i] = 1(开关本身的操作必然会影响自己的当前状态);

X[i]取值为0或1，这样就是N个N元一次方程组，利用高斯消元求解即可。将增广矩阵化简为上三角的形式后，剩余全为0的行的个数为自由变元的个数F（自由变元就是它在取值范围内可以取任意值，这题是个方阵，所以自由变元的个数等于全为0的行的个数），所以，由于开关一共两种状态，取值为0和1，所以总的解的个数为2^F。特殊的，如果某一行系数全为零，而增广矩阵最后一列对应行的值不为0，则表示无解。

以上内容摘自：http://www.cppblog.com/menjitianya/archive/2014/06/08/207226.html

#include"cstdlib"#include"cstdio"#include"cstring"#include"cmath"#include"queue"#include"algorithm"#include"iostream"using namespace std;struct matrix{    int mat[55][55];} a,b;int equ,var;__int64 poww(int x){    __int64 ans=1;    while(x--) ans*=2;    return ans;}void debug(){    for(int i=0; i<equ; i++)    {        for(int j=0; j<=var; j++) printf("%d ",b.mat[i][j]);        puts("");    }    puts("");}int gauss(){    int i,j;    int row,col;    for(row=0,col=0; row<equ&&col<var; col++,row++)    {        for(i=row; i<equ; i++) if(b.mat[i][col]) break;        if(i==equ)        {            row--;            continue;        }        for(j=0; j<=var; j++) swap(b.mat[i][j],b.mat[row][j]);        for(i=row+1; i<equ; i++)        {            if(b.mat[i][col])                for(j=col; j<=var; j++) b.mat[i][j]^=b.mat[row][j];        }    }    for(i=row; i<equ; i++) if(b.mat[i][var]) return -1;    return var-row;}int main(){    int t;    cin>>t;    while(t--)    {        int n,i;        int s[33];        cin>>n;        equ=var=n;        memset(a.mat,0,sizeof(a.mat));        for(i=0;i<n;i++) cin>>s[i];        for(i=0;i<n;i++)        {            int xx;cin>>xx;            s[i]^=xx;            a.mat[i][var]=s[i];            a.mat[i][i]=1;        }        int xx,yy;        while(scanf("%d%d",&xx,&yy),(xx+yy))        {            a.mat[--yy][--xx]=1;        }        b=a;        //debug();        int ans=gauss();        if(ans==-1) puts("Oh,it's impossible~!!");        else printf("%I64d\n",poww(ans));    }}

0 0