OFFSET
1,1
PROG
(PARI) lista(nn) = {forprime(p=3, nn, if ((((p-1) % 4) == 0) && isprime(p+2), print1(p*(p+2), ", ")); ); } \\ Michel Marcus, Feb 24 2015
(Magma) [16*n^2+16*n+3: n in [1..700]| IsPrime(4*n+1) and IsPrime(4*n+3)]; // Vincenzo Librandi, Feb 24 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jun 04 2002
EXTENSIONS
More terms from Michel Marcus, Feb 24 2015
STATUS
approved