OFFSET
1,1
COMMENTS
Indices of the primes listed in A031934.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..9790
MATHEMATICA
Select[Range[2500], Prime[#] + 16 == Prime[# + 1] &] (* Vincenzo Librandi, Mar 19 2019 *)
Position[Partition[Prime[Range[2100]], 2, 1], _?(#[[2]]-#[[1]]== 16&), 1, Heads-> False]//Flatten (* Harvey P. Dale, Nov 01 2020 *)
PROG
(PARI) A(N=100, g=16, p=2, i=primepi(p)-1, L=List())={forprime(q=1+p, , i++; if(p+g==p=q, listput(L, i); N--||break)); Vec(L)} \\ returns the list of first N terms of the sequence
(Magma) [n: n in [1..2100] | NthPrime(n+1) - NthPrime(n) eq 16]; // Vincenzo Librandi, Mar 19 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Oct 19 2018
STATUS
approved