OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..4000
MAPLE
spf:= proc(n) local t; add(t[1]*t[2], t=ifactors(n)[2]) end proc:
f:= proc(n) local a, b;
a:= spf(n);
b:= numtheory:-sigma(n);
a+b+spf(b)+numtheory:-sigma(a)
end proc:
S:= select(t -> t < 1000 and isprime(t), map(f, {$1..1000})):
sort(convert(S, list));
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Apr 14 2021
STATUS
approved