高斯消元-poj1222熄灯问题状态压缩解法

最新推荐文章于 2020-02-28 17:52:34 发布

weixin_30491641

最新推荐文章于 2020-02-28 17:52:34 发布

阅读量153

点赞数

原文链接：http://www.cnblogs.com/zsben991126/p/10660498.html

版权

有点自闭的。。为什么我最后的答案是倒着来的啊。。

搞明白了：因为一开始构造的系数就是反着的，，所以主元也倒过来了。。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int main(){
    int t;
    cin>>t;
    for(int tt=1;tt<=t;tt++){
        long long tmp,ans[100]={},A[100]={};
        for(int i=1;i<=30;i++){
            cin>>tmp;
            A[i]^=tmp;
        }
        for(int i=1;i<=5;i++)
            for(int j=1;j<=6;j++){
                int c=(i-1)*6+j;
                A[c]|=1<<c;
                if(j<=5)A[c]|=1<<(c+1);
                if(j>=2)A[c]|=1<<(c-1);
                if(i<=4)A[c]|=1<<(c+6);
                if(i>=2)A[c]|=1<<(c-6);
            }
        for(int i=1;i<=30;i++){
            for(int j=i;j<=30;j++) 
                if(A[j]>A[i])swap(A[j],A[i]);
            if(A[i]==0||A[i]==1)break;
            for(int k=30;k;k--)
                if(A[i]>>k & 1){ 
                    for(int j=1;j<=30;j++)
                        if(j!=i && (A[j]>>k & 1))A[j]^=A[i]; 
                    break;
                }
        }    
        printf("PUZZLE #%d\n",tt);
        for(int i=30;i>=1;i--){//为什么答案是倒的啊。。 
            cout<<(A[i] & 1);
            if(i%6==1)puts("");
            else cout<<" ";
        }
    }
}