https://www.luogu.org/problemnew/show/P1074

这道题单独以每个位置递归纯暴力搜索的话,复杂度9^81,考虑剪枝,和八皇后类似,在同一行同一列同一宫则不能放。另外,想象解答树,先搜索情况少的位置和先搜索情况多的位置总结点数是一样的,不一样的地方在于先搜情况少的的话,靠近树根的分叉少,靠近枝叶稠密,若先搜情况多的,则树根处稠密,枝叶稀疏。考虑剪枝,显然先搜索情况少的,容易很快剪掉很大的一棵棵子树,先搜情况多的则达不到这样的效果。因此按每行的空位排序,空位少的行优先搜索,显然在实际玩数独游戏时,也是这样的方法。

这题写的遇到了很多的曲折。

先是想着以一行为单位,递归9次,但显然,这样的思路是没有办法写出代码的。

改变思路以一格为单位,递归81次,又想到了一个贪心,越靠近中间分数越高,那就优先填大数字,等到搜索到第一组可行解时,就把它作为最优解,但显然,以及一般情况下,这样的贪心都是大错特错的。

然后用先搜索冲突少的地方的方法来写,对于整个9*9方阵按每个点行列宫已有数排序,设置好搜索完每一个点后的下一个点,写了120+行,TLE得很惨。原因是:假如每一行有5个空位,按上面的方法,先搜一下这里,下次跳到很远的地方搜一下,右跳到别处...而这样几次没有相互干扰,这样不利于剪枝。

因此,最后策略:按每行的空位排序,空位少的行优先搜索这一行,这一行搜索完才搜索其他行。在实际玩数独游戏时也是这样的方法。

交上去RE,结果是因为在函数里面定义的数组并不赋初值0,使用里面的值越界了,把它当全局变量用了...

改了以后终于AC了。

#include<bits/stdc++.h>
using namespace std;

int n=9,a[10][10],num[10][10],ans=-1;
int row[10][10],col[10][10],block[10][10];
int nxt[10],hang[10];

const int score[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}};

inline int Block(int i,int j){return (i-1)/3*3+(j-1)/3+1;}

bool cmp(int i,int j)
{
	return hang[i]>hang[j];
}

void init()
{
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)
	{
		scanf("%d",&a[i][j]);
		if(a[i][j])
		{
			hang[i]++;
			row[i][a[i][j]]=1;
			col[j][a[i][j]]=1;
			block[Block(i,j)][a[i][j]]=1;
		}
	}
	
	int d[10]={0};for(int i=1;i<=n;i++)d[i]=i;
	sort(d+1,d+1+n,cmp);
	for(int i=0;i<n;i++)nxt[d[i]]=d[i+1];
}

void search(int r,int c,int tot)
{
	if(!r){ans=max(ans,tot);}
	else if(a[r][c])
	{
		if(c+1<=n)search(r,c+1,tot+a[r][c]*score[r][c]);
		else search(nxt[r],1,tot+a[r][c]*score[r][c]);
	}
	else for(int i=9;i>=1;i--)
	{
		if(row[r][i]||col[c][i]||block[Block(r,c)][i])continue;
		row[r][i]=col[c][i]=block[Block(r,c)][i]=1;
		num[r][c]=i;
		if(c+1<=n)search(r,c+1,tot+num[r][c]*score[r][c]);
		else search(nxt[r],1,tot+num[r][c]*score[r][c]);
		row[r][i]=col[c][i]=block[Block(r,c)][i]=0;
	}
}

int main()
{
	freopen("input.in","r",stdin);
	init();
	search(nxt[0],1,0);
	printf("%d\n",ans);
	//debug();
	return 0;
}