%I #18 Sep 08 2022 08:46:15
%S 7,13,21,24,25,30,33,44,51,53,55,60,65,68,71,79,80,84,87,88,98,104,
%T 106,108,110,113,116,122,135,136,140,148,152,158,159,162,165,169,174,
%U 176,184,189,191,196,197,199,204,209,211,216,218,223,227,234,237,245
%N Numbers k such that prime(k) == 1 (mod 8).
%C The asymptotic density of this sequence is 1/4 (by Dirichlet's theorem). - _Amiram Eldar_, Mar 01 2021
%H Robert Israel, <a href="/A269704/b269704.txt">Table of n, a(n) for n = 1..10000</a>
%e a(1) = 7 because prime(7) = 17 and 17 == 1 (mod 8).
%p Res:= NULL: count:= 0:
%p p:= 2:
%p for n from 2 while count < 100 do
%p p:= nextprime(p);
%p if p mod 8 = 1 then count:= count+1; Res:= Res, n fi
%p od:
%p Res; # _Robert Israel_, May 06 2019
%t Select[Range[300], Mod[Prime[#], 8] == 1 &]
%o (Magma) [n: n in [1..500] | NthPrime(n) mod 8 eq 1];
%o (PARI) lista(nn) = for(n=1, nn, if(Mod(prime(n),8)==1, print1(n, ", "))); \\ _Altug Alkan_, Mar 04 2016
%Y The associated primes are in A007519.
%Y Cf. similar sequences listed in A269703.
%K nonn
%O 1,1
%A _Vincenzo Librandi_, Mar 04 2016