任意阶幻方(魔方矩阵)C语言实现
------------------来自百度
先标出引用地址:
http://blog.ddedu.com.cn/user1/88/archives/2007/2007420143329.shtml //任意阶幻方构造方法
http://blog.csdn.net/cmutoo/article/details/5487157 //任意阶幻方C语言代码实现(有些许错误)
基础知识这里看:http://blog.csdn.net/oowgsoo/article/details/1567910
任意阶幻方的构造方法有很多种,所以要选定一种易于代码实现的一种
在上篇博客中说道:
/************************************************************************************
幻方的数量:
与我们大多数人的常识不同,幻方的数量不是唯一的,而且也不是一个简单的问题
3阶幻方只有1种
4阶幻方有880种,通过旋转和反射,总共可以有7040个幻方
5阶幻方有275 305 224个,这是用计算机算的
6阶幻方,大概是1.7743*10**19~1.7766*10**19之间,这是用统计学方法计算的,居然计算机也算不出来,更不要说6阶以上的幻方数量了
************************************************************************************/
所以代码实现的就有很大的局限性,只能实现某一种构造方法的幻方
幻方构造分为
1、奇数阶
2、双偶阶
3、单偶阶
三种。
对于奇数阶的幻方:
/*******************************************************************
把1(或最小的数)放在第一行正中; 按以下规律排列剩下的n×n-1个数:
(1)每一个数放在前一个数的右上一格;
(2)如果这个数所要放的格已经超出了顶行那么就把它放在底行,仍然要放在右一列;
(3)如果这个数所要放的格已经超出了最右列那么就把它放在最左列,仍然要放在上一行;
(4)如果这个数所要放的格已经超出了顶行且超出了最右列,那么就把它放在前一个数的下一行同一列的格内;
(5)如果这个数所要放的格已经有数填入,处理方法同(4)。
这种写法总是先向“右上”的方向,象是在爬楼梯。
void Odd(int n,int index)
{
while(st != index)
{
cube[ox+stx][oy+sty] = st++; if((st-) % n != )
{
stx--;
sty++;
}
else
{
stx++;
} stx = ((stx-)%n+n)%n+;
sty = (sty%n == ? n : sty%n); }
}
对于双偶数阶的
就是一个中心对称
void DouEven(int n)
{
int i,j;
int st = ;
for(i=; i<=n; i++)
{
for(j=; j<=n; j++)
{
cube[i][j] = st++;
}
} int zx,zy,fx,fy;
for(i=; i<=n*n; i+=)
{ for(j=; j<=n*n; j+=)
{
zx=i-,zy=j-,fx=i,fy=j-;
cube[zx][zy]=n*n-cube[zx][zy]+;
cube[zx+][zy+]=n*n-cube[zx+][zy+]+;
cube[zx+][zy+]=n*n-cube[zx+][zy+]+;
cube[zx+][zy+]=n*n-cube[zx+][zy+]+; cube[fx][fy]=n*n-cube[fx][fy]+;
cube[fx-][fy+]=n*n-cube[fx-][fy+]+;
cube[fx-][fy+]=n*n-cube[fx-][fy+]+;
cube[fx-][fy+]=n*n-cube[fx-][fy+]+;
}
}
}
对于单偶数阶的,麻烦许多
因为要用到奇数阶的构造方法
void SingleEven(int n)
{
int i,j; ox=oy=;st=;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //A ox=oy=n/;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //D ox=,oy=n/;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //B ox=n/,oy=;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //C int k=(n-)/,tmp;
for(j=; j<=n/; j++)
{
for(i=; i<=k; i++)
{
if(j != (n/+)/)
{
tmp = cube[j][i];
cube[j][i] = cube[j+n/][i];
cube[j+n/][i] = tmp;
}
else
{
tmp = cube[j][i+(n/+)/-];
cube[j][i+(n/+)/-] = cube[j+n/][i+(n/+)/-];
cube[j+n/][i+(n/+)/-] = tmp;
}
}
}
if(k-)
{
for(i=; i<=n/; i++)
{
int tmpp = (*n+)/-;
for(j=; j<=k-; j++)
{
tmp = cube[i][j+tmpp];
cube[i][j+tmpp] = cube[i+n/][j+tmpp];
cube[i+n/][j+tmpp] = tmp; }
}
}
}
最后贴一个完整的代码:
#include <stdio.h>
#include <string.h> int cube[][];
int stx,sty;
int st;
int num;
int ox,oy; void Odd(int n,int index)
{
while(st != index)
{
cube[ox+stx][oy+sty] = st++; if((st-) % n != )
{
stx--;
sty++;
}
else
{
stx++;
} stx = ((stx-)%n+n)%n+;
sty = (sty%n == ? n : sty%n); }
} void DouEven(int n)
{
int i,j;
int st = ;
for(i=; i<=n; i++)
{
for(j=; j<=n; j++)
{
cube[i][j] = st++;
}
} int zx,zy,fx,fy;
for(i=; i<=n*n; i+=)
{ for(j=; j<=n*n; j+=)
{
zx=i-,zy=j-,fx=i,fy=j-;
cube[zx][zy]=n*n-cube[zx][zy]+;
cube[zx+][zy+]=n*n-cube[zx+][zy+]+;
cube[zx+][zy+]=n*n-cube[zx+][zy+]+;
cube[zx+][zy+]=n*n-cube[zx+][zy+]+; cube[fx][fy]=n*n-cube[fx][fy]+;
cube[fx-][fy+]=n*n-cube[fx-][fy+]+;
cube[fx-][fy+]=n*n-cube[fx-][fy+]+;
cube[fx-][fy+]=n*n-cube[fx-][fy+]+;
}
}
} void SingleEven(int n)
{
int i,j; ox=oy=;st=;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //A ox=oy=n/;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //D ox=,oy=n/;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //B ox=n/,oy=;
stx=,sty=(n/+)/;
Odd(n/,n*n*/+); //C int k=(n-)/,tmp;
for(j=; j<=n/; j++)
{
for(i=; i<=k; i++)
{
if(j != (n/+)/)
{
tmp = cube[j][i];
cube[j][i] = cube[j+n/][i];
cube[j+n/][i] = tmp;
}
else
{
tmp = cube[j][i+(n/+)/-];
cube[j][i+(n/+)/-] = cube[j+n/][i+(n/+)/-];
cube[j+n/][i+(n/+)/-] = tmp;
}
}
}
if(k-)
{
for(i=; i<=n/; i++)
{
int tmpp = (*n+)/-;
for(j=; j<=k-; j++)
{
tmp = cube[i][j+tmpp];
cube[i][j+tmpp] = cube[i+n/][j+tmpp];
cube[i+n/][j+tmpp] = tmp; }
}
}
} void print(int n)
{
int i,j; for(i=; i<=n; i++)
{
int sum = ;
for(j=; j<=n ; j++)
{
sum += cube[i][j];
printf("%4d",cube[i][j]);
}
printf(" sum = %d",sum);
printf("\n");
}
} int main()
{ int i,j,k,t,n,m;
do
{
printf("Please Input n(3-17): ");
scanf("%d",&n);
if(n<)continue; memset(cube,,sizeof(cube));
if(n % == )
{
stx=,sty=(n+)/;
ox=oy=;
st=;
Odd(n,n*n+);
print(n);
}
else if(n % == )
{
DouEven(n);
print(n);
}
else if(n % == && n % != )
{ SingleEven(n);
print(n);
}
}while(); return ;
}
再上一个专门关于介绍幻方的博客:http://blog.sina.com.cn/u/1225071715
任意阶幻方(魔方矩阵)C语言实现的更多相关文章
- 任意阶魔方阵(幻方)的算法及C语言实现
写于2012.10: 本来这是谭浩强那本<C程序设计(第四版)>的一道课后习题,刚开始做得时候去网上找最优的算法,结果发现奇数和双偶数(4的倍数)的情况下算法都比较简单,但是单偶数(2的倍 ...
- Java 实现任意N阶幻方的构造
一.关于单偶数阶幻方和双偶数阶幻方 (一)单偶数阶幻方(即当n=4k+2时) 任何4k+2 阶幻方都可由2k+1阶幻方与2×2方块复合而成,6是此类型的最小阶. 以6阶为例,可由3阶幻方与由0,1,2 ...
- 【C++小白成长撸】--(续)单偶数N阶魔方矩阵
1 /*程序的版权和版本声明部分: **Copyright(c) 2016,电子科技大学本科生 **All rights reserved. **文件名:单偶数N阶魔方矩阵 **程序作用:单偶数N阶魔 ...
- 【C++小白成长撸】--N阶幻方(魔阵)矩阵
解决方法:1.第一个元素放在第一行中间一列 2.下一个元素存放在当前元素的上一行.下一列. 3.如果上一行.下一列已经有内容,则下一个元素的存放位置为当前列的下一行. 在找上一行.下一行或者下一列的时 ...
- Java 实现奇数阶幻方的构造
一.设计的流程图如下所示 二.Java 语言的代码实现 package MagicSquare; //奇数幻方的实现 public class Magic_Odd { //n 为幻方的阶数 publi ...
- Codeforces 710C. Magic Odd Square n阶幻方
C. Magic Odd Square time limit per test:1 second memory limit per test:256 megabytes input:standard ...
- n阶幻方
前序 最近在学习一些经典的算法,搞得头昏脑涨,就想换换脑子.在家里的旧书堆里面乱翻,无意中将一本具有十多年历史的小学数学奥林匹克竞赛的书发掘了出来,能放到现在挺不容易的,就拿起来随便翻翻.看了看目录, ...
- codeforces 710C Magic Odd Square(构造或者n阶幻方)
Find an n × n matrix with different numbers from 1 to n2, so the sum in each row, column and both ma ...
- n阶幻方问题
转载自:http://blog.csdn.net/fengchaokobe/article/details/7437767 目录 第一节 n阶幻方问题 第二节 由n阶幻方引发 ...
随机推荐
- IOS第八天(3:UITableViewController团购, 点击底部代码调整)
****代理者的方法中 // 通知页脚视图调整视图显示状态 [footerView endRefresh]; //发送代理通知的类中 /** 视图控制器刷新完成调用方法 */ - (void)endR ...
- ServletContext读取Web应用中的资源文件
package cn.itcast; import java.io.FileInputStream; import java.io.IOException; import java.io.InputS ...
- String类StringBuffer类与StringBuilder类gc垃圾回收
String类的特点 直接赋值和new调用构造方法两种, 直接赋值时会将字符串常量入内存池,当其他变量再赋相同值时,不再在堆空间开辟内存 new构造方法会开辟两块堆内存空间,可以使用intern手工入 ...
- Codevs 1022 覆盖
时间限制: 1 s 空间限制: 128000 KB 题目等级 : 大师 Master 题目描述 Description 有一个N×M的单位方格中,其中有些方格是水塘,其他方格是陆地.如果要用1×2的矩 ...
- DNS视图以及日志压力测试
1 访问控制列表 配置在/etc/named.conf文件的最顶端 acl innct { 192.168.1.0/24; 127.0.0.0/8; }; ...
- Ruby 类案例
#!/user/bin/ruby # -*-coding:UTF-8-*- class Customer @@no_of_customers=0 def initialize(id,name,addr ...
- [BS-31]导航控制器的interactivePopGestureRecognizer属性
导航控制器的interactivePopGestureRecognizer属性 如何自定义导航控制器push出来的视图控制器的左侧返回按钮? 首先需要知道的是:如果要自定义push出来的VC的左侧返回 ...
- tcpdump教程入门
tcpdump是一个最基本重要的网络分析工具, 掌握好这, 对于学习tcp/ip协议也是很有帮助的. 理解了tcp/ip协议栈的知识, 分析调优网络的能力才会更高. 所以使用tcpdump相比其它的工 ...
- SpringMVC源码剖析(二)- DispatcherServlet的前世今生
上一篇文章<SpringMVC源码剖析(一)- 从抽象和接口说起>中,我介绍了一次典型的SpringMVC请求处理过程中,相继粉墨登场的各种核心类和接口.我刻意忽略了源码中的处理细节,只列 ...
- Apache Commons 工具集
一.Commons BeanUtils http://jakarta.apache.org/commons/beanutils/index.html 说明:针对Bean的一个工具集.由于Bean往往是 ...