login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A106328 Numbers j such that 8*(j^2) + 9 = k^2 for some positive number k. 13

%I #58 Nov 26 2022 02:45:49

%S 0,3,18,105,612,3567,20790,121173,706248,4116315,23991642,139833537,

%T 815009580,4750223943,27686334078,161367780525,940520349072,

%U 5481754313907,31950005534370,186218278892313,1085359667819508,6325939728024735,36870278700328902,214895732473948677

%N Numbers j such that 8*(j^2) + 9 = k^2 for some positive number k.

%C The ratio k(n) /(2*j(n)) tends to sqrt(2) as n increases.

%C The squares of the numbers in this sequence are one less than a triangular number: a(n)^2 = A164080(n). For example, 18^2 is 324, and 325 is a triangular number. a(n)^2 + 1 = A164055(n). a(n)^2 = A072221(n)(A072221(n)+1)/2 - 1. - _Tanya Khovanova_ & Alexey Radul, Aug 09 2009

%C For n > 0, a(n+1) is the n-th almost balancing number of first type (see Tekcan and Erdem). - _Stefano Spezia_, Nov 25 2022

%H Reinhard Zumkeller, <a href="/A106328/b106328.txt">Table of n, a(n) for n = 1..1000</a>

%H Tanya Khovanova, <a href="http://www.tanyakhovanova.com/RecursiveSequences/RecursiveSequences.html">Recursive Sequences</a>

%H Soumeya M. Tebtoub, Hacène Belbachir, and László Németh, <a href="https://hal.archives-ouvertes.fr/hal-02918958/document#page=18">Integer sequences and ellipse chains inside a hyperbola</a>, Proceedings of the 1st International Conference on Algebras, Graphs and Ordered Sets (ALGOS 2020), hal-02918958 [math.cs], 17-18.

%H Ahmet Tekcan and Alper Erdem, <a href="https://arxiv.org/abs/2211.08907">General Terms of All Almost Balancing Numbers of First and Second Type</a>, arXiv:2211.08907 [math.NT], 2022.

%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (6,-1).

%F a(1)=0, a(2)=3 then a(n) = 6*a(n-1) - a(n-2).

%F a(n) = ((3+2*sqrt(2))^(n-1) - (3-2*sqrt(2))^(n-1))*3/4/sqrt(2). - _Max Alekseyev_, Jan 11 2007

%F a(n) = 3*A001109(n). - _M. F. Hasler_, _R. J. Mathar_, Jun 03 2009

%F a(n) = (3/4)*A005319(n-1).

%F G.f.: 3*x^2/(1 - 6*x + x^2). - _Philippe Deléham_, Nov 17 2008

%F E.g.f.: 3 - 3*exp(3*x)*(4*cosh(2*sqrt(2)*x) - 3*sqrt(2)*sinh(2*sqrt(2)*x))/4. - _Stefano Spezia_, Nov 25 2022

%t s=0;lst={};Do[s+=n;If[Sqrt[s-1]==Floor[Sqrt[s-1]],AppendTo[lst,Sqrt[s-1]]],{n,8!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Apr 02 2009 *)

%t Rest@ CoefficientList[Series[3 x^2/(1 - 6 x + x^2), {x, 0, 24}], x] (* _Michael De Vlieger_, Nov 02 2020 *)

%o (Haskell)

%o a106328 n = a106328_list !! (n-1)

%o a106328_list = 0 : 3 : zipWith (-) (map (* 6) (tail a106328_list)) a106328_list

%o -- _Reinhard Zumkeller_, Jan 10 2012

%o (PARI) concat(0, Vec(3*x^2/(1-6*x+x^2) + O(x^40))) \\ _Michel Marcus_, Sep 07 2016

%o (PARI) a(n)=([0,1;-1,6]^n*[-3;0])[1,1] \\ _Charles R Greathouse IV_, Sep 07 2016

%Y Cf. A001109, A005319, A072221, A103328, A164080, A164055.

%K nonn,easy

%O 1,2

%A _Pierre CAMI_, Apr 29 2005

%E More terms from _Max Alekseyev_, Jan 11 2007

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 18 22:56 EDT 2024. Contains 370952 sequences. (Running on oeis4.)