UVA 11464 - Even Parity 状态压缩，分析 难度: 2

题目

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2459

题意

N*N 的01方阵，可用操作为把任意0变为1，求操作的最小次数，使得任意位置的上下左右之和（不包含自身）为偶数

思路

如刘书，关键在于状态只有第一行的2^15个。

感想

1. 忘了memset

代码

#include <algorithm>
#include <cassert>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <tuple>
#include <cassert>

using namespace std;
#define LOCAL_DEBUG
#define NOW_QUEUE (ques[queId])
#define LAST_QUEUE (ques[1 - queId])
const int MAXN = ;
struct Node{
    int lineSta;
    int lineSum;
    int changed;
    Node(int _lineSta, int _lineSum, int _changed) {
        lineSta = _lineSta;
        lineSum = _lineSum;
        changed = _changed;
    }
};
queue<Node> ques[];
int queId;
int orgStatus[MAXN][MAXN];
int stackedStatus[MAXN];
int statusLimit;

int getLineSum(int formerlineSta, int lineSta) {
    return (formerlineSta ^ (lineSta >> ) ^ (lineSta << )) & statusLimit;
}

int getChanged(int lineno, int lineSta) {
    int changedSta = stackedStatus[lineno] ^ lineSta;
    int changed = ;
    while (changedSta > ) {
        changed++;
        changedSta -= (changedSta & (-changedSta));
    }
    return changed;
}

int main() {
#ifdef LOCAL_DEBUG
    freopen("input.txt", "r", stdin);
    //freopen("output2.txt", "w", stdout);
#endif // LOCAL_DEBUG
    int T;
    scanf("%d", &T);
    for (int ti = ; ti <= T; ti++) {
        queId = ;

        int n;
        scanf("%d", &n);
        for (int i = ; i < n; i++) {
            stackedStatus[i] = ;
            for (int j = ; j < n; j++) {
                scanf("%d", orgStatus[i] + j);
                stackedStatus[i] = stackedStatus[i] *  + orgStatus[i][j];
            }
        }

        statusLimit = ( << n) - ;
        for (int sta = ; sta <= statusLimit; sta++) {
            if((sta & stackedStatus[]) == stackedStatus[])NOW_QUEUE.push(Node(sta, getLineSum(, sta), getChanged(, sta)));
        }
        for (int i = ; i < n; i++) {
            queId =  - queId;
            while (!LAST_QUEUE.empty()) {
                Node node = LAST_QUEUE.front(); LAST_QUEUE.pop();
                Node newNode(node);
                newNode.lineSta = stackedStatus[i];
                int lastLineSum = node.lineSum ^ stackedStatus[i];
                bool fl = true;
                for (int j = ; j < n && fl; j++) {
                    if (lastLineSum & ( << j)) {
                        if (stackedStatus[i] & ( << j)) { fl = false; break; }
                        else {
                            newNode.lineSta |=  << j;
                        }
                    }
                }
                if (fl) {
                    newNode.lineSum = getLineSum(node.lineSta, newNode.lineSta);
                    newNode.changed = node.changed + getChanged(i, newNode.lineSta);
                    NOW_QUEUE.push(newNode);
                }
            }
        }
        int ans = n * n + ;
        while (!NOW_QUEUE.empty()) {
            Node node = NOW_QUEUE.front(); NOW_QUEUE.pop();
            ans = min(ans, node.changed);
        }
        if (ans > n * n)printf("Case %d: -1\n", ti);
        else printf("Case %d: %d\n", ti, ans);
    }

    return ;
}