- prime <- function(x){
- if(x==1) return ("NULL")
- y <- c()
- a <- c(2:x)
- b <- min(a[which(x%%a==0)])
- x<-x/b
- if(x==1) return ("NULL")
- y <-c(y,b)
- while(!x==1){
- a <- c(2:x)
- b <- min(a[which(x%%a==0)])
- x<-x/b
- y <- c(y,b)
- }
- return(max(y))
- }
复制代码
参考之前的帖子,把union改成max
https://bbs.pinggu.org/forum.php?mod=viewthread&tid=5986371&page=1#pid46704368