require 'prime' def fact(n) (1..n).inject(:*); end fact42 = fact(42) p Prime.prime_division(fact42)