poj 1222 高斯消元（改改模板）

#include <cstdio>
#include <cstring>
#include <iostream>
#include <stack>
#include <queue>
#include <map>
#include <algorithm>
#include <vector>
#include <cstdlib>
#include <cmath>

using namespace std;

const int maxn = 1000005;

typedef long long LL;

int A[35][35];
int ans[35];

void gauss(int m,int n) //m个方程,n个变量
{
int i,j,k,r,u;
for(i=0,j=0;i<m&&j<n;j++){ //当前处理i个方程j个变量
r = i;
for(k = i;k < m;k++){
if(A[k][j]){
r = k;
break;
}
}
if(A[r][j]){
for(k = 0;k <= n;k++) swap(A[r][k],A[i][k]);
for(u = 0;u < m;u++){
if(u!=i && A[u][j]) for(k = i;k <= n;k++) A[u][k] ^= A[i][k];
}
i++;
}

}
for(int i=0;i<30;i++){
ans[i] = A[i][30];
}

}

int main()
{
int t;
scanf("%d",&t);
int cas = 0;
while(t--){
cas++;
memset(A,0,sizeof(A));
for(int i=0;i<30;i++){
scanf("%d",&A[i][30]);
//ans[i] = 0;
}
int x,xx,y,yy; //构造矩阵
for(int i=0;i<30;i++){
x = i/6;
y = i%6;
for(int j=0;j<30;j++){
xx = j/6;
yy = j%6;
if((abs(x - xx)+abs(y - yy)) <= 1) A[i][j] = 1;
else A[i][j] = 0;
}
}
gauss(30,30);
printf("PUZZLE #%d\n",cas);
for(int i=0;i<30;i++){
printf("%d",ans[i]);
if((i+1)%6==0) printf("\n");
else printf(" ");
}
}
return 0;
}