洛谷 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 提高组 第四题

分析：这题很坑，调了好久才弄出来。首先是数独，数独是什么我就不说了，一般人都懂，我们用三个数组来判断每个格能填什么。对于3*3的小格麻烦一点，然后就是最优填，先从可以填的数少的格填（从可以填很多数的地方填出来的基本都是错的，数独都做不到），然后就是统计分数。

代码：

var
n,ans,l,time:longint;
a:array[1..81,1..3]of longint;
c:array[1..27,0..9]of boolean;
v:array[1..81]of longint;
d,d2,d3:array[1..81]of longint;
o:array[1..81]of longint;

function max(a,b:longint):longint;
 begin
  if a>b then exit(a) else exit(b);
 end;

function min(a,b:longint):longint;
 begin
  if a<b then exit(a) else exit(b);
 end;

procedure init;
var
 i,j,k:longint;
begin
 for i:=1 to 9 do
  for j:=1 to 9 do
   begin
    k:=(i-1)*9+j;
    a[k,1]:=i;
    a[k,2]:=j+9;
    a[k,3]:=(i-1)div 3*3+(j-1)div 3+1+18;
    v[k]:=10-max(abs(i-5),abs(j-5));
   end;
 fillchar(c,sizeof(c),1);
 for i:=1 to 81 do
  begin
   read(d[i]);
   c[a[i,1],d[i]]:=false;
   c[a[i,2],d[i]]:=false;
   c[a[i,3],d[i]]:=false;
  end;
end;

procedure check;
var
 i,t:longint;
begin
 t:=0;
 for i:=1 to 81 do
  inc(t,d[i]*v[i]);
 if t>ans then ans:=t;
end;

procedure dfs(dep:longint);
var
  i,k:longint;
begin
 if dep>l then begin check; exit;end;
 inc(time);
 k:=o[dep];
 for i:=1 to 9 do
  if c[a[k,1],i] and c[a[k,2],i] and c[a[k,3],i] then
   begin
    c[a[k,1],i]:=false;
    c[a[k,2],i]:=false;
    c[a[k,3],i]:=false;
    d[k]:=i;
    dfs(dep+1);
    d[k]:=0;
    c[a[k,1],i]:=true;
    c[a[k,2],i]:=true;
    c[a[k,3],i]:=true;
   end;
end;

procedure main;
 var
  i,j,k,t:longint;
begin
 d2:=d;
 for i:=1 to 81 do
  d3[i]:=9;
 for i:=1 to 81 do
  if d2[i]>0 then
   begin
    for j:=1 to 81 do
     if (a[i,1]=a[j,1])or(a[i,2]=a[j,2])or(a[i,3]=a[j,3]) then
    dec(d3[j]);
   end;
 l:=0;
 while 1=1 do
  begin
   k:=maxlongint;
   for i:=1 to 81 do
    if (d2[i]=0)and(d3[i]<k) then
     begin
      k:=d3[i]*11+v[i];
      j:=i;
     end;
   if k=maxlongint then break;
   inc(l);
   o[l]:=j;
   d2[j]:=10;
   for i:=1 to 81 do
    if (a[i,1]=a[j,1])or(a[i,2]=a[j,2])or(a[i,3]=a[j,3]) then
     dec(d3[j]);
  end;
 time:=0;
 ans:=-1;
 dfs(1);
 writeln(ans);
end;

begin
init;
main;
end.