靶形数独

本文介绍了一种针对靶形数独问题的优化算法。通过使用数组记录行列及九宫格状态来剪枝不可行解,并调整搜索顺序以减少搜索树规模。此外还探讨了更高效的位运算方法。

靶形数独

原题地址:洛谷


(听说数独是一个NP完全问题)
本题使用搜索,且需要一定的优化,从网上百度到了以下优化:


剪掉不可行解

为了提高效率,可以像八皇后一样,开三个数组,记录行列以及九宫格的状态,只有在三个数组中都未出现的值才可以被填入。
可以参考以下实现:

//定义数组,以及判断格子是否可行,设置格子已被占用及回溯的方法(使用宏定义)
#define get_block(a, b) ((a/3*3)+b/3)
#define check(x, y, i) (col[x][i]&&row[y][i]&&block[get_block(x,y)][i])
#define set(x, y, i) (col[x][i]=false,row[y][i]=false,block[get_block(x,y)][i]=false)
#define reset(x, y, i) (col[x][i]=true,row[y][i]=true,block[get_block(x,y)][i]=true)
#define get_score(x, y, i) (score[x][y]*i)

int col[10][10], row[10][10], block[10][10], score[10][10], ns[10][10];
int main() {
    memset(col, -1, sizeof(col));
    memset(row, -1, sizeof(line));
    memset(block, -1, sizeof(block));
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            scanf("%d", &ns[i][j]);
            col[i][ns[i][j]] = false;
            line[j][ns[i][j]] = false;
            block[get_block(i, j)][ns[i][j]] = false;
            score[i][j] = 10 - max(abs(i - 4), abs(j - 4));
        }
    }

    dfs(0, 0);
    //…………
}

这是使用数组的方法,更好更快的方法是使用位运算。但是笔者尚未能熟练掌握使用位运算,因此还是使用了简单的数组解法。


改变搜索顺序

但是这题由于给出的格子比较少,因此只靠剪去不可行解只能过掉部分数据,因此还有其他一些优化方法,比如改变搜索的顺序。
要加快程序运行速度,可以先搜索可拓展节点数较少的节点,这样搜索树的节点数会少一些,可以自己实验一个例子。
具体的实现方法可以参考以下代码:

inline
void get_next(int &x, int &y) {
    int cnt = 10;
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < n; ++j) {
            if (vis[i][j]) continue;
            if (ns[i][j]) {
                x = i;
                y = j;
                return;
            }
            int tmp = 0;
            for (int k = 1; k <= 9; ++k) {
                if (check(i, j, k)) ++tmp;
            }
            if (tmp < cnt) {
                cnt = tmp;
                x = i;
                y = j;
            }
        }
    }
}

vis数组用来标记哪些点已经被访问过了。虽然每一次都要访问整个数组才能确定下一步要访问的节点,但是加上以上优化后代码顺利通过,因此这是一个不错的剪枝,反映出了搜索顺序不同所带来的时间上的巨大差异。


高级专题

除了以上两种方法之外,还可以使用Dancing Link(跳舞链),具体的可以参考最后的列表,从16*16且有多组数据的数独问题中得到启发,但是笔者没打过跳舞链只能给出链接。


参考来源:

c++:# P1074 [NOIP 2009 提高组] 数独 ## 题目背景 **此为远古题,不保证存在可以通过任意符合要求的输入数据的程序**。 ## 题目描述 小城和小华都是热爱数学的好学生,最近,他们不约而同地迷上了数独游戏,好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了,于是他们向 Z 博士请教,Z 博士拿出了他最近发明的“数独”,作为这两个孩子比试的题目。 数独的方格同普通数独一样,在 $9$ 格宽且 $9$ 格高的大九宫格中有 $9$ 个 $3$ 格宽且 $3$ 格高的小九宫格(用粗黑色线隔开的)。在这个大九宫格中,有一些数字是已知的,根据这些数字,利用逻辑推理,在其他的空格上填入 $1$ 到 $9$ 的数字。每个数字在每个小九宫格内不能重复出现,每个数字在每行、每列也不能重复出现。但数独有一点和普通数独不同,即每一个方格都有一个分值,而且如同一个子一样,离中心越近则分值越高。(如图) ![](https://cdn.luogu.com.cn/upload/image_hosting/dhvuc32i.png) 上图具体的分值分布是:最里面一格(黄色区域)为 $10$ 分,黄色区域外面的一圈(红色区域)每个格子为 $9$ 分,再外面一圈(蓝色区域)每个格子为 $8$ 分,蓝色区域外面一圈(棕色区域)每个格子为 $7$ 分,最外面一圈(白色区域)每个格子为 $6$ 分,如上图所示。比赛的要求是:每个人必须完成一个给定的数独(每个给定数独可能有不同的填法),而且要争取更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和 总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和。如图,在以下的这个已经填完数字的数独游戏中,总分数为 $2829$。游戏规定,将以总分数的高低决出胜负。 ![](https://cdn.luogu.com.cn/upload/image_hosting/yx82mmnc.png) 由于求胜心切,小城找到了善于编程的你,让你帮他求出,对于给定的数独,能够得到的最高分数。 ## 输入格式 一共 $9$ 行。每行 $9$ 个整数(每个数都在 $0 \sim 9$ 的范围内),表示一个尚未填满的数独方格,未填的空格用“$0$”表示。每两个数字之间用一个空格隔开。 ## 输出格式 输出共 $1$ 行。输出可以得到的数独的最高分数。如果这个数独无解,则输出整数 $-1$。 ## 输入输出样例 #1 ### 输入 #1 ``` 7 0 0 9 0 0 0 0 1 1 0 0 0 0 5 9 0 0 0 0 0 2 0 0 0 8 0 0 0 5 0 2 0 0 0 3 0 0 0 0 0 0 6 4 8 4 1 3 0 0 0 0 0 0 0 0 7 0 0 2 0 9 0 2 0 1 0 6 0 8 0 4 0 8 0 5 0 4 0 1 2 ``` ### 输出 #1 ``` 2829 ``` ## 输入输出样例 #2 ### 输入 #2 ``` 0 0 0 7 0 2 4 5 3 9 0 0 0 0 8 0 0 0 7 4 0 0 0 5 0 1 0 1 9 5 0 8 0 0 0 0 0 7 0 0 0 0 0 2 5 0 3 0 5 7 9 1 0 8 0 0 0 6 0 1 0 0 0 0 6 0 9 0 0 0 0 1 0 0 0 0 0 0 0 0 6 ``` ### 输出 #2 ``` 2852 ``` ## 说明/提示 #### 数据规模与约定 - 对于 $40\%$ 的数据,数独中非 $0$ 数的个数不少于 $30$; - 对于 $80\%$ 的数据,数独中非 $0$ 数的个数不少于 $26$; - 对于 $100\%$ 的数据,数独中非 $0$ 数的个数不少于 $24$。 NOIP 2009 提高组 第三题
最新发布
07-16
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值