cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A364540 Numbers k for which A005940(k) < k.

This page as a plain text file.
%I A364540 #8 Jul 28 2023 10:34:01
%S A364540 9,17,18,33,34,35,36,65,66,67,68,69,70,72,129,130,131,132,133,134,135,
%T A364540 136,137,138,140,144,257,258,259,260,261,262,263,264,265,266,267,268,
%U A364540 270,272,273,274,276,280,288,289,385,513,514,515,516,517,518,519,520,521,522,523,524,525,526,527,528,529,530,531
%N A364540 Numbers k for which A005940(k) < k.
%F A364540 a(n) = 1 + A356450(n).
%t A364540 nn = 540; 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 *)
%o A364540 (PARI)
%o A364540 A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
%o A364540 isA364540(n) = (A005940(n)<n);
%Y A364540 Positions of negative terms in A364499.
%Y A364540 One more than A356450.
%Y A364540 Subsequence of A364541, complement of A364542.
%Y A364540 Cf. A005940.
%K A364540 nonn
%O A364540 1,1
%A A364540 _Antti Karttunen_, Jul 28 2023