R程序分享2：prime

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

用筛法求素数，最简单的方法
用我4g，32位机试了下，能算到1e6，到1e7就不行了
写着玩的，跟计算pi的代码差不多可以用于评测机器性能，呵呵

1. # 本函数用筛法计算素数。
   
2. 
   
3. # x是取值范围，用c(min, max)表示，若仅有一个值表示最大值。
   
4. 
   
5. # 最后返回的是指定范围内的素数向量。
   
6. 
   
7. # Created on 2014.5.30
   
8. # Modified on 2014.5.30
   
9. 
   
10. prime <- function(x){
    
11.   stopifnot(is.numeric(x))
    
12.   
    
13.   # 检查x的合法性
    
14.   if (length(x) == 1){
    
15.     x[2] <- x
    
16.     x[1] <- 2
    
17.   }
    
18.   
    
19.   res <- c(x[1]:x[2])
    
20.   flag <- rep(TRUE, x[2]-x[1]+1)
    
21.   
    
22.   for (i in 2:sqrt(x[2])){
    
23.     flag[which((res %% i) == 0 & res > i)] <- FALSE
    
24.   }
    
25.   
    
26.   return(res[flag])
    
27. }
    
28. 
    
29. Demo <- function(){
    
30.   p <- prime(1000)
    
31.   p <- prime(c(800,1000))
    
32. }

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：prime IME Rim R程序 Modified 程序

感谢分享。。。。。