poj3279

本文介绍了一个通过矩阵操作来解决特定问题的算法实现过程。该算法的目标是最小化矩阵操作步数,以达到期望的状态。文章详细展示了如何使用位运算、数组复制等手段进行矩阵的转换,并最终找到最优解。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include <iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 20;
const int inf = 500;
int a[maxn][maxn]={0};
int ans[maxn][maxn];
int step;
int main()
{
    int m,n;
    cin>>m>>n;
    step=inf;
    for(int i=0;i<m;i++)
        for(int j=0;j<n;j++) cin>>a[i][j];
    for(int i=0;i<(1<<n);i++)
    {
        bool y=true;
        int b[maxn][maxn];
        int c[maxn][maxn]={0};
        memcpy(b,a,sizeof b);
        for(int j=0;j<n;j++)
            c[0][n-1-j] = (i>>j)&1;
        for(int k=1;k<=m;k++)
        {
            b[k-1][0]=b[k-1][0]^c[k-1][0];
            if(n>1)
            {
                b[k-1][0]=b[k-1][0]^c[k-1][1];
                for(int j=1;j<n-1;j++) b[k-1][j]=b[k-1][j]^c[k-1][j]^c[k-1][j-1]^c[k-1][j+1];
                b[k-1][n-1]=b[k-1][n-1]^c[k-1][n-1]^c[k-1][n-2];
            }
            for(int j=0;j<n;j++) b[k][j]=b[k][j]^c[k-1][j];
            for(int j=0;j<n;j++) c[k][j]=b[k-1][j];
        }
        for(int j=0;j<n;j++)
            if(b[m-1][j]!=0)
            {
                y=false;
                break;
            }
        if(y)
        {
            int tmp=0;
            for(int k=0;k<m;k++)
                for(int j=0;j<n;j++) tmp+=c[k][j];
            if(tmp<step)
            {
                step=tmp;
                memcpy(ans,c,sizeof ans);
            }
        }
    }
    if(step!=inf)
    {
        for(int i=0;i<m;i++)
            {
                for(int j=0;j<n-1;j++) cout<<ans[i][j]<<' ';
                cout<<ans[i][n-1]<<endl;
            }
    }
    else cout<<"IMPOSSIBLE\n";
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值