NOIP 考前 高斯消元练习

POJ 1830

列出n个方程右边为最后的情况

每一行代表第几个灯的情况，每一行代表是否按第几个按钮写出方程即可。

#include <cstdio>
#include <cstring>
const int Maxn=100;
int M[Maxn][Maxn],Ans,a[Maxn],b[Maxn],p,q,KASE,n;
inline void Swap(int &x,int &y) {int t=x;x=y;y=t;}
inline bool Check(int x)
{
    for (int i=1;i<=n;i++) if (M[x][i]) return false;
    return true;
}
bool Gauss()
{
    int Pos=1;
    for (int i=1;i<=n;i++)
    {
        int k=0;
        for (int j=Pos;j<=n;j++) if (M[j][i]) {k=j; break;}
        if (k==0) continue;
        for (int j=1;j<=n+1;j++) Swap(M[Pos][j],M[k][j]);
        for (int j=1;j<=n;j++)
            if (M[j][i] && j!=Pos)
                for (int l=1;l<=n+1;l++) M[j][l]^=M[Pos][l];
        Pos++;
    }
    int cnt=0;
    for (int i=1;i<=n;i++)
        if (Check(i))
        {
            if (M[i][n+1]) return false;
            cnt++;
        }
    return Ans=(1<<cnt);
}
int main()
{
    scanf("%d",&KASE);
    for (int Kase=1;Kase<=KASE;Kase++)
    {
        scanf("%d",&n);
        memset(M,0,sizeof(M));
        for (int i=1;i<=n;i++) scanf("%d",&a[i]);
        for (int i=1;i<=n;i++) scanf("%d",&b[i]);
        for (int i=1;i<=n;i++) M[i][n+1]=a[i]^b[i];
        for (int i=1;i<=n;i++) M[i][i]=1;
        while (scanf("%d%d",&p,&q)!=EOF)
        {
            if (p==0 && q==0) break;
            M[q][p]=1;
        }
        if (Gauss()) printf("%d
",Ans); else puts("Oh,it's impossible~!!");
    }
    return 0;
}