Fliptile (DFS)

题目链接: Fliptile

大致题意:

给定一个01矩阵, 大小为N*M, 每次可以选定一个格子, 将该格子与其上下左右四个格子进行反转, (0变1, 1变0), 问至少操作多少次可以将矩阵中的所有数值变为0. 若有解则输出翻转次数最少的翻转情况(若翻转次数最少的情况有多种, 则输出翻转字典序最小的.), 否则输出 IMPOSSIBLE.

解题思路:

我们很容易发现, 如果对1个格子进行偶数次翻转, 那将是没有意义的, 所以每个格子的情况只有翻转和不翻转两种情况.
一共N*M个格子, 每一个格子都有两种情况, 这样我们发现复杂度最坏为 2^nm.

考虑到优化, 其实暴力的方法有许多的情况是不可达的, 一个格子的数值改变只能由自己及相邻4个格子来进行改变. 如果我们以行的角度去看, 若第i行已经确定了其翻转情况, 如果在该行的某个位置为1, 则只能通过i+1行来进行改变.

不妨枚举出第一行所有翻转情况(最多2^m), 然后分别判断每一个点是否为0, 若不为0, 则翻转该位置下方的格子, 直至第n-1行. 完成所有翻转后, 检查第n行是否符合要求, 若符合要求则为一种可行解. 由于要求翻转次数少, 且字典序小, 所以我们应当从第一行都不翻转到第一行都翻转(从右向左)来继续(遵循字典序小->大), 然后遍历2-n-1行, 检查每一个[i, j]位置是否需要翻转. 最后检查第n行是否符合要求即可.

在上述解决方案中, 如果我们想知道某个位置[i, j]是否需要翻转, 则需要检查[i - 1, j], 为1则翻转. 而某个位置是否为1我们可以通过如下判断: 如果该点初始为1, 我们不妨也认为是经历了一次翻转得到的, 然后检测当前位置及当前位置左右上三个方向即可. 若总翻转次数为偶数, 则表明为0, 反之为1. 而对于最后一行而言, 不需要检测当前位置是否翻转(但是检测了也没问题, 这样可以省函数.)

AC代码:

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int N = 20;
bool g[N][N], vis[N][N], res[N][N]; //g存图, vis表示当前位置是否需要翻转, res存答案翻转情况
int fx[][2] = { 0, 0, 0, 1, -1, 0, 0, -1 };
int n, m;
bool judge(int x, int y) { //判断当前位置数值
    int cou = g[x][y];
    for (int i = 0; i < 4; ++i) {
        int dx = x + fx[i][0], dy = y + fx[i][1];
        if (dx < 1 || dy < 1 || dx > n || dy > m) continue;
        cou += vis[dx][dy];
    }
    return cou % 2;
}
int fact() {
    int cou = 0;
    for (int i = 2; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            if (judge(i - 1, j)) vis[i][j] = 1, cou++;
        }
    }

    for (int j = 1; j <= m; ++j)
        if (judge(n, j)) return 0x3f3f3f3f;

    return cou;
}
int main()
{
    cin >> n >> m;
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) scanf("%d", &g[i][j]);
    }

    int fmin = 0x3f3f3f3f; //初始不可达
    for (int i = 0; i <= 1 << m; ++i) {
        int num = 0; memset(vis, 0, sizeof vis); //初始化当前情况
        for (int j = 1; j <= m; ++j) vis[1][m - j + 1] = i >> (j - 1) & 1, num += vis[1][m - j + 1];
        num += fact();

        if (num < fmin) { //更优
            fmin = num;
            memcpy(res, vis, sizeof res); //复制操作情况
        }
    }

    if (fmin == 0x3f3f3f3f) printf("IMPOSSIBLE\n");
    else {
        for (int i = 1; i <= n; ++i) {
            for (int j = 1; j <= m; ++j) printf("%d%c", res[i][j], j == m ? '\n' : ' ');
        }
    }
    return 0;
}

END