Revision as of 01:55, 5 November 2022 (view source) GordonCharlton (talk \| contribs) (Added Quackery.) ← Older edit		Revision as of 12:17, 2 January 2023 (view source) Davgot (talk \| contribs) m (Fixed indentation) Newer edit →
Line 528: sieve <- function(n) { if (n < 2) return (NULL) primes <- rep(TRUE, n) Line 534: for (i in 1:floor(sqrt(n))) if (primes[i]) primes[seq(i*i, n, by = i)] <- FALSE which(primes) Line 544: factors <- primes[n %% primes == 0] if (length(factors) == 0) n n else { for (p in factors) { # add all elements of n/p that are also prime d <- n / p if (d != p && all(d %% primes[primes <= floor(sqrt(d))] != 0)) factors <- c(factors, d) } } factors } } cat("The prime factors of 600,851,475,143 are", paste(prime.factors(600851475143), collapse = ", "), "\n") </syntaxhighlight> {{Out}}

Largest prime factor: Difference between revisions

Largest prime factor (view source)