OFFSET
1,1
COMMENTS
Indices of the primes listed in A061779.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..10040
MATHEMATICA
Select[Range[3000], Prime[#] + 22 == Prime[# + 1] &] (* Vincenzo Librandi, Mar 22 2019 *)
PROG
(PARI) A(N=100, g=22, 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..3000] | NthPrime(n+1) - NthPrime(n) eq 22]; // Vincenzo Librandi, Mar 22 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Oct 19 2018
STATUS
approved