OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 83 = 3*3*3*3+2.
a(2) = 137 = 3*3*3*5+2.
a(3) = 191 = 3*3*3*7+2.
a(4) = 227 = 3*3*5*5+2.
MATHEMATICA
With[{nn=50}, Take[Select[Union[Times@@@Tuples[Prime[Range[2, nn]], 4]+2], PrimeQ], nn]] (* Harvey P. Dale, Oct 18 2013 *)
PROG
(Sage) is_A126711 = lambda n: is_prime(n) and sum(m for p, m in factor(n-2)) == 4 # [D. S. McNeil, Dec 11 2010]
(PARI) list(lim)=my(t, tt, v=List()); forprime(p=3, lim\27, forprime(q=p, lim\p\9, forprime(r=q, lim\p\q\3, t=p*q*r; forprime(s=r, lim\t, if(ispseudoprime(tt=t*s+2), listput(v, tt)))))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Feb 17 2011
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jonathan Vos Post, Feb 12 2007
EXTENSIONS
Corrected (299 removed) by D. S. McNeil, Dec 10 2010
STATUS
approved