OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 8 is a term because 8 is composite, sopfr(8) = 3*2 = 6, 8 mod 6 = 2, and 8 is divisible by 2+6.
MAPLE
filter:= proc(n) local s, m, t;
if isprime(n) then return false fi;
s:= add(t[1]*t[2], t=ifactors(n)[2]);
m:= n mod s;
n mod (m+s) = 0
end proc:
select(filter, [$4..10000]);
MATHEMATICA
sopfr[n_] := Total[Times @@@ FactorInteger[n]];
okQ[n_] := CompositeQ[n] && With[{s = sopfr[n]}, Divisible[n, Mod[n, s]+s]];
Select[Range[4, 1000], okQ] (* Jean-François Alcover, Nov 08 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 08 2021
STATUS
approved