A322923 Primes of the form 3*p + 4, where p is a prime.
13, 19, 37, 43, 61, 73, 97, 127, 163, 181, 223, 241, 271, 307, 313, 331, 397, 421, 457, 523, 541, 547, 577, 601, 673, 691, 727, 757, 811, 853, 883, 937, 997, 1051, 1063, 1123, 1153, 1171, 1231, 1297, 1303, 1321, 1531, 1567, 1627, 1693, 1783, 1801
Offset: 1
Links
- Muniru A Asiru, Table of n, a(n) for n = 1..10000
Programs
-
GAP
P:=Filtered([1..1000],IsPrime);; a:=Filtered(List(P,i->3*i+4),k->IsPrime(k)); # Muniru A Asiru, Mar 23 2019
-
Magma
[a: p in PrimesUpTo(600) | IsPrime(a) where a is 3*p+4];
-
Maple
select(isprime,[3*ithprime(p)+4$p=1..120]); # Muniru A Asiru, Mar 23 2019
-
Mathematica
Select[Table[p=Prime[n];3p+4,{n,85}],PrimeQ] -
PARI
terms(n) = my(x=0, i=0); forprime(p=1, , if(i >= n, break); x=3*p+4; if(ispseudoprime(x), print1(x, ", "); i++)) /* Print initial 50 terms as follows: */ terms(50) \\ Felix Fröhlich, Mar 23 2019