SDUTRound #4 - 2018新春大作战SDUT [4153]A 某电气鼠的幻方【n阶幻方的构建板子】

某电气鼠的幻方

Time Limit: 1000 ms Memory Limit: 65536 KiB

Submit Statistic Discuss

Problem Description

作为日服第一老鼠，皮卡丘还是有必要出来给大家拜个年滴。顺便说一下，大家去KFC的时候别忘了去尝一下皮卡丘套餐哦。

皮卡丘现在抛给你一个难题，可要接好了。皮卡丘喜欢奇数，又喜欢幻方，请你给他构造一个他喜欢的奇数阶幻方哦。

什么！你没听说过幻方？！

幻方是一个n*n的方阵，且是一个满足每行之和，每列之和还有对角线之和均相等的方阵哦。

Input

先输入一个整数t，代表数据组数。

每组数据输入一个奇数n（1<=n<=365），代表幻方阶数。

Output

对于每组数据，用1至n*n这n*n个正整数构造n阶幻方并输出，数字之间用空格隔开。

Sample Input

2
1
3

Sample Output

1
4 9 2
3 5 7
8 1 6

Hint

Source

xuanhuang

题意： 就是让你构建一个 n*n 矩阵，使得矩阵的横竖斜相加都相同，n为奇数

分析： 这个很有规律，就当做板子吧，很全面，奇偶的都有 QAQ，具体请看这里.

参考代码

#include <bits/stdc++.h>

using namespace std;
#define maxn 500

int m[maxn][maxn];

void OddMagic(int n)                         //奇数阶幻方
{
    memset(m,0,sizeof(m));                   //幻方清零
    int x=0, y=n/2;
    for(int i=1; i<=n*n; i++)
    {
        m[x][y]=i;
        x--;                                 //依次沿右上角填充
        y++;
        if(x<0 && y>n-1) {x=x+2; y=y-1;}     //沿对角线超出
        else if(x<0) x=x+n;                  //沿上边界超出
        else if(y>n-1) y=y-n;                //沿右边界超出
        else if(m[x][y]!=0) {x=x+2; y=y-1;}  //右上角已填充
    }
}

void DoubleEvenMagic(int n)                  //双偶数阶幻方
{
    memset(m,0,sizeof(m));                   //幻方清零
    for(int i=1, x=0, y=0; i<=n*n; i++)      //依次按顺序赋初值
    {
        m[x][y]=i;
        y++;
        if(y>n-1) {x++; y-=n;}
    }
    for(int i=0; i<n; i++)                   //将幻方分解成m*m个4阶幻方，并将每个4阶幻方的对角线元素换成其互补数
        for(int j=0; j<n; j++)
            if(i%4==0 && j%4==0)             //左对角线
                for(int k=0; k<4; k++)
                    m[i+k][j+k]=(n*n+1)-m[i+k][j+k];
            else if(i%4==3 &&j%4==0)         //右对角线
                for(int k=0; k<4; k++)
                    m[i-k][j+k]=(n*n+1)-m[i-k][j+k];
}

void SingleEvenMagic(int n)                  //单偶数阶幻方
{
    memset(m,0,sizeof(m));                   //幻方清零
    int n0=n/2;
    OddMagic(n0);                            //将幻方分解成2*2个奇数阶幻方，调用奇数阶幻方函数填充左上角奇数阶幻方
    for(int i=0; i<n0; i++)
        for(int j=0; j<n0; j++)
        {
            m[i+n0][j+n0]=m[i][j]+n0*n0;       //填充右下角奇数阶幻方
            m[i][j+n0]=m[i+n0][j+n0]+n0*n0;     //填充右上角奇数阶幻方
            m[i+n0][j]=m[i][j+n0]+n0*n0;       //填充左下角奇数阶幻方
        }
    int k=(n-2)/4;                             //满足公式n=4*k+2
    for(int i=0; i<n0; i++)
        for(int j=0; j<k; j++)
            if(i==n0/2) swap(m[i][i+j],m[i+n0][i+j]);   //将左上角幻方最中间元素从左向右的k个元素与左下角幻方相应位置元素交换
            else swap(m[i][j],m[i+n0][j]);              //将左上角幻方除最中间行外的每行前k个元素与左下角幻方相应位置元素交换
    for(int i=0; i<n0; i++)
        for(int j=n0+n0/2; j>n0+n0/2-(k-1); j--)
            swap(m[i][j],m[i+n0][j]);                   //将右上角幻方最中间列起从右向左的k-1列元素与右下角幻方相应位置元素交换
}

bool Check(int n)
{
    int cnt=n*(n*n+1)/2;                                //每行每列以及对角线之和
    for(int i=0; i<n; i++)
    {
        int sum_row=0,sum_line=0;
        for(int j=0; j<n; j++)
            {
                sum_row+=m[i][j];                       //检查各行
                sum_line+=m[j][i];                      //检查各列
            }
        if(sum_row!=cnt || sum_line!=cnt) return false;
    }
    int sum_left=0,sum_right=0;
    for(int i=0; i<n; i++)
    {
        sum_left+=m[i][i];                              //检查左对角线
        sum_right+=m[n-i-1][i];                         //检查右对角线
    }
    if(sum_left!=cnt || sum_right!=cnt) return false;
    return true;
}

void print(int n)                                       //按格式输出
{
    for(int i = 0;i < n;i++) {
        for(int j = 0;j < n;j++) {
            printf("%d%c",m[i][j],j == n - 1?'\n':' ');
        }
    }
}

int main()
{
   // ios::sync_with_stdio(false);
    int n;
    int T;cin>>T;
    while(T--)
    {
        cin>>n;
        if(n == 1) {
            cout<<1<<endl;continue;
        }
        if(n&1) {OddMagic(n); if(Check(n)) print(n);}                    //输出奇数阶幻方
        else if(!(n%4)) {DoubleEvenMagic(n); if(Check(n)) print(n);}     //输出双偶数阶幻方
        else {SingleEvenMagic(n); if(Check(n)) print(n);}                //输出单偶数阶幻方
    }
    return 0;
}

• 如有错误或遗漏，请私聊下UP，thx