login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A124363 a(n) = n^3 + 71*n + 1 4
1, 73, 151, 241, 349, 481, 643, 841, 1081, 1369, 1711, 2113, 2581, 3121, 3739, 4441, 5233, 6121, 7111, 8209, 9421, 10753, 12211, 13801, 15529, 17401, 19423, 21601, 23941, 26449, 29131, 31993, 35041, 38281, 41719, 45361, 49213, 53281, 57571 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

A prime-producing polynomial from a paper by Kermit Rose.

Produces 51 primes between 0 and 99: Length[Flatten[Table[If[PrimeQ[n^3 + 71*n + 1], n^3 + 71*n + 1, {}], {n, 0, 99}]]]

REFERENCES

Prime-Producing Cubic Polynomials (Make Corrections) Joe L. Mott, KERMIT ROSE http://citeseer.ist.psu.edu/403593.html

LINKS

Table of n, a(n) for n=1..39.

Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

FORMULA

G.f.: x*(1+69*x-135*x^2+71*x^3) / (x-1)^4 . - R. J. Mathar, Jan 25 2016

MATHEMATICA

Table[n^3 + 71*n + 1, {n, 0, 99}]

LinearRecurrence[{4, -6, 4, -1}, {1, 73, 151, 241}, 50] (* Harvey P. Dale, Oct 01 2021 *)

PROG

(Maxima) A124363(n):=n^3 + 71*n + 1$

makelist(A124363(n), n, 0, 30); /* Martin Ettl, Nov 08 2012 */

CROSSREFS

Sequence in context: A037994 A138642 A142179 * A044324 A044705 A248399

Adjacent sequences:  A124360 A124361 A124362 * A124364 A124365 A124366

KEYWORD

nonn,easy

AUTHOR

Roger L. Bagula, May 27 2007

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 21 08:47 EST 2022. Contains 350475 sequences. (Running on oeis4.)