OFFSET
1,1
COMMENTS
The numbers prime(n)^5 are in the sequence.
Also contains 4^p p and 2^p p^3 for any prime p>2, and 3^p p^2 for any prime > 3. - Robert Israel, Jun 20 2017
From David A. Corneth, Apr 24 2020: (Start)
Suppose we look for terms below (inclusive) u. Let maxp be the largest prime factor that is the multiple of at least one of those terms. Then maxp is the largest prime below (inclusive) u^(1/5).
Proof: The sum of prime factors counted with multiplicity of a term t divisible by maxp is 5*maxp. The smallest product of primes summing to 5*maxp where the largest prime factor of t is maxp is maxp^5 which must be <= u. Solving this gives maxp is the largest prime below (inclusive) u^(1/5).
This enables us to search through the positive integers via a tree starting at 1. (End)
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000 (first 200 terms from Robert Israel. Following 4800 terms from Amiram Eldar)
David A. Corneth, PARI program
David A. Corneth, Terms <= 10^12
EXAMPLE
192 is in the sequence because 192 = 2^6 * 3 => sum of prime divisors = 2*6 + 3 = 15 = 5*3 where 3 is the greatest prime divisor.
MAPLE
with(numtheory):A:= proc(n) local e, j; e := ifactors(n)[2]: add (e[j][1]*e[j][2], j=1..nops(e)) end: for m from 2 to 500000 do: x:=factorset(m):n1:=nops(x):if A(m)=5*x[n1] then printf(`%d, `, m):else fi:od:
MATHEMATICA
Select[Range[2, 10^5], Plus @@ Times @@@ (f = FactorInteger[#]) == 5 * f[[-1, 1]] &] (* Amiram Eldar, Apr 24 2020 *)
PROG
(PARI) See Corneth link \\ David A. Corneth, Apr 24 2020
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michel Lagneau, May 29 2012
STATUS
approved