OFFSET
1,2
MATHEMATICA
nn = 2^20; Array[Set[a[#], #] &, 2]; Do[If[EvenQ[n], Set[a[n], 2 a[n/2]], Set[a[n], Times @@ Power @@@ Map[{Prime[PrimePi[#1] + 1], #2} & @@ # &, FactorInteger[a[(n + 1)/2]]]]], {n, 3, nn}]; Select[Range[1, nn, 2], Divisible[a[#], #] &] (* Michael De Vlieger, Jul 28 2023 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 28 2023
STATUS
approved