login
A373807
Indices of prime terms in A373805.
4
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,1
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)
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 *)
PROG
(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
CROSSREFS
Sequence in context: A154788 A219617 A100278 * A075333 A297250 A014190
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Aug 11 2024.
STATUS
approved