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

#include<cmath>
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;
#define maxn 10
int mark[maxn][maxn],a[maxn][maxn],b[maxn][maxn],in[maxn][maxn];
int ans=-1;
int heng[maxn],shu[maxn],xie[maxn];
void DFS(int now){
    int next=0x3f3f3f3f,nx=0,ny=0;
    for(int i=1;i<=9;++i)
        for(int j=1;j<=9;++j)
            if(!b[i][j]&&!a[i][j]){
                int cnt=0,tmp=heng[i]|shu[j]|xie[in[i][j]];
                for(int k=1;k<=9;++k)
                    if(!(tmp&(1<<k))) ++cnt;
                if(cnt==0) return;
                if(cnt<next) next=cnt,nx=i,ny=j;

            }
    if(next==0x3f3f3f3f){
        ans=max(ans,now);
        return ;
    }
    int tmp=heng[nx]|shu[ny]|xie[in[nx][ny]];
    for(int i=9;i>=1;--i)
        if(!(tmp&(1<<i))){
            heng[nx]^=1<<i,shu[ny]^=1<<i,xie[in[nx][ny]]^=1<<i;
            b[nx][ny]=1,a[nx][ny]=i;
            DFS(now+i*mark[nx][ny]);
            heng[nx]^=1<<i,shu[ny]^=1<<i,xie[in[nx][ny]]^=1<<i;
            b[nx][ny]=0,a[nx][ny]=0;
        }
}
int main(){
    for(int i=1;i<=9;++i)
        for(int j=1;j<=9;++j)
            mark[i][j]=min(min(i,10-i),min(j,10-j))+5;
    for(int i=1;i<=9;++i)
        for(int j=1;j<=9;++j)
            in[i][j]=(i-1)/3*3+(j+2)/3;
    int now=0;
    for(int i=1;i<=9;++i)
        for(int j=1;j<=9;++j){
            scanf("%d",&a[i][j]);
            if(a[i][j]){
                if( (heng[i]&(1<<a[i][j])) || (shu[j]&(1<<a[i][j])) || (xie[in[i][j]]&(1<<a[i][j])) ){
                    printf("-1\n");
                    return 0;
                }
                heng[i]+=(1<<a[i][j]);
                shu[j]+=(1<<a[i][j]);
                xie[in[i][j]]+=(1<<a[i][j]);
                b[i][j]=1;
                now+=a[i][j]*mark[i][j];
            }
        }
    DFS(now);
    printf("%d\n",ans);
    return 0;
}

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#define Ll long long
using namespace std;
int v[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 g[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},
};
bool h[10][10],l[10][10],G[10][10],c[90];
int map[10][10],x[90],y[90],n,ans;
void UpDate_Ans(){
    int sum=0;
    for(int i=1;i<=9;++i)
        for(int j=1;j<=9;++j)
            sum+=map[i][j]*v[i][j];
    ans=max(ans,sum);
}
void DFS(int k){
    if(k>n){ UpDate_Ans(); return ; }
    int w=1e9,p;
    for(int i=1;i<=n;++i)
        if(!c[i]){
            int ww=0;
            for(int j=1;j<=9;++j)
                if(!h[x[i]][j]&&!l[y[i]][j]&&!G[ g[x[i]][y[i]] ][j])
                    if(++ww==w) break;
            if(ww<w) w=ww,p=i;
        }
    int xx=x[p],yy=y[p];
    c[p]=1;
    for(int i=1;i<=9;++i)
        if(!h[xx][i]&&!l[yy][i]&&!G[g[xx][yy]][i]){
            map[xx][yy]=i,h[xx][i]=1,l[yy][i]=1,G[g[xx][yy]][i]=1;
            DFS(k+1);
            h[xx][i]=0,l[yy][i]=0,G[g[xx][yy]][i]=0;
        }
    c[p]=0;
}
int main(){
    for(int i=1;i<=9;++i)
        for(int j=1;j<=9;++j){
            scanf("%d",&map[i][j]);
            if(map[i][j]) h[i][map[i][j]]=1,l[j][map[i][j]]=1,G[g[i][j]][map[i][j]]=1;
            else x[++n]=i,y[n]=j;
        }
    ans=-1;
    DFS(1);
    printf("%d\n",ans);
    return 0;
}

这道题并不是什么A*，纯粹的启发式搜索，就是暴搜啊！