HDU 5163 状压DP

HDU 5163
题目链接：
http://acm.hdu.edu.cn/showproblem.php?pid=5613
题意：
给一个矩阵，矩阵的值为原来矩阵此位置所在的九宫格里所有数字和（均为01）
还原出原来的矩阵。
思路：
状压枚举第一列的状态，然后递推后面状态就可以，很简单的模型题。
然而并不知道怎么合法，看到别人有用左上角的g值来推右下角，瞬间感觉智商再一次被碾压……详见代码
原来想合着一起发的，突然发现这题还没人发代码赶紧来一发赚一赚访问量[京都脸]
源码：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <string>
#include <algorithm>
#include <iostream>
using namespace std;
int g[105][105];
int out[105][105];
int output[105][105];
int n, m;
int cal(int i, int j)
{
    int ans = g[i][j];
    for(int t1 = -1 ; t1 <= 1 ; t1++){
        for(int t2 = -1 ; t2 <= 1 ; t2++){
            if(t1 == 1 && t2 == 1) continue;
            ans -= out[i + t1][j + t2];
        }
    }
    return ans;
}
bool check(int i, int j)
{
    int temp = 0;
    for(int t1 = -1 ; t1 <= 1 ; t1++){
        for(int t2 = -1 ; t2 <= 1 ; t2++){
            temp += out[i + t1][j + t2];
        }
    }
    if(temp == g[i][j]) return true;
    return false;
}
bool valid(int s)
{
    if((1 << n) & s) return false;
    for(int i = 1 ; i <= m ; i++) out[1][i] = out[n][i] = 0;
    for(int i = 2 ; i <= n - 1 ; i++){
        if((1 << i) & s) out[i][1] = 1;
        else out[i][1] = 0;
    }
// printf("s = %d\n");
// for(int i = 0 ; i < n ; i++) printf("%d ", out[i][0]);
// printf("\n");
    for(int j = 2 ; j <= m ; j++){
        for(int i = 2 ; i <= n - 1 ; i++){
            out[i][j] = cal(i - 1, j - 1);
            if(out[i][j] > 1 || out[i][j] < 0) return false;
        }
        for(int i = 1 ; i <= n ; i++){
            if(!check(i, j - 1)) return false;
        }
    }
    for(int i = 1 ; i <= n ; i++) if(!check(i, m)) return false;
    for(int i = 1 ; i <= n ; i++) for(int j = 1 ; j <= m ; j++) output[i][j] = out[i][j];
// for(int i)
    return true;
}
int main()
{
// freopen("HDU 5163.in", "r", stdin);
    int T;
    scanf("%d", &T);
    while(T--){
        scanf("%d%d", &n, &m);
        memset(g, 0, sizeof(g));
        memset(out, 0, sizeof(out));
        for(int i = 1 ; i <= n ; i++) for(int j = 1 ; j <= m ; j++) scanf("%d", &g[i][j]);
        int ok = 0;
        for(int s = 0 ; s < (1 << (n + 1)); s += 4){
            if(valid(s)){
// printf("s = %d\n", s);
                ok++;
            }
        }
        if(ok == 0) printf("Impossible\n");
        else if(ok >= 2) printf("Multiple\n");
        else{
            for(int i = 1 ; i <= n ; i++){
                for(int j = 1 ; j <= m ; j++){
                    printf("%d", output[i][j]);
                    if(j == m) printf("\n");
                    else printf(" ");
                }
            }
        }
    }
    return 0;
}