NOIP提高组 2009 靶形数独 题解

Posted 2020-10-06

此文为博主原创题解，转载时请通知博主，并把原文链接放在正文醒目位置。

此题解目前在洛谷上无法AC，只能拿到95分，使用洛谷的朋友请注意。

题目链接：  洛谷   　　CODEVS 　　 Vijos

描述

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他
们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 Z博士请教，
Z 博士拿出了他最近发明的“靶形数独” ，作为这两个孩子比试的题目。
靶形数独的方格同普通数独一样，在 9 格宽×9 格高的大九宫格中有 9 个 3 格宽×3 格
高的小九宫格（用粗黑色线隔开的） 。在这个大九宫格中，有一些数字是已知的，根据这些
数字，利用逻辑推理，在其他的空格上填入 1到 9 的数字。每个数字在每个小九宫格内不能
重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即
每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。 （如图）

上图具体的分值分布是：最里面一格（黄色区域）为 10 分，黄色区域外面的一圈（红
色区域）每个格子为 9 分，再外面一圈（蓝色区域）每个格子为 8分，蓝色区域外面一圈（棕
色区域）每个格子为 7分，最外面一圈（白色区域）每个格子为 6 分，如上图所示。比赛的
要求是：每个人必须完成一个给定的数独（每个给定数独可能有不同的填法） ，而且要争取
更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字
的乘积的总和。如图，在以下的这个已经填完数字的靶形数独游戏中，总分数为 2829。游
戏规定，将以总分数的高低决出胜负。

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

格式

输入格式

一共 9 行。每行 9 个整数（每个数都在 0—9 的范围内） ，表示一个尚未填满的数独方
格，未填的空格用“0”表示。每两个数字之间用一个空格隔开。

输出格式

输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-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

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


分析：
思路很容易想到，按照正常的靶形数独来做，求出所有可能的情况并计算得分，最后输出最高的得分。
但是裸的DFS交上去之后就会疯狂TLE。
第一个优化：在搜索填入的数字时，从9往1搜。
效果：CODEVS - 2个点TLE   Vijos - 3个点TLE
第二个优化：搜索时从整张数独的右下角往左上角搜。
效果：CODEVS - AC　　　Vijos - AC　　　　洛谷：1个点TLE
然后无论怎么改都在洛谷上TLE，心态崩溃，不改了qwq搜索苦手
等我什么时候AC了再回来填这个坑qwq

其实洛谷的数据可能并不是最强的，然而洛谷时限1s，Vijos时限2s，CODEVS时限4s...于是洛谷的题变得极其恶心。

AC（TLE）代码：

 1 #include<algorithm>
 2 #include<cmath>
 3 #include<cstdio>
 4 #include<queue>
 5 
 6 inline void read(int &x)
 7 {
 8     char ch = getchar(),c = ch;x = 0;
 9     while(ch < ‘0‘ || ch > ‘9‘) c = ch,ch = getchar();
10     while(ch <= ‘9‘ && ch >= ‘0‘) x = (x<<1)+(x<<3)+ch-‘0‘,ch = getchar();
11     if(c == ‘-‘) x = -x;
12 }
13 
14 int mp[15][15],notel[15][15],notec[15][15],noteg[15][15];
15 int sum,ans = -1,base,tmp;
16 int cnt[15][15] =  {{0,0,0,0,0,0,0,0,0,0},
17                     {0,6,6,6,6,6,6,6,6,6},
18                     {0,6,7,7,7,7,7,7,7,6},
19                     {0,6,7,8,8,8,8,8,7,6},
20                     {0,6,7,8,9,9,9,8,7,6},
21                     {0,6,7,8,9,10,9,8,7,6},
22                     {0,6,7,8,9,9,9,8,7,6},
23                     {0,6,7,8,8,8,8,8,7,6},
24                     {0,6,7,7,7,7,7,7,7,6},
25                     {0,6,6,6,6,6,6,6,6,6}};
26 
27 struct SPACE
28 {
29     int x,y;
30 }q[300];
31 
32 void dfs(int x,int y,int step)
33 {
34     if(step == 0)
35     {
36         /*for(int i = 1;i <= 9;++ i)
37         {
38             for(int j = 1;j <= 9;++ j)
39                 printf("%d ",mp[i][j]);
40             puts("");
41         }*/
42         if(tmp && tmp > ans)    ans = tmp;
43         return;
44     }
45     else
46     {
47         int box = 3*((x-1)/3)+(y-1)/3+1;
48         for(register int i = 9;i >= 1;-- i)
49         {
50             if(notel[x][i] || notec[y][i] || noteg[box][i])
51                 continue;
52             else
53             {
54                 notel[x][i] = 1;
55                 notec[y][i] = 1;
56                 noteg[box][i] = 1;
57                 mp[x][y] = i;
58                 tmp += mp[x][y]*cnt[x][y];
59                 dfs(q[step-1].x,q[step-1].y,step-1);
60                 tmp -= mp[x][y]*cnt[x][y];
61                 notel[x][i] = 0;
62                 notec[y][i] = 0;
63                 noteg[box][i] = 0;
64             }
65         }
66         return;
67     }
68 }
69 
70 int main()
71 {
72     for(register int i = 1;i <= 9;++ i)
73         for(register int j = 1;j <= 9;++ j)
74         {
75             read(mp[i][j]);
76             int box = 3*((i-1)/3)+(j-1)/3+1;
77             if(mp[i][j])
78             {
79                 notel[i][mp[i][j]] = 1;
80                 notec[j][mp[i][j]] = 1;
81                 noteg[box][mp[i][j]] = 1;
82                 base += mp[i][j]*cnt[i][j];
83             }
84             else
85             {
86                 sum ++;
87                 q[sum].x = i,q[sum].y = j;
88             }
89         }
90     tmp = base;
91     dfs(q[sum].x,q[sum].y,sum);
92     printf("%d\n",ans);
93     return 0;
94 }