洛谷—— P1074 靶形数独

https://www.luogu.org/problem/show?pid=1074

题目描述

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

们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 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 提高组 第四题

记录需要搜索的点，即未填数的点，挨个搜索

从第一个未填数点开始搜

 #include <cstdio>

 int op[][]={{,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,}};
 int pos[][]={{,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,},
               {,,,,,,,,}};
 int map[][],ans=-;
 int vx[],vy[],cnt;
 bool h[][],l[][],g[][];

 int where(int i,int j)
 {
     if(i<&&j<) return ;
     else if(i<&&j<) return ;
     else if(i<&&j<) return ;
     else if(i<&&j<) return ;
     else if(i<&&j<) return ;
     else if(i<&&j<) return ;
     else if(j<) return ;
     else if(j<) return ;
     else if(j<) return ;
 }

 void DFS(int num)
 {
     if(num>cnt)
     {
         int tmp=;
         for(int i=; i<; ++i)
           for(int j=; j<; ++j)
               tmp+=map[i][j]*op[i][j];
         ans=ans>tmp?ans:tmp;
         return ;
     }
     int tox=vx[num],toy=vy[num],topos=pos[tox][toy];
     for(int i=; i<=; ++i)
     {
         if(!h[tox][i]&&!l[toy][i]&&!g[topos][i])
         {
             h[tox][i]=;
             l[toy][i]=;
             g[topos][i]=;
             map[tox][toy]=i;
             DFS(num+);
             h[tox][i]=;
             l[toy][i]=;
             g[topos][i]=;
             map[tox][toy]=;
         }
     }
 }

 int AC()
 {
     for(int i=; i<; ++i)
         for(int j=; j<; ++j)
         {
             scanf("%d",&map[i][j]);
             if(!map[i][j])
             {
                 vx[++cnt]=i;
                 vy[cnt]=j;
             }
             else
             {
                 h[i][map[i][j]]=;
                 l[j][map[i][j]]=;
                 g[pos[i][j]][map[i][j]]=;
             }
         }
     DFS(); printf("%d\n",ans);
     return ;
 }

 int Aptal=AC();
 int main(){;}

TLE6个点，70分

改变搜索顺序，从最后一个未填数点开始搜

 #include <cstdio>

 const int op[][]={{,,,,,,,,},
                     {,,,,,,,,},
                     {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,}};
 const int pos[][]={{,,,,,,,,},
                     {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,}};
 int map[][],ans=-;
 int vx[],vy[],cnt;
 bool h[][],l[][],g[][];

 void DFS(int num)
 {
     if(!num)
     {
         int tmp=;
         for(int i=; i<; ++i)
           for(int j=; j<; ++j)
               tmp+=map[i][j]*op[i][j];
         ans=ans>tmp?ans:tmp;
         return ;
     }
     int tox=vx[num],toy=vy[num],topos=pos[tox][toy];
     for(int i=; i<=; ++i)
     {
         if(!h[tox][i]&&!l[toy][i]&&!g[topos][i])
         {
             h[tox][i]=;
             l[toy][i]=;
             g[topos][i]=;
             map[tox][toy]=i;
             DFS(num-);
             h[tox][i]=;
             l[toy][i]=;
             g[topos][i]=;
             map[tox][toy]=;
         }
     }
 }

 int AC()
 {
     for(int i=; i<; ++i)
         for(int j=; j<; ++j)
         {
             scanf("%d",&map[i][j]);
             if(!map[i][j])
             {
                 vx[++cnt]=i;
                 vy[cnt]=j;
             }
             else
             {
                 h[i][map[i][j]]=;
                 l[j][map[i][j]]=;
                 g[pos[i][j]][map[i][j]]=;
             }
         }
     DFS(cnt); printf("%d\n",ans);
     return ;
 }

 int Aptal=AC();
 int main(){;}

TLE1个点，95分

增加搜索状态，记录当前的得分，如果当前得分+9*10<ans，说明这样填数不可能会搜到最优解

判断搜索层数，快超时的时候return、、、

 #include <cstdio>

 const int op[][]={{,,,,,,,,},
                     {,,,,,,,,},
                     {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,}};
 const int pos[][]={{,,,,,,,,},
                     {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,},
                       {,,,,,,,,}};
 int vx[],vy[],cnt,deep;
 int map[][],ans=-,fraction;
 bool h[][],l[][],g[][];

 void DFS(int num,int now_fraction)
 {
     if(!num)
     {
         ans=ans>now_fraction?ans:now_fraction;
         return ;
     }
     if(++deep>1e7) return ;
     if(num**+now_fraction<ans) return ;
     int tox=vx[num],toy=vy[num],topos=pos[tox][toy],tmp;
     for(int i=; i<=; ++i)
     {
         if(!h[tox][i]&&!l[toy][i]&&!g[topos][i])
         {
             h[tox][i]=;
             l[toy][i]=;
             g[topos][i]=;
             map[tox][toy]=i;
             DFS(num-,now_fraction+map[tox][toy]*op[tox][toy]);
             h[tox][i]=;
             l[toy][i]=;
             g[topos][i]=;
             map[tox][toy]=;
         }
     }
 }

 int AC()
 {
     for(int i=; i<; ++i)
         for(int j=; j<; ++j)
         {
             scanf("%d",&map[i][j]);
             if(!map[i][j])
             {
                 vx[++cnt]=i;
                 vy[cnt]=j;
             }
             else
             {
                 h[i][map[i][j]]=;
                 l[j][map[i][j]]=;
                 g[pos[i][j]][map[i][j]]=;
                 fraction+=op[i][j]*map[i][j];
             }
         }
     DFS(cnt,fraction);
     printf("%d\n",ans);
     return ;
 }

 int Aptal=AC();
 int main(){;}

AC