This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A087125 #91 Feb 16 2025 08:32:50
%S A087125 0,5,54,539,5340,52865,523314,5180279,51279480,507614525,5024865774,
%T A087125 49741043219,492385566420,4874114620985,48248760643434,
%U A087125 477613491813359,4727886157490160,46801248083088245,463284594673392294,4586044698650834699,45397162391834954700
%N A087125 Indices k of hex numbers H(k) that are also triangular.
%C A087125 From the law of cosines, the non-Pythagorean triple {a(n), a(n)+1=A253475(n+1), A072256(n+1)} forms a near-isosceles triangle with the angle bounded by the consecutive sides equal to the regular tetrahedron's central angle (see A156546 and A247412). This implies also that a(n) are those numbers k such that (16/3)*A000217(k)+1 is a perfect square. - _Federico Provvedi_, Apr 04 2023
%H A087125 Colin Barker, <a href="/A087125/b087125.txt">Table of n, a(n) for n = 0..1000</a>
%H A087125 Editors, L'Intermédiaire des Mathématiciens, <a href="/A072256/a072256.pdf">Query 4500: The equation x(x+1)/2 = y*(y+1)/3</a>, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (I).
%H A087125 Editors, L'Intermédiaire des Mathématiciens, <a href="/A072256/a072256_1.pdf">Query 4500: The equation x(x+1)/2 = y*(y+1)/3</a>, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (II).
%H A087125 Editors, L'Intermédiaire des Mathématiciens, <a href="/A072256/a072256_2.pdf">Query 4500: The equation x(x+1)/2 = y*(y+1)/3</a>, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (III).
%H A087125 Editors, L'Intermédiaire des Mathématiciens, <a href="/A072256/a072256_3.pdf">Query 4500: The equation x(x+1)/2 = y*(y+1)/3</a>, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (IV).
%H A087125 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/HexNumber.html">Hex Number</a>
%H A087125 <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (11,-11,1).
%F A087125 G.f.: (-x^2+5*x)/((1-x)*(1-10*x+x^2)).
%F A087125 a(n) = 11*a(n-1) - 11*a(n-2) + a(n-3) for n > 2. - _Colin Barker_, Jun 23 2015
%F A087125 a(n) = (-4 - (5-2*sqrt(6))^n*(-2 + sqrt(6)) + (2+sqrt(6))*(5+2*sqrt(6))^n)/8. - _Colin Barker_, Mar 05 2016
%F A087125 a(n) = 10*a(n-1) - a(n-2) + 4 for n > 1. - _Charlie Marion_, Feb 14 2023
%F A087125 a(n) = ((x^(n+1)+1)*(x^n-1))/(2*x^n*(x-1)), with x=5+2*sqrt(6). - _Federico Provvedi_, Apr 04 2023
%F A087125 a(n) = sqrt(3*A161680(A054318(n+1)) + 1/4) - 1/2 = floor(sqrt(3*A000217(A054318(n+1)-1) + 1/4)). - _Federico Provvedi_, Apr 16 2023
%t A087125 CoefficientList[Series[(-x^2+5*x)/((1-x)*(1-10*x+x^2)), {x, 0, 25}], x] (* _G. C. Greubel_, Nov 04 2017 *)
%t A087125 LinearRecurrence[{11,-11,1},{0,5,54},30] (* _Harvey P. Dale_, Jun 14 2022 *)
%t A087125 Table[(x Sqrt[z^(2 n + 1) + z^-(2 n + 1) - 2] - 4) / 8 //. {x -> Sqrt[2], y -> Sqrt[3], z -> (5 + 2 x y)}, {n, 0, 100}] // Round (* _Federico Provvedi_, Apr 16 2023 *)
%o A087125 (PARI) concat(0, Vec(x*(x-5)/((x-1)*(x^2-10*x+1)) + O(x^50))) \\ _Colin Barker_, Jun 23 2015
%o A087125 (Magma) [Round((-4-(5-2*Sqrt(6))^n*(-2+Sqrt(6)) + (2+Sqrt(6))*(5 + 2*Sqrt(6))^n)/8): n in [0..25]]; // _G. C. Greubel_, Nov 04 2017
%Y A087125 Cf. A006244, A031138, A072256, A054318.
%K A087125 nonn,easy
%O A087125 0,2
%A A087125 _Eric W. Weisstein_, Aug 14 2003