login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A046176 Indices of square numbers which are also hexagonal. 13
1, 35, 1189, 40391, 1372105, 46611179, 1583407981, 53789260175, 1827251437969, 62072759630771, 2108646576008245, 71631910824649559, 2433376321462076761, 82663163018885960315, 2808114166320660573949 (list; graph; refs; listen; history; internal format)
OFFSET

1,2

COMMENTS

Bisection (even part) of Chebyshev sequence with Diophantine property.

(3*b(n))^2 - 2*(2*a(n+1))^2 = 1 with companion sequence b(n)=A077420(n),n>=0.

Sequence also refers to inradius of primitive Pythagorean triangles with consecutive legs, odd followed by even. - Lekraj Beedassy, Apr 23 2003

As n increases, this sequence is approximately geometric with common ratio r = lim(n -> Infinity, a(n)/a(n-1)) = ( 1+ sqrt(2))^4 = 17 + 12 * sqrt(2). - Ant King, Nov 08 2011

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 1..200

Tanya Khovanova, Recursive Sequences

Eric Weisstein's World of Mathematics, Link to a section of The World of Mathematics.

Index entries for sequences related to Chebyshev polynomials.

FORMULA

a(n) = 34*a(n-1)-a(n-2), a(0)=-1, a(1)=1.

a(n+1) = S(2*n, 6)= S(n, 34) + S(n-1, 34), n>=1, with S(n, x) := U(n, x/2) Chebyshev's polynomials of the second kind. See A049310. S(n, 34)=A029547(n).

G.f.: x*(1+x)/(1-34*x+x^2).

a(n+1) = sum((-1)^k*binomial(2*n-k, k)*6^(2*(n-k)), k=0..n), n>=0.

a(n) = A001109(2n+1). - Lekraj Beedassy, Apr 23 2003

Define f[x,s] = s x + Sqrt[(s^2-1)x^2+1]; f[0,s]=0. a(n) = f[f[a(n-1),3],3]. - Marcos Carreira, Dec 27 2006

Contribution from Antonio Alberto Olivares, Mar 22 2008: (Start)

a(n) = (sqrt(2)/8)(3+2*sqrt(2))*(17+12*sqrt(2))^(n-1) -(sqrt(2)/8)(3-2*sqrt(2))*(17-12*sqrt(2))^(n-1).

a(n) = (sqrt(2)/8)*(17+12*sqrt(2))^(n-1/2) -(sqrt(2)/8)*(17-12*sqrt(2))^(n-1/2).

a(n) = (sqrt(2)/8)*(3+2*sqrt(2))^(2n-1) -(sqrt(2)/8)*(3-2*sqrt(2))^(2n-1).

a(n) = (sqrt(2)/8)*(1+sqrt(2))^(4n-2) -(sqrt(2)/8)*(1-sqrt(2))^(4n-2).

a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3). (End)

a(n+1) = 17*a(n)+6*(8*a(n)^2+1)^0.5 for n>=0. - Richard Choulet (richardchoulet(AT)yahoo.fr), May 01 2009

a(n) = b such that (-1)^(n+1) * Integral_{x=-Pi/2..Pi/2} cos((2*n-1)*x)/(3-sin(x)) dx = c + b*ln(2). - Francesco Daddi, Aug 01 2011

MATHEMATICA

q=9; s=0; lst={}; Do[s+=n; If[Sqrt[q*s+1] == Floor[Sqrt[q*s+1]], AppendTo[lst, Sqrt[q*s+1]]], {n, 0, 9!}]; lst [From Vladimir Orlovsky, Apr 02 2009]

LinearRecurrence[{34, -1}, {1, 35}, 15] (* Ant King, Nov 08 2011 *)

PROG

(MAGMA) I:=[1, 35]; [n le 2 select I[n] else 34*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Nov 22 2011

CROSSREFS

Cf. A008844, A046177.

Cf. A001109.

Sequence in context: A049395 A115492 A187538 * A162847 A029546 A163218

Adjacent sequences:  A046173 A046174 A046175 * A046177 A046178 A046179

KEYWORD

nonn,easy

AUTHOR

Eric Weisstein (eric(AT)weisstein.com)

EXTENSIONS

Chebyshev comments from Wolfdieter Lang, Nov 29 2002

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 14 00:47 EST 2012. Contains 205567 sequences.