login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A253477
Indices of centered heptagonal numbers (A069099) which are also centered triangular numbers (A005448).
3
1, 10, 46, 1045, 5005, 114886, 550450, 12636361, 60544441, 1389884770, 6659338006, 152874688285, 732466636165, 16814825826526, 80564670640090, 1849477966229521, 8861381303773681, 203425761459420730, 974671378744464766, 22374984282570050725
OFFSET
1,2
COMMENTS
Also positive integers y in the solutions to 3*x^2 - 7*y^2 - 3*x + 7*y = 0, the corresponding values of x being A253476.
FORMULA
a(n) = a(n-1)+110*a(n-2)-110*a(n-3)-a(n-4)+a(n-5).
G.f.: -x*(x^4+9*x^3-74*x^2+9*x+1) / ((x-1)*(x^4-110*x^2+1)).
EXAMPLE
10 is in the sequence because the 10th centered heptagonal number is 316, which is also the 15th centered triangular number.
MATHEMATICA
LinearRecurrence[{1, 110, -110, -1, 1}, {1, 10, 46, 1045, 5005}, 30] (* Harvey P. Dale, Aug 13 2018 *)
PROG
(PARI) Vec(-x*(x^4+9*x^3-74*x^2+9*x+1)/((x-1)*(x^4-110*x^2+1)) + O(x^100))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Colin Barker, Jan 02 2015
STATUS
approved