c语言怎么编辑幻方程序,幻方的程序

满意答案

cohqxo

2014.02.08

采纳率：57%    等级：12

已帮助：15822人

function M = magic(n)

n=input('input N please,n=')

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+2*J-2,n);

M = n*A + 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:n*n,n,n)';

M(K) = n*n+1 - M(K);

% Singly even order.

else

p = n/2;

M = magic(p);

M = [M M+2*p^2; M+3*p^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

00分享举报