[NOIP2009] 提高组洛谷P1074 靶形数独

题目描述

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

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

首先需要搜索出所有的数独解，然后每次算出该解对应的靶形数独值，记录最大值即可。

剪枝：存每行、每列、每正方块出现的数的集合h[],r[],block[]，对于每个0，求其对应h,r,block的并，就可以很快地知道该位置可以填的数。

剪枝2：↑用位运算处理以上操作。

剪枝3:计算每个0可能可以填的数的数量，每次填可填数量最小的位置。

刚开始没有用位运算，最高卡到90分……

之后加上了剪枝2和3，然而因为剪枝3中最小位置的信息记录错误，无限T……

↑又花了一个小时才查出来错，悲伤。

↑最近怎么尽是犯这些低级错误……

 /*by SilverN*/

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 using namespace std;

 int sc[][]={{,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,},

                 {,,,,,,,,,}};

 int mp[][];

 int ans=-;

 int block[],h[],r[];

 int in[][];

 int sx[],sy[],cnt=;

 inline int blcnt(int x){

     int res=;

     while(x){

         res-=(x&);

         x>>=;

     }

     return res;

 }

 void init(){

     for(int i=;i<=;i++)

      for(int j=;j<=;j++)

         in[i][j]=(i+)/*+(j+)/-;

 }

 void dfs(int sum){

     int i,j;

     int mini=,pos=;

     int tmp,now;

     int ctmp;

     for(i=;i<=cnt;i++){

         if(!mp[sx[i]][sy[i]]){

             tmp=h[sx[i]]|r[sy[i]]|block[in[sx[i]][sy[i]]];

             now=blcnt(tmp);

             if(!now)return;

             if(now<mini){

                 ctmp=tmp;//记录最小值对应的tmp

                 //先前由于没有记录ctmp，而在后面计算时用tmp算，出现了bug，费了一个小时查出来错

                 mini=now;

                 pos=i;

             }

         }

     }

     if(mini==){

         ans=max(ans,sum);

         return;

     }

     i=sx[pos];

     j=sy[pos];

              for(int k=;k<=;k++){

                  if(!(ctmp&(<<k))){

                      mp[i][j]=k;

                      h[i]^=(<<k);

                      r[j]^=(<<k);

                      block[in[i][j]]^=(<<k);

                      dfs(sum+k*sc[i][j]);

                      //回溯

                      block[in[i][j]]^=(<<k);

                      h[i]^=(<<k);

                      r[j]^=(<<k);

                      mp[i][j]=;

                  }

              }

     return;

 }

 int main(){

     int i,j;

     cnt=;

     int su=;

     init();

     for(i=;i<=;i++)//行

      for(j=;j<=;j++){//列

          scanf("%d",&mp[i][j]);

          if(mp[i][j]){

              r[j]^=<<mp[i][j];

              h[i]^=<<mp[i][j];

              block[in[i][j]]^=<<mp[i][j];

              su+=sc[i][j]*mp[i][j];

         }

         else{

             cnt++;

             sx[cnt]=i;

             sy[cnt]=j;

         }

     }

     dfs(su);

     printf("%d\n",ans);

     return ;

 }