A373807 Indices of prime terms in A373805.

2, 4, 8, 10, 13, 15, 17, 23, 25, 27, 37, 41, 46, 48, 50, 52, 55, 57, 60, 62, 64, 74, 76, 79, 81, 83, 88, 90, 92, 94, 97, 99, 115, 117, 119, 124, 129, 133, 138, 141, 153, 156, 163, 166, 171, 173, 175, 179, 181, 183, 205, 207, 219, 224, 226, 236, 240, 245, 247, 250, 254, 258, 276, 278, 281, 283

Offset: 1

Author

N. J. A. Sloane, Aug 11 2024

Keywords

• nonn

Comments

Please see the comments on A373805 and A373806 for further information.

a(608) > 65536. - Michael De Vlieger, Aug 12 2024

Links

• Michael De Vlieger, Table of n, a(n) for n = 1..607 (First 319 terms from N. J. A. Sloane)

Crossrefs

Cf. A373805, A373806.

Programs

Mathematica

nn = 2^14; s = j = 1; Reap[Monitor[Do[If[PrimeQ[j], Sow[n - 1]; s = -s; k = Prime[n] + s, k = 2 j + s]; j = k, {n, 2, nn}], n] ][[-1, 1]] (* Michael De Vlieger, Aug 12 2024 *)

Python

# uses imports and A373805_gen in A373805
def agen(): # generator of terms
    yield from (i for i, ai in enumerate(A373805_gen(), 1) if isprime(ai))
print(list(islice(agen(), 66))) # Michael S. Branicky, Aug 12 2024