|
%I #16 Sep 16 2022 09:02:13
%S 1,73,151,241,349,481,643,841,1081,1369,1711,2113,2581,3121,3739,4441,
%T 5233,6121,7111,8209,9421,10753,12211,13801,15529,17401,19423,21601,
%U 23941,26449,29131,31993,35041,38281,41719,45361,49213,53281,57571
%N a(n) = n^3 + 71*n + 1
%C Produces 51 primes between 0 and 99.
%H Joe L. Mott and Kermit Rose, <a href="http://www.math.fsu.edu/~aluffi/archive/paper134.ps">Prime-Producing Cubic Polynomials</a>.
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-6,4,-1).
%F G.f.: x*(1+69*x-135*x^2+71*x^3) / (x-1)^4 . - _R. J. Mathar_, Jan 25 2016
%t Table[n^3 + 71*n + 1, {n, 0, 99}]
%t LinearRecurrence[{4,-6,4,-1},{1,73,151,241},50] (* _Harvey P. Dale_, Oct 01 2021 *)
%o (Maxima) A124363(n):=n^3 + 71*n + 1$
%o makelist(A124363(n),n,0,30); /* _Martin Ettl_, Nov 08 2012 */
%K nonn,easy
%O 1,2
%A _Roger L. Bagula_, May 27 2007
|
