|
|
A166524
|
|
a(n) = 9*n - a(n-1), with n>1, a(1)=1.
|
|
1
|
|
|
1, 17, 10, 26, 19, 35, 28, 44, 37, 53, 46, 62, 55, 71, 64, 80, 73, 89, 82, 98, 91, 107, 100, 116, 109, 125, 118, 134, 127, 143, 136, 152, 145, 161, 154, 170, 163, 179, 172, 188, 181, 197, 190, 206, 199, 215, 208, 224, 217, 233, 226, 242, 235, 251, 244, 260, 253
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
G.f.: -x*(-1-16*x+8*x^2) / ( (1+x)*(x-1)^2 ). - R. J. Mathar, Mar 08 2011
E.g.f.: (1/4)*(23*exp(-x) + 9*(1 + 2*x)*exp(x) - 32).
a(n) = a(n-1) + a(n-2) - a(n-3). (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = 1/8 + cot(Pi/9)*Pi/9. - Amiram Eldar, Feb 24 2023
|
|
MATHEMATICA
|
CoefficientList[Series[-(- 1 - 16 x + 8 x^2)/((1 + x) (x - 1)^2), {x, 0, 80}], x] (* Vincenzo Librandi, Sep 13 2013 *)
LinearRecurrence[{1, 1, -1}, {1, 17, 10}, 60] (* Harvey P. Dale, Dec 24 2014 *)
|
|
PROG
|
(Magma) [n eq 1 select 1 else 9*n-Self(n-1): n in [1..80]]; // Vincenzo Librandi, Sep 13 2013
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|