P1687 - 【NOIP2009】靶形数独
Description
小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 Z博士请教，Z 博士拿出了他最近发明的“靶形数独” ，作为这两个孩子比试的题目。
靶形数独的方格同普通数独一样，在 9 格宽×9 格高的大九宫格中有 9 个 3 格宽×3 格高的小九宫格（用粗黑色线隔开的）
。在这个大九宫格中，有一些数字是已知的，根据这些数字，利用逻辑推理，在其他的空格上填入 1到 9
的数字。每个数字在每个小九宫格内不能重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即每一个方格都有一个分值，而且
如同一个靶子一样，离中心越近则分值越高。 （如图）
上图具体的分值分布是：最里面一格（黄色区域）为 10 分，黄色区域外面的一圈（红色区域）每个格子为 9
分，再外面一圈（蓝色区域）每个格子为 8分，蓝色区域外面一圈（棕色区域）每个格子为 7分，最外面一圈（白色区域）每个格子为 6
分，如上图所示。比赛的要求是：每个人必须完成一个给定的数独（每个给定数独可能有不同的填法）
，而且要争取更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和。如图，在以下的这个已经填完数字的靶形数独
游戏中，总分数为 2829。游戏规定，将以总分数的高低决出胜负。
由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能够得到的最高分数。
Input
一共 9 行。每行 9 个整数（每个数都在 0—9 的范围内） ，表示一个尚未填满的数独方格，未填的空格用“0”表示。每两个数字之间用一个空格隔开。
Output
输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。
Sample Input
样例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
样例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
Sample Output
样例1：
2829
样例2：
2852
我们先找到一个可填的数的数量最小的那个填数,这样更容易得出答案.
所以每次都找到可填的数的数量最小的那个地方搜索,这样可以减少很多状态量.
1 #include<set>
2 #include<map>
3 #include<queue>
4 #include<stack>
5 #include<ctime>
6 #include<cmath>
7 #include<string>
8 #include<vector>
9 #include<cstdio>
10 #include<cstdlib>
11 #include<cstring>
12 #include<iostream>
13 #include<algorithm>
14 #define LL long long
15 #define RG register
16 using namespace std;
17 int mp[10][10],vis1[10][10],vis2[10][10],vis3[10][10];
18 int w[10][10]={{0},{0,6,6,6,6,6,6,6,6,6},{0,6,7,7,7,7,7,7,7,6},{0,6,7,8,8,8,8,8,7,6},{0,6,7,8,9,9,9,8,7,6},{0,6,7,8,9,10,9,8,7,6},{0,6,7,8,9,9,9,8,7,6},{0,6,7,8,8,8,8,8,7,6},{0,6,7,7,7,7,7,7,7,6},{0,6,6,6,6,6,6,6,6,6}};
19 int kt[10][10]={{0},{0,1,1,1,2,2,2,3,3,3},{0,1,1,1,2,2,2,3,3,3},{0,1,1,1,2,2,2,3,3,3},{0,4,4,4,5,5,5,6,6,6},{0,4,4,4,5,5,5,6,6,6},{0,4,4,4,5,5,5,6,6,6},{0,7,7,7,8,8,8,9,9,9},{0,7,7,7,8,8,8,9,9,9},{0,7,7,7,8,8,8,9,9,9}};
20 LL ans=-1;int sum=0;
21 struct data{
22 int x,y,s;
23 }po[100];
24 bool bj[100];
25 inline void answer()
26 {
27 LL kd=0;
28 for(RG int i=1;i<=9;i++)
29 for(RG int j=1;j<=9;j++)
30 kd+=w[i][j]*mp[i][j];
31 ans=kd>ans?kd:ans;
32 }
33 inline void DFS(int now)
34 {
35 if(now==sum+1){answer();return;}
36 int zd=100,kp;
37 for(int i=1;i<=sum;i++)
38 {
39 if(!bj[i] && po[i].s<zd) zd=po[i].s,kp=i;
40 }
41 bj[kp]=1;
42 int x=po[kp].x,y=po[kp].y;
43 for(RG int k=1;k<=9;k++)
44 if(!vis1[x][k] && !vis2[y][k] && !vis3[kt[x][y]][k]) {
45 mp[x][y]=k;
46 vis1[x][k]=1,vis2[y][k]=1,vis3[kt[x][y]][k]=1;
47 for(RG int i=1;i<=sum;i++)
48 {
49 if(bj[i]) continue;
50 int xx=po[i].x,yy=po[i].y,sumsum=0;
51 for(RG int ki=1;ki<=9;ki++)
52 if(!vis1[xx][ki] && !vis2[yy][ki] && !vis3[kt[xx][yy]][ki]) sumsum++;
53 po[i].s=sumsum;
54 }
55 DFS(now+1);
56 vis1[x][k]=0,vis2[y][k]=0,vis3[kt[x][y]][k]=0;
57 mp[x][y]=0;
58 for(RG int i=1;i<=sum;i++)
59 {
60 if(bj[i]) continue;
61 int xx=po[i].x,yy=po[i].y,sumsum=0;
62 for(RG int ki=1;ki<=9;ki++)
63 if(!vis1[xx][ki] && !vis2[yy][ki] && !vis3[kt[xx][yy]][ki]) sumsum++;
64 po[i].s=sumsum;
65 }
66 }
67 bj[kp]=0;
68 }
69 int main()
70 {
71 freopen("!.in","r",stdin);
72 freopen("!.out","w",stdout);
73 for(int i=1;i<=9;i++)
74 for(int j=1;j<=9;j++)
75 scanf("%d",&mp[i][j]);
76 for(int i=1;i<=9;i++)
77 for(int j=1;j<=9;j++)
78 if(mp[i][j]) vis1[i][mp[i][j]]=1,vis2[j][mp[i][j]]=1,vis3[kt[i][j]][mp[i][j]]=1;
79 for(RG int i=1;i<=9;i++)
80 for(RG int j=1;j<=9;j++)
81 {
82 if(mp[i][j]) continue;
83 int sumsum=0;
84 for(RG int k=1;k<=9;k++)
85 if(!vis1[i][k] && !vis2[j][k] && !vis3[kt[i][j]][k]) sumsum++;
86 po[++sum].x=i,po[sum].y=j,po[sum].s=sumsum;
87 }
88 DFS(1);
89 printf("%lld",ans);
90 return 0;
91 }