poj3279

gdymind

于 2016-07-13 11:08:01 发布

阅读量342

点赞数

CC 4.0 BY-SA版权

分类专栏： acm poj dp 状态压缩文章标签： poj acm 算法

本文链接：https://blog.youkuaiyun.com/gdymind/article/details/51896226

acm 同时被 3 个专栏收录

72 篇文章

订阅专栏

31 篇文章

订阅专栏

poj

25 篇文章

订阅专栏

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

摘要生成于 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;
}