OFFSET
1,1
COMMENTS
If n = p_1^c_1 * p_2^c_2 * p_3^c_3 * ... * p_k^c_k, where c's are positive integers and p's are distinct primes, then product{j=1 to k}[c_j] = sum{j=1 to k}[p_j+c_j].
EXAMPLE
9000 = 3^2 * 2^3 * 5^3. Product of exponents is 2*3*3=18, sum of prime factors and exponents is 3+2+2+3+5+3=18, hence 9000 is in the sequence.
MATHEMATICA
t = {}; n = 1; While[Length[t] < 38, n++; f = FactorInteger[n]; sm = Total[Flatten[f]]; pr = Times @@ Transpose[f][[2]]; If[sm == pr, AppendTo[t, n]]]; t (* T. D. Noe, Nov 08 2013 *)
peQ[n_]:=Module[{fi=FactorInteger[n]}, Times@@fi[[All, 2]]==Total[ Flatten[ fi]]]; Select[Range[400000], peQ] (* Harvey P. Dale, May 21 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Nov 06 2013
STATUS
approved