【题】【搜索(预处理剪枝)】NOIP2009 靶形数独

NKOJ1693 【NOIP2009】靶形数独
时间限制 : 40000 MS 空间限制 : 128000 KB

问题描述
小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向Z博士请教，Z博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。
　　靶形数独的方格同普通数独一样，在9格宽×9格高的大九宫格中有9个3格宽×3格高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些数字，利用逻辑推理，在其他的空格上填入1到9的数字。每个数字在每个小九宫格内不能重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。（如图）

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

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

输入格式
输入文件名为sudoku.in。
　　　一共9行，每行9个整数（每个数都在0—9的范围内），表示一个尚未填满的数独方格，未填满的空格用“0”表示。每两个数字之间用一个空格隔开。

输出格式
输出文件sudoku.out共1行。
　　输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。

样例输入
【输入样例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

样例输出
【输出样例1】
2829
【输出样例2】
2852

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

思路：
只说一个超强剪枝…
考虑到每次填一个数会让后面填的数的选择减少，那每次就选最少选择的位置填。
定义最少选择为同一行已填数、同一列已填数、同一区域已填数中最大值 最大的一点。
/ * 但是，定义最少选择为同一行已填数+同一列已填数+同一区域已填数 最大 可以快3倍….,不是很懂 */

#include<cstdio>
#include<iostream>
using namespace std;
const int n=9;
typedef int int_[n+1][n+1];
typedef bool bool_[n+1][n+1];

int_ ma;
int_ so={   {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_ jj={   {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_ hh,ll,qq,use;
int ch[n],cl[n],cq[n];
int ans=-1;
int sum[n*n];
int tot,xx[n*n],yy[n*n];

inline int max(int &a,int &b,int &c)
{
    return max(a,max(b,c));
}

void dfs(int pos,int now)
{
    if(pos==tot)
    {
        ans=max(ans,now);
        return ;
    }
    if(now+sum[tot]-sum[pos-1]<ans) return ;
    int i=xx[pos+1],j=yy[pos+1];
    for(int k=1;k<=n;k++)
    {
        if(hh[i][k]||ll[j][k]||qq[jj[i][j]][k]) continue;
        ma[i][j]=k;
        hh[i][k]=ll[j][k]=qq[jj[i][j]][k]=true;
        dfs(pos+1,now+k*so[i][j]);
        ma[i][j]=0;
        hh[i][k]=ll[j][k]=qq[jj[i][j]][k]=false;
    }
}

int main()
{   //in 
    for(int i=1,j;i<=n;i++)
     for(j=1;j<=n;j++)
      scanf("%d",&ma[i][j]);
    //ini 1 
    int tots=0;
    for(int i=1,j;i<=n;i++)
     for(j=1;j<=n;j++)
      if(ma[i][j])
      {
        hh[i][ma[i][j]]=true;
        ll[j][ma[i][j]]=true;
        qq[jj[i][j]][ma[i][j]]=true;
        tots+=ma[i][j]*so[i][j];
        ch[i]++,ch[j]++,cq[jj[i][j]]++;
        use[i][j]=1;
      }
      else tot++;
    //ini 2 
    for(int k=1,i,j,temp,tm;k<=tot;++k)
    {
        tm=0;
        for(i=1;i<=n;i++)
         for(j=1;j<=n;j++)
         {
            if(use[i][j]||(temp=max(ch[i],cl[j],cq[jj[i][j]]))<=tm) continue;
            tm=temp;
            xx[k]=i,yy[k]=j;
         }
        ch[xx[k]]++,ch[yy[k]]++,cq[jj[xx[k]][yy[k]]]++;
        use[xx[k]][yy[k]]=1;
        for(i=9;i>0;--i)
        {
            if(hh[xx[k]][i]||ll[yy[k]][i]||qq[jj[xx[k]][yy[k]]][i]) continue;
            break;
        }
        sum[k]=sum[k-1]+so[xx[k]][yy[k]]*i;
    }

    dfs(0,tots);
    cout<<ans;
}