24小时热门版块排行榜    

查看: 1659  |  回复: 1

Simewe

新虫 (初入文坛)

[求助] Matlab中自定义的函数怎么能支持向量化操作?

CODE:
function r=primeQ(n)
r=1;
for i=2:sqrt(n);
    if mod(n,i)==0;
        r=0;
        return
    end
end;

上面的primeQ是个判断素数的函数,调用的时候只会这样用
sum( arrayfun(@primeQ, 2:1e6) )
听说arraygfun的速度不快,我希望可以支持向量化操作,就像下面这样
sum(primeQ(2:1e5))
应该怎么做?
注意:我知道判断素数有更好的算法,我关心的不是算法问题
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

prs1022

铁虫 (初入文坛)

【答案】应助回帖

你首先测试是否为偶数,如果不是,后边的偶数因子可以全部跳过,这样可以减少更多的时间
当然可以也对3做同样的操作,判断数不是2或3的倍数后只要判断那些6k+1和6k+5的因子
可以继续下去,但是越往后提升的效果就越差,所以自行斟酌

对于比较大的数可以考虑用Fermat小定理或者Miller-Rabin检验去判断;双精度最大连续整数上限是2^53

另外可以用matlab coder把m文件转换成mex,这样会更快
stayhungry,stayfoolish
2楼2013-10-24 17:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Simewe 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +5 于海海海海 2026-03-19 5/250 2026-03-19 11:43 by laoshidan
[考研] 085601材料工程专硕求调剂 +7 慕寒mio 2026-03-16 7/350 2026-03-19 10:57 by w沐阳w
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
信息提示
请填处理意见