Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A002315
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A002315 NSW numbers: a(n) = 6*a(n-1) - a(n-2); also a(n)^2 - 2*b(n)^2 = -1 with b(n)=A001653(n).
(Formerly M4423 N1869)
+0
74
1, 7, 41, 239, 1393, 8119, 47321, 275807, 1607521, 9369319, 54608393, 318281039, 1855077841, 10812186007, 63018038201, 367296043199, 2140758220993, 12477253282759, 72722761475561, 423859315570607, 2470433131948081 (list; graph; listen)
OFFSET

0,2

COMMENT

Named after the Newman-Shanks-Williams reference.

Also numbers n such that A125650[ 3*n^2 ] is an odd perfect square, where A124650(n) is a numerator of n(n+3)/(4(n+1)(n+2)) = Sum[ 1/(k(k+1)(k+2)), {k,1,n} ]. Sequence of numbers 3*n^2 is a bisection of A125651(n). - Alexander Adamchuk (alex(AT)kolmogorov.com), Nov 30 2006

For positive n, a(n) corresponds to the sum of legs of near-isosceles primitive Pythagorean triangles (with consecutive legs). - Lekraj Beedassy (blekraj(AT)yahoo.com), Feb 06 2007

Also numbers n such that n^2 is a centered 16-gonal number; or a number of the form 8k(k+1)+1, where k = A053141(n) = {0, 2, 14, 84, 492, 2870, ...}. - Alexander Adamchuk (alex(AT)kolmogorov.com), Apr 21 2007

A002315(n)=A001333(2*n+1) [From Ctibor O. ZIZKA (ctibor.zizka(AT)seznam.cz), Aug 13 2008]

The lower principal convergents to 2^(1/2), beginning with 1/1, 7/5, 41/29, 239/169, comprise a strictly increasing sequence; numerators=A002315 and denominators=A001653. - Clark Kimberling (ck6(AT)evansville.edu), Aug 27 2008

The upper intermediate convergents to 2^(1/2) beginning with 10/7, 58/41, 338/239, 1970/1393 form a strictly decreasing sequence; essentially, numerators=A075870, denominators=A002315. - Clark Kimberling (ck6(AT)evansville.edu), Aug 27 2008

General recurrence is a(n)=(a(1)-1)*a(n-1)-a(n-2), a(1)>=4, lim n->infinity a(n)= x*(k*x+1)^n, k =(a(1)-3), x=(1+sqrt((a(1)+1)/(a(1)-3)))/2. Examples in OEIS: a(1)=4 gives A002878, primes in it A121534. a(1)=5 gives A001834, primes in it A086386. a(1)=6 gives A030221, primes in it not in OEIS {29,139,3191,...}. a(1)=7 gives A002315, primes in it A088165. a(1)=8 gives A033890, primes in it not in OEIS (does there exist any ?). a(1)=9 gives A057080, primes in it not in OEIS {71,34649,16908641,...}. a(1)=10 gives A057081, primes in it not in OEIS {389806471,192097408520951,...}. [From Ctibor O. Zizka (ctibor.zizka(AT)seznam.cz), Sep 02 2008]

Numbers n such that (ceiling(sqrt(n*n/2)))^2 = (1+n*n)/2 [From Ctibor O. Zizka (c.zizka(AT)email.cz), Nov 09 2009]

REFERENCES

E. Barcucci et al., A combinatorial interpretation of the recurrence f_{n+1} = 6 f_n - f_{n-1}, Discrete Math., 190 (1998), 235-240.

A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 256.

J. Bonin, L. Shapiro and R. Simion, Some q-analogues of the Schroeder numbers arising from combinatorial statistics on lattice paths, H. Statistical Planning and Inference, 16,1993,35-55 (p. 50).

A. S. Fraenkel, Recent results and questions in combinatorial game complexities, Theoretical Computer Science, vol. 249, no. 2 (2000), 265-288.

D. H. Lehmer, Lacunary recurrence formulas for the numbers of Bernoulli and Euler, Annals Math., 36 (1935), 637-649.

