|
|
A253458
|
|
Indices of centered heptagonal numbers (A069099) which are also centered hexagonal numbers (A003215).
|
|
3
|
|
|
1, 13, 325, 8425, 218713, 5678101, 147411901, 3827031313, 99355402225, 2579413426525, 66965393687413, 1738520822446201, 45134575989913801, 1171760454915312613, 30420637251808214125, 789764808092098254625, 20503464373142746406113, 532300308893619308304301
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Also positive integers y in the solutions to 6*x^2 - 7*y^2 - 6*x + 7*y = 0, the corresponding values of x being A253457.
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 27*a(n-1)-27*a(n-2)+a(n-3).
G.f.: -x*(x^2-14*x+1) / ((x-1)*(x^2-26*x+1)).
a(n) = 1/2+(13+2*sqrt(42))^(-n)*(7+sqrt(42)-(-7+sqrt(42))*(13+2*sqrt(42))^(2*n))/28. - Colin Barker, Mar 03 2016
|
|
EXAMPLE
|
13 is in the sequence because the 13th centered heptagonal number is 547, which is also the 14th centered hexagonal number.
|
|
MATHEMATICA
|
LinearRecurrence[{27, -27, 1}, {1, 13, 325}, 20] (* Harvey P. Dale, Oct 13 2022 *)
|
|
PROG
|
(PARI) Vec(-x*(x^2-14*x+1)/((x-1)*(x^2-26*x+1)) + O(x^100))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|