login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

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

  NODES
COMMUNITY 1
INTERN 1