OFFSET
1,1
COMMENTS
Except for the first, all the terms above generate brilliant numbers.
Numbers n such that n - 1 + i and n + 1 + i are (twin) Gaussian primes, see Shanks. - Charles R Greathouse IV, Apr 20 2011
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Daniel Shanks, A Note on Gaussian Twin Primes, Mathematics of Computation 14:70 (1960), pp. 201-203.
FORMULA
a(k) = A096012(k) + 1. (Because n^4+4 = ((n-1)^2+1)((n+1)^2+1).) - Jeppe Stig Nielsen, Feb 26 2016
MATHEMATICA
Select[Range[5000], PrimeOmega[#^4+4]==2&] (* Harvey P. Dale, Sep 07 2017 *)
PROG
(PARI) forstep(n=1, 1e5, 2, if(isprime(n^2-2*n+2) && isprime(n^2+2*n+2), print1(n", "))) \\ Charles R Greathouse IV, Apr 20 2011
(Magma) IsSemiprime:=func< n | &+[ k[2]: k in Factorization(n) ] eq 2 >; [ n: n in [1..5000] | IsSemiprime(n^4+4)]; // Vincenzo Librandi, Apr 20 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Jul 10 2005
STATUS
approved