%I #22 Sep 08 2022 08:44:52
%S 3,5,7,10,12,16,22,23,25,26,32,35,36,38,40,46,50,51,57,58,61,63,65,67,
%T 68,71,73,77,78,80,85,87,91,98,100,102,105,107,123,137,138,140,142,
%U 143,150,155,158,161,165,166,172,177,178,180,187,190,191,192,197,200,205
%N Numbers n such that 666*n + 1 is prime.
%H Vincenzo Librandi, <a href="/A037030/b037030.txt">Table of n, a(n) for n = 1..1000</a>
%t f[n_]:=666*n+1; lst={};Do[p=f[n];If[PrimeQ[p],AppendTo[lst,n]],{n,0,6!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Jul 02 2009 *)
%t Select[Range[3, 1000], PrimeQ[666 # + 1] &] (* _Vincenzo Librandi_, May 19 2014 *)
%o (Magma) [n: n in [1..1000] | IsPrime(666*n + 1)]; // _Vincenzo Librandi_, Nov 17 2010
%o (PARI) is(n)=isprime(666*n+1) \\ _Charles R Greathouse IV_, Feb 17 2017
%Y Cf. A037029.
%K nonn,easy
%O 1,1
%A _G. L. Honaker, Jr._
%E More terms from _Michel ten Voorde_