noip2009靶形数独

题目描述

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向 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 提高组 第四题

解题思路：启发式搜索，先预处理分数，并根据每个未填的空所在行列（优先为列，排除行列已知数相同的情况）已知数的个数入队，搜索

#include <cstdio>
#define max(a,b) a>b?a:b
//启发式搜索，周围数字多的先搜索， 预处理分数 
int mp[9][9][11], x[10][11], y[10][11], t = 6, a[11][11], m, f[11][11], x1, y1;
int s[11][11];
int que[1235][3], cnt, ans;
int read(){
	int x, f = 1;
	char ch;
	while(ch = getchar(), ch<'0'||ch>'9') if(ch == '-') f = -1;
	x = ch - 48;
	while(ch = getchar(), '0'<=ch&&ch<='9') x = x * 10 + ch - 48;
	return x * f;
}
void dfs(int k, int sum){
	int i = que[k][1], j = que[k][2], t;
	if(k > cnt){
		ans = max(ans, sum);
	}
	else{
		for(t = 1; t <= 9; t++){
			if(!mp[i/3][j/3][t]&&!x[i][t]&&!y[j][t]){
				mp[i/3][j/3][t] = x[i][t] = y[j][t] = 1;
				dfs(k+1,sum+a[i][j]*t);
				mp[i/3][j/3][t] = x[i][t] = y[j][t] = 0;
			}
		}
	}
}
int main(){
	int k, i, j;
	for(k = 0; k <= 4; k++){
		for(i = k; i < 9-k; i++){
			for(j = k; j < 9-k; j++){
				a[i][j] = t;
			}
		}
		t++;
	}
	for(i = 0; i < 9; i++){
		for(j = 0; j < 9; j++){
			t = read();
			if(!t) continue;
			if(!mp[i/3][j/3][t]++&&!x[i][t]++&&!y[j][t]++){
				m += t*a[i][j];
				f[i][j] = k;
			}
			else{
				printf("-1");
				return 0;
			}
		}
	}
	while(1){//看那一个空周围数字多 
		x1 = y1 = 9;
		for(i = 0; i < 9; i++){
			for(j = 0; j < 9; j++){
				if(!f[i][j]){
					s[i][j] = 0;//重置 
					for(k = 0; k < 9; k++){
						if(f[i][k]) s[i][j] += 123;
						if(f[k][j]) s[i][j]++;
					}
					if(s[i][j] > s[x1][y1]){
						x1 = i;
						y1 = j;
					}
				}
			}
		}
		if(s[x1][y1]){
			f[x1][y1] = -1; que[++cnt][1] = x1; que[cnt][2] = y1;//标记为已搞过 , 入队 
		}
		else break;
	}
	dfs(1,m);
	printf("%d", ans?ans:-1);
	return 0;
}