POJ 2965 The Pilots Brothers' refrigerator（高斯消元三种方法）

最新推荐文章于 2019-10-10 20:27:13 发布
最新推荐文章于 2019-10-10 20:27:13 发布
阅读量985 收藏 2
点赞数 1
分类专栏： ACM_数论 ACM_矩阵+高斯 ACM_PKU ITAK的ACM之路 文章标签： 高斯消元
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qingshui23/article/details/52490608
版权
ITAK的ACM之路 同时被 3 个专栏收录
295 篇文章 26 订阅
订阅专栏
ACM_数论
97 篇文章 0 订阅
订阅专栏
ACM_PKU
47 篇文章 0 订阅
订阅专栏

传送门

The game “The Pilots Brothers: following the stripy elephant” has a quest where a player needs to open a refrigerator.

There are 16 handles on the refrigerator door. Every handle can be in one of two states: open or closed. The refrigerator is open only when all handles are open. The handles are represented as a matrix 4х4. You can change the state of a handle in any location [i, j] (1 ≤ i, j ≤ 4). However, this also changes states of all handles in row i and all handles in column j.

The task is to determine the minimum number of handle switching necessary to open the refrigerator.

Input

The input contains four lines. Each of the four lines contains four characters describing the initial state of appropriate handles. A symbol “+” means that the handle is in closed state, whereas the symbol “−” means “open”. At least one of the handles is initially closed.

Output

The first line of the input contains N – the minimum number of switching. The rest N lines describe switching sequence. Each of the lines contains a row number and a column number of the matrix separated by one or more spaces. If there are several solutions, you may give any one of them.

Sample Input

-+--
----
----
-+--

Sample Output

6
1 1
1 3
1 4
4 1
4 3
4 4

Source

题目大意：

给你一个 4∗4 的矩阵，每个格子只有两个状态 “+” 和 "−" ，现在要将初始状态转化为全部是 "−" 的，当翻转一个格子的时候，每个格子

对应的行和列都得进行翻转，输出最小步数以及如何翻转的步骤。

解题思路：

第一种方法：还是高斯消元，首先找出系数矩阵，然后进行高斯消元就好啦
My Code： 

/**
2016 - 09 - 09  晚上
Author: ITAK

Motto:

今日的我要超越昨日的我，明日的我要胜过今日的我，
以创作出更好的代码为目标，不断地超越自己。
**/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <queue>
#include <algorithm>
#include <set>
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
const int INF = 1e9+5;
const int MAXN = 400;
const int MOD = 1e9+7;
const double eps = 1e-7;
const double PI = acos(-1);
using namespace std;
int equ, var;///equ个方程 var个变量
int a[MAXN][MAXN];///增广矩阵
int x[MAXN];///解集
int x_i[MAXN];
bool free_x[MAXN];///判断是不是自由变元
int free_num;///自由变元的个数
int Gauss()
{
    int Max_r;///当前列绝对值最大的存在的行
    ///col：处理当前的列
    int row,col = 0;
    int free_x_num;
    int free_index;
    free_num = 0;
    for(int i=0; i<=var; i++)
    {
        x[i] = 0;
        free_x[i] = 1;
    }
    for(row=0; row<equ&&col<var; row++,col++)
    {
        Max_r = row;
        for(int i=row+1; i<equ; i++)
            if(abs(a[i][col]) > abs(a[Max_r][col]))
                Max_r = i;
        if(a[Max_r][col] == 0)
        {
            free_x[col] = 1;
            x_i[free_num++] = col;
            row--;
            continue;
        }
        if(Max_r != row)
            for(int i=col; i<var+1; i++)
                swap(a[row][i], a[Max_r][i]);
        ///消元
        for(int i=row+1; i<equ; i++)
            if(a[i][col])
                for(int j=col; j<var+1; j++)
                    a[i][j] ^= a[row][j];
    }
    for(int i=row; i<equ; i++)
        if(a[i][col])
            return -1;///无解
    ///保证对角线主元非 0
    for(int i=0;  i<equ; i++)
    {
        if(!a[i][i])
        {
            int j;
            for(j=i+1; j<var; j++)
                if(a[i][j])
                    break;
            if(j == var)
                break;
            for(int k=0; k<equ; k++)
                swap(a[k][i], a[k][j]);
        }
    }
    if(row < var)
        return var - row;///自由变元的个数
    ///回代，得到解集
    for(int i=var-1; i>=0; i--)
    {
        x[i] = a[i][var];
        for(int j=i+1; j<var; j++)
            x[i] ^= (a[i][j] && x[j]);
    }
    return 0;///唯一解
}
void Debug()
{
    puts("");
    cout<<"+++++++++++++++++++++++++++分界线++++++++++++++++++++++++++++++"<<endl;
    for(int i=0; i<equ; i++)
    {
        for(int j=0; j<var+1; j++)
        {
            cout<<a[i][j]<<" ";
        }
        cout<<endl;
    }
    cout<<"+++++++++++++++++++++++++++分界线++++++++++++++++++++++++++++++"<<endl;
    puts("");
}
void Init()
{
    memset(a, 0, sizeof(a));
    memset(x, 0, sizeof(x));
    for(int i=0; i<equ; i++)
    {
        a[i][i] = 1;
        int tmp = i;
        while(tmp-4 >= 0)
        {
            a[tmp-4][i] = 1;
            tmp -= 4;
        }
        tmp = i;
        while(tmp+4 < 16)
        {
            a[tmp+4][i] = 1;
            tmp += 4;
        }
        int ta = i % 4, tb = i / 4;
        tmp = i;
        for(int j=0; j<ta; j++)
        {
            tmp--;
            a[tmp][i] = 1;
        }
        tmp = i;
        for(int j=ta; j<3; j++)
        {
            tmp++;
            a[tmp][i] = 1;
        }
    }
}
char str[MAXN][MAXN];
int main()
{
    equ = var = 16;
    /**Init();
    Gauss();
    Debug();**/
    while(cin>>str[0])
    {
        for(int i=1; i<4; i++)
            cin>>str[i];
        Init();
        for(int i=0; i<4; i++)
            for(int j=0; j<4; j++)
                if(str[i][j] == '+')
                    a[i*4+j][var] = 1;
                else
                    a[i*4+j][var] = 0;
        int tmp = Gauss(), sum = 0;
        for(int i=0; i<var; i++)
            sum += x[i];
        cout<<sum<<endl;
        for(int i=0; i<var; i++)
            if(x[i])
                cout<<i/4+1<<" "<<i%4+1<<endl;
    }
    return 0;
}

