#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[33][33];
int Matrix[666][666];
int m,n,d,id;
bool inmap(int x,int y)
{
if(0<=x&&x<n&&0<=y&&y<m) return true;
else return false;
}
void f(int nn,int x,int y)
{
Matrix[nn][x*m+y]=1;
for(int l=0;l<=d;l++)
{
if(inmap(x+l,y+d-l))
{
Matrix[nn][(x+l)*m+y+d-l]=1;
}
if(inmap(x+l,y+l-d))
{
Matrix[nn][(x+l)*m+y+l-d]=1;
}
if(inmap(x-l,y+d-l))
{
Matrix[nn][(x-l)*m+y+d-l]=1;
}
if(inmap(x-l,y+l-d))
{
Matrix[nn][(x-l)*m+y+l-d]=1;
}
}
Matrix[nn][n*m]=a[x][y];
}
void build()
{
//对每个点建立一个方程
memset(Matrix,0,sizeof(Matrix));
id=0;
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
{
f(id++,i,j);
}
}
}
void output()
{
for(int i=0;i<id;i++)
{
for(int j=0;j<=id;j++)
{
cout<<Matrix[i][j]<<' ';
}
cout<<endl;
}
}
int gaosi_rank(int nn,int mm)
{
int i=0,j=0;
while(i<nn&&j<mm)
{
//找一行第j列不为0的
int r=i;
for(int k=i;k<nn;k++)
{
if(Matrix[k][j])
{
r=k;
break;
}
}
if(Matrix[r][j])//如果找到一行第j列不为0
{
if(r!=i)//不相等就替换
{
for(int k=0;k<=mm;k++) swap(Matrix[r][k],Matrix[i][k]);
}
for(int u=i+1;u<nn;u++)//
{
if(Matrix[u][j])
for(int k=j;k<=mm;k++)
{
Matrix[u][k]^=Matrix[i][k];
}
}
i++;
}
j++;
}
for(int k=i;k<nn;k++)
{
if(Matrix[k][mm])
{
return false;
}
}
return true;
}
int main()
{
while(scanf("%d%d%d",&m,&n,&d)!=EOF)
{
if(!n&&!m&&!d) break;
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
{
scanf("%d",&a[i][j]);
}
}
build();
if(gaosi_rank(id,id))
{
cout<<1<<endl;
}
else
{
cout<<0<<endl;
}
}
return 0;
}
bnu 11952 高斯消元法求秩
最新推荐文章于 2020-01-16 19:09:34 发布