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!)
A001080 a(n) = 16*a(n-1) - a(n-2) with a(0) = 0, a(1) = 3.
(Formerly M3155 N1278)
9

%I M3155 N1278 #95 Dec 27 2023 10:00:49

%S 0,3,48,765,12192,194307,3096720,49353213,786554688,12535521795,

%T 199781794032,3183973182717,50743789129440,808716652888323,

%U 12888722657083728,205410845860451325,3273684811110137472,52173546131901748227,831503053299317834160

%N a(n) = 16*a(n-1) - a(n-2) with a(0) = 0, a(1) = 3.

%C Also 7*x^2 + 1 is a square; n=7 in PARI script below. - _Cino Hilliard_, Mar 08 2003

%C That is, the terms are solutions y of the Pell-Fermat equation x^2 - 7 * y^2 = 1. The corresponding values of x are in A001081. (x,y) = (1,0), (8,3), (127,48), ... - _Bernard Schott_, Feb 23 2019

%C The first solution to the equation x^2 - 7*y^2 = 1 is (X(0); Y(0)) = (1; 0) and the other solutions are defined by: (X(n); Y(n))= (8*X(n-1) + 21*Y(n-1); 3*X(n-1) + 8*Y(n-1)), with n >= 1. - _Mohamed Bouhamida_, Jan 16 2020

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

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

%D V. Thébault, Les Récréations Mathématiques. Gauthier-Villars, Paris, 1952, p. 281.

%H T. D. Noe, <a href="/A001080/b001080.txt">Table of n, a(n) for n = 0..200</a>

%H H. Brocard, <a href="http://resolver.sub.uni-goettingen.de/purl?PPN598948236_0004/DMDLOG_0053">Notes élémentaires sur le problème de Peel [sic]</a>, Nouvelle Correspondance Mathématique, 4 (1878), 337-343.

%H M. Davis, <a href="http://www.rand.org/pubs/research_memoranda/RM5494.html">One equation to rule them all</a>, Trans. New York Acad. Sci. Ser. II, 30 (1968), 766-773.

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

%H Simon Plouffe, <a href="https://arxiv.org/abs/0911.4975">Approximations de séries génératrices et quelques conjectures</a>, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.

%H Simon Plouffe, <a href="/A000051/a000051_2.pdf">1031 Generating Functions</a>, Appendix to Thesis, Montreal, 1992

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

%F G.f.: 3*x/(1-16*x+x^2).

%F From _Mohamed Bouhamida_, Sep 20 2006: (Start)

%F a(n) = 15*(a(n-1) + a(n-2)) - a(n-3).

%F a(n) = 17*(a(n-1) - a(n-2)) + a(n-3). (End)

%F a(n) = 16*a(n-1) - a(n-2) with a(1)=0 and a(2)=3. - _Sture Sjöstedt_, Nov 18 2011

%F E.g.f.: exp(8*x)*sinh(3*sqrt(7)*x)/sqrt(7). - _G. C. Greubel_, Feb 23 2019

%p A001080:=3*z/(1-16*z+z**2); # conjectured (correctly) by _Simon Plouffe_ in his 1992 dissertation

%t LinearRecurrence[{16,-1},{0,3},30] (* _Harvey P. Dale_, Nov 01 2011 *)

%t CoefficientList[Series[3*x/(1-16*x+x^2), {x, 0, 30}], x] (* _G. C. Greubel_, Dec 20 2017 *)

%o (PARI) nxsqp1(m,n) = { for(x=1,m, y = n*x*x+1; if(issquare(y),print1(x" ")) ) }

%o (PARI) x='x+O('x^30); concat([0], Vec(3*x/(1-16*x+x^2))) \\ _G. C. Greubel_, Dec 20 2017

%o (Magma) I:=[0,3]; [n le 2 select I[n] else 16*Self(n-1) - Self(n-2): n in [1..30]]; (* _G. C. Greubel_, Dec 20 2017 *)

%o (Sage) (3*x/(1-16*x+x^2)).series(x, 30).coefficients(x, sparse=False) # _G. C. Greubel_, Feb 23 2019

%o (GAP) a:=[0,3];; for n in [3..30] do a[n]:=16*a[n-1]-a[n-2]; od; a; # _G. C. Greubel_, Feb 23 2019

%Y Equals 3 * A077412. Bisection of A084069.

%Y Cf. A048907.

%Y Cf. A001081, A010727. - _Vincenzo Librandi_, Feb 16 2009

%K nonn,easy

%O 0,2

%A _N. J. A. Sloane_

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 28 04:13 EDT 2024. Contains 371235 sequences. (Running on oeis4.)