第二种方法：

首先我们通过系数矩阵，先进行消元发现只有唯一解，所以根据求逆的那个公式，做出来了，这就是求逆的过程，也已经讲过了。

My Code： 

/**
2016 - 09 - 09  晚上
Author: ITAK

Motto:

今日的我要超越昨日的我，明日的我要胜过今日的我，
以创作出更好的代码为目标，不断地超越自己。
**/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <queue>
#include <algorithm>
#include <set>
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
const int INF = 1e9+5;
const int MAXN = 400;
const int MOD = 1e9+7;
const double eps = 1e-7;
const double PI = acos(-1);
using namespace std;
int equ, var;///equ个方程 var个变量
int a[MAXN][MAXN];///增广矩阵
int c[MAXN][MAXN];
void Gauss()
{
    int row = 0;
    for(int col=0; col<var; col++)///代表列，主元的位置
    {
        int k = col;///k 代表行，从对角线的行开始
        for(; k<equ; k++)
            if(a[k][col])///找到当前列不为 0 的数作为主元
                break;
        if(k == equ)///含有自由变量，row表示开始自由变量的行数
        {
            row = col;
            break;
        }
        for(int i=0; i<var+1; i++)///交换行
            swap(a[col][i], a[k][i]);
        for(int i=0; i<equ; i++)///开始消元,i 代表行
        {
            if(col!=i && a[i][col])///不是主元行，且要消元的行必须是1
                for(int j=0; j<var+1; j++)///j 代表列
                    a[i][j] ^= a[col][j];
        }
    }
    for(int i=0; i<16; i++)
        for(int j=0; j<16; j++)
            c[i][j] = a[i][j+16];
}
void Debug()
{
    puts("");
    cout<<"+++++++++++++++++++++++++++分界线++++++++++++++++++++++++++++++"<<endl;
    for(int i=0; i<equ; i++)
    {
        for(int j=0; j<16; j++)
        {
            cout<<c[i][j]<<" ";
        }
        cout<<endl;
    }
    cout<<"+++++++++++++++++++++++++++分界线++++++++++++++++++++++++++++++"<<endl;
    puts("");
}
void Init()
{
    memset(a, 0, sizeof(a));
    for(int i=0; i<equ; i++)
    {
        a[i][i] = 1;
        int tmp = i;
        while(tmp-4 >= 0)
        {
            a[tmp-4][i] = 1;
            tmp -= 4;
        }
        tmp = i;
        while(tmp+4 < 16)
        {
            a[tmp+4][i] = 1;
            tmp += 4;
        }
        int ta = i % 4, tb = i / 4;
        tmp = i;
        for(int j=0; j<ta; j++)
        {
            tmp--;
            a[tmp][i] = 1;
        }
        tmp = i;
        for(int j=ta; j<3; j++)
        {
            tmp++;
            a[tmp][i] = 1;
        }
    }
    for(int i=0; i<16; i++)
        a[i][i+16] = 1;
}
int x[MAXN], b[MAXN];
char str[MAXN][MAXN];
int main()
{
    equ = 16;
    var = 32;
    Init();
    Gauss();
    ///Debug();
    while(cin>>str[0])
    {
        for(int i=1; i<4; i++)
            cin>>str[i];
        ///!!!注意初始化
        memset(x, 0, sizeof(x));
        memset(b, 0, sizeof(b));
        for(int i=0; i<4; i++)
            for(int j=0; j<4; j++)
                if(str[i][j] == '+')
                    b[i*4+j] = 1;
        for(int i=0; i<16; i++)
            for(int j=0; j<16; j++)
                x[i] ^= (c[i][j]*b[j]);
        int sum = 0;
        for(int i=0; i<16; i++)
            sum += x[i];
        cout<<sum<<endl;
        for(int i=0; i<16; i++)
            if(x[i])
                cout<<i/4+1<<" "<<i%4+1<<endl;
    }
    return 0;
}

