|
|
A157716
|
|
One-eighth of triangular numbers (integers only).
|
|
2
|
|
|
0, 15, 17, 62, 66, 141, 147, 252, 260, 395, 405, 570, 582, 777, 791, 1016, 1032, 1287, 1305, 1590, 1610, 1925, 1947, 2292, 2316, 2691, 2717, 3122, 3150, 3585, 3615, 4080, 4112, 4607, 4641, 5166, 5202, 5757, 5795, 6380, 6420, 7035, 7077, 7722, 7766, 8441
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Also numbers k such that 1+64*k is a square. (End)
|
|
LINKS
|
|
|
FORMULA
|
G.f.: x^2*(15+2*x+15*x^2)/((1+x)^2*(1-x)^3 ). [Maksym Voznyy (voznyy(AT)mail.ru), Jul 26 2009; checked and corrected by R. J. Mathar, Sep 16 2009]
a(n) = (2*n-1 + 7/8*(-1)^n)^2 -1/64. - Robert Israel, Apr 20 2014
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5). - Wesley Ivan Hurt, Nov 10 2020
Sum_{n>=2} 1/a(n) = 16 - (sqrt(2*(2+sqrt(2))) + sqrt(2) + 1)*Pi. - Amiram Eldar, Mar 17 2022
|
|
EXAMPLE
|
|
|
MAPLE
|
seq((2*n-1 + 7/8*(-1)^n)^2 - 1/64, n = 1 .. 1000); # Robert Israel, Apr 20 2014
|
|
MATHEMATICA
|
Array[(2 # - 1 + 7/8*(-1)^#)^2 - 1/64 &, 46] (* or *)
Rest@ CoefficientList[Series[x^2*(15 + 2 x + 15 x^2)/((1 + x)^2*(1 - x)^3), {x, 0, 46}], x] (* Michael De Vlieger, Nov 05 2020 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|