|
|
A364541
|
|
Numbers k for which A005940(k) <= k.
|
|
5
|
|
|
1, 2, 3, 4, 5, 6, 8, 9, 10, 12, 16, 17, 18, 20, 24, 32, 33, 34, 35, 36, 40, 48, 64, 65, 66, 67, 68, 69, 70, 72, 80, 96, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 140, 144, 160, 192, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 270, 272, 273, 274, 276, 280, 288, 289, 320, 384, 385, 512
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
A029747 is included as a subsequence, because it gives the known fixed points of map n -> A005940(n).
|
|
LINKS
|
|
|
MATHEMATICA
|
nn = 512; 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[nn], a[#] <= # &] (* Michael De Vlieger, Jul 28 2023 *)
|
|
PROG
|
(PARI)
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
|
|
CROSSREFS
|
Positions of nonpositive terms in A364499.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|