cp's OEIS Frontend

A309844 Primes of the form n^4 + n^2 + 3.

%I A309844 #20 Sep 28 2019 22:37:02
%S A309844 3,5,23,653,10103,83813,160403,234743,280373,1049603,3420653,6252503,
%T A309844 11319863,52207853,92246423,146422103,174913853,221548343,442071653,
%U A309844 479807123,577224653,607597853,655385603,937921253,1222865933,1249233683,1387525253,1506177293
%N A309844 Primes of the form n^4 + n^2 + 3.
%C A309844 Digital root of all values > 3 is 5, compare A017221.
%t A309844 f[n_] := n^4 + n^2 + 3; Select[f /@ Range[0, 200], PrimeQ] (* _Amiram Eldar_, Aug 24 2019 *)
%o A309844 (MATLAB)
%o A309844 a = [];
%o A309844 for n = 0:1e3
%o A309844     x = n.^4+n.^2+3;
%o A309844     if isprime(x); a = [a,x]; end;
%o A309844 end
%o A309844 (Python)
%o A309844 from sympy import isprime
%o A309844 a = []
%o A309844 for n in range(0,1000):
%o A309844     x = n**4+n**2+3
%o A309844     if isprime(x):
%o A309844         a.append(x)
%Y A309844 Subset of A027753. Subset of A017221.
%Y A309844 Cf. A037896, A182344, A182345, A049423.
%K A309844 nonn,easy
%O A309844 1,1
%A A309844 _Christopher R. Madan_, Aug 19 2019