OFFSET
0,1
COMMENTS
a(n) = largest odd divisor of A059100(n+1). Proof: Observe that a(2n) = A059100(2n+1) and a(2n+1) = (A059100(2n+2))/2 and note that (A059100(m))/2 is odd for even m. - Jeremy Gardiner, Aug 25 2013
a(n) is also the denominator of the (n+1)-st largest circle in a special case of the Pappus chain inspired by the Yin-Yang symbol. See illustration in the links. - Kival Ngaokrajang, Jun 20 2015
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Kival Ngaokrajang, Illustration of initial terms.
Index entries for linear recurrences with constant coefficients, signature (0,3,0,-3,0,1).
FORMULA
G.f.: ( -3-3*x-2*x^2-3*x^4-x^5 ) / ( (x-1)^3*(1+x)^3 ). - R. J. Mathar, Jan 21 2011
a(n) = ((-1)^n+3)*(n^2+2*n+3)/4. - Bruno Berselli, Jan 21 2011
From Amiram Eldar, Aug 09 2022: (Start)
a(n) = numerator(((n+1)^2 + 2)/2).
Sum_{n>=0} 1/a(n) = (coth(Pi/sqrt(2))*Pi/sqrt(2) + tanh(Pi/sqrt(2))*Pi/(2*sqrt(2)) - 1)/2. (End)
E.g.f.: ((6 + 3*x + 2*x^2)*cosh(x) + (3 + 6*x + x^2)*sinh(x))/2. - Stefano Spezia, Oct 19 2024
MATHEMATICA
LinearRecurrence[{0, 3, 0, -3, 0, 1}, {3, 3, 11, 9, 27, 19}, 50] (* Amiram Eldar, Aug 09 2022 *)
PROG
(Magma) [((-1)^n+3)*(n^2+2*n+3)/4: n in [0..50]]; // Vincenzo Librandi, Aug 07 2011
(PARI) vector(100, n, n--; (1/4)*((-1)^n+3)*(n^2+2*n+3)) \\ Derek Orr, Jun 27 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paul Curtz, Aug 30 2009
STATUS
approved