完美正方形-深度遍历

题目：如果一些边长互不相同的正方形，可以恰好拼出一个更大的正方形，则称其为完美正方形。

22阶完美正方形一共有8种。下面的组合是另一种： 2 5 9 11 16 17 19 21 22 24 26 30 31 33 35 36  41 46 47 50 52 61 如果告诉你该方案紧贴着上边沿的是从左到右依次为：47 46 61， 你能计算出紧贴着下边沿的是哪几个正方形吗？

思路：深度搜索，将完美正方形看成154*154的方阵，将47，46，61分别放置在方阵各位并填写数值。根据行遍历，发现第一个为0的位置，开始按照a数组顺序填充数值，若judge判断成功且该数值未被访问过，则在方阵中填充此数并且标以访问。接下来以填充后的方阵开始深度遍历，继续发现第一个不为0的位置，如上，若judge判断失败则将之前填充的数值全部改为0，并且标记为未被访问。直到全部填充完整，深度遍历结束。

#include<iostream>
#include <cstdio>
#include <cstring>
#include<algorithm>
#include<map>
#define maxn 100010
using namespace std;
typedef long long LL;
int a[19]= { 2,5,9,11,16,17,19,21,22,24,26,30,31,33,35,36,41,50,52};
int mat[154][154];
int visited[100];
void fill(int x, int y, int n, int s) {
    for (int i = x; i < x + n; i++) {
        for (int j = y; j < y + n; j++) {
            mat[i][j] = s;
        }
    }
}
int solve() {
    for (int i = 0; i < 154; i++) {
        for (int j = 0; j < 154; j++) {
            if (!mat[i][j]) {
                return 0;
            }
        }
    }
    return 1;
}
int judge(int x, int y, int n) {
    if (x + n > 154 || y + n > 154) {
        return 0;
    }
    for (int i = x; i < x + n; i++) {
        for (int j = y; j < y + n; j++) {
            if (mat[i][j]) {
                return 0;
           }
        }
    }
    return 1;
}
int dfs(int x, int y) {
    if (solve()) {
        return 1;
    }
    else {
        int flag = 1;
        for (int i = 0; i < 154 && flag; i++) {
            for (int j = 0; j < 154 && flag; j++) {
                if (!mat[i][j]) {
                    x = i, y = j, flag = 0;
                }
            }

        }
        for (int k = 0; k < 19; k++) {
            if (judge(x, y, a[k])) {
                if (!visited[k]) {
                    fill(x, y, a[k], a[k]);
                    visited[k] = 1;
                    cout << a[k] << " ";
                    if (dfs(x, y + a[k])) {
                        return 1;
                    }
                    fill(x, y, a[k], 0);
                    visited[k] = 0;
                }
            }
            else{
                break;
            }
        }
    }
    return 0;
}
int main() {
    memset(mat, 0, sizeof(mat));
    memset(visited, 0, sizeof(visited));
    fill(0, 0, 47, 47);
    fill(0, 47, 46, 46);
    fill(0, 93, 61, 61);
    dfs(0, 0);
    int ans = 0;
    for (int i = 0; i < 154; i++) {
        if (mat[153][i] != ans) {
          //  cout << mat[153][i] << " ";
            ans = mat[153][i];
        }
    }
    return 0; 
}