Noip2009 靶形数独 - 搜索 - 减少搜索树规模

当搜索数较浅的一层分支特别多时，这个搜索数规模比较大，因为搜索树的规模主要和“搜索中间节点”有关，因为搜素树底层节点是固定的，想要减少规模就应该削减搜索中的那些节点。
因此我们让浅层分支少，深层分支多，其规模一定是比原来小的

假设第一个节点有5条分支，第二个节点有2条分支
如果先搜索第一个节点，5条分支，每个分支再连到第二个节点，又有两条分支
整个搜索树有15个（5 + 10）“中间节点”

如果先搜索第二个节点，两条分支，每个分支又有5条分支
整个搜索树有12个（2 + 10）“中间节点”

因此我们把分支多的点放在后面搜索可以有效减少搜索树规模

对于这题，一个填数的位置，其可填数越少，分支就越少，先算出一个要填的点，其所在行列宫内有多少0，然后由0数少的点先开始搜索就可以了

另外，这个题好像卡常？？我在洛谷上开O2勉强过了。。。

#include <algorithm>
#include <iostream>
#include <cstdio>
using namespace std;
#define debug(x) cerr << #x << "=" << x << endl;
const int MAXN = 90;
int gra[12][12], ans, cnt, gong[10][10], line[10][10], col[10][10];
struct ppos{
    int x, y, k;
}pos[MAXN];
int calc() {
    int sum = 0;
    for(int i=1; i<=9; i++) {
        for(int j=1; j<=9; j++) {
            if(i == 1 || j == 1 || i == 9 || j == 9) {
                sum += gra[i][j] * 6;
            } else if(i == 2 || j == 2 || i == 8 || j == 8) {
                sum += gra[i][j] * 7;
            } else if(i == 3 || j == 3 || i == 7 || j == 7) {
                sum += gra[i][j] * 8;
            } else if(i == 4 || j == 4 || i == 6 || j == 6) {
                sum += gra[i][j] * 9;
            } else if(i == 5 && j == 5) {
                sum += gra[i][j] * 10;
            }
        }
    }
    return sum;
}
int find(int x, int y) {
    int num;
    if(x <= 3) {
        if(y <= 3) num = 1;
        else if(y <= 6) num = 2;
        else if(y <= 9) num = 3;
    } else if(x <= 6) {
        if(y <= 3) num = 4;
        else if(y <= 6) num = 5;
        else if(y <= 9) num = 6;
    } else if(x <= 9) {
        if(y <= 3) num = 7;
        else if(y <= 6) num = 8;
        else if(y <= 9) num = 9;
    }
    return num;
}
bool judge(int x, int y, int k, int num) {
    if(gong[num][k] || line[x][k] || col[y][k])
        return false;
    return true;
}
void dfs(int n) {
    if(n == cnt + 1) {
        ans = max(ans, calc());
        return;
    }
    int x = pos[n].x;
    int y = pos[n].y;
    int num = find(x, y);
    for(int i=1; i<=9; i++) {
        if(judge(x, y, i, num)) {
            gong[num][i] = line[x][i] = col[y][i] = 1;
            gra[x][y] = i;
            dfs(n + 1);
            gra[x][y] = 0;
            gong[num][i] = line[x][i] = col[y][i] = 0;
        }
    }
}
void init() {
    for(int i=1; i<=cnt; i++) {
        int sum = 0;
        int stx, sty, edx, edy;
        int x = pos[i].x, y = pos[i].y;
        if(x<=3) stx = 1, edx = 3;
        else if (x<=6) stx = 4, edx = 6;
        else if (x<=9) stx = 7, edx = 9;
        if(y<=3) sty = 1, edy = 3;
        else if(y<=6) sty = 4, edy = 6;
        else if(y<=9) sty = 7, edy = 9;
        for(int i=stx; i<=edx; i++) {
            for(int j=sty; j<=edy; j++) {
                if(!gra[i][j]) {
                    sum++;
                }
            }
        }
        for(int i=1; i<=9; i++) {
            if(!gra[x][i]) sum++;
            if(!gra[i][y]) sum++;
        }
        pos[i].k = sum;
    }
    for(int i=1; i<=9; i++) {
        for(int j=1; j<=9; j++) {
            if(gra[i][j]) {
                int num = find(i, j);    
                gong[num][gra[i][j]] = 1;
                line[i][gra[i][j]] = 1;
                col[j][gra[i][j]] = 1;
            }
        }
    }
}
bool cmp(ppos a, ppos b) {
    return a.k < b.k;
}
int main() {
    for(int i=1; i<=9; i++) {
        for(int j=1; j<=9; j++) {
            cin >> gra[i][j];
            if(!gra[i][j]) {
                pos[++cnt].x = i, pos[cnt].y = j;
            }
        }
    }
    init();
    sort(pos+1,pos+cnt+1,cmp);
    dfs(1);
    ans = ans == 0 ? -1 : ans;
    cout << ans << endl;
    return 0;
}