|
%I #13 Oct 25 2018 11:09:01
%S 738,1315,3022,3042,3607,4112,4300,4362,4555,4619,4761,4893,5204,5358,
%T 5615,5637,6215,6265,6479,6610,6706,6933,7295,7829,7884,8049,8198,
%U 8548,9085,9155,9524,9588,9641,9826,9924,10463,10824,11367,11590,11701,11729,11869,12159,12258,12275,12327
%N Indices of primes followed by a gap (distance to next larger prime) of 32.
%C Indices of the primes listed in A126784.
%H Robert Israel, <a href="/A320714/b320714.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Pri#gaps">Index entries for primes, gaps between</a>.
%F a(n) = A000720(A126784(n)).
%F A320714 = { i>0 | prime(i+1) = prime(i) + 32 }.
%p p:= 2: Res:= NULL: count:= 0:
%p for k from 1 while count < 100 do
%p q:= nextprime(p);
%p if q - p = 32 then
%p count:= count+1;
%p Res:= Res, k;
%p fi;
%p p:= q;
%p od:
%p Res; # _Robert Israel_, Oct 25 2018
%o (PARI) A(N=100,g=32,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
%Y Equals A000720 o A126784.
%Y Indices of 32's in A001223.
%Y Row 16 of A174349.
%Y Cf. A029707, A029709, A320701, A320702, ..., A320720 (analog for gaps 2, 4, 6, 8, ..., 44), A116493 (gap 70), A116496 (gap 100), A116497 (gap 200), A116495 (gap 210).
%K nonn
%O 1,1
%A _M. F. Hasler_, Oct 19 2018
|
