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输入样例#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 数的个数不少于30。
80%的数据，数独中非 0 数的个数不少于26。
100%的数据，数独中非0数的个数不少于24。
NOIP 2009 提高组 第四题
这真是一道恶心的搜索题，首先你得先熟悉数独。
玩家需要根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫（3*3）内的数字均含1-9，不重复。
————百度百科
所以我们有三个数组
分别代表行，列，九宫格中数字使用情况。
行和列都还好，循环判断即可，但九宫格，就需要一个函数
1 int ninth( int i , int j ) { 2 if( i <= 3 && j <= 3 ) return 1 ; 3 if( i <= 3 && j <= 6 ) return 2 ; 4 if( i <= 3 ) return 3 ; 5 if( i <= 6 && j <= 3 ) return 4 ; 6 if( i <= 6 && j <= 6 ) return 5 ; 7 if( i <= 6 ) return 6 ; 8 if( j <= 3 ) return 7 ; 9 if( j <= 6 ) return 8 ; 10 return 9 ; 11 }
接着我们处理每个格子的分数，我是用函数判断，而ly用数组，最后时间来看，我还是太蒟了
1 inline int check(int x,int y) 2 { 3 if(x==1||y==1||x==9||y==9)return 6; 4 if(x==2||y==2||x==8||y==8)return 7; 5 if(x==3||y==3||x==7||y==7)return 8; 6 if(x==4||y==4||x==6||y==6)return 9; 7 if(x==5&&y==5)return 10; 8 } ————by cx
1 int point[ 6 ] = { 0 , 6 , 7 , 8 , 9 , 10} ; 2 ————by ly
预处理完后便可以开始考虑如何搜索，这里我与ly有了不同的搜法。
1.ly
ly选择的是一行一行搜下去，全部答案搜出来后，再求最大值。
1 int dfs( int h , int x , int y ) { 2 if( h == 10 ) { 3 print( ); 4 return 0 ; 5 } 6 if( y == 10 ) { 7 dfs( h + 1 , sss[ h + 1 ].line , 1 ); 8 return 0 ; 9 } 10 if( !map[ x ][ y ] ) { 11 for( int i = 1 ; i <= 9 ; ++i ) { 12 if( line[ x ][ i ] == 0 && list[ y ][ i ] == 0 && nine[ ninth( x , y ) ][ i ] == 0 ){ 13 line[ x ][ i ] = 1 , list[ y ][ i ] = 1 , nine[ ninth( x , y ) ][ i ] = 1 ; 14 map[ x ][ y ] = i ; 15 dfs( h , x , y + 1 ); 16 map[ x ][ y ] = 0 ; 17 line[ x ][ i ] = 0 , list[ y ][ i ] = 0 , nine[ ninth( x , y ) ][ i ] = 0 ; 18 } 19 } 20 } 21 else dfs( h , x , y + 1 ); 22 }
2.cx
我选择的是用一个数组存要填的的点，一个一个搜，我一开始以为我的方法会快一点，结果我被打脸了
1 void dfs(int x,int y) 2 { 3 if(tot==pos-1){maxn=max(ans,maxn);return;} 4 for(register int i=1;i<=9;++i) 5 { 6 if(!line[x][i]&&!list[y][i]&&!nine[ninth(x,y)][i]) 7 { 8 pos++;ans+=i*check(x,y); 9 line[x][i]=list[y][i]=1; 10 nine[ninth(x,y)][i]=1; 11 dfs(b[0][pos].w,b[1][pos].w); 12 pos--;ans-=i*check(x,y); 13 line[x][i]=list[y][i]=0; 14 nine[ninth(x,y)][i]=0; 15 } 16 } 17 }
接下来便是重点，搜索剪枝
像我这种没玩过数独的乡里人，不知道玩数独有这样一个方法：
从数多的一行开始填，这样要选择的数就少了，不合法的情况就可以省掉一些
所以我们定义一个（
struck），用一个来存每行
的个数。
用它作为关键字一遍后，再从最少的开始搜。
ly程序：(用时: 2678ms / 内存: 920KB)
1 #include<iostream> 2 #include<cstdio> 3 #include<algorithm> 4 using namespace std ; 5 6 inline int read( ) {//日常快读 7 int x = 0 , f = 1 ; 8 char c = getchar( ) ; 9 while( c > '9' || c < '0' ) { 10 if( c == '-' ) f = -1 ; 11 c = getchar( ); 12 } 13 while( c >= '0' && c <= '9' ) x = x * 10 + c - '0' , c = getchar( ) ; 14 return f == 1 ? x : -x ; 15 } 16 17 int ninth( int i , int j ) {//判断i行j列在第几个九宫格里，因为九个九宫格不重复，所以某些地方可以少判断一些条件。 18 if( i <= 3 && j <= 3 ) return 1 ; 19 if( i <= 3 && j <= 6 ) return 2 ; 20 if( i <= 3 ) return 3 ; 21 if( i <= 6 && j <= 3 ) return 4 ; 22 if( i <= 6 && j <= 6 ) return 5 ; 23 if( i <= 6 ) return 6 ; 24 if( j <= 3 ) return 7 ; 25 if( j <= 6 ) return 8 ; 26 return 9 ; 27 } 28 29 int point[ 6 ] = { 0 , 6 , 7 , 8 , 9 , 10} ;//存放分数; 30 int map[ 10 ][ 10 ] ;//记录某个位置上的数; 31 long long ans = -1 ;//无解则输出-1; 32 struct node { 33 int line , sum ; 34 }sss[ 11 ] ;//记录每行需要填的零的个数; 35 bool cmp( node i , node j ) { 36 return i.sum < j.sum ; 37 } 38 bool line[ 10 ][ 10 ] , nine[ 10 ][ 10 ] , list[ 10 ][ 10 ] ;//进行数独游戏的判断; 39 //为了方便观察，函数都扔下去; 40 int dfs( int , int , int ); 41 int print( ); 42 43 int main( ) 44 { 45 for( int i = 1 ; i <= 9 ; ++i ) { 46 int k = 0 ; 47 for( int j = 1 ; j <= 9 ; ++j ) { 48 map[ i ][ j ] = read( ) ; 49 if( !map[ i ][ j ] ) ++k; 50 line[ i ][ map[ i ][ j ] ] = 1 ; 51 nine[ ninth( i , j ) ][ map[ i ][ j ] ] = 1 ; 52 list[ j ][ map[ i ][ j ] ] = 1 ; 53 } 54 sss[ i ].sum = k , sss[ i ].line = i ; 55 } 56 sort( sss + 1 , sss + 10 , cmp ); 57 dfs( 1 , sss[ 1 ].line , 1 ) ; 58 printf( "%lld" , ans ); 59 return 0 ; 60 } 61 int dfs( int h , int x , int y ) { 62 if( h == 10 ) {//全部搜完了并成立，进行算分 63 print( ); 64 return 0 ; 65 } 66 if( y == 10 ) {//为避免特判过多而加的中转; 67 dfs( h + 1 , sss[ h + 1 ].line , 1 ); 68 return 0 ; 69 } 70 if( !map[ x ][ y ] ) { 71 for( int i = 1 ; i <= 9 ; ++i ) { 72 if( line[ x ][ i ] == 0 && list[ y ][ i ] == 0 && nine[ ninth( x , y ) ][ i ] == 0 ){ 73 line[ x ][ i ] = 1 , list[ y ][ i ] = 1 , nine[ ninth( x , y ) ][ i ] = 1 ; 74 map[ x ][ y ] = i ; 75 dfs( h , x , y + 1 ); 76 //记得回溯 77 map[ x ][ y ] = 0 ; 78 line[ x ][ i ] = 0 , list[ y ][ i ] = 0 , nine[ ninth( x , y ) ][ i ] = 0 ; 79 } 80 } 81 } 82 else dfs( h , x , y + 1 ); 83 } 84 int print( ) {//统计当前方案的分数 85 long long sum = 0 ; 86 for( int i = 1 ; i <= 9 ; ++i ) { 87 for( int j = 1 ; j <= 9 ; ++j ) { 88 sum += ( map[ i ][ j ] * point[ min( min( i , 10 - i ) , min ( j , 10 - j ) ) ] ); 89 //越靠近中心，x与y越接近5,否则远离5 90 //所以min( min( i , 10 - i ) , min ( j , 10 - j ) )与 91 //point数组搭配即可算出当前位置的分值 92 } 93 } 94 ans = max( ans , sum );//更新最大值 95 }
cx程序：(用时: 3361ms / 内存: 928KB)
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm> 5 using namespace std; 6 inline int max(int x,int y){return(x)<(y)?(y):(x);} 7 bool line[10][10],list[10][10],nine[10][10]; 8 int ans=0,tot=0,maxn=-1,pos=1,Line=0,a[10]; 9 //a数组记录每行0的个数，maxn记录最大值，tot记录总共要填的点数 10 struct node 11 {int w,list;}b[2][82];//list记录那行0的个数 12 inline void read(int &x) 13 { 14 x=0;int f=1; 15 char ch=getchar(); 16 while(ch<'0'||ch>'9') 17 {if(ch=='-') f=-1; ch=getchar();} 18 while(ch>='0'&&ch<='9') 19 {x=x*10+ch-'0';ch=getchar();} 20 x*=f; 21 } 22 inline void write(int x) 23 { 24 if(x<0){putchar('-');write(~x+1);} 25 else{if(x>9)write(x/10);putchar(x%10+'0');} 26 } 27 inline int ninth(int i,int j) 28 { 29 if(i<=3&&j<=3)return 1;if(i<=3&&j<=6)return 2;if(i<=3)return 3; 30 if(i<=6&&j<=3)return 4;if(i<=6&&j<=6)return 5;if(i<=6)return 6; 31 if(j<=3)return 7;if(j<=6)return 8;return 9; 32 } 33 inline int check(int x,int y) 34 { 35 if(x==1||y==1||x==9||y==9)return 6; 36 if(x==2||y==2||x==8||y==8)return 7; 37 if(x==3||y==3||x==7||y==7)return 8; 38 if(x==4||y==4||x==6||y==6)return 9; 39 if(x==5&&y==5)return 10; 40 } 41 void dfs(int x,int y) 42 { 43 if(tot==pos-1){maxn=max(ans,maxn);return;}//搜完了 44 for(register int i=1;i<=9;++i)//判断数字1-9 45 { 46 if(!line[x][i]&&!list[y][i]&&!nine[ninth(x,y)][i])//判断数字是否被填过 47 { 48 pos++;ans+=i*check(x,y); 49 line[x][i]=list[y][i]=1; 50 nine[ninth(x,y)][i]=1; 51 dfs(b[0][pos].w,b[1][pos].w); 52 pos--;ans-=i*check(x,y); 53 line[x][i]=list[y][i]=0; 54 nine[ninth(x,y)][i]=0;//回溯 55 } 56 } 57 } 58 bool cmp(node i,node j) 59 {return i.list<j.list;} 60 int main() 61 { 62 bool flag1=0,flag2=0; 63 for(register int i=1;i<=9;++i) 64 { 65 for(register int j=1;j<=9;++j) 66 { 67 register int k;read(k); 68 if(!k)b[0][tot+1].w=i,b[1][tot+1].w=j,tot++,Line++;//代表要填 69 else//代表填过 70 { 71 if(i==1||i==2)if(j==4)if(k==9||k==2)flag1=flag2=1; 72 if((i==3&&j==4&&k==1))flag1=1;if(i==4&&j==4&&k==9)flag2=1; 73 ans+=k*check(i,j);nine[ninth(i,j)][k]=1; 74 line[i][k]=1,list[j][k]=1; 75 } 76 } 77 a[i]=Line;Line=0;//记录每行0的个数 78 } 79 if(!flag1&&!flag2)//特判，为什么后面讲了 80 { 81 for(register int i=1;i<=tot;++i)b[0][i].list=b[1][i].list=a[b[0][i].w]; 82 std::sort(b[0]+1,b[0]+tot+1,cmp); 83 std::sort(b[1]+1,b[1]+tot+1,cmp);//按行中0的个数排序 84 } 85 dfs(b[0][pos].w,b[1][pos].w);//搜索 86 write(maxn); 87 }
因为我的搜法不同，在之后仍有两个点过不去，所以我特判了一下不
的情况，仔细反思一下，我好像懂了。
ly是一行行搜，所以他一直搜同一行，而我存的是点，又因为不稳定，我的程序可能搜完一个，又去搜另一行，所以不行，兴许可以改成
或桶排，相较之下，很明显还是ly程序更优。
