OFFSET
1,2
LINKS
Colin Barker, Table of n, a(n) for n = 1..1000
Ana Rechtman, Mai 2017, 2e défi, Images des Mathématiques, CNRS, 2017 (in French).
Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).
FORMULA
a(n) = 72*(n-1)^6. - David A. Corneth, May 16 2017
O.g.f.: 72*x^2*(1 + x)*(1 + 56*x + 246*x^2 + 56*x^3 + x^4) / (1 - x)^7. - Colin Barker, May 17 2017
E.g.f.: 72*(-1 + (1 - x + x^2 + 10*x^3 + 20*x^4 + 9*x^5 + x^6)*exp(x)). - Bruno Berselli, May 17 2017
EXAMPLE
From Michael De Vlieger, May 16 2017: (Start)
72 is a term because 2*72 = 144 = 12^2 and 3*72 = 216 = 6^3.4608 is a term because 2*4608 = 96^2 and 3*4608 = 24^3. (End)
MATHEMATICA
Array[72 (# - 1)^6 &, 26] (* Michael De Vlieger, May 16 2017 *)
LinearRecurrence[{7, -21, 35, -35, 21, -7, 1}, {0, 72, 4608, 52488, 294912, 1125000, 3359232}, 30] (* Harvey P. Dale, May 07 2022 *)
PROG
(PARI) isok(x) = issquare(2*x) && ispower(3*x, 3);
(PARI) concat(0, Vec(72*x^2*(1 + x)*(1 + 56*x + 246*x^2 + 56*x^3 + x^4) / (1 - x)^7 + O(x^30))) \\ Colin Barker, May 17 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michel Marcus, May 16 2017
EXTENSIONS
More terms from Michael De Vlieger, May 16 2017
STATUS
approved