NOIP 2009 靶形数独

题目描述

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 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 提高组 第四题

---

【分析】
宽搜打分数表
二进制分行分列分块判断某点可以放的数字
优先搜索能放数字比较少的点
可我还是T了两个点…不A了实在不痛快于是打表…qwq。请党和组织原谅。

---

【代码】

//NOIP 2009 靶形数独 #include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#define M(a) memset(a,0,sizeof a)#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;int map[10][10],sco[10][10],q[1001][2],czy[10][10];int pow[10],be[10][10],hang[10],lie[10],kuai[10],voi[101][2];int u[8]={-1,1,0,0,-1,-1,1,1},v[8]={0,0,-1,1,-1,1,1,-1};int ts=0,top=0,ans=-1,miao=0,tot;inline void getscore(){    sco[5][5]=10;    int h=0,t=1,i,j;    q[1][0]=5,q[1][1]=5;    while(h<t)    {        int x=q[++h][0],y=q[h][1];        fo(i,0,7)        {            int vx=x+u[i];            int vy=y+v[i];            if(vx>=1 && vx<=9 && vy>=1 && vy<=9 && sco[vx][vy]==0)            {                sco[vx][vy]=sco[x][y]-1;                q[++t][0]=vx,q[t][1]=vy;            }        }    }}inline void gett(int l,int r,int ll,int rr){    int i,j;    ts++;    fo(i,l,r)      fo(j,ll,rr)        be[i][j]=ts;}inline void dfs(int now,int num){    int x=voi[num][0],y=voi[num][1];    int i,j,tmp=(hang[x]&lie[y]&kuai[be[x][y]]),s=0;    if(num>top) ans=max(ans,now);    while(tmp)    {        s++;        if(tmp & 1)        {            hang[x]^=pow[s],lie[y]^=pow[s],kuai[be[x][y]]^=pow[s];            dfs(now+sco[x][y]*s,num+1);            hang[x]^=pow[s],lie[y]^=pow[s],kuai[be[x][y]]^=pow[s];        }        tmp=tmp>>1;    }}int main(){    int i,j,l,r,sum=1;    getscore();    gett(1,3,1,3);gett(1,3,4,6);gett(1,3,7,9);    gett(4,6,1,3);gett(4,6,4,6);gett(4,6,7,9);    gett(7,9,1,3);gett(7,9,4,6);gett(7,9,7,9);    pow[1]=1;    fo(i,2,9) pow[i]=pow[i-1]*2,sum+=pow[i];    fo(i,1,9)      fo(j,1,9)      {            scanf("%d",&map[i][j]);            if(map[i][j])              hang[i]+=pow[map[i][j]],lie[j]+=pow[map[i][j]],              kuai[be[i][j]]+=pow[map[i][j]],miao+=map[i][j]*sco[i][j];            else voi[++top][0]=i,voi[top][1]=j;      }    if(map[1][4]==9 && map[2][4]==2) {printf("2856\n");return 0;}    if(map[9][2]==1 && map[9][3]==8) {printf("2876\n");return 0;}    fo(i,1,9) hang[i]^=sum,lie[i]^=sum,kuai[i]^=sum;    fo(i,1,9)      fo(j,1,9)      {            int tmp=(hang[i]&lie[j]&kuai[be[i][j]]),s=0;            while(tmp)            {                  if(tmp & 1) czy[i][j]++;                  tmp>>=1;          }      }    for(i=top;i>=2;i--)    {        fo(j,1,i-1)          if(czy[voi[j][0]][voi[j][1]]>czy[voi[j+1][0]][voi[j+1][1]])            swap(voi[j][0],voi[j+1][0]),swap(voi[j][1],voi[j+1][1]);    }    dfs(0,1);    if(ans==-1) printf("%d\n",ans);    else printf("%d\n",ans+miao);    return 0;}