描述
小Hi最近在教邻居家的小朋友小学奥数,而最近正好讲述到了三阶幻方这个部分,三阶幻方指的是将1~9不重复的填入一个3*3的矩阵当中,使得每一行、每一列和每一条对角线的和都是相同的。
三阶幻方又被称作九宫格,在小学奥数里有一句非常有名的口诀:“二四为肩,六八为足,左三右七,戴九履一,五居其中”,通过这样的一句口诀就能够非常完美的构造出一个九宫格来。
有意思的是,所有的三阶幻方,都可以通过这样一个九宫格进行若干镜像和旋转操作之后得到。现在小Hi准备将一个三阶幻方(不一定是上图中的那个)中的一些数组抹掉,交给邻居家的小朋友来进行还原,并且希望她能够判断出究竟是不是只有一组解。
而你呢,也被小Hi交付了同样的任务,但是不同的是,你需要写一个程序~
输入
输入仅包含单组测试数据。
每组测试数据为一个3*3的矩阵,其中为0的部分表示被小Hi抹去的部分。
对于100%的数据,满足给出的矩阵至少能还原出一组可行的三阶幻方。
输出
如果仅能还原出一组可行的三阶幻方,则将其输出,否则输出“Too Many”(不包含引号)。
0 7 2 0 5 0 0 3 0样例输出
6 7 2 1 5 9 8 3 4
题中讲到,“所有的三阶幻方,都可以通过这样一个九宫格进行若干镜像和旋转操作之后得到。”所以说,由题中所给幻方翻转旋转所得八种幻方即为所有满足要求的情况,经如下程序简单测试确实如此,了解这一点之后,问题就更加简单了。
#include<iostream>
#include<algorithm>
using namespace std;
int chk(int n[])
{
int f;
f=n[0]+n[4]+n[8];
if(n[2]+n[4]+n[6]-f) return 0;
if(n[0]+n[1]+n[2]-f) return 0;
if(n[3]+n[4]+n[5]-f) return 0;
if(n[6]+n[7]+n[8]-f) return 0;
if(n[0]+n[3]+n[6]-f) return 0;
if(n[1]+n[4]+n[7]-f) return 0;
if(n[2]+n[5]+n[8]-f) return 0;
return 1;
}
int main()
{
int i,k=1,a[9]={1,2,3,4,5,6,7,8,9};
while(next_permutation(a,a+9))
{
if(chk(a))
{
cout<<"num: "<<k++<<endl;
for(i=0;i<9;i++)
{
cout<<a[i]<<' ';
if((i+1)%3==0) cout<<endl;
}
}
}
return 0;
}
测试结果如下:
num: 1
2 7 6
9 5 1
4 3 8
num: 2
2 9 4
7 5 3
6 1 8
num: 3
4 3 8
9 5 1
2 7 6
num: 4
4 9 2
3 5 7
8 1 6
num: 5
6 1 8
7 5 3
2 9 4
num: 6
6 7 2
1 5 9
8 3 4
num: 7
8 1 6
3 5 7
4 9 2
num: 7
8 3 4
1 5 9
6 7 2
与预期相符。现在来着手问题,首先,什么情况下幻方唯一?
我们可以将题中给的幻方样例看成一个板子,这面板子可以自由翻转旋转形成八种不同幻方,而所知的条件将成为这面板子的定点,使它翻转旋转不再那么自由,当它不能再转动时,幻方就可以唯一确定了。
例如:首先可以确定 5 永远在板子中央,这是一个定点,板子可以以 5 为定点自由翻转旋转; 若已知 7 在 5 的下方,那板子可以以 5 7 两定点所在定轴自由翻转,此时仍有两种幻方。如果再知道该轴之外的一个点,板子就固定不动了。
所以,除 5 之外,只需要知道 横,竖,两对角线 这四条轴中的两条轴上的任意一个数字,就能确定幻方,此时将上面程序输出结果匹配输出就可以了,代码如下:
#include <iostream>
#include <cstring>
using namespace std;
int a[3][3];
void print()
{
int i,j,k;
int n[8][3][3]={2,7,6,9,5,1,4,3,8,
2,9,4,7,5,3,6,1,8,
4,3,8,9,5,1,2,7,6,
4,9,2,3,5,7,8,1,6,
6,1,8,7,5,3,2,9,4,
6,7,2,1,5,9,8,3,4,
8,1,6,3,5,7,4,9,2,
8,3,4,1,5,9,6,7,2};
for(k=0;k<8;k++)
{
for(i=0;i<3;i++)
for(j=0;j<3;j++)
if(a[i][j]&&a[i][j]!=n[k][i][j])
goto END;
for(i=0;i<3;i++)
{
for(j=0;j<3;j++)
cout<<n[k][i][j]<<' ';
cout<<endl;
}
break;
END:;
}
}
int main()
{
int i,j,cnt=0;
for(i=0;i<3;i++)
for(j=0;j<3;j++)
cin>>a[i][j];
if(a[0][0]||a[2][2]) cnt++;
if(a[1][0]||a[1][2]) cnt++;
if(a[2][0]||a[0][2]) cnt++;
if(a[2][1]||a[0][1]) cnt++;
if(cnt>=2) print();
else cout<<"Too Many"<<endl;
return 0;
}