A230307 a(n) = 139*n^2 - 2307*n + 3331.

3331, 1163, -727, -2339, -3673, -4729, -5507, -6007, -6229, -6173, -5839, -5227, -4337, -3169, -1723, 1, 2003, 4283, 6841, 9677, 12791, 16183, 19853, 23801, 28027, 32531, 37313, 42373, 47711, 53327, 59221, 65393, 71843, 78571, 85577, 92861, 100423, 108263

Offset

0,1

Comments

|a(n)| are distinct noncomposite numbers for n = 0 to 35.

The values of this polynomial are never divisible by a prime less than 53.

Links

Arkadiusz Wesolowski, Table of n, a(n) for n = 0..1000

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

Formula

G.f.: (3331 - 8830*x + 5777*x^2)/(1 - x)^3.

Maple

seq(139*n^2-2307*n+3331, n=0..37);

Mathematica

Table[139*n^2 - 2307*n + 3331, {n, 0, 37}]
LinearRecurrence[{3, -3, 1}, {3331, 1163, -727}, 40] (* or *) CoefficientList[ Series[(3331-8830x+5777x^2)/(1-x)^3, {x, 0, 40}], x] (* Harvey P. Dale, Jul 06 2021 *)

Prog

(Magma) [139*n^2-2307*n+3331 : n in [0..37]];
(PARI) for(n=0, 37, print1(139*n^2-2307*n+3331, ", "));

Crossrefs

Sequence in context: A253865 A253538 A147881 * A043504 A251331 A261656

Adjacent sequences: A230304 A230305 A230306 * A230308 A230309 A230310

Keyword

sign,easy

Author

Arkadiusz Wesolowski, Oct 29 2013

Status

approved