poj1222 EXTENDED LIGHTS OUT

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
int dir[5][2]={{0,0},{0,1},{0,-1},{1,0},{-1,0}};
int map[5][6];
int A[30][31];
#define in(x,y) (x>=0&&x<5&&y>=0&&y<6)
void change(int a[31],int b[31])
{
int c[31];
int i;
for(i=0;i<31;i++)
{
c[i]=a[i];
}
for(i=0;i<31;i++)
{
a[i]=b[i];
}
for(i=0;i<31;i++)
{
b[i]=c[i];
}
}
void jian(int a[31],int b[31],int k)
{
int i;
int n=(-1)*a[k]/b[k];
for(i=0;i<31;i++)
{
a[i]=(a[i]+(n*b[i]))%2;
}
}
int main()
{
//freopen("out.txt","w+",stdout);
int total;
int counter=1;
scanf("%d",&total);
while(total--)
{
printf("PUZZLE #%d\n",counter++);
int i,j,k;
for(i=0;i<5;i++)
{
for(j=0;j<6;j++)
{
scanf("%d",&map[i][j]);
}
}
int x[30];
int temp1,temp2;
memset(A,0,sizeof(A));
int np=0;
for(i=0;i<5;i++)
{
for(j=0;j<6;j++)
{
for(k=0;k<5;k++)
{
temp1=i+dir[k][0];
temp2=j+dir[k][1];
if(!in(temp1,temp2))
{
continue;
}
A[np][temp1*6+temp2]=1;
}
A[np][30]=map[i][j];
np++;
}
}
int maxi;
for(i=0,k=0;i<30&&k<30;i++,k++)//确定第i行，并将i+1,i+2,,,,29行的第k列清0
{
maxi=i;
for(j=i+1;j<30;j++)
{
if(abs(A[j][k])>abs(A[maxi][k]))
{
maxi=j;
}
}
//第j行的第k列为最大
if(A[maxi][k]==0)
{
i--;
continue;
}
j=maxi;
if(i!=j)
{
change(A[i],A[j]);
}
//已确定第i行
for(j=i+1;j<30;j++)
{
if(A[j][k]!=0)
{
jian(A[j],A[i],k);//用A[i]将A[j]的第k列清0
}
}
/*
int ii,jj;
for(ii=0;ii<30;ii++)
{
for(jj=0;jj<31;jj++)
{
printf("%d ",A[ii][jj]);
}
printf("\n");
}
printf("\n\n");
*/
}
for(i=29;i>=0;i--)
{
//A[i][i],A[i][i+1],
int temp;
temp=A[i][30];
for(j=i+1;j<=29;j++)
{
temp-=x[j]*A[i][j];
}
if(A[i][i]==0)
{
x[i]=0;
}
else
{
x[i]=temp/A[i][i];
}
}
for(i=0;i<5;i++)
{
for(j=0;j<5;j++)
{
printf("%d ",abs((x[i*6+j])%2));
}
printf("%d\n",abs((x[i*6+5])%2));
}
}
return 0;
}