bnu 11952 高斯消元法求秩

#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 高斯消元法求秩)