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!)
A106525 Values of x in x^2 - 49 = 2*y^2. 6

%I #35 Sep 08 2022 08:45:18

%S 9,11,21,43,57,119,249,331,693,1451,1929,4039,8457,11243,23541,49291,

%T 65529,137207,287289,381931,799701,1674443,2226057,4660999,9759369,

%U 12974411,27166293,56881771,75620409,158336759,331531257,440748043

%N Values of x in x^2 - 49 = 2*y^2.

%C From _Wolfdieter Lang_, Sep 27 2016: (Start)

%C These are the x members of all positive solutions (x(n), y(n)), proper and improper, of the Pell equation x^2 - 2*y^2 = 7^2.

%C The y(n) members are given in 2*A276600(n+2).

%C This sequence is composed of the two y members of the two proper classes of solutions of the Pell equation x^2 - 2*y^2 = 7^2 and of 7 times the proper solutions X of the Pell equation X^2 - 2*Y^2 = +1. See A275793, A275795 and 7*A001541. See A275793 for further information, and the Nagell reference. (End)

%C The sums of the consecutive integers in the following sequences will be squares: for n, i >= 1, if mod(i,3)=0 then 7*n+1, 7*n+2, ..., a(i)*n + (A001541(i/3)-1)/2; otherwise, if mod(i,3)=1 or 2 then 7*n+4, 7*n+5, ..., a(i)*n + (a(i)-1)/2.

%H Colin Barker, <a href="/A106525/b106525.txt">Table of n, a(n) for n = 1..1001</a> (offset adapted by _Georg Fischer_, Jan 31 2019).

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

%F a(3*k) = A001541(k)*7 for k >= 2.

%F a(3*k+1) = (A001541(k+2) + A001541(k-1) + A001541(k) - A001541(k+1))/2;

%F a(3*k+2) = (A001541(k+2) + A001541(k-1) - A001541(k) + A001541(k+1))/2.

%F a(3*n) = A275793(n), a(3*n+1) = A275795(n), a(3*n+2) = 7*A001541(n+1), n >= 0. - _Wolfdieter Lang_, Sep 27 2016

%F From _Colin Barker_, Mar 29 2012: (Start)

%F a(n) = 6*a(n-3) - a(n-6).

%F G.f.: x*(9 + 11*x + 21*x^2 - 11*x^3 - 9*x^4 - 7*x^5)/(1 - 6*x^3 + x^6). (End)

%e In the following, aa(n) denotes A001541(n):

%e a(9)=693; as mod(9,3)=0, a(9)=aa(3)*7=99*7=693, also 693^2-49=2*490^2

%e a(10)=1451; as mod(10,3)=1, a(10)=(aa(5)+aa(2)+aa(3)-aa(4))/2 =(3363+17+99-577)/2=1451, also 1451^2-49=2*1026^2.

%e The solutions (proper and every third pair improper) of x^2 - 2*y^2 = +49 begin [9, 4], [11, 6], [21, 14], [43, 30], [57, 40], [119, 84], [249, 176], [331, 234], [693, 490], [1451, 1026], [1929, 1364], [4039, 2856], [8457, 5980], [11243, 7950], [23541, 16646], ... - _Wolfdieter Lang_, Sep 27 2016

%t LinearRecurrence[{0,0,6,0,0,-1}, {9,11,21,43,57,119}, 50] (* _Vincenzo Librandi_, Oct 26 2018 *)

%o (PARI) Vec((9+11*x+21*x^2-11*x^3-9*x^4-7*x^5)/(1-6*x^3+x^6) + O(x^50)) \\ _Colin Barker_, Sep 28 2016

%o (Magma) I:=[9,11,21,43,57,119]; [n le 6 select I[n] else 6*Self(n-3)-Self(n-6): n in [1..40]]; // _Vincenzo Librandi_, Oct 26 2018

%o (Sage)

%o def A106525_list(prec):

%o P.<x> = PowerSeriesRing(ZZ, prec)

%o return P( x*(9+11*x+21*x^2-11*x^3-9*x^4-7*x^5)/(1-6*x^3+x^6) ).list()

%o a=A106525_list(41); a[1:] # _G. C. Greubel_, Sep 15 2021

%Y Cf. 14*A001109, 7*A001541, A106525, A106526.

%Y Cf. A275793, A275794, A275795, A275796, 2*A276600.

%K nonn,easy

%O 1,1

%A Andras Erszegi (erszegi.andras(AT)chello.hu), May 07 2005

%E Entry revised by _N. J. A. Sloane_, Oct 26 2018 at the suggestion of _Georg Fischer_.

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 29 04:23 EDT 2024. Contains 371264 sequences. (Running on oeis4.)