靶形数独

靶形数独

时间限制: 2 Sec  内存限制: 128 MB

题目描述

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

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

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

输入

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

输出

输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-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

样例输出

2829

提示

40%的数据，数独中非 0数的个数不少于 30。 

80%的数据，数独中非 0数的个数不少于 26。 

100%的数据，数独中非 0 数的个数不少于 24。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
int p[10][10][10];
bool l[10][10];
int o[10][10],c;
int s[10][10],ans,mmax,mmin,i2,j2;
int  a[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
                };
bool find()
{
    mmin=8888;
    int i,j;
    for(i=1; i<=9; i++)
        for(j=1; j<=9; j++)
            if(o[i][j]==0&&l[i][j]==false)
                return false;
            else if(l[i][j]==false&&o[i][j]<mmin)
            {
                mmin=o[i][j];
                i2=i;
                j2=j;
            }
    return true;
}
void putin(int x,int y,int z)
{
    int i,j,x2=(x-1)/3+1,y2=(y-1)/3+1;
    for(i=1; i<=9; i++)
    {
        if(p[x][i][z]==0)o[x][i]--;
        p[x][i][z]++;
    }
    for(i=1; i<=9; i++)
        {
            if(p[i][y][z]==0)o[i][y]--;
            p[i][y][z]++;
        }
    for(i=x2*3-2; i<=x2*3; i++)
        for(j=y2*3-2; j<=y2*3; j++)
        {
            if(p[i][j][z]==0)o[i][j]--;
            p[i][j][z]++;
        }
}
void putout(int x,int y,int z)
{
    int i,j,x2=(x-1)/3+1,y2=(y-1)/3+1;
    for(i=1; i<=9; i++)
    {
        p[x][i][z]--;
        if(p[x][i][z]==0)o[x][i]++;
    }
    for(i=1; i<=9; i++)
        {
            p[i][y][z]--;
            if(p[i][y][z]==0)o[i][y]++;
        }
    for(i=x2*3-2; i<=x2*3; i++)
        for(j=y2*3-2; j<=y2*3; j++)
        {
            p[i][j][z]--;
            if(p[i][j][z]==0)o[i][j]++;
        }
}
void dfs(int k)
{
    int i,j;
    if(k==82)
    {
        if(mmax<ans)mmax=ans;
    }
    else
    {
        if(!find())
        return;
        int d=i2,f=j2;
        for(i=1; i<=9; i++)
        {
            if(p[d][f][i]==0)
            {
                putin(d,f,i);
                l[d][f]=true;
                ans=ans+i*a[d][f];
                dfs(k+1);
                ans=ans-i*a[d][f];
                l[d][f]=false;
                putout(d,f,i);
            }
        }
    }
}
int main()
{
    int i,j,k;
    for(i=0; i<=9; i++)
        for(j=0; j<=9; j++)
            o[i][j]=9;
    for(i=1; i<=9; i++)
        for(j=1; j<=9; j++)
        {
            cin>>s[i][j];
            if(s[i][j]!=0)
            {
                putin(i,j,s[i][j]);
                l[i][j]=true;
                c++;
                ans=ans+s[i][j]*a[i][j];
            }
        }
    mmax=-1;
    dfs(c+1);
    cout<<mmax;
    return 0;
}