NOIP2009 靶形数独

https://daniu.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。

玩过数独的朋友都知道，填数先找只能填一个数的格子，或者可以选的数最少的格子

所以记录所有没有数的格子，找出可以填数最少的那一个，dfs

很重要的细节优化：如果出现一个格子可以填的数=1，那么break

没有这个第9和13个点一直TLE

这启示我们不要轻视dfs中任何一个剪枝

#include<cstdio>
#include<algorithm>
using namespace std;
int sum,ans,tmp,can;
int x[100],y[100];
bool v[100],h[10][10],l[10][10],bl[10][10];
int k[10][10],f[10][10];
void dfs(int now,int tot)
{
    if(now==sum+1) 
    {
        ans=max(ans,tot);
        return;
    }
    int minn=10,dx=0;
    for(int i=1;i<=sum;i++)
     if(!v[i])
     {
          can=0;
         for(int j=1;j<=9;j++)
           if(!h[x[i]][j]  && !l[y[i]][j] && !bl[k[x[i]][y[i]]][j]) can++;    
              if(can<minn) 
              { 
                 minn=can; dx=i;
                 if(minn==1) break; 
              }
     }
    v[dx]=true;
    for(int P=9;P>=1;P--)
    {
        if(!h[x[dx]][P] && !l[y[dx]][P] && !bl[k[x[dx]][y[dx]]][P])
          {
                h[x[dx]][P]=true; l[y[dx]][P]=true;  bl[k[x[dx]][y[dx]]][P]=true;
                dfs(now+1,tot+P*f[x[dx]][y[dx]]);
                h[x[dx]][P]=false; l[y[dx]][P]=false;  bl[k[x[dx]][y[dx]]][P]=false;
          }
    }
    v[dx]=false;
}
void pre()
{
    for(int i=1;i<=9;i++)
     for(int j=1;j<=9;j++)
     {
          k[i][j]=(i-1)/3*3+(j-1)/3+1;
          if(i==1 || i==9 || j==1 || j==9) f[i][j]=6;
          else if(i==2 || i==8 || j==2 || j==8) f[i][j]=7;
          else if(i==3 || i==7 || j==3 || j==7) f[i][j]=8;
          else if(i==4 || i==6 || j==4 || j==6) f[i][j]=9;
          else f[i][j]=10;     
     }
}
int main()
{
    pre();
    int u;
    for(int i=1;i<=9;i++)
     for(int j=1;j<=9;j++)
      {
          scanf("%d",&u);
          if(!u) x[++sum]=i,y[sum]=j;
          else h[i][u]=true,l[j][u]=true,bl[k[i][j]][u]=true,tmp+=u*f[i][j];
      }
    dfs(1,0);
    if(!ans) printf("-1");
    else printf("%d",ans+tmp);
}