第三种方法：

直接打表，直接通过高斯消元得到 A−1 ，然后把 A−1 的矩阵保存在一个数组中，直接进行矩阵乘法就 OK 了

My Code： 

/**
2016 - 09 - 09  晚上
Author: ITAK

Motto:

今日的我要超越昨日的我，明日的我要胜过今日的我，
以创作出更好的代码为目标，不断地超越自己。
**/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <queue>
#include <algorithm>
#include <set>
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
const int INF = 1e9+5;
const int MAXN = 400;
const int MOD = 1e9+7;
const double eps = 1e-7;
const double PI = acos(-1);
using namespace std;
int x[MAXN], b[MAXN];
char str[MAXN][MAXN];
int cc[16][16]={1 ,1 ,1 ,1 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,
1 ,1 ,1 ,1 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,
1 ,1 ,1 ,1 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,
1 ,1 ,1 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,
1 ,0 ,0 ,0 ,1 ,1 ,1 ,1 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,
0 ,1 ,0 ,0 ,1 ,1 ,1 ,1 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,
0 ,0 ,1 ,0 ,1 ,1 ,1 ,1 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,
0 ,0 ,0 ,1 ,1 ,1 ,1 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,
1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,1 ,1 ,1 ,1 ,0 ,0 ,0 ,
0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,1 ,1 ,1 ,1 ,0 ,1 ,0 ,0 ,
0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,1 ,1 ,1 ,1 ,0 ,0 ,1 ,0 ,
0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,1 ,1 ,1 ,1 ,0 ,0 ,0 ,1 ,
1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,1 ,1 ,1 ,
0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,1 ,1 ,1 ,1 ,
0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,1 ,1 ,1 ,1 ,
0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,0 ,0 ,0 ,1 ,1 ,1 ,1 ,1 ,
};
int main()
{
    /**equ = 16;
    var = 32;
    Init();
    Gauss();
    Debug();*/
    while(cin>>str[0])
    {
        for(int i=1; i<4; i++)
            cin>>str[i];
        ///!!!注意初始化
        memset(x, 0, sizeof(x));
        memset(b, 0, sizeof(b));
        for(int i=0; i<4; i++)
            for(int j=0; j<4; j++)
                if(str[i][j] == '+')
                    b[i*4+j] = 1;
        for(int i=0; i<16; i++)
            for(int j=0; j<16; j++)
                x[i] ^= (cc[i][j]*b[j]);
        int sum = 0;
        for(int i=0; i<16; i++)
            sum += x[i];
        cout<<sum<<endl;
        for(int i=0; i<16; i++)
            if(x[i])
                cout<<i/4+1<<" "<<i%4+1<<endl;
    }
    return 0;
}
ITAK
关注
专栏目录
The Pilots Brothers' refrigerator&&开灯问题
54188的博客
04-08 779
同样，一个连蒙再猜翻译了半个多小时还翻译错了好多的让我自闭一天的题T^T 这个题的意思是给出四行四列的一个字符串，要将其中的“+”全部变为“-”，输出最少的变换方式。 一开始想到的是dfs，但做了半天除了一堆bug什么都没写出来，后来看题解才明白是什么意思。 ...
POJ 2965 The Pilots Brothers' refrigerator （DFS + 枚举）
LiQiyao的博客
04-18 631
The Pilots Brothers' refrigerator Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 26183   Accepted: 10089   Special Judge Description The game “The Pil
POJ2965
Kizuki
02-02 409
这道题有多种方法 1.暴力枚举 #include"stdio.h" int i,j,num=0,min=17,turns[16]={0xf888,0xf444,0xf222,0xf111,0x8f88,0x4f44,0x2f22,0x1f11,0x88f8,0x44f4, 0x22f2,0x11f1,0x888f,0x444f,0x222f,0x111f}; int temchge[16
poj 2965 The Pilots Brothers' refrigerator-1
weixin_30371469的博客
05-07 195
poj 2965 The Pilots Brothers' refrigerator-1 枚举 #include <iostream> using namespace std; void flip(int* a,int i) { int row,col; row=i/4+1; col=i...
POJ - 2965 The Pilots Brothers' refrigerator （三种方法）（状压BFS 状压DFS 思维）
QDU_ZXF的博客
08-29 4024
题目链接： http://poj.org/problem?id=2965 题目： The game “The Pilots Brothers: following the stripy elephant” has a quest where a player needs to open a refrigerator. There are 16 handles on the refriger...
【POJ - 2965】The Pilots Brothers' refrigerator（暴力枚举，思维）
xuanweiace的博客
10-10 335
题干： The game “The Pilots Brothers: following the stripy elephant” has a quest where a player needs to open a refrigerator. There are 16 handles on the refrigerator door. Every handle can be in one o...
POJ2965-The Pilots Brothers' refrigerator
07-29
在提供的代码文件中，"POJ2965-The Pilots Brothers' refrigerator(DFS+enum).cpp" 应该是使用DFS和枚举实现的解决方案，而 "POJ2965-The Pilots Brothers' refrigerator(DFS+Bit).cpp" 是使用DFS和位运算的版本。...
POJ 2965 The Pilots Brothers' refrigerator（技巧）
da_kao_la的博客
09-07 197
题目来源：http://poj.org/problem?id=2965 问题描述 The Pilots Brothers' refrigerator Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 29748   Accepted: 1...
POJ 2965 The Pilots Brothers' refrigerator【枚举+dfs】
free斩
07-22 6144
题目:http://poj.org/problem?id=2965 来源:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=26732#problem/B 题意: 就是把题目中给出的状态图,全部翻转成 ---- ---- ---- ----状态 翻转: 每次翻转一个,那么它所在的行和列都要翻转 问最小翻转次数,
POJ---3185：The Water Bowls【高斯消元】
KobeDuu
03-06 494
题意： 一排有20个碗，给出每个碗的状态，改变一个碗的状态，其左右两边的碗的状态也随之改变(0 -> 1，1 -> 0)，求全部碗的状态为0时的最少操作数 分析: （1）每个碗的状态只有0或1（所以mod 2），设第i个碗旋转了 xi 次，得到以下方程： 说明：如果转动第j个碗影响第i个碗，aij = 1，否则aij = 0；vi 为第 i 碗的初末状态的异...
POJ 2695 The Pilots Brothers' refrigerator(神奇的规律)
N诺 - 计算机考研保研必备神器 noobdream.com
09-15 1221
转载请注明出处：http://blog.csdn.net/a1dark 分析：如果想要将一个“+”翻转成“-”，那么必然会把对应的行和列上的所有点翻转一次、由于一个点翻偶数次就相当于不翻转、所以我需要统计“+”、然后将对应行和列的值+1、最后统计奇数值的个数、便是要翻转的点、 #include int mpt[5][5]; int main(){ char ch; for(
poj2965 The Pilots Brothers' refrigerator
qq_36386435的博客
03-09 255
he Pilots Brothers’ refrigerator Time Limit: 1000 ms / Memory Limit: 65536 kb Description The game “The Pilots Brothers: following the stripy elephant” has a quest where a player needs to open a ...
The Pilots Brothers' refrigerator
weixin_30408309的博客
07-16 64
#include<iostream> #include<cstdio> #include<string> using namespace std; int mem[4][4];//用于储存翻转次数，如果是偶数，则相当于没有翻转 int main(){ char ch; for (int x = 0; x<4; x++)...
POJ2965 The Pilots Brothers' refrigerator
sunacmer的专栏
03-14 636
考查：BFS提交情况：一次AC时间:313MS 内存:1672K代码：#include typedef struct _n { int n; int step; }node; typedef struct _p { int row; int col; int pre; }path; node queue[65536]; int visit[65536]; path paths[65536]; int num; int map[16]= { 0x111f,0x
POJ 2965-The Pilots Brothers' refrigerator(贪心+枚举)
Rocky0429
01-09 1372
The Pilots Brothers' refrigerator Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 19464   Accepted: 7462   Special Judge Description The game “The Pilo
北大POJ 1002题解题技巧与实现方法
资源摘要信息:"poj_1002_487.rar_poj 1002"是一道编程题目，主要考察数据结构中的映射关系，特别是哈希表的应用。该题目来源于北京大学在线评测系统POJ（Programming Online Judge），题目编号为1002。 描述中提到...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值