2021-02-09 洛谷 P2615 神奇的幻方（NOIP 2015 提高组 d1t1）题解

洛谷 P2615 神奇的幻方（NOIP 2015 提高组 d1t1）题解

题目描述

幻方是一种很神奇的 N∗N 矩阵：它由数字1,2,3,⋯⋯,N×N 构成，且每行、每列及两条对角线上的数字之和都相同。

当 N 为奇数时，我们可以通过下方法构建一个幻方：

首先将 1 写在第一行的中间。

之后，按如下方式从小到大依次填写每个数 K(K=2,3,⋯,N×N) ：

若 (K-1) 在第一行但不在最后一列，则将 K 填在最后一行， (K−1) 所在列的右一列；
若 (K-1) 在最后一列但不在第一行，则将 K 填在第一列， (K−1) 所在行的上一行；
若 (K−1) 在第一行最后一列，则将 K 填在 (K−1) 的正下方；
若 (K−1) 既不在第一行，也不在最后一列，如果 (K−1) 的右上方还未填数，则将 K 填在 (K-1) 的右上方，否则将 K 填在 (K−1) 的正下方。
现给定 N ，请按上述方法构造 N×N 的幻方。

输入格式

一个正整数 N ，即幻方的大小。

输出格式

共 N 行 ，每行 N 个整数，即按上述方法构造出的 N×N 的幻方，相邻两个整数之间用单空格隔开。

输入输出样例

Sample Input:

3

Sample Output:

8 1 6
3 5 7
4 9 2

说明/提示

对于 100% 的数据，对于全部数据， 1 ≤ N ≤ 39 且 N 为奇数。
（NOIp2015 提高组 d1t1）

这道题拿到一看，数据规模很小，甚至可以打表。 但是要是打表的话就太暴力，体现不出这道题的价值了

看一眼洛谷的题解，还是用模拟的比较多。模拟算是这道题比较正规（也是题面引导你用的）方法。

首先开一组循环输入。然后将第一行中间的数赋值为1，定义一个x,y表示当前位置，把x赋值为1，y赋值为(n+1)/2。

其次再开一组循环，根据x,y的值（即K-1的坐标)情况判断怎么填写数字，在填写数字完毕后将x,y的值赋为当前坐标。

最后当循环数到达n*n后结束并输出~

轻松AC~

下面就是模拟的代码：

#include <iostream>
#include <cstdio>
using namespace std;
int n,a[10010][10010],x=1,y,cnt=2;//cnt赋值为2，因为从2开始填数
int main()
{
	scanf("%d",&n);//输入不解释
	y=(n+1)/2;//将x,y的值（k-1)坐标赋值为1,(n+1)/2;
	a[x][y]=1;//将第一行最中间的数赋值为1.
	int now=n*n;
	for(int i=2;i<=now;++i)
	{
	    if(x==1&&y!=n)//模拟第一种情况
	    {
	    	a[n][y+1]=cnt;
	    	x=n;
	    	y+=1;
	    	cnt++;
		}
		else if(x!=1&&y==n)//注意有else模拟第二种情况
		{
			a[x-1][1]=cnt;
			x-=1;
			y=1;
			cnt++;
		}
		else if(x==1&&y==n)//模拟第三种情况
		{
			a[x+1][y]=cnt;
			x+=1;
			cnt++;
		}
		else //第四种情况
		{
			if(a[x-1][y+1]==0)//如果k-1右上方还未填数
			{
				a[x-1][y+1]=cnt;
				x-=1;
				y+=1;
				cnt++;
			}
			else//反之~
			{
				a[x+1][y]=cnt;
				x+=1;
				cnt++;
			}
		}
	}
	for(int i=1;i<=n;++i)//完美输出
	{
		for(int e=1;e<=n;++e)
		printf("%d ",a[i][e]);
		printf("\n");
	}
    return 0;//好习惯；
}

其实这种模拟还是很费时费力（if 打多了容易晕），凡事多想一步，能不能再进行一点优化（主要是让自己舒服）呢？

还真有，这道题可以从纯数学的角度考虑，很好地化繁为简（也是逛洛谷题解的时候看见的）

下面照搬大佬的题解：

其实这题完全没必要很多的if语句，把幻方扩展开来，

题目的意思就是从第一行中点开始，把每一个下面的数

放到“右上角”，若右上角有数，则放到正下方。

（有点像NOIp初赛的一道题）

#include<cstdio>
using namespace std;
int n,a[40][40],x,y;
int main(){
	scanf("%d",&n);
	x=1,y=(n+1)/2;
	for(int i=1;i<=n*n;i++){
		a[x][y]=i;
		if(!a[(x-2+n)%n+1][y%n+1]) x=(x-2+n)%n+1,y=y%n+1;
		else x=x%n+1;//数学运算
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			printf("%d ",a[i][j]);
		}
		printf("\n");
	}
} 

这样代码也更短， 简洁明了， 又是一个数学推导与思维碰撞的的神奇产物