码迷,mamicode.com
首页 > 其他好文 > 详细

POJ_1222_高斯消元

时间:2016-09-03 14:52:13      阅读:115      评论:0      收藏:0      [点我收藏+]

标签:

题目描述:

  每组数据给出一个5*6的0 1矩阵,每次操作可以把某个位置及其四周的位置0 1置换,要求输出操作位置的矩阵。

思路:

  每个位置操作2次则等于没有操作,所以每个位置有操作和不操作两种选择,爆搜应该会超时。

  在网上看到了高斯消元的做法,按照每个操作位置影响的位置构造系数矩阵,然后读入题目的数据构成增广矩阵。

  接下来的做法便和高斯消元一样,只是把原来的-变成了^。

  30条方程,30个未知量,所以最终的解也是唯一。

 

#include<cstdio>
#include<iostream>
using namespace std;

int a[31][31] = {0};

int main()
{
    int n;
    cin >> n;
    for(int x = 1;x <= n;x++)
    {
        for(int i = 0;i < 30;i++)
        {
            a[i][i] = 1;
            if(i > 5)   a[i][i-6] = 1;
            if(i < 24)  a[i][i+6] = 1;
            if(i%6)     a[i][i-1] = 1;
            if((i+1)%6) a[i][i+1] = 1;
        }
        for(int i = 0;i < 30;i++)   cin >> a[i][30];
        for(int i = 0;i < 30;i++)
        {
            int temp = i;
            for(;temp < 30;temp++)
            {
                if(a[temp][i])
                {
                    for(int j = 0;j <= 30;j++)  swap(a[temp][j],a[i][j]);
                    break;
                }
            }
            for(int j = 0;j < 30;j++)
            {
                if(j != i && a[j][i])
                {
                    for(int k = i;k <= 30;k++)
                    {
                        a[j][k] ^= a[i][k];
                    }
                }
            }
        }
        cout << "PUZZLE #" << x << endl;
        for(int i = 0;i < 30;i++)
        {
            if((i+1)%6) cout << a[i][30] << " ";
            else        cout << a[i][30] << endl;
        }
    }
    return 0;
}

 

POJ_1222_高斯消元

标签:

原文地址:http://www.cnblogs.com/zhurb/p/5836868.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!