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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第廿八题：旋转矩阵对角线的和
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 1474  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +7 c297914 2026-04-05 8/400 2026-04-05 22:13 by hemengdong
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 306分材料与化工求调剂 +7 黎吧啦啦你很有� 2026-04-03 7/350 2026-04-05 17:18 by Hdyxbekcb
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 315求调剂 +13 小羊小羊_ 2026-04-02 14/700 2026-04-04 20:30 by 蓝云思雨
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 292求调剂 +11 2022080213 2026-04-04 13/650 2026-04-04 18:38 by macy2011
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 一志愿重庆大学085404，总分314分，求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 312求调剂 +6 小小墨123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 0817化工学硕调剂 +11 努力上岸中！ 2026-03-31 11/550 2026-04-01 20:30 by 赖春艳
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
信息提示
关闭
请填处理意见
关闭