Luogu1074靶形数独【启发式搜索】

Luogu1074靶形数独

题目描述

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他

们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 Z 博士请教，

Z 博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。

靶形数独的方格同普通数独一样，在 9 格宽×9 格高的大九宫格中有 9 个 3 格宽×3 格

高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些数字，利用逻辑推理，在其他的空格上填入 1 到 9 的数字。每个数字在每个小九宫格内不能

重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即

每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。（如图）

上图具体的分值分布是：最里面一格（黄色区域）为 10 分，黄色区域外面的一圈（红

色区域）每个格子为 9 分，再外面一圈（蓝色区域）每个格子为 8 分，蓝色区域外面一圈（棕

色区域）每个格子为 7 分，最外面一圈（白色区域）每个格子为 6 分，如上图所示。比赛的

要求是：每个人必须完成一个给定的数独（每个给定数独可能有不同的填法），而且要争取

更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字

的乘积的总和

总分数即每个方格上的分值和完成这个数独时填在相应格上的数字

的乘积的总和。如图，在以下的这个已经填完数字的靶形数独游戏中，总分数为 2829。游戏规定，将以总分数的高低决出胜负。

由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能

够得到的最高分数。

输入输出格式

输入格式：

一共 9 行。每行 9 个整数（每个数都在 0―9 的范围内），表示一个尚未填满的数独方

格，未填的空格用“0”表示。每两个数字之间用一个空格隔开。

输出格式：

输出文件 sudoku.out 共 1 行。

输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。

输入输出样例

输入样例#1：

sudoku1
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

sudoku2
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

输出样例#1：

sudoku1
2829

sudoku2
2852

说明

【数据范围】

40%的数据，数独中非 0 数的个数不少于 30。

80%的数据，数独中非 0 数的个数不少于 26。

100%的数据，数独中非 0 数的个数不少于 24。

NOIP 2009 提高组 第四题

一道非常好的启发式搜索题，是noip2009提高组第四题，其实这题的剪枝并不难，用类似八皇后的方法，用heng[]与zong[]还有ge[]三个数组分别表示行／列／大格中一个数字是否出现过，然后ｄｆｓ。先每个空格扫过去，找剩下的可能性最小的一格，再对这一格搜索就够了。数独进行到最后可能性会越来越少，不用担心搜索量过大。

程序头上我悄悄地开了一个O2开关，可以减少大约40ms。

 #pragma GCC optimize(2)
 #include<cstdio>
 #include<iostream>
 #include<cstring>
 using namespace std;
 ][]=
 {
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,},
 {,,,,,,,,}
 };
 ;
 struct node{
     bool heng[n][n],zong[n][n],ge[n][n];
     int matrix[n][n],filled;
 }tmp;
 ;
 inline int judge(int x,int y)
 {
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
     &&y<=) ;
 }
 inline void dfs()
 {
     )
     {
         ans=max(ans,value);
         return;
     }
     ];
     ,unknown[],point,size=;
     ;i<n;i++)
         ;j<n;j++) )
         {
             size=;
             ;k<n;k++)
             {
                 t[k]=(tmp.heng[i][k]||tmp.zong[j][k]||tmp.ge[judge(i,j)][k]);
                 ) size++;
             }
             ) return;
             if(size<flag)
             {
                 flag=size;
                 p=i;
                 q=j;
                 point=;
                 ;k<n;k++) ) unknown[++point]=k+;
             }
         }
     tmp.filled++;
     ;i<=point;i++)
     {
         tmp.matrix[p][q]=unknown[i];
         tmp.heng[p][unknown[i]-]=;
         tmp.zong[q][unknown[i]-]=;
         tmp.ge[judge(p,q)][unknown[i]-]=;
         value+=unknown[i]*a[p][q];
         dfs();
         value-=unknown[i]*a[p][q];
         tmp.heng[p][unknown[i]-]=;
         tmp.zong[q][unknown[i]-]=;
         tmp.ge[judge(p,q)][unknown[i]-]=;
     }
     tmp.matrix[p][q]=;
     tmp.filled--;
 }
 int main()
 {
     ;i<n;i++)
         ;j<n;j++)
         {
             scanf("%d",&tmp.matrix[i][j]);
             )
             {
                 tmp.heng[i][tmp.matrix[i][j]-]=;
                 tmp.zong[j][tmp.matrix[i][j]-]=;
                 tmp.ge[judge(i,j)][tmp.matrix[i][j]-]=;
                 value+=tmp.matrix[i][j]*a[i][j];
                 tmp.filled++;
             }
         }
     dfs();
     ) printf("-1");
     else printf("%d",ans);
     ;
 }