NOIP2009 靶形数独

Posted 2020-10-11 欢迎来到Frank滑稽至上主义的教室

题目描述

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

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

 

思路：暴力做法基本上与数独挑战一样，都是直接枚举+判断，然而这样有一个问题。数据给的空位可能会很多，这样搜索复杂度也会很高，怎么办呢？这时候我们就需要一个玄学优化：优先搜索可能情况少的，这样先局部减少复杂度，再限制后面的复杂度的增长，可以减少一些总体的复杂度，当然还有其他的优化，例如：倒着搜、倒着枚举数字、O2优化、inline等等。我主要是用的优先搜索少的。

搜索很简单，枚举即可。

PS:这题排序有坑，我把sort换成stable_sort就玄学过了。。。（stable_sort大法好）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
const int maxn=12;
inline int read()
{
    int ret=0,f=1;
    char c=getchar();
    while(c<‘0‘||c>‘9‘)
    {if(c==‘-‘) f=-1;c=getchar();}
    while(c>=‘0‘&&c<=‘9‘)
    {ret=ret*10+c-‘0‘;c=getchar();}
    return ret*f;
}
int mmp[maxn][maxn];
struct node{
    int x,y,k;
}p[maxn*maxn];
bool row[maxn][maxn],col[maxn][maxn],gz[4][4][maxn];
int tot=0;
int f[maxn][maxn]={{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},
};
inline void place(int x,int y,int v,bool t)
{
    row[x][v]=t,col[y][v]=t;
    gz[(x-1)/3][(y-1)/3][v]=t;
}
bool cmp(node A,node B)
{
    return A.k>B.k;
}
void make_pre()
{
    for(int i=1;i<=tot;i++)
    {
        int x=p[i].x,y=p[i].y;
        int xx=(x-1)/3,yy=(y-1)/3;
        for(int j=1;j<=9;j++)
        {
            if(row[x][j]||col[y][j]||gz[xx][yy][j]) p[i].k++;
        }
    }
    stable_sort(p+1,p+1+tot,cmp);
}
int ans=0;
void dfs(int step,int sum)
{
    if(step==tot+1)
    {
        ans=max(ans,sum);
        return;
    }
    if(sum+(tot+1-step)*90<=ans) return;
    int x=p[step].x,y=p[step].y;
    int xx=(x-1)/3,yy=(y-1)/3;
    for(int i=1;i<=9;i++)
    {
        if(row[x][i]||col[y][i]||gz[xx][yy][i]) continue;
        mmp[x][y]=i;
        place(x,y,i,1);
        dfs(step+1,sum+i*f[x][y]);
        mmp[x][y]=0;
        place(x,y,i,0);
    }
    return;
}
int main()
{
    int ansp=0;
    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
        {
            mmp[i][j]=read();
            ansp+=mmp[i][j]*f[i][j];
            if(!mmp[i][j])
            {
                tot++;
                p[tot].x=i,p[tot].y=j;
                
            }
            else place(i,j,mmp[i][j],1);
        }
    make_pre();
    dfs(1,ansp);
    if(!ans) printf("-1\n");
    else printf("%d\n",ans);
    return 0;
}