bzoj 3503 和谐矩阵

答案有空格!
答案有空格!
答案有空格!
重要的事情说三遍。
然后就是各种打错变量名

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<bitset>
#define md
#define ll long long
#define inf (int) 1e9
#define eps 1e-8
#define N 50
using namespace std;
int lx[4]={-1,0,0,0},ly[4]={0,-1,0,1};
ll ju[N][N];
bool ans[N][N];
bitset<N> a[N];
int m,n;
void gaosi()
{
for (int i=1;i<=m;i++)
{
int k=i;
while (k<=m&&!a[k][i]) k++;
if (k>m) continue;
swap(a[i],a[k]);
for (int j=i+1;j<=m;j++)
if (a[j][i]) a[j]^=a[i];
}
for (int i=m;i;i--)
{
if (a[i][i]==0) ans[1][i]=a[i][0]=1; else ans[1][i]=a[i][0];
for (int j=1;j<i;j++) if (a[j][i]) a[j][0]=a[i][0]^a[j][0];
}
}


int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=m;i++) ju[1][i]=(ll)1<<i;
for (int i=2;i<=n+1;i++)
for (int j=1;j<=m;j++)
{
ll now=0;
for (int k=0;k<=3;k++)
{
int x=i-1+lx[k],y=j+ly[k];
if (x<1||y<1||x>n||y>m) continue;
now^=ju[x][y];
}
ju[i][j]=now;
}
for (int i=1;i<=m;i++)
for (int j=1;j<=m;j++)
a[i][j]=(ju[n+1][i]>>j)&1;
gaosi();
for (int i=2;i<=n;i++)
for (int j=1;j<=m;j++)
{
ans[i][j]=0;
for (int k=0;k<=3;k++)
{
int x=i-1+lx[k],y=j+ly[k];
if (x<1||y<1||x>n||y>m) continue;
ans[i][j]^=ans[x][y];
}
}
for (int i=1;i<=n;i++)
{
printf("%d",ans[i][1]);
for (int j=2;j<=m;j++)
printf(" %d",ans[i][j]);
printf("\n");
}
return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值