NOIP2009靶形数独

题目描述：

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 Z 博士请教，Z 博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。

靶形数独的方格同普通数独一样，在 9 格宽× 9 格高的大九宫格中有 9 个 3 格宽× 3 格高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些数字，利用逻辑推理，在其他的空格上填入 1到 9 的数字。每个数字在每个小九宫格内不能重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。（如图）

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

总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和。如图，在以下的这个已经填完数字的靶形数独游戏中，总分数为 2829。游戏规定，将以总分数的高低决出胜负。

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

输入输出格式

输入格式：

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

输出格式：

输出共 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

输出样例#1：

2829

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

输出样例#2：

2852

说明

【数据范围】

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

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

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

NOIP 2009 提高组 第四题

（转自洛谷P1074）

这道题做法较多，用拓扑排序或状压dp也行

这里是dfs+剪枝（优化搜索顺序）

若是每遍都遍历所有没填的位置，那程序运行太慢了

所以可以每遍优先搜索可填的数最少的点

预处理时将每行每列已知数的个数存起来，每遍搜索时找行列已知数最多的格子填数，这样可以大大优化程序运行的时间。

如输入样例#1，开始时应找（8,8）这个格子填数

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;

inline int read()
{
    int f=,x=;
    char ch=getchar();
    while(ch<'' || ch>'') {if(ch=='-') f=-; ch=getchar();}
    while(ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
    return x*f;
}

int cnt,ans;
int a[][],rowcnt[],colcnt[];
bool row[][],col[][],sqr[][];
int sco[][]={{,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,},
                 {,,,,,,,,,}};

int g(int x,int y)
{
    return (x-)/*++(y-)/;
}

void dfs(int r,int c,int step)
{
    int i,j,k;
    if(step==)
    {
        int score=;
        for(i=;i<=;i++)
            for(j=;j<=;j++)
                score+=a[i][j]*sco[i][j];
        ans=max(ans,score);
        return ;
    }
    for(k=;k<=;k++)
    {
        if(row[r][k] || col[c][k] || sqr[g(r,c)][k]) continue;
        row[r][k]=;
        col[c][k]=;
        sqr[g(r,c)][k]=;
        a[r][c]=k;
        rowcnt[r]++;
        colcnt[c]++;
        int maxrow=-,r1=,maxcol=-,c1=;
        for(i=;i<=;i++)
            if(rowcnt[i]>maxrow && rowcnt[i]<)
            {
                maxrow=rowcnt[i];
                r1=i;
            }
        for(i=;i<=;i++)
            if(colcnt[i]>maxcol && !a[r1][i])
            {
                maxcol=colcnt[i];
                c1=i;
            }
        dfs(r1,c1,step+);
        row[r][k]=;
        col[c][k]=;
        sqr[g(r,c)][k]=;
        a[r][c]=;
        rowcnt[r]--;
        colcnt[c]--;
    }
}

int main()
{
    int i,j;
    for(i=;i<=;i++)
    {
        for(j=;j<=;j++)
        {
            a[i][j]=read();
            if(a[i][j])
            {
                row[i][a[i][j]]=;
                col[j][a[i][j]]=;
                sqr[g(i,j)][a[i][j]]=;
                rowcnt[i]++;
                colcnt[j]++;
                cnt++;
            }
        }
    }
    int maxrow=-,r=,maxcol=-,c=;
    for(i=;i<=;i++)
        if(rowcnt[i]>maxrow && rowcnt[i]<)
        {
            maxrow=rowcnt[i];
            r=i;
        }
    for(i=;i<=;i++)
        if(colcnt[i]>maxcol && !a[r][i])
        {
            maxcol=colcnt[i];
            c=i;
        }
    dfs(r,c,cnt);
    if(ans==) ans=-;
    printf("%d",ans);
    return ;
}

这是代码