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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第廿八题：旋转矩阵对角线的和
51/1返回列表
查看: 1714  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 17:34 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 14:03 by GkCVyODlx4mR
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[基金申请] 国社科又开始会评了，不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
信息提示
关闭
请填处理意见
关闭