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 A156772 #22 Sep 08 2022 08:45:41 %S A156772 531,1260,1989,2718,3447,4176,4905,5634,6363,7092,7821,8550,9279, %T A156772 10008,10737,11466,12195,12924,13653,14382,15111,15840,16569,17298, %U A156772 18027,18756,19485,20214,20943,21672,22401,23130,23859,24588,25317,26046 %N A156772 a(n) = 729*n - 198. %C A156772 The identity (6561*n^2 - 3564*n + 485)^2 - (81*n^2 - 44*n + 6)*(729*n - 198)^2 = 1 can be written as A156774(n)^2 - A156676(n)*a(n)^2 = 1. %H A156772 Vincenzo Librandi, <a href="/A156772/b156772.txt">Table of n, a(n) for n = 1..10000</a> %H A156772 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (2,-1). %F A156772 a(n) = 2*a(n-1) - a(n-2). %F A156772 G.f.: x*(531 + 198*x)/(1-x)^2. %F A156772 E.g.f.: 9*(22 - (22 - 81*x)*exp(x)). - _G. C. Greubel_, Jun 19 2021 %t A156772 LinearRecurrence[{2,-1}, {531,1260}, 40] %o A156772 (Magma) I:=[531, 1260]; [n le 2 select I[n] else 2*Self(n-1)-Self(n-2): n in [1..50]]; %o A156772 (PARI) a(n)=729*n-198 \\ _Charles R Greathouse IV_, Dec 23 2011 %o A156772 (Sage) [9*(81*n -22) for n in [1..50]] # _G. C. Greubel_, Jun 19 2021 %Y A156772 Cf. A156676, A156774. %K A156772 nonn,easy %O A156772 1,1 %A A156772 _Vincenzo Librandi_, Feb 15 2009