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

---

二进制神题：只需要记录这行这列这个方块那些数字不能填，可以填的数就是并集的补集。
一开始：枚举格子，枚举每一个可以填的数。（75分）
考虑剪枝：先搜索可能填的数较少的格子。（90分）
运用位运算优化，可以ac。

---

#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int ans,blc[15],r[15],c[15],a[15][15],cnt;
int score[15][15]=
{
    {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 clc()
{
    int sum=0;
    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
            sum+=score[i][j]*a[i][j];
    return sum;
}
int count(int x)
{
    int res=9;
    for(int i=1;i<=9;i++)
    {
        res-=x&1;
        x>>=1;
    }
    return res;
}
void dfs(int u)
{
    if(u>cnt)
    {
        ans=max(ans,clc());
        return ;
    }
    //从剩下的可以填的格中选一个可能性最少的 
    int mn=1e9+7,tx,ty,tres,bt;
    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
            if(a[i][j]==0)
            {
                int t,res,bcnt;
                t=(ceil((double)i/3)-1)*3+ceil((double)j/3);
                res=r[i]|c[j]|blc[t];
                res>>=1;
                bcnt=count(res);
                if(bcnt<mn)
                {
                    mn=bcnt;
                    tx=i;
                    ty=j;
                    tres=res;
                    bt=t;
                }
            }
    for(int i=1;i<=9;i++)
    {
        int tmp=tres&1;//取最后一位 
        if(tmp==0)
        {
            a[tx][ty]=i;
            r[tx]^=1<<i;
            c[ty]^=1<<i;
            blc[bt]^=1<<i;
            dfs(u+1);
            r[tx]^=1<<i;
            c[ty]^=1<<i;
            blc[bt]^=1<<i;
            a[tx][ty]=0;
        }
        tres>>=1;
    }
}
int main()
{
    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
        {
            scanf("%d",&a[i][j]);
            if(a[i][j]>0)
            {
                int t=(ceil((double)i/3)-1)*3+ceil((double)j/3);
                r[i]^=1<<a[i][j];
                c[j]^=1<<a[i][j];
                blc[t]^=1<<a[i][j];
            }
            else
                cnt++;
        }
    dfs(1);
    if(ans!=0)
        printf("%d\n",ans);
    else
        printf("-1\n");
    return 0;
}