【NOIP 2009】靶形数独

题目描述

小城和小华都是热爱数学的好学生,最近,他们不约而同地迷上了数独游戏,好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了,于是他们向 Z 博士请教,Z 博士拿出了他最近发明的“靶形数独”,作为这两个孩子比试的题目。

靶形数独的方格同普通数独一样,在 9 格宽× 9 格高的大九宫格中有 个 3 格宽× 3 格高的小九宫格(用粗黑色线隔开的)。在这个大九宫格中,有一些数字是已知的,根据这些数字,利用逻辑推理,在其他的空格上填入 1 到 9 的数字。每个数字在每个小九宫格内不能重复出现,每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同,即每一个方格都有一个分值,而且如同一个靶子一样,离中心越近则分值越高。(如图)

上图具体的分值分布是:最里面一格(黄色区域)为 10 分,黄色区域外面的一圈(红色区域)每个格子为 9 分,再外面一圈(蓝色区域)每个格子为 8 分,蓝色区域外面一圈(棕色区域)每个格子为 7 分,最外面一圈(白色区域)每个格子为 6 分,如上图所示。比赛的要求是:每个人必须完成一个给定的数独(每个给定数独可能有不同的填法),而且要争取更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和

总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和。如图,在以下的这个已经填完数字的靶形数独游戏中,总分数为 2829。游戏规定,将以总分数的高低决出胜负。

   

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

输入输出格式

输入格式:

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

输出格式:

输出共 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: 
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

题解

状态的存储

数独类问题比较难的是状态的存储,所以我们需要写一些数组去存储状态。

bool a[10][10]; //a[i][j]判断i行有没有j这个数;
bool b[10][10]; //b[i][j]判断i列有没有j这个数 ;
bool c[10][10]; //c[i][j]判断第i个九宫格是否有j这个数;

本题特殊性

对于本题特殊带有的权值,我们可以用一种极其简单的方法去解决  打表!!

int score[10][10]{                                  //权值
    {0,0,0,0,0,0,0,0,0,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},
};
int gong[10][10]{                                  //所在宫的位置
    {0,0,0,0,0,0,0,0,0,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},
};

这样简化完程序的复杂性,下面我们要做的就是优化时间;

优化时间

想一下我们平时是如何玩数独的(虽然我没有玩过)

我们一定会先去找数已经填很多的地方去填,因为那个位置可以填的数少,选择的少;

在搜索中,这样做可以减少搜索树的分支,从而提高效率;

所以我们可以记录每一行 0 的个数,然后 sort 再存储搜索 0 的次序;

最后DFS 更新 ans;即可

code

下面是代码实现

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 using namespace std;
 6 bool a[10][10]; //a[i][j]判断i行有没有j这个数
 7 bool b[10][10]; //b[i][j]判断i列有没有j这个数 
 8 bool c[10][10]; //c[i][j]判断第i个九宫格是否有j这个数
 9 
10 int score[10][10]{
11     {0,0,0,0,0,0,0,0,0,0},
12     {0,6,6,6,6,6,6,6,6,6},
13     {0,6,7,7,7,7,7,7,7,6},
14     {0,6,7,8,8,8,8,8,7,6},
15     {0,6,7,8,9,9,9,8,7,6},
16     {0,6,7,8,9,10,9,8,7,6},
17     {0,6,7,8,9,9,9,8,7,6},
18     {0,6,7,8,8,8,8,8,7,6},
19     {0,6,7,7,7,7,7,7,7,6},
20     {0,6,6,6,6,6,6,6,6,6},
21 };
22 
23 int gong[10][10]{               
24     {0,0,0,0,0,0,0,0,0,0},
25     {0,1,1,1,2,2,2,3,3,3},
26     {0,1,1,1,2,2,2,3,3,3},
27     {0,1,1,1,2,2,2,3,3,3},
28     {0,4,4,4,5,5,5,6,6,6},
29     {0,4,4,4,5,5,5,6,6,6},
30     {0,4,4,4,5,5,5,6,6,6},
31     {0,7,7,7,8,8,8,9,9,9},
32     {0,7,7,7,8,8,8,9,9,9},
33     {0,7,7,7,8,8,8,9,9,9},
34 };
35 int map[10][10],s[100][5]; 
36 
37 int ans,cnt,w;
38 struct node{
39     int x,ord;
40 }zero[10];
41 
42 bool cmp(node a,node b){
43     return a.x<b.x;
44 }
45 
46 void dfs(int dep,int sum){           
47     if(dep>cnt){
48         ans=max(ans,sum);
49         return;
50     }
51     for(int i=1;i<=9;i++){
52         if(!a[s[dep][1]][i]&&!b[s[dep][2]][i]&&!c[s[dep][3]][i]){
53             a[s[dep][1]][i]=1;
54             b[s[dep][2]][i]=1;
55             c[s[dep][3]][i]=1;
56             dfs(dep+1,sum+i*score[s[dep][1]][s[dep][2]]);
57             a[s[dep][1]][i]=0;
58             b[s[dep][2]][i]=0;
59             c[s[dep][3]][i]=0;
60         }
61     }
62 }
63 
64 int main()
65 {
66 //    freopen("sudo.in","r",stdin);
67 //    freopen("sudo.out","w",stdout);
68     for(int i=1;i<=9;i++)
69         for(int j=1;j<=9;j++){
70             int x;
71             scanf("%d",&x);
72             map[i][j]=x;
73             if(x){
74                 a[i][x]=1;
75                 b[j][x]=1;
76                 c[gong[i][j]][x]=1;
77                 w+=map[i][j]*score[i][j];
78             }
79             else{           
80                 zero[i].x++;
81                 zero[i].ord=i;
82             }
83         }
84     sort(zero+1,zero+10,cmp);//排序
85     for(int i=1;i<=9;i++){
86         for(int j=1;j<=9;j++){
87             if(map[zero[i].ord][j]==0){
88                 s[++cnt][1]=zero[i].ord;       //
89                 s[cnt][2]=j;                   //
90                 s[cnt][3]=gong[zero[i].ord][j];// 0的位置 
91             }
92         }
93     }
94     dfs(1,w);
95     if(ans!=0) printf("%d",ans); 
96     else printf("-1");
97     return 0;
98 }

还有什么不透彻可以随时询问。

 

原文地址:https://www.cnblogs.com/ZAGER/p/9428498.html