靶形数独[NOIP2009]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

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

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

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

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

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

思考

直接$\mathrm{DLX}$。
跑到最大即可。
但有一点要注意:

for(int j=r[i];j!=i;j=r[j])remove(col[j]);
for(int j=l[i];j!=i;j=l[j])resume(col[j]);

执行这两个操作时，一定要不同方向枚举，不然会超时。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
const int g[10][10]=
{
	{0,0,0,1,1,1,2,2,2},
	{0,0,0,1,1,1,2,2,2},
	{0,0,0,1,1,1,2,2,2},
	{3,3,3,4,4,4,5,5,5},
	{3,3,3,4,4,4,5,5,5},
	{3,3,3,4,4,4,5,5,5},
	{6,6,6,7,7,7,8,8,8},
	{6,6,6,7,7,7,8,8,8},
	{6,6,6,7,7,7,8,8,8},
};
const int w[10][10]=
{
	{6,6,6,6,6,6,6,6,6},
	{6,7,7,7,7,7,7,7,6},
	{6,7,8,8,8,8,8,7,6},
	{6,7,8,9,9,9,8,7,6},
	{6,7,8,9,10,9,8,7,6},
	{6,7,8,9,9,9,8,7,6},
	{6,7,8,8,8,8,8,7,6},
	{6,7,7,7,7,7,7,7,6},
	{6,6,6,6,6,6,6,6,6},
};
const int N=330;
const int M=3300;
int l[M],r[M],col[M],row[M],u[M],d[M],sz[N],num,head;
int ans,a[N];
void remove(int p)
{
	r[l[p]]=r[p];l[r[p]]=l[p];
	for(int i=u[p];i!=p;i=u[i])
		for(int j=r[i];j!=i;j=r[j])
			u[d[j]]=u[j],d[u[j]]=d[j],--sz[col[j]];
}
void resume(int p)
{
	r[l[p]]=l[r[p]]=p;
	for(int i=d[p];i!=p;i=d[i])
		for(int j=l[i];j!=i;j=l[j])
			u[d[j]]=d[u[j]]=j,++sz[col[j]];
}
void link(int i,int j)
{
	++num;
	u[num]=u[j];d[num]=j;col[num]=j;row[num]=i;
	u[d[num]]=d[u[num]]=num;++sz[j];
	if(!head)head=num;
	else l[num]=num-1,r[num-1]=num;
}
void make(int i,int j,int k)
{
	int len=(i*9+j)*9+k;head=0;
	link(len,i*9+j+1);link(len,81+i*9+k);link(len,162+j*9+k);link(len,243+g[i][j]*9+k);
	l[head]=num;r[num]=head;
}
bool flag;
bool dance(int val)
{
	if(!r[0])
	{
		if(ans<val)ans=val;
		return true;
	}
	int p=r[0];
	for(int i=r[p];i;i=r[i])if(sz[i]<sz[p])p=i;
	remove(p);
	for(int i=u[p];i!=p;i=u[i])
	{ 
		int x=(row[i]-1)/9/9;
		int y=(row[i]-1)/9%9;
		int c=(row[i]-1)%9;++c;
		for(int j=r[i];j!=i;j=r[j])remove(col[j]);
		if(dance(val+w[x][y]*c))flag=true;
		for(int j=l[i];j!=i;j=l[j])resume(col[j]);
	}
	resume(p);
	return flag;
}
int main()
{
	num=324;memset(sz,0,sizeof(sz));ans=0;flag=false;
	for(int i=0;i<=num;i++)r[i]=i+1,l[i]=i-1,u[i]=d[i]=i;
	l[0]=num;r[num]=0;
	for(int i=0;i<9;i++)
		for(int j=0;j<9;j++)
		{
			int x;scanf("%d",&x);
			if(x)make(i,j,x);
			else for(int k=1;k<=9;k++)make(i,j,k);
		}
	if(dance(0))printf("%d\n",ans);
	else puts("-1");
	return 0;
}