OFFSET
1,1
COMMENTS
Composite numbers k such that k + sopfr(k) is a square, where sopfr(k) is the sum of prime factors of k with multiplicity.
Contains no semiprimes.
Includes 9*p if p is a prime of the form (x^2-6)/10 where x == 4 or 6 (mod 10).
Is this sequence disjoint from A386245?
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 171 is a term because 171 = 3^2 * 19 is composite and 171 + 3 + 3 + 19 = 196 = 14^2 is a square.
MAPLE
filter:= proc(n) local t;
if isprime(n) then return false fi;
issqr(n + add(t[1]*t[2], t=ifactors(n)[2]))
end proc:
select(filter, [$4..10000]);
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Jul 16 2025
STATUS
approved
