[交流] Euler 工程 第廿八题：旋转矩阵对角线的和已有4人参与

l = 2;
s = 1;
b = 3;
e = 9;
while 2*l-1 <= 1001
    s = s + sum(linspace(b, e, 4));
    b = e + 2*l;
    e = b + 2*l*3;
    l = l + 1;
end
fprintf('sum=%d\n', s);

Originally posted by fan6cy at 2011-06-19 16:48:28:
sum =

   669171001

Elapsed time is 0.000000 seconds.



clear,clc
tic
N0=1001;
i=3;sum=1;
while i<=N0
    sum=sum+i*(4*i-6)+6;
    i=i+2;
end
sum
toc