A361502 - OEIS

A361502

Index of n-th prime in A359804.

2, 3, 4, 8, 13, 42, 347, 3466, 49012, 528231, 717126, 63056215, 1375559400, 7038527851

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,1

COMMENTS

Theorem: Every prime appears in A359804. For proof see A359804.

It appears that the primes in A359804 appear in order.

LINKS

Rémy Sigrist, C++ program

MATHEMATICA

nn = 2^20; c[_] = False; q[_] = 1;

i = 1; j = 2; c[1] = c[2] = True; u = 3;

{2}~Join~Reap[Monitor[Do[

(k = q[#]; While[c[k #], k++]; k *= #;

While[c[# q[#]], q[#]++]) &[(p = 2;

While[Divisible[i j, p], p = NextPrime[p]]; p)];

If[PrimeQ[k], Sow[n]; Print[n]];

Set[{c[k], i, j}, {True, j, k}];

If[k == u, While[c[u], u++]], {n, 3, nn}], n]][[-1, -1]] ] (* Michael De Vlieger, Mar 19 2023 *)

PROG

(C++) See Links section.

CROSSREFS

KEYWORD

nonn,more

AUTHOR

N. J. A. Sloane, Mar 18 2023, based on a comment made by Michael De Vlieger in A359804 in which he gave the values of a(1) to a(12).

EXTENSIONS

a(13)-a(14) from Rémy Sigrist, Mar 19 2023

STATUS

approved