OFFSET
0,2
LINKS
Michael De Vlieger, Table of n, a(n) for n = 0..10000
Michael De Vlieger, Annotated plot of prime factors p | a(n) at (n, pi(p)) for n = 0..50.
MATHEMATICA
nn = 58; r = c[_] = 0; m = 1; Array[Set[{a[#], c[#]}, {#, #}] &, 2, 0]; Do[k = SelectFirst[Union@ Map[Total, Rest@ Subsets[2^Reverse[Length[#] - Position[#, 1][[All, 1]]] &@ IntegerDigits[2^(r + 2) - m - 1, 2]]], c[#] == 0 &]; Set[{a[n], c[k]}, {k, n}]; m += a[n]; If[And[IntegerQ[#], # > 0], m -= a[#]] &[n/2]; If[And[EvenQ[k], PrimePowerQ[k], k > 2^r], r++], {n, 2, nn}]; Table[Times @@ Map[Prime, 1 + Length[#] - Position[#, 1][[All, 1]]] &@ IntegerDigits[a[n], 2], {n, 0, nn}] (* Michael De Vlieger, Jul 14 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, Jul 14 2022
STATUS
approved