login
A198017
a(n) = n*(7*n + 11)/2 + 1.
11
1, 10, 26, 49, 79, 116, 160, 211, 269, 334, 406, 485, 571, 664, 764, 871, 985, 1106, 1234, 1369, 1511, 1660, 1816, 1979, 2149, 2326, 2510, 2701, 2899, 3104, 3316, 3535, 3761, 3994, 4234, 4481, 4735, 4996, 5264, 5539, 5821, 6110, 6406, 6709, 7019, 7336, 7660, 7991
OFFSET
0,2
COMMENTS
First bisection of A193053 (see also the numerical spiral illustrated in the Links section).
The inverse binomial transform yields 1, 9, 7, 0, 0 (0 continued).
The subsequence of primes begins: 79, 211, 269, 571, 1511, 1979, 3761, 4481, 5821, 6709, 7019. [Jonathan Vos Post, Oct 21 2011]
FORMULA
G.f.: (1+7*x-x^2)/(1-x)^3.
a(n) = A195020(2n) + 2n + 1.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) = 2*a(n-1) - a(n-2) + 7.
MATHEMATICA
Table[(n(7n+11))/2+1, {n, 0, 60}] (* or *) LinearRecurrence[{3, -3, 1}, {1, 10, 26}, 60] (* Harvey P. Dale, Mar 03 2013 *)
PROG
(PARI) for(n=0, 47, print1(n*(7*n+11)/2+1", "));
(Magma) [n*(7*n+11)/2+1: n in [0..47]];
CROSSREFS
Cf. A195020 (vertices of the numerical spiral in link).
Cf. A017005 (first differences).
Sequence in context: A229309 A044452 A299409 * A137351 A134406 A099978
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Oct 21 2011 - based on remarks and sequences by Omar E. Pol
STATUS
approved