%I #15 Mar 24 2017 04:07:13
%S 19,89,271,4159,10099,20879,28729,38611,50849,130681,391249,457651,
%T 1049599,1187009,1501849,1875529,3113459,3420649,3750031,4102649,
%U 6767801,7893289,9837631,10559249,11319859,14780179,17854849,21385999,31646249
%N Primes of form n^4 + n^2 - 1.
%C Primes of the form A071253(n) - 1. - _Altug Alkan_, Mar 24 2017
%H Robert Israel, <a href="/A174819/b174819.txt">Table of n, a(n) for n = 1..10000</a>
%e a(2) = 89 is in the sequence because 3^4 + 3^2 - 1 = 89 is prime.
%p select(isprime, [seq(n^4+n^2-1, n=1..1000)]); # _Robert Israel_, Mar 24 2017
%t a={}; Do[p=n^4+n^2-1; If[PrimeQ[p], AppendTo[a, p]], {n, 10^2}]; Print[a];
%o (PARI) for(n=1, 1e3, if(isprime(p=n^4+n^2-1), print1(p ", "))) \\ _Altug Alkan_, Mar 24 2017
%Y Cf. A071253.
%K nonn
%O 1,1
%A _Michel Lagneau_, Dec 01 2010
|