M. Newman, D. Shanks and H. C. Williams, Simple groups of square order and an interesting sequence of primes, Acta Arith. 38 (1980/81), no. 2, 129-140. MR82b:20022

Problem 47, Amer. Math. Monthly, 4 (1897), 25-28.

P. Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 288.

Santana, S. F. and Diaz-Barrero, J. L. (2006). Some properties of sums involving Pell numbers. Missouri Journal of Mathematical Sciences 18(1).

N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.

R. A. Sulanke, Bijective recurrences concerning Schroeder paths, Electron. J. Combin. 5 (1998), Research Paper 47, 11 pp.

P.-F. Teilhet, Reply to Query 2094, L'Interm\'{e}diaire des Math\'{e}maticiens, 10 (1903), 235-238.

LINKS

T. D. Noe, Table of n, a(n) for n = 0..200

Index entries for two-way infinite sequences

Index entries for sequences related to linear recurrences with constant coefficients

S. Plouffe, Approximations de S\'{e}ries G\'{e}n\'{e}ratrices et Quelques Conjectures, Dissertation, Universit\'{e} du Qu\'{e}bec \`{a} Montr\'{e}al, 1992.

S. Plouffe, 1031 Generating Functions and Conjectures, Universit\'{e} du Qu\'{e}bec \`{a} Montr\'{e}al, 1992.

Enrica Duchi, Andrea Frosini, Renzo Pinzani and Simone Rinaldi, A Note on Rational Succession Rules, J. Integer Seqs., Vol. 6, 2003.

A. S. Fraenkel, Arrays, numeration systems and Frankenstein games, Theoret. Comput. Sci. 282 (2002), 271-284.

Tanya Khovanova, Recursive Sequences

The Prime Glossary, NSW number.

R. A. Sulanke, Moments of generalized Motzkin paths, J. Integer Sequences, Vol. 3 (2000), #00.1.

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

Eric Weisstein, Link to a section of The World of Mathematics, Centered Polygonal Number.

Index entries for sequences related to Chebyshev polynomials.

Missouri J. Mathematical Sciences (New web address)

FORMULA

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

a(n)=(1+sqrt(2))/2*(3+sqrt(8))^n+(1-sqrt(2))/2*(3-sqrt(8))^n. - Ralf Stephan, Feb 23 2003

a(n) = sqrt(2*(A001653(n))^2-1)

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

a(n) = S(n, 6)+S(n-1, 6) = S(2*n, sqrt(8)), S(n, x) = U(n, x/2) are Chebyshev's polynomials of the 2nd kind. Cf. A049310. S(n, 6)= A001109(n+1).

a(n) ~ 1/2*(sqrt(2) + 1)^(2*n+1) - Joe Keane (jgk(AT)jgk.org), May 15 2002

Lim n -> inf. a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson (omomom(AT)hotmail.com), Oct 06 2002

Let q(n, x)=sum(i=0, n, x^(n-i)*binomial(2*n-i, i)); then (-1)^n*q(n, -8)=a(n) - Benoit Cloitre (benoit7848c(AT)orange.fr), Nov 10 2002

With a=3+2sqrt(2), b=3-2sqrt(2): a(n)=(a^((2n+1)/2)-b^((2n+1)/2))/2. a(n)=A077444(n)/2. - Mario Catalani (mario.catalani(AT)unito.it), Mar 31 2003

a(n)=sum(k=0, n, 2^k*binomial(2*n+1, 2*k)) - Zoltan Zachar (zachar(AT)fellner.sulinet.hu), Oct 08 2003

Same as: i such that Mod(sigma(i^2+1, 2), 2)=1 - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 26 2004

a(n) = L(n, -6)*(-1)^n, where L is defined as in A108299; see also A001653 for L(n, +6). - Reinhard Zumkeller (reinhard.zumkeller(AT)gmail.com), Jun 01 2005

a(n)=A001652(n)+A046090(n); e.g. 239=119+120 - Charlie Marion (charliem(AT)bestweb.net), Nov 20 2003

A001541(n)*a(n+k)=A001652(2n+k)+A001652(k)+1; e.g. 3*1393=4069+119+1; for k>0, A001541(n+k)*a(n)=A001652(2n+k)-A001652(k-1); e.g. 99*7=696-3 - Charlie Marion (charliemath(AT)verizon.net), Mar 17 2003

a(n)=Jacobi_P(n,1/2,-1/2,3)/Jacobi_P(n,-1/2,1/2,1); - Paul Barry (pbarry(AT)wit.ie), Feb 03 2006

P_{2n}+P_{2n+1} where P_i are the Pell numbers (A000129). Also the square root of the partial sums of Pell numbers: P_{2n}+P_{2n+1} = sqrt(sum_{i=0}^{4n+1} P_i) (Santana and Diaz-Barrero, 2006). - David Eppstein (eppstein(AT)ics.uci.edu), Jan 28 2007

a(n) = 2*A001652(n) + 1 = 2*A046729(n) + (-1)^n. - Lekraj Beedassy (blekraj(AT)yahoo.com), Feb 06 2007

a(n) = sqrt[A001108(2*n+1)] - Anton Vrba (antonvrba(AT)yahoo.com), Feb 14 2007

a(n) = Sqrt[ 8*A053141(n)*(A053141(n) + 1) + 1 ]. - Alexander Adamchuk (alex(AT)kolmogorov.com), Apr 21 2007

a(n+1) = 3*a(n)+(8*a(n)^2+8)^0.5, a(1)=1. - Richard Choulet (richardchoulet(AT)yahoo.fr), Sep 18 2007

a(n)= third binomial transform of 1,4,8,32,64,256,512 [From Al Hakanson (hawkuu(AT)gmail.com), Aug 15 2009]

MAPLE

a[0]:=1: a[1]:=7: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..20); - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Jul 26 2006

A002315:=(1+z)/(1-6*z+z**2); [S. Plouffe in his 1992 dissertation.]

MATHEMATICA

a[0] = 1; a[1] = 7; a[n_] := a[n] = 6a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 20}] (from Robert G. Wilson v Jun 09 2004)

q=16; s=0; lst={}; Do[s+=n; If[Sqrt[q*s+1]==Floor[Sqrt[q*s+1]], AppendTo[lst, Sqrt[q*s+1]]], {n, 0, 8!}]; lst [From Vladimir Orlovsky (4vladimir(AT)gmail.com), Apr 02 2009]

PROGRAM

(PARI) a(n)=subst(poltchebi(abs(n+1))-poltchebi(abs(n)), x, 3)/2

(PARI) a(n)=if(n<0, -a(-1-n), polsym(x^2-2*x-1, 2*n+1)[2*n+2]/2)

(PARI) a(n)=local(w=3+quadgen(32)); imag((1+w)*w^n)

(PARI) for (i=1, 10000, if(Mod(sigma(i^2+1, 2), 2)==1, print1(i, ", ")))

(Other) sage: [(lucas_number2(n, 6, 1)-lucas_number2(n-1, 6, 1))/4 for n in xrange(1, 22)]# [From Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Nov 10 2009]

CROSSREFS

Bisection of A001333. Cf. A001109, A001653. A065513(n)=a(n)-1.

First differences of A001108 and A055997. Bisection of A084068 and A088014. Pairwise sums of A001109. Cf. A077444.

Cf. A125650, A125651, A125652.

Row sums of unsigned triangle A127675.

Cf. A053141.

Sequence in context: A026002 A057009 A140480 this_sequence A141813 A088165 A108983

Adjacent sequences: A002312 A002313 A002314 this_sequence A002316 A002317 A002318

KEYWORD

nonn,easy,nice,new

AUTHOR

N. J. A. Sloane (njas(AT)research.att.com).

EXTENSIONS

More terms from James A. Sellers (sellersj(AT)math.psu.edu), Feb 16 2000

Removed attribute "conjectured" from Plouffe g.f R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Mar 11 2009

page 1

Search completed in 0.007 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 20:09 EST 2009. Contains 167514 sequences.


AT&T Labs Research