login
Numbers k such that 6*k + 1 is prime.
33

%I #59 Feb 04 2024 01:11:05

%S 1,2,3,5,6,7,10,11,12,13,16,17,18,21,23,25,26,27,30,32,33,35,37,38,40,

%T 45,46,47,51,52,55,56,58,61,62,63,66,68,70,72,73,76,77,81,83,87,90,91,

%U 95,96,100,101,102,103,105,107,110,112,115,118,121,122,123,125,126,128,131,135,137

%N Numbers k such that 6*k + 1 is prime.

%C For all elements of this sequence there are no (x,y) positive integers such that a(n)=6*x*y+x+y or a(n)=6*x*y-x-y. - _Pedro Caceres_, Apr 19 2019

%H Ray Chandler, <a href="/A024899/b024899.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from T. D. Noe)

%F a(n) = A024892(n)/2 = (A034936(n)+1)/2. - _Ray Chandler_, Dec 26 2003

%F a(n) = (A002476(n)-1)/6. - _Zak Seidov_, Aug 31 2016

%p a := [ ]: for n from 0 to 400 do if isprime(6*n+1) then a := [ op(a), n ]; fi; od: A002476 := n->a[n];

%t Select[Range@ 140, PrimeQ[6 # + 1] &] (* _Michael De Vlieger_, Jan 23 2018 *)

%o (Magma) [n: n in [0..200]| IsPrime(6*n+1)] // _Vincenzo Librandi_, Nov 20 2010

%o (PARI) select(n->n%6==1,primes(100))\6 \\ _Charles R Greathouse IV_, Apr 28 2015

%Y A002476 gives primes, A091178 gives prime index.

%Y Complement of A046954 relative to A000027.

%K nonn,easy

%O 1,2

%A _Clark Kimberling_