OFFSET
1,1
COMMENTS
Indices of the primes listed in A124595.
LINKS
MAPLE
P:= select(isprime, [2, seq(i, i=3..10^5, 2)]): nP:= nops(P):
DP:= P[2..-1]-P[1..-2]:
select(t -> DP[t]=28, [$1..nP-1]); # Robert Israel, Mar 18 2026
PROG
(PARI) A(N=100, g=28, 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
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Oct 19 2018
STATUS
approved
