rmstrong <- function(x) {
tmp <- strsplit(as.character(x), split="")
cubic <- sapply(tmp, function(y)sum(as.numeric(y)^3))
return(cubic == x)
}

s <- 1:1000000
s[armstrong(s)]


Elimde bu şekilde bir kod var ama sadece 3 basamaklı armstrong sayılarını gösteriyor.