魔方阵算法及C语言实现
1 魔方阵概念
魔方阵是指由1,2,3……n2填充的,每一行、每一列、对角线之和均相等的方阵,阶数n = 3,4,5…。魔方阵也称为幻方阵。
例如三阶魔方阵为:
魔方阵有什么的规律呢?
魔方阵分为奇幻方和偶幻方。而偶幻方又分为是4的倍数(如4,8,12……)和不是4的倍数(如6,10,14……)两种。下面分别进行介绍。
2 奇魔方的算法
2.1 奇魔方的规律与算法
奇魔方(阶数n = 2 * m + 1,m =1,2,3……)规律如下:
- 数字1位于方阵中的第一行中间一列;
- 数字a(1 < a ≤ n2)所在行数比a-1行数少1,若a-1的行数为1,则a的行数为n;
- 数字a(1 < a ≤ n2)所在列数比a-1列数大1,若a-1的列数为n,则a的列数为1;
- 如果a-1是n的倍数,则a(1 < a ≤ n2)的行数比a-1行数大1,列数与a-1相同。
2.2 奇魔方算法的C语言实现
#include <stdio.h>
// Author: http://furzoom.com/
// N为魔方阶数
#define N 11 int main()
{
int a[N][N];
int i;
int col,row; col = (N-)/;
row = ; a[row][col] = ; for(i = ; i <= N*N; i++)
{
if((i-)%N == )
{
row++;
}
else
{
// if row = 0, then row = N-1, or row = row - 1
row--;
row = (row+N)%N; // if col = N, then col = 0, or col = col + 1
col ++;
col %= N;
}
a[row][col] = i;
}
for(row = ;row<N;row++)
{
for(col = ;col < N; col ++)
{
printf("%6d",a[row][col]);
}
printf("\n");
}
return ;
}
3 偶魔方的算法
偶魔方的情况比较特殊,分为阶数n = 4 * m(m =1,2,3……)的情况和阶数n = 4 * m + 2(m = 1,2,3……)情况两种。
3.1 阶数n = 4 * m(m =1,2,3……)的魔方(双偶魔方)
算法1:阶数n = 4 * m(m =1,2,3……)的偶魔方的规律如下:
- 按数字从小到大,即1,2,3……n2顺序对魔方阵从左到右,从上到下进行填充;
- 将魔方中间n/2列的元素上、下进行翻转;
- 将魔方中间n/2行的元素左、右进行翻转。
C语言实现
#include <stdio.h>
// Author: http://furzoom.com/
// N为魔方阶数,
#define N 12 int main()
{
int a[N][N];//存储魔方
int i, temp;//临时变量
int col, row;//col 列,row 行 //初始化
i = ;
for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
a[row][col] = i;
i++;
}
} //翻转中间列
for(row = ; row < N/; row ++)
{
for(col = N/;col < N/*;col ++)
{
temp = a[row][col];
a[row][col] = a[N-row-][col];
a[N-row-][col] = temp;
}
} //翻转中间行
for(col = ; col < N/; col ++)
{
for(row = N/;row < N/ * ;row ++)
{
temp = a[row][col];
a[row][col] = a[row][N-col-];
a[row][N-col-] = temp;
}
} for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
printf("%5d",a[row][col]);
}
printf("\n");
}
return ;
}
算法2:阶数n = 4 * m(m =1,2,3……)的偶魔方的规律如下:
- 按数字从小到大,即1,2,3……n2顺序对魔方阵从左到右,从上到下进行填充;
- 将魔方阵分成若干个4×4子方阵,将子方阵对角线上的元素取出;
- 将取出的元素按从大到小的顺序依次填充到n×n方阵的空缺处。
C语言实现
#include <stdio.h>
// Author: http://furzoom.com/
// N为魔方阶数
#define N 12 int main()
{
int a[N][N];//存储魔方
int temparray[N*N/];//存储取出的元素
int i;//循环变量
int col, row;// col 列,row 行 //初始化
i = ;
for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
a[row][col] = i;
i++;
}
}
//取出子方阵中对角线上的元素,且恰好按从小到大的顺序排放
i = ;
for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
if((col % == row % ) || ( == ( col % + row % )))
{
temparray[i] = a[row][col];
i++;
}
}
}
//将取出的元素按照从大到小的顺序填充到n×n方阵中
i = N*N/ -;
for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
if((col % == row % ) || ( == ( col % + row % )))
{
a[row][col] = temparray[i];
i--;
}
}
}
//输出方阵
for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
printf("%5d",a[row][col]);
}
printf("\n");
}
return ;
}
3.2 阶数n = 4 * m + 2(m =1,2,3……)的魔方(单偶魔方)
算法
设k = 2 * m + 1;单偶魔方是魔方中比较复杂的一个。
- 将魔方分成A、B、C、D四个k阶方阵,如下图
这四个方阵都为奇方阵,利用上面讲到的方法依次将A、D、B、C填充为奇魔方。
- 交换A、C魔方元素,对魔方的中间行,交换从中间列向右的m列各对应元素;对其他行,交换从左向右m列各对应元素。
- 交换B、D魔方元素,交换从中间列向左m – 1列各对应元素。
C语言实现
#include <stdio.h>
// Author: http://furzoom.com/
// N为魔方阶数
#define N 10 int main()
{
int a[N][N] = { {} };//存储魔方
int i,k,temp;
int col,row;// col 列,row 行 //初始化
k = N / ;
col = (k-)/;
row = ;
a[row][col] = ;
//生成奇魔方A
for(i = ; i <= k*k; i++)
{
if((i-)%k == )//前一个数是3的倍数
{
row++;
}
else
{
// if row = 0, then row = N-1, or row = row - 1
row--;
row = (row+k)%k; // if col = N, then col = 0, or col = col + 1
col ++;
col %= k;
}
a[row][col] = i;
} //根据A生成B、C、D魔方
for(row = ;row < k; row++)
{
for(col = ;col < k; col ++)
{
a[row+k][col+k] = a[row][col] + k*k;
a[row][col+k] = a[row][col] + *k*k;
a[row+k][col] = a[row][col] + *k*k;
}
} // Swap A and C
for(row = ;row < k;row++)
{
if(row == k / )//中间行,交换从中间列向右的m列,N = 2*(2m+1)
{
for(col = k / ; col < k - ; col++)
{
temp = a[row][col];
a[row][col] = a[row + k][col];
a[row + k][col] = temp;
}
}
else//其他行,交换从左向右m列,N = 2*(2m+1)
{
for(col = ;col < k / ;col++)
{
temp = a[row][col];
a[row][col] = a[row + k][col];
a[row + k][col] = temp;
}
}
} // Swap B and D
for(row = ; row < k;row++)//交换中间列向左m-1列,N = 2*(2m+1)
{
for(i = ;i < (k - )/ - ;i++)
{
temp = a[row][k+ k/ - i];
a[row][k+ k / -i] = a[row + k][k+k/ -i];
a[row + k][k+k/ -i] = temp;
}
} //输出魔方阵
for(row = ;row < N; row++)
{
for(col = ;col < N; col ++)
{
printf("%5d",a[row][col]);
}
printf("\n");
} return ;
}
==========================================》
你若要放弃,那我就可以安慰自己了,因为这样我不再孤独了,可以参考能够不断往前走,成功难道不是早晚的是事吗?
魔方阵算法及C语言实现的更多相关文章
- 任意阶魔方阵(幻方)的算法及C语言实现
写于2012.10: 本来这是谭浩强那本<C程序设计(第四版)>的一道课后习题,刚开始做得时候去网上找最优的算法,结果发现奇数和双偶数(4的倍数)的情况下算法都比较简单,但是单偶数(2的倍 ...
- C语言——打印魔方阵(每一行,每一列,对角线之和相等)
<一>魔方阵说明: 魔方阵是一个N*N的矩阵: 该矩阵每一行,每一列,对角线之和都相等: <二>魔方阵示例: 三阶魔方阵: 8 1 6 3 5 7 4 9 ...
- n阶魔方阵(奇数阵)的输出
需求 要求输出1~n²的自然数构成的魔方阵. STEP 1 什么是魔方阵? 魔方阵,古代又称“纵横图”,是指组成元素为自然数1.2…n2的平方的n×n的方阵,其中每个元素值都不相等,且每行.每列以及主 ...
- 一个UUID生成算法的C语言实现 --- WIN32版本 .
一个UUID生成算法的C语言实现——WIN32版本 cheungmine 2007-9-16 根据定义,UUID(Universally Unique IDentifier,也称GUID)在时 ...
- 无限大整数相加算法的C语言源代码
忙里偷闲,终于完成了无限大整数相加算法的C语言代码,无限大整数相加算法的算法分析在这里. 500位的加法运行1000次,不打印结果的情况下耗时0.036秒,打印结果的情况下耗时16.285秒. 下面是 ...
- 数据结构算法集---C++语言实现
//数据结构算法集---C++语言实现 //各种类都使用模版设计,可以对各种数据类型操作(整形,字符,浮点) /////////////////////////// // // // 堆栈数据结构 s ...
- 【2(2N+1)魔方阵 】
/* 2(2N+1)魔方阵 */ #include<stdio.h> #include<stdlib.h> #define N 6 #define SWAP(x, y) {in ...
- 1164: 零起点学算法71——C语言合法标识符(存在问题)
1164: 零起点学算法71——C语言合法标识符 Time Limit: 1 Sec Memory Limit: 64 MB 64bit IO Format: %lldSubmitted: 10 ...
- 【最全】经典排序算法(C语言)
算法复杂度比较: 算法分类 一.直接插入排序 一个插入排序是另一种简单排序,它的思路是:每次从未排好的序列中选出第一个元素插入到已排好的序列中. 它的算法步骤可以大致归纳如下: 从未排好的序列中拿出首 ...
随机推荐
- Nodejs与Net 和SQL 交互利器Edge.js
https://github.com/tjanczuk http://www.cnblogs.com/joylee/archive/2013/02/05/msnodesql.html edge.js这 ...
- cloud maintenance of OpenNebula
OpenNebula 4.4.1 maintenance release,官方建议当前的生产环境使用3.x or 4.x的其它版本; php调用curl工具伪造ip Upgrading from Op ...
- BestCoder Round #66 (div.2) hdu5592
GTW likes math Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others) ...
- [置顶] 很荣幸被选为2013年度 CSDN博客之星评选,如果觉得我的文章可以,请投我一票!
亲爱的小伙伴们,很荣幸我被选为<2013年度CSDN博客之星候选人>,希望大家多多支持,geekguy会继续努力,为大家奉献更好的文章. 投票地址:http://vote.blog.csd ...
- C#多线程(上) 分类: C# 线程 2015-03-09 10:35 174人阅读 评论(0) 收藏
一.多线程的相关概念 什么是进程? 当一个程序开始运行时,它就是一个进程,进程包括运行中的程序和程序所使用到的内存和系统资源. 而一个进程又是由多个线程所组成的. 什么是线程? 线程是程序中的一个执行 ...
- shell脚本的入参
shell脚本参数可以任意多,但只有前9个可以被访问,使用shift命令可以改变这个限制.参数从第一个开始,在第九个结束.$0 程序名字$n 第n个参数值,n=1..9 $* 所有命令行参数$@ ...
- 如何克隆路由器MAC地址,怎么操作?
路由器的“MAC地址克隆”的意思是: 不克隆时,从外网访问你的电脑,获得的MAC地址是路由器的mac地址. 克隆后,从外网访问你的电脑,获得的MAC地址是你电脑网卡的mac地址. 实用举例如下: 中国 ...
- [Selenium]中使用css选择器进行元素定位
参考:http://www.cnblogs.com/webblog/archive/2009/07/07/1518274.html 常见语法 * 通用元素选择器,匹配任何元素 E 标签选择器,匹配所有 ...
- HTML输出 二 控制行背景颜色
$Infors = Get-Content ports01.txt$Temp_PortStatustxt = "C:\Windows\Temp\PortStatustxt.txt" ...
- 【转】使用GDB调试Coredump文件
来自:http://blog.ddup.us/?p=176 写C/C++程序经常要直接和内存打交道,一不小心就会造成程序执行时产生Segment Fault而挂掉.一般这种情况都是因为数组越界访问,空 ...