应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第廿八题：旋转矩阵对角线的和
71/1返回列表
查看: 1439  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

holmescn

金虫 (正式写手)

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

这是旋转矩阵是这样生成的:  中心是1, 然后向右移动, 并按顺时针方向旋转, 这样生成的一个5x5的矩阵如下所示

21 22 23 24 25
20   7   8   9 10
19   6   1   2 11
18   5   4   3 12
17 16 15 14 13

其对角线上的和为101

如果这个把这个矩阵扩大到1001x1001,那它的对角元素的和是多少呢?

[ Last edited by holmescn on 2011-6-17 at 08:58 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-06-17 08:55:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★ ★
余泽成(金币+3): 辛苦了！ 2011-06-18 16:11:07
Matlab版:
CODE:
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);

结果: 669171001
赞一下(1人) 回复此楼
2楼2011-06-17 09:18:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎常来程序语言版讨论 2011-06-17 18:22:26
余泽成(金币+3): 鼓励讨论！ 2011-06-18 16:11:21
这个问题直接用数列算可以算出来。
这个蜗牛壳的每圈圈数是一个等差数列：
第一圈是1
第二圈是4*(2*2-1-1) = 8
第三圈是4*(2*3-1-1) = 16
第k圈是4*(2*k-1-1) = 8k-8
第k圈的边长是2k-1，第k圈的右上角那个数是前面k圈求和所得
第一圈是1
第二圈是1+8 = 9
第三圈是1+8+16=25
第k圈是4k^2-4k+1
根据边长可以前推3次，为了数学美，当然也可以从k-1圈后推3次，就可以依次求出四个角的通项公式：
第一圈是1
第二圈是1+2=3,1+2*2=5,1+2*3=7,1+2*4=9
第三圈是9+4=13,9+4*2=17,9+4*3=21,9+4*4=25
第k圈的边长是2k-1，步长是2k-2，计算可得：
右下角的：4k^2-10k+7
左下角的：4k^2-8k+5
左上角的：4k^2-6k+3
右上角的：4k^2-4k+1
进而直接按照通项公式可以求解。
赞一下(3人) 回复此楼
漩涡的中心有一块空地，空空的。
3楼2011-06-17 11:32:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎常来程序语言版讨论 2011-06-17 18:22:33
余泽成(金币+2): 鼓励交流！ 2011-06-18 16:11:33
matlab的,
奇数周长, 不计中间的1,每圈4个数从大到小是: len^2开始,间隔len-1
len^2
len^2-(len-1)
len^2-2*(len-1)
len^2-3*(len-1)
CODE:
function result = euler28()
tic;
result = 1;
% faster than using SUM(...)
for i=3:2:1001
    d = i*i;
    result = result+d;
    d = d+1-i;
    result = result+d;
    d = d+1-i;
    result = result+d;
    d = d+1-i;
    result = result+d;
end
toc;
end

结果
CODE:
% Elapsed time is 0.000014 seconds.
% ans =
%    669171001

[ Last edited by libralibra on 2011-6-19 at 16:55 ]
赞一下(3人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2011-06-17 15:04:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 鼓励交流，欢迎常来！ 2011-06-19 17:22:20
本帖仅楼主可见
赞一下
5楼2011-06-19 16:48:28
已阅   申请程序强帖   回复此楼   编辑   查看我的主页

holmescn

金虫 (正式写手)
引用回帖:
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

发代码要学会用BBCode喔
赞一下 回复此楼
6楼2011-06-19 17:00:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+2): 感谢参与 2011-06-20 11:59:52
n=1001;
switch mod(n,2)
case 1
result=(sum([1:n].^2)+(n-1)/2)*2-sum(2:2:n-1)*2-1
case 0
result=(sum([1:n].^2)+n/2)*2-sum(1:2:n-1)*2
end

result2=(sum([1:n].^2)+(n-mod(n,2))/2)*2-sum(1+mod(n,2):2:n-1)*2-mod(n,2)
赞一下(2人) 回复此楼
7楼2011-06-20 08:35:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 holmescn 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分，求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 304求调剂（085602一志愿985） +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 研究生招生 +3 徐海涛11 2026-03-10 7/350 2026-03-12 14:26 by 徐海涛11
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
信息提示
关闭
请填处理意见
关闭