login
Primes q such that q = floor(p/3) for some prime p.
1

%I #14 Sep 08 2022 08:45:46

%S 2,3,5,7,13,17,19,23,29,37,43,59,79,83,89,97,103,127,139,149,163,167,

%T 173,197,199,227,233,239,257,269,293,313,317,337,349,353,367,383,397,

%U 409,419,433,439,457,479,499,503,523,569,577,607,643,659,709,757,769

%N Primes q such that q = floor(p/3) for some prime p.

%C Primes q such that 3*q+1 or 3*q+2 is prime. Agrees with A023208 except for initial term 2.

%C Essentially the same as A023208. - _R. J. Mathar_, Jul 05 2009

%e 3 is in the sequence since 11 is prime and floor(11/3) = 3; 11 is not in the sequence since 11 = floor(34/3) = floor(35/3) and neither 34 nor 35 is prime.

%t lst={};Do[r=Prime[n];If[PrimeQ[p=Floor[r/3]],AppendTo[lst,p]],{n,6!}];lst

%t Select[Floor[Prime[Range[350]]/3],PrimeQ] (* _Harvey P. Dale_, Aug 26 2013 *)

%t Select[Prime[Range[200]],AnyTrue[3#+{1,2},PrimeQ]&] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Aug 07 2019 *)

%o (Magma) [ q: q in PrimesUpTo(800) | IsPrime(3*q+1) or IsPrime(3*q+2) ]; // _Klaus Brockhaus_, Jul 06 2009

%o (PARI) isA162338(n) = isprime(n) && (isprime(3*n+1) || isprime(3*n+2)) \\ _Michael B. Porter_, Jan 30 2010

%Y Cf. A162337. Essentially the same as A023208 (n and 3n+2 are both prime).

%K nonn

%O 1,1

%A _Vladimir Joseph Stephan Orlovsky_, Jul 01 2009

%E Edited and listed terms verified by _Klaus Brockhaus_, Jul 06 2009