奇数阶魔方阵(幻方阵)的生成,方法较多,Yidie教授于2007-6-15在算法设计栏目中,发表的“[原创]奇数幻方阵生成打印程序”,是属于右下方斜行法。我也转载过一个右上方斜行法的程序。这里将刚刚学会的控件数组的添加,用在奇数魔方阵生成程序上,加上个Timer,让它动态显示:如何右上方斜行布阵的。右上方斜行法布阵规则,参见算法设计中我转载的程序。敬请大侠高手给出改进意见!程序代码如如下: Dim n As Integer, i As Integer, j As Integer, k As Integer Private Sub Command2_Click() Dim m As Integer For m = 0 To n n - 1 Label2(m)Visible = False Next m k = 0: i = 0: j = n \ 2 Label2(k)Top = i (11000 / n) Label2(k)Left = j (12000 / n) Label2(k)Visible = True Timer1Enabled = True End Sub Private Sub Form_Load() List1Selected(0) = True Command2Enabled = False '布阵命令失效 End Sub Private Sub list1_click() '列表框中选择魔方阵的阶数,布阵命令失效 Command2Enabled = False End Sub Private Sub command1_click() '“生成”命令过程代码 For i = 1 To n n - 1 Unload Label2(i) Next i n = List1Text For i = 1 To n n - 1 Load Label2(i) Next i For i = 0 To n n - 1 With Label2(i) Top = (i \ n) (11000 / n) Left = (i Mod n) (12000 / n) + 100 Width = 11000 / n Height = 10000 / n Caption = i + 1 BackColor = RGB(255, 0, 0) FontSize = 360 / n Visible = True End With Next i Command2Enabled = True End Sub Private Sub timer1_timer() k = k + 1 If k < n n Then If k Mod n = 0 Then i = i + 1 Else i = IIf(i = 0, n - 1, i - 1) j = IIf(j = n - 1, 0, j + 1) End If Label2(k)Top = i (11000 / n) Label2(k)Left = j (12000 / n) + 100 Label2(k)Visible = True Else Timer1Enabled = fase End If End Sub Private Sub command3_click() Unload Me End Sub
求采纳
转的:
下面是两种方法的代码:
#include <stdioh>
void main()
{
int a[16][16];
int i,j,k,p,m,n;
/初始化/
p=1;
while(p==1)
{
printf("请输入n(0<n<=15,n是奇数)\n");
scanf("%d",&n);
if((n!=0)&&(n<=15)&&(n%2!=0))
{
printf("矩阵阶数是:%d\n",n);
p=0;
}
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
a[i][j]=0;
}
/建立魔方阵/
j=n/2+1;
a[1][j]=1;
for(k=2;k<=nn;k++)
{
i=i-1;
j=j+1;
if((i<1)&&(j>n))
{
i=i+2;
j=j-1;
}
else
{
if(i<1) i=n;
if(j>n) j=1;
}
if(a[i][j]==0)
a[i][j]=k;
else
{
i=i+2;
j=j-1;
a[i][j]=k;
}
}
/输出/
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
printf("%5d",a[i][j]);
printf("\n");
}
}
#include <stdioh>
int i,j,n,a[16][16]; /全局定义方阵行、列、阶数以及最大方阵容量/
void creat() /建立魔方阵/
{
int k;
j=n/2+1;
a[1][j]=1;
for(k=2;k<=nn;k++)
{
i=i-1;
j=j+1;
if((i<1)&&(j>n))
{
i=i+2;
j=j-1;
}
else
{
if(i<1) i=n;
if(j>n) j=1;
}
if(a[i][j]==0)
a[i][j]=k;
else
{
i=i+2;
j=j-1;
a[i][j]=k;
}
}
}
void print() /输出魔方阵/
{
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
printf("%5d",a[i][j]);
printf("\n");
}
}
void main()
{
int k,p,m;
while(1)
{
p=1;
while(p==1)
{
printf("请输入n(0<n<=15,n是奇数)\n");
scanf("%d",&n);
if((n!=0)&&(n<=15)&&(n%2!=0))
{
printf("矩阵阶数是:%d\n",n);
p=0;
}
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
a[i][j]=0;
}
creat();
print();
}
}
>
这个简单的很啊,你既然知道是magic函数可以实现这功能就好办了。
直接在matlab里输入命令open magic就可以打开magic函数了,里面自然就是代码咯 呵呵
如果你要的是C代码的话推荐你到图书馆看本书《C语言算法速查手册》里面有滴,下面的是MATLAB里的具体函数了
function M = magic(n)
%MAGIC Magic square
% MAGIC(N) is an N-by-N matrix constructed from the integers
% 1 through N^2 with equal row, column, and diagonal sums
% Produces valid magic squares for all N > 0 except N = 2
% Copyright 1984-2002 The MathWorks, Inc
% $Revision: 515 $ $Date: 2002/04/15 03:44:23 $
% Historically, MATLAB's magic was a built-in function
% This M-file uses a new algorithm to generate the same matrices
n = floor(real(double(n(1))));
% Odd order
if mod(n,2) == 1
[J,I] = meshgrid(1:n);
A = mod(I+J-(n+3)/2,n);
B = mod(I+2J-2,n);
M = nA + B + 1;
% Doubly even order
elseif mod(n,4) == 0
[J,I] = meshgrid(1:n);
K = fix(mod(I,4)/2) == fix(mod(J,4)/2);
M = reshape(1:nn,n,n)';
M(K) = nn+1 - M(K);
% Singly even order
else
p = n/2;
M = magic(p);
M = [M M+2p^2; M+3p^2 M+p^2];
if n == 2, return, end
i = (1:p)';
k = (n-2)/4;
j = [1:k (n-k+2):n];
M([i; i+p],j) = M([i+p; i],j);
i = k+1;
j = [1 i];
M([i; i+p],j) = M([i+p; i],j);
end
你的问题在于没有回退机制:
即,如果尝试的某一个数不符合要求,在试下一个数的时候,你没有把原先这个数的b[i]重置(而仍然是0)
整个思路上还是比较繁复,你这样写的代码很不清晰,建议改用递归来解决,你稍微查下递归的思路,参考几个代码就知道怎么用了。
对平面魔方的构造,分为三种情况:N为奇数、N为4的倍数、N为其它偶数(4n+2的形式)
⑴ N 为奇数时,最简单
(1) 将1放在第一行中间一列;
(2) 从2开始直到n×n止各数依次按下列规则存放:
按 45°方向行走,如向右上
每一个数存放的行比前一个数的行数减1,列数加1
(3) 如果行列范围超出矩阵范围,则回绕。
例如1在第1行,则2应放在最下一行,列数同样加1;
(4) 如果按上面规则确定的位置上已有数,或上一个数是第1行第n列时,
则把下一个数放在上一个数的下面。
⑵ N为4的倍数时
采用对称元素交换法。
首先把数1到n×n按从上至下,从左到右顺序填入矩阵
然后将方阵的所有4×4子方阵中的两对角线上位置的数关于方阵中心作对
称交换,即a(i,j)与a(n+1-i,n+1-j)交换,所有其它位置上的数不变。
(或者将对角线不变,其它位置对称交换也可)
⑶ N 为其它偶数时
当n为非4倍数的偶数(即4n+2形)时:首先把大方阵分解为4个奇数(2m+1阶)子方阵。
按上述奇数阶魔方给分解的4个子方阵对应赋值
上左子阵最小(i),下右子阵次小(i+v),下左子阵最大(i+3v),上右子阵次大(i+2v)
即4个子方阵对应元素相差v,其中v=nn/4
四个子矩阵由小到大排列方式为 ① ③
④ ②
然后作相应的元素交换:a(i,j)与a(i+u,j)在同一列做对应交换(j<t或j>n-t+2),
a(t-1,0)与a(t+u-1,0);a(t-1,t-1)与a(t+u-1,t-1)两对元素交换
其中u=n/2,t=(n+2)/4 上述交换使每行每列与两对角线上元素之和相等。
snjsj 我的程序算法:
这个魔方阵的算法可以对除2以外的任意阶数的方阵进行输出,结果保存在运行程序的目录下面的Magictxt文件中,用ie或者写字板打开以保持格式的一致(主要是回车符在记事本中为黑方框,认不出来)。当然具体的程序中,有内存空间以及变量范围的约束,我试过了,100以内的是可以的。
偶数阶的算法都是建立在奇数阶的基础之上,设方阵的阶数为n,则魔方阵常数(即每列每行以及对角线元素之和)为n(nn+1)/2。
请对照程序代码看,否则可能看不懂,可以一边看一边用笔对小阶的进行演算。
先说奇数阶的算法,这是最容易的算法:
n=2m+1,m为自然数
1)将数字1填在(0,(n+1)/2) ;要注意c中是从下标0开始
2)从左上往右下依次填。
3)由2),列的下标出界(超过n-1)时,行加1,以n为摸的余数为应填的列数;
4)由2),行的下标出界(超过n-1)时,列加1,以n为摸的余数为应填的行数;
5)由2),行列都未出界,但已添上其他数,应在当前位置左横移一个位置进行填数。
然后是偶数阶:
分两种情况,一种是n%4==2,一种是n%4==0
前一种:n=2(2m+1),m为自然数
1)将n阶方阵分为四个小魔方阵ABCD如下排列:
B C
D A
因为nn=4(2m+1)(2m+1),
记u=n/2=2m+1,分为1~uu,uu+1~2uu,2uu+1~3uu,3uu+1~4uu
即在调用子函数的时候分别如下面传递参数:
A(0),B(uu),C(2uu),D(3uu)
分别在ABCD中按照前面的填法把奇数阶填好(注意加上所传参数作为基数,每一个元素都要加上这个值),最后做如下交换:
(1)B中第0~(m-1)-1行中元素与C中相对应元素交换
(2)D中第(n-1)-m+1~(n-1)共m行的每行中的元素与A中相对应元素交换
(3)交换D:(u+m,m)与A中对应元素(矩阵中心值)
(4)交换D:(n-1,m)与A中对应元素(实际为矩阵最大值nn)
所谓对应位置,指相对于小魔方阵的左顶角的相对的行列位置
上面的这些你可以用数学进行证明,利用魔方阵常数(注意n阶的和u阶的关系)
后一种:n=4m,m为自然数
因为行列都是4的倍数,因而可以将整个矩阵分为每44的小矩阵。
先判断一个数是否在划为44小矩阵的对角线上,
如果在,则填该位置的数为nn-i+1(i为该元素的相对位置,从1开始,比如n阶的第s行第t个元素则其i=sn+t)
如果不在,则填上i。
很简单,初衷是这样:你while循环里面是输入奇数,退出循环,那么i=n+1是为了和后面的第一个for循环中i=i-1对应起来,把第二个数填在最下一行,假设你输入n=3,那么i=n+1的值就是4,那么执行第一个for循环后,i=i-1,i的值变成3……
中间的没有继续看下去了,因为逻辑不是很清晰……
只看最后的输出,貌似是定义了N x N的大数组,然后只根据n的值在其中部分二维数组中填写魔方阵,最后只输出部分,这样不是很清楚,给你一个也是初学者写的奇阶魔方阵程序示例吧……
#define N 3 //此处修改N值为其他奇数,即可求其他的奇阶魔方阵;
#include<stdioh>
int main()
{
int array[N][N],i,j,k; //定义array数组用于存放魔方阵,i和j用于循环,k是要填入的数;
for(i=0;i<N;i++) //array数组初始化,全赋值为0;
{
for(j=0;j<N;j++)
{
array[i][j]=0;
}
}
i=0;
j=N/2;
array[i][j]=1; //先将第一行中间一列填上1;
for(k=2;k<NN+1;k++) //开始填数,每次填一个,共循环NN-1次,刚好全部填满;
{
if((i==0)&&(j==N-1)) //规则1:若当前格子在右上顶角,则将下一个数填在当前格子正下方一行;
{
i=i+1;
}
else if(i==0) //规则2:若在最上一行(不包含右上顶角),则下一个数填在底行右侧;
{
i=N-1;
j=j+1;
}
else if(j==N-1) //规则3:若在最右一列(不包含右上顶角),则下一个数填在第一列上一行;
{
i=i-1;
j=0;
}
else //不是上述情况下,下一个数填在当前格子右上方;
{
i=i-1;
j=j+1;
}
if(array[i][j]==0) //若符合条件的格子里数值为0,即代表还未填数,则将k值填入;
{
array[i][j]=k;
}
else //若符合条件的格子已填入数,则填入当前格子的正下方一行;
{
i=i+2;
j=j-1;
array[i][j]=k;
}
}
for(i=0;i<N;i++) //输出魔方阵;
{
for(j=0;j<N;j++)
{
printf("%d",array[i][j]);
}
printf("\n");
}
return(0);
}
以上就是关于求制作VB奇数阶幻方阵的动画布阵方法和详细代码全部的内容,包括:求制作VB奇数阶幻方阵的动画布阵方法和详细代码、求程序代码打印魔方阵,所谓的魔方阵是nxn方针,其中n为奇数,它由nxn个正整数组成,它的每一行,每一列、用matlab做一个产生魔方阵的程序等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)