OFFSET
0,1
EXAMPLE
MAPLE
sopfr:= proc(n) local F, t;
F:= ifactors(n)[2];
add(t[1]*t[2], t=F)
end proc:
M:= 2*10^6: Mh:=M/2-2:
g:= proc(t) local v; v:= t - sopfr(t); if v > 0 and v < Mh then v fi end proc:
R:= sort(map(g, [$1..M])):
V:= Array(0..20):
V[0]:= 4:
state:= 1: r:= R[1]:
for i from 2 to nops(R) do
if R[i] = r then
state:= state+1
else
if state <= 20 and V[state] = 0 then V[state]:= r fi;
state:= 1: r:= R[i]
fi
od:
convert(V, list);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Will Gosnell and Robert Israel, Sep 14 2025
STATUS
approved
