divisors(1) # 1
divisors(2) # 1 2
divisors(2^5) # 1 2 4 8 16 32
divisors(1000) # 1 2 4 5 8 10 ... 100 125 200 250 500 1000
divisors(1001) # 1 7 11 13 77 91 143 1001
# unitary divisors function
unitary_divisors <- function(n) {
divs <- divisors(n)
adiv <- apply(cbind(divs, rev(divs)), 1, mGCD) == 1
return(divs[adiv])
}
divisors(120)
## [1] 1 2 3 4 5 6 8 10 12 15 20 24 30 40 60 120
unitary_divisors(120)
## 1 3 5 8 15 24 40 120
Run the code above in your browser using DataLab