作者haolihy (好厉害)
看板R_Language
标题[问题] 质因数分解
时间Wed Nov 4 17:17:24 2015
[问题类型]:
效能谘询(我想让R 跑更快)
[软体熟悉度]:
入门(写过其他程式，只是对语法不熟悉)
[问题叙述]:
质因数分解
[程式范例]:
想写一个程式对正整数n作质因数分解
目标是n<10^8 可以1min跑完
我的程式大概要6~7min 改进空间应该很大
而太慢的主因应该是列出<n质数表的程式太慢
请大家赐教
质数表的程式:
primeY <- function(num){
sqnum <- floor(sqrt(num))
x <- c(2:num)
y <- c()
repeat{
if(x[1] > sqnum) break
y <- c(y, x[1])
x <- x[x %% x[1] != 0]
}
return(c(y,x))
}
质因数分解的程式:
primefactor <- function(num){
p <- primeY(num)[which(num %% primeY(num) == 0)]
q <- c()
e <- c()
for(i in 1:length(p)){
e[i] <- 1
q[i] <- num / p[i]
repeat{
if(q[i] %% p[i] != 0) break
e[i] <- e[i] + 1
q[i] <- q[i] / p[i]
}
}
return(list('pfactor'=p,
'multi'=e))
}
[环境叙述]:
R 3.2.2
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 210.61.134.14
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1446628646.A.D6C.html