OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
98420 is in the sequence because the prime divisors are 2, 5, 7, 19, 37 and the sum 2 + 5 + 7 + 19 + 37 = 70 = 2*(37 - 2).
MAPLE
filter:= proc(n) local P; P:= numtheory:-factorset(n);
convert(P, `+`) = 2*(max(P)-min(P))
end proc:
select(filter, [$1..10000]); # Robert Israel, Apr 09 2019
MATHEMATICA
Select[Range[5500], Plus@@((pl=First/@FactorInteger[#])/2)==pl[[-1]]-pl[[1]]&]
PROG
(PARI) isok(n) = if (n>1, my(f=factor(n)[, 1]); 2*(vecmax(f) - vecmin(f)) == vecsum(f)); \\ Michel Marcus, Apr 10 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Nov 13 2011
STATUS
approved