A157649 -id:A157649 - cp's OEIS frontend

A155923 Positive numbers y such that y^2 is of the form x^2+(x+17)^2 with integer x.

13, 17, 25, 53, 85, 137, 305, 493, 797, 1777, 2873, 4645, 10357, 16745, 27073, 60365, 97597, 157793, 351833, 568837, 919685, 2050633, 3315425, 5360317, 11951965, 19323713, 31242217, 69661157, 112626853, 182092985, 406014977, 656437405

Offset: 1

Klaus Brockhaus, Feb 09 2009

(-5,a(1)) and (A118120(n), a(n+1)) are solutions (x, y) to the Diophantine equation x^2+(x+17)^2 = y^2. (Offset 1 is assumed for A118120.)
lim_{n -> infinity} a(n)/a(n-3) = 3+2*sqrt(2).
lim_{n -> infinity} a(n)/a(n-1) = (19+6*sqrt(2))/17 for n mod 3 = {0, 2}.
lim_{n -> infinity} a(n)/a(n-1) = (387+182*sqrt(2))/17^2 for n mod 3 = 1.
For the generic case x^2+(x+p)^2=y^2 with p=2*m^2-1 a prime number in A066436, m>=2, the x values are given by the sequence defined by: a(n)=6*a(n-3)-a(n-6)+2p with a(1)=0, a(2)=2m+1, a(3)=6m^2-10m+4, a(4)=3p, a(5)=6m^2+10m+4, a(6)=40m^2-58m+21.Y values are given by the sequence defined by: b(n)=6*b(n-3)-b(n-6) with b(1)=p, b(2)=2*m^2+2m+1, b(3)=10m^2-14m+5, b(4)=5p, b(5)=10m^2+14m+5, b(6)=58m^2-82m+29. [From Mohamed Bouhamida, Sep 09 2009]

			(-5,a(1)) = (-5,13) is a solution: (-5)^2+(-5+17)^2 = 25+144 = 169 = 13^2;
(A118120(1), a(2)) = (0, 17) is a solution: 0^2+(0+17)^2 = 289 = 17^2;
(A118120(2), a(3)) = (7, 25) is a solution: 7^2+(7+17)^2 = 49+576 = 625 = 25^2.

Index entries for linear recurrences with constant coefficients, signature (0,0,6,0,0,-1).

Cf. A118120, A156035 (decimal expansion of 3+2*sqrt(2)), A156163 (decimal expansion of (19+6*sqrt(2))/17), A157649 (decimal expansion of (387+182*sqrt(2))/17^2).

Cf. A156156 (first trisection), A156157 (second trisection), A156158 (third trisection).

LinearRecurrence[{0,0,6,0,0,-1},{13,17,25,53,85,137},50] (* Harvey P. Dale, Feb 11 2015 *)

{forstep(n=-5, 660000000, [1,3], if(issquare(2*n*(n+17)+289, &k), print1(k, ",")))}

a(n) = 6*a(n-3)-a(n-6) for n > 6; a(1) = 13, a(2) = 17, a(3) = 25, a(4) = 53, a(5) = 85, a(6) = 137.

G.f.: x*(1-x)*(13+30*x+55*x^2+30*x^3+13*x^4)/(1-6*x^3+x^6).

G.f. corrected, first and fourth comment and examples edited, cross-reference added by Klaus Brockhaus, Sep 22 2009

A156159 Squares of the form k^2+(k+17)^2 with integer k.

Original entry on oeis.org

169, 289, 625, 2809, 7225, 18769, 93025, 243049, 635209, 3157729, 8254129, 21576025, 107267449, 280395025, 732947329, 3643933225, 9525174409, 24898630849, 123786459889, 323575532569, 845820499225, 4205095700689

Offset: 1

Klaus Brockhaus, Feb 09 2009

Square roots of k^2+(k+17)^2 are in A155923, values k (except for -5) are in A118120.

			625 = 25^2 is of the form k^2+(k+17)^2 with k = 7: 7^2+24^2 = 625. Hence 625 is in the sequence.

Index entries for linear recurrences with constant coefficients, signature (1,0,34,-34,0,-1,1).

Equals A155923^2. Cf. A156160 (first trisection), A156161 (second trisection), A156162 (third trisection).

Cf. A118120, A156035 (decimal expansion of 3+2*sqrt(2)), A156164 (decimal expansion of 17+12*sqrt(2)), A156163 (decimal expansion of (19+6*sqrt(2))/17), A157649 (decimal expansion of (387+182*sqrt(2))/17^2).

LinearRecurrence[{1,0,34,-34,0,-1,1},{169,289,625,2809,7225,18769,93025},30] (* Harvey P. Dale, Apr 22 2022 *)

{forstep(n=-5, 1600000, [1, 3], if(issquare(a=2*n*(n+17)+289), print1(a, ",")))}

a(n) = 34*a(n-3)-a(n-6)-2312 for n > 6; a(1)=169, a(2)=289, a(3)=625, a(4)=2809, a(5)=7225, a(6)=18769.
G.f.: x*(169+120*x+336*x^2-3562*x^3+336*x^4+120*x^5+169*x^6)/((1-x)*(1-34*x^3+x^6)).
Limit_{n -> oo} a(n)/a(n-3) = (17+12*sqrt(2)).
Limit_{n -> oo} a(n)/a(n-1) = ((19+6*sqrt(2))/17)^2 for n mod 3 = {0, 2}.
Limit_{n -> oo} a(n)/a(n-1) = ((387+182*sqrt(2))/17^2)^2 for n mod 3 = 1.

G.f. corrected, fourth comment and cross-references edited by Klaus Brockhaus, Sep 23 2009

cp's OEIS Frontend

A155923 Positive numbers y such that y^2 is of the form x^2+(x+17)^2 with integer x.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions