[WIKIOI 1174]靶形数独(深度优先搜索)

题目描述 Description

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他
们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向Z 博士请教，
Z 博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。
靶形数独的方格同普通数独一样，在 9 格宽×9 格高的大九宫格中有9 个3 格宽×3 格
高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些

数字，利用逻辑推理，在其他的空格上填入1 到9 的数字。每个数字在每个小九宫格内不能
重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即
每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。

 

上图具体的分值分布是：最里面一格（黄色区域）为 10 分，黄色区域外面的一圈（红
色区域）每个格子为9 分，再外面一圈（蓝色区域）每个格子为8 分，蓝色区域外面一圈（棕
色区域）每个格子为7 分，最外面一圈（白色区域）每个格子为6 分，如上图所示。比赛的
要求是：每个人必须完成一个给定的数独（每个给定数独可能有不同的填法），而且要争取
更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字
的乘积的总和。如图，在以下的这个已经填完数字的靶形数独游戏中，总分数为2829。游
戏规定，将以总分数的高低决出胜负。

由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能
够得到的最高分数。

输入描述 Input Description

一共 9 行。每行9 个整数（每个数都在0—9 的范围内），表示一个尚未填满的数独方
格，未填的空格用“0”表示。每两个数字之间用一个空格隔开。

输出描述 Output Description

输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。

样例输入 Sample Input

【输入输出样例 1】

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

【输入输出样例 2】

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

样例输出 Sample Output

【输入输出样例 1】

2829

【输入输出样例 1】

2852

数据范围及提示 Data Size & Hint

【数据范围】
40%的数据，数独中非0 数的个数不少于30。
80%的数据，数独中非0 数的个数不少于26。
100%的数据，数独中非0 数的个数不少于24。

好恶心的一道题...优化过的深度优先搜索还是可以过的.

#include<stdio.h>
#include<math.h>
int h[10]={},hs[10]={},zs[10]={},xj[5][5]={},hq[10]={};
int ans=-1,st[10],a[10][10];
void make()
{
    int sum=0,i,j;
    for (i=1;i<5;i++)
      {for (j=i;j<11-i;j++)
        sum+=(a[i][j]+a[10-i][j])*(5+i);
       for (j=i+1;j<10-i;j++)
        sum+=(a[j][i]+a[j][10-i])*(5+i);
        }
    sum+=a[5][5]*10;   
    if (sum>ans) ans=sum;
}
void dfs(int k)
{
     if (k==10) make();
     else
     {
      int x,y,j,pos,p,i=st[k];
      x=511-h[i];
      y=x&-x; //
      h[i]|=y;///
      j=(int)log2(y)+1; ///
      pos=511-(hs[i]|zs[j]|xj[(i-1)/3][(j-1)/3]); ///
      while (pos>0)
      {p=pos&-pos;  ///
       pos-=p;       ///
       a[i][j]=(int)log2(p)+1; ///
       hs[i]|=p;    //
       zs[j]|=p;
       xj[(i-1)/3][(j-1)/3]|=p;
       if (x==y) dfs(k+1);///
       else dfs(k); //
       hs[i]-=p; //
       zs[j]-=p;
       xj[(i-1)/3][(j-1)/3]-=p;
       };
       h[i]-=y;  
       };
}
int main()
{
    int i,j,p0;
    for (i=1;i<10;i++)
      for (j=1;j<10;j++)
        {scanf("%d",&a[i][j]); 
         if (a[i][j]>0)
           {h[i]|=1<<(j-1);  
                       
            p0=1<<(a[i][j]-1);
            if (((hs[i]&p0)!=0)||((zs[j]&p0)!=0)||
                       ((xj[(i-1)/3][(j-1)/3]&p0)!=0))
            {printf("-1\n");return 0;};  //
            hs[i]|=p0;  //
                       //
            zs[j]|=p0;    //
 
            xj[(i-1)/3][(j-1)/3]|=p0; //
            }
          else hq[i]++;}  //0
    for (i=1;i<10;i++) st[i]=i;        
    for (i=1;i<9;i++)                 //
      for (j=i+1;j<10;j++)      //
       if (hq[st[i]]>hq[st[j]])
         {st[i]^=st[j];               //0
          st[j]^=st[i];
          st[i]^=st[j];}
         
    for (i=1;hq[st[i]]==0;i++);         //
    dfs(i);                          
    printf("%d\n",ans);              
    return 0;
}