洛谷 P1074 靶形数独

题目描述

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他

们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 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

思路：

瞎搞，从后往前搜（改一下搜索顺序就能A太神奇了），具体看给的注释。

#include<cstdio>
#include<cstring>

#define inf 0x3fffffff

int a[10][10], ans = 0;		//二维数组存靶形数独，ans为最终答案。

int direct( int x, int y ) {
	return  (x - 1) / 3 * 3 + (y + 2) / 3;
}	//获取（x，y）在第几个九宫格内。

int hahamax( int x, int y ) {
	return x > y ? x : y; 
}

int hahamin( int x, int y ) {
	return x < y ? x : y;
}

int haha_abs( int x ) {
	return x > 0 ? x : -x;
}

int score( int x, int y) {
	return 10 - hahamax( haha_abs( x - 5 ), haha_abs( y - 5 ) );
}	//（x，y）可以得到的分数。

bool h[10][10], l[10][10], nine[10][10];	//h[第几行][某个数]是否被用过，l[第几列][某个数]是否被用过，nine[第几个九宫格][某个数]是否被用过。

void haharead() {
	for( int i = 1; i <= 9; i++ ) {
		for( int j = 1; j <= 9; j++ ) {
			scanf( "%d", &a[i][j] );
			ans += a[i][j] * score( i, j );		//初始分数。
		}
	}
	return;
}

bool judge = 0;

void hahaput() {
	for( int i = 1; i <= 9; i++ ) {
		for( int j = 1; j <= 9; j++ ) {
			h[i][j] = l[i][j] = nine[i][j] = 1;	//预处理。
		}
	}
	for( int i = 1; i <= 9; i++ ) {
		for( int j = 1; j <= 9; j++ ) {
			if( a[i][j] ) {
				if( h[i][a[i][j]] && l[j][a[i][j]] && nine[direct( i, j )][a[i][j]] ) {
					h[i][a[i][j]] = l[j][a[i][j]] = nine[direct( i, j )][a[i][j]] = 0;	//先把已知的标记。
				}
				else {
					judge = 1;	//用于特判，给出的数独本身就不符合条件。
					return;
				}
			}
		}
	}
	return;
}

bool flag = 0;

void dfs( int mark ) {
	int _min = inf, px = 0, py = 0, cnt, goal = 0, K;
	for( int i = 1; i <= 9; i++ ) {
      		for( int j = 1; j <= 9; j++) {
       			if( !a[i][j] ) {
				cnt = 0;
				K = 0;
				for( int k = 1; k <= 9; k++ ) {//枚举可以用的数。
					if( h[i][k] && l[j][k] && nine[direct( i, j )][k] ) {
						cnt++;//得到该点可用的数的数量。
						K = k;
					}
				}
				goal += K * score( i, j );//此状况下能得到的最大得分。
				if( cnt < _min ) {
					_min = cnt;
					px = i;
					py = j;
				}//得到需要填的最后一个点的坐标。
			}
		}
       }
	//剪枝。
	if( goal + mark <= ans ) {
		return;
	}
	if( goal == 0 ) {
		flag = 1;
		ans = hahamax( ans, mark );
		return;
	}
	if( cnt == 0 ) {
		return;
	}
	//搜索。
	int scr = score( px, py ), _area = direct( px, py );
	for( int i = 1; i <= 9 ; i++ ) {
		if( h[px][i] && l[py][i] && nine[_area][i] ) {
			h[px][i] = l[py][i] = nine[_area][i] = 0;
			a[px][py] = i;
			dfs( mark + i * scr );
			a[px][py] = 0;
			h[px][i] = l[py][i] = nine[_area][i] = 1;
		}
	}
	return;       
}

int main() {
	haharead();
	hahaput();
	if( judge ) {
		printf( "-1" );
		return 0;
	}
	dfs(ans);
	if(flag) {
		printf( "%d", ans );
	}
	else {
		printf( "-1" );
	}
	return 0;
}