A129642 -id:A129642 - cp's OEIS frontend

A129725 Nonnegative values x of solutions (x, y) to the Diophantine equation x^2+(x+521)^2 = y^2.

0, 100, 1159, 1563, 2079, 8080, 10420, 13416, 48363, 61999, 79459, 283140, 362616, 464380, 1651519, 2114739, 2707863, 9627016, 12326860, 15783840, 56111619, 71847463, 91996219, 327043740, 418758960, 536194516, 1906151863, 2440707339, 3125171919, 11109868480

Offset: 1

Mohamed Bouhamida, Jun 02 2007

Also values x of Pythagorean triples (x, x+521, y).
Corresponding values y of solutions (x, y) are in A160583.
lim_{n -> infinity} a(n)/a(n-3) = 3+2*sqrt(2).
lim_{n -> infinity} a(n)/a(n-1) = (537+92*sqrt(2))/521 for n mod 3 = {1, 2}.
lim_{n -> infinity} a(n)/a(n-1) = (520659+314170*sqrt(2))/521^2 for n mod 3 = 0.

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

Cf. A160583, A001652, A129642, A156035 (decimal expansion of 3+2*sqrt(2)), A160584 (decimal expansion of (537+92*sqrt(2))/521), A160585 (decimal expansion of (520659+314170*sqrt(2))/521^2).

LinearRecurrence[{1, 0, 6, -6, 0, -1, 1}, {0, 100, 1159, 1563, 2079, 8080, 10420}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 13 2012 *)

{forstep(n=0, 10000000, [3, 1], if(issquare(2*n^2+1042*n+271441), print1(n, ",")))}

a(n) = 6*a(n-3)-a(n-6)+1042 for n > 6; a(1)=0, a(2)=100, a(3)=1159, a(4)=1563, a(5)=2079, a(6)=8080.
G.f.: x*(100+1059*x+404*x^2-84*x^3-353*x^4-84*x^5) / ((1-x)*(1-6*x^3+x^6)).
a(3*k+1) = 521*A001652(k) for k >= 0.

Edited and two terms added by Klaus Brockhaus, Jun 08 2009

A160581 Decimal expansion of (601+276*sqrt(2))/457.

Original entry on oeis.org

2, 1, 6, 9, 1, 9, 6, 8, 1, 2, 2, 8, 6, 5, 9, 5, 6, 9, 6, 8, 6, 9, 2, 9, 1, 2, 2, 0, 7, 4, 8, 0, 8, 8, 9, 9, 2, 7, 4, 6, 6, 7, 2, 7, 2, 9, 9, 8, 6, 9, 5, 3, 3, 2, 2, 0, 9, 3, 5, 7, 4, 6, 3, 3, 8, 8, 4, 9, 9, 3, 7, 9, 4, 4, 3, 2, 2, 5, 7, 2, 0, 5, 0, 8, 2, 5, 0, 9, 7, 5, 8, 1, 4, 9, 7, 3, 5, 3, 9, 1, 8, 4, 8, 7, 1

Offset: 1

Klaus Brockhaus, Jun 08 2009

Equals lim_{n -> oo} b(n)/b(n-1) for n mod 3 = {1, 2}, b = A129642.

Equals lim_{n -> oo} b(n)/b(n-1) for n mod 3 = {0, 2}, b = A160580.

			(601+276*sqrt(2))/457 = 2.16919681228659569686...

G. C. Greubel, Table of n, a(n) for n = 1..10000

Cf. A002193, A129642, A160580, A160582.

(601 +276*Sqrt(2))/457; // G. C. Greubel, Apr 08 2018

RealDigits[(601 +276*Sqrt[2])/457, 10, 100][[1]] (* G. C. Greubel, Apr 08 2018 *)

(601 +276*sqrt(2))/457 \\ G. C. Greubel, Apr 08 2018

Equals (23+6*sqrt(2))/(23-6*sqrt(2)).

A160582 Decimal expansion of (213651 +31850*sqrt(2))/457^2.

Original entry on oeis.org

1, 2, 3, 8, 6, 6, 3, 8, 2, 8, 7, 0, 6, 7, 8, 3, 7, 3, 9, 9, 4, 7, 6, 8, 3, 6, 6, 5, 5, 4, 8, 2, 1, 3, 7, 0, 3, 6, 9, 2, 3, 5, 2, 1, 2, 6, 3, 2, 4, 5, 4, 9, 2, 9, 7, 0, 0, 7, 2, 9, 8, 3, 3, 3, 5, 0, 8, 9, 2, 1, 8, 9, 7, 8, 2, 5, 7, 8, 8, 5, 3, 2, 3, 2, 4, 2, 3, 8, 2, 9, 1, 6, 2, 8, 5, 7, 0, 8, 0, 5, 1, 8, 2, 7, 4

Offset: 1

Klaus Brockhaus, Jun 08 2009

Equals lim_{n -> infinity} b(n)/b(n-1) for n mod 3 = 0, b = A129642.
Equals lim_{n -> infinity} b(n)/b(n-1) for n mod 3 = 1, b = A160580.
A quadratic number with minimal polynomial 208849*x^2 - 427302*x + 208849. - Charles R Greathouse IV, Dec 06 2016

			(213651 +31850*sqrt(2))/457^2 = 1.23866382870678373994...

G. C. Greubel, Table of n, a(n) for n = 1..10000
Index entries for algebraic numbers, degree 2

Cf. A129642, A160580, A002193 (decimal expansion of sqrt(2)), A160581 (decimal expansion of (601+276*sqrt(2))/457).

(213651 +31850*Sqrt(2))/457^2; // G. C. Greubel, Apr 07 2018

with(MmaTranslator[Mma]): Digits:=150:
RealDigits(evalf((213651+31850*sqrt(2))/457^2))[1]; # Muniru A Asiru, Apr 08 2018

RealDigits[(213651+31850Sqrt[2])/457^2,10,120][[1]] (* Harvey P. Dale, Jan 06 2013 *)

(213651+31850*sqrt(2))/457^2 \\ Charles R Greathouse IV, Dec 06 2016

Equals (650 +49*sqrt(2))/(650 -49*sqrt(2)).

Equals (3 +2*sqrt(2))*(23 -6*sqrt(2))^2/(23 +6*sqrt(2))^2.

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

Original entry on oeis.org

325, 457, 877, 1073, 2285, 4937, 6113, 13253, 28745, 35605, 77233, 167533, 207517, 450145, 976453, 1209497, 2623637, 5691185, 7049465, 15291677, 33170657, 41087293, 89126425, 193332757, 239474293, 519466873, 1126825885, 1395758465

Offset: 1

Klaus Brockhaus, Jun 08 2009

nonn

(-204, a(1)) and (A129642(n), a(n+1)) are solutions (x, y) to the Diophantine equation x^2+(x+457)^2 = y^2.
lim_{n -> infinity} a(n)/a(n-3) = 3+2*sqrt(2).
lim_{n -> infinity} a(n)/a(n-1) = (601+276*sqrt(2))/457 for n mod 3 = {0, 2}.
lim_{n -> infinity} a(n)/a(n-1) = (213651+31850*sqrt(2))/457^2 for n mod 3 = 1.

			(-204, a(1)) = (-204, 325) is a solution: (-204)^2+(-204+457)^2 = 41616+64009 = 105625 = 325^2.
(A129642(1), a(2)) = (0, 457) is a solution: 0^2+(0+457)^2 = 208849 = 457^2.
(A129642(3), a(4)) = (495, 1073) is a solution: 495^2+(495+457)^2 = 245025+906304 = 1151329 = 1073^2.

Cf. A129642, A001653, A156035 (decimal expansion of 3+2*sqrt(2)), A160581 (decimal expansion of (601+276*sqrt(2))/457), A160582 (decimal expansion of (213651+31850*sqrt(2))/457^2).

{forstep(n=-204, 10000000, [3, 1], if(issquare(2*n^2+914*n+208849, &k), print1(k, ",")))}

a(n) = 6*a(n-3)-a(n-6) for n > 6; a(1)=325, a(2)=457, a(3)=877, a(4)=1073, a(5)=2285, a(6)=4937.
G.f.: (1-x)*(325+782*x+1659*x^2+782*x^3+325*x^4) / (1-6*x^3+x^6).
a(3*k-1) = 457*A001653(k) for k >= 1.

A201916 Nonnegative values x of solutions (x, y) to the Diophantine equation x^2+(x+2737)^2 = y^2.

Original entry on oeis.org

0, 75, 203, 323, 552, 708, 1020, 1127, 1311, 1428, 1608, 1820, 1955, 2336, 2675, 3128, 3311, 3627, 3927, 4140, 4508, 4743, 5535, 6003, 6800, 7280, 7848, 8211, 8588, 9240, 9860, 11063, 11895, 13583, 14168, 15180, 15827, 16827, 18011, 18768, 20915, 22836

Offset: 1

T. D. Noe, Feb 09 2012

Note that 2737 = 7 * 17 * 23, the product of the first three distinct primes in A058529 (and A001132) and hence the smallest such number. This sequence satisfies a linear difference equation of order 55 whose 55 initial terms can be found by running the Mathematica program.
There are many sequences like this one. What determines the order of the linear difference equation? All primes p have order 7. For those p, it appears that p^2 has order 11, p^3 order 15, and p^i order 3+4*i. It appears that for semiprimes p*q (with p > q), the order is 19. What is the next term of the sequence beginning 3, 7, 19, 55, 163? This could be sequence A052919, which is 1 + 2*3^f, where f is the number of primes.
The crossref list is thought to be complete up to Feb 14 2012.

Colin Barker, Table of n, a(n) for n = 1..1000

Cf. A001652 (1), A076296 (7), A118120 (17), A118337 (23), A118674 (31).
Cf. A129288 (41), A118675 (47), A118554 (49), A118673 (71), A129289 (73).
Cf. A118676 (79), A129298 (89), A129836 (97), A157119 (103), A161478 (113).
Cf. A129837 (119), A129992 (127), A129544 (137), A161482 (151).
Cf. A206426 (161), A130608 (167), A161486 (191), A185394 (193).
Cf. A129993 (199), A198294 (217), A130609 (223), A129625 (233).
Cf. A204765 (239), A129991 (241), A207058 (263), A129626 (281).
Cf. A205644 (287), A207059 (289), A129640 (313), A205672 (329).
Cf. A129999 (337), A118611 (343), A130610 (359), A207060 (401).
Cf. A129641 (409), A207061 (433), A130645 (439), A130004 (449).
Cf. A129642 (457), A129725 (521), A101152 (569), A130005 (577).
Cf. A207075 (479), A207076 (487), A207077 (497), A207078 (511).
Cf. A111258 (601), A115135 (617), A130013 (647), A130646 (727).
Cf. A122694 (761), A123654 (809), A129010 (833), A130647 (839).
Cf. A129857 (857), A130014 (881), A129974 (937), A129975 (953).
Cf. A130017 (967), A118630 (2401), A118576 (16807).

d = 2737; terms = 100; t = Select[Range[0, 55000], IntegerQ[Sqrt[#^2 + (#+d)^2]] &]; Do[AppendTo[t, t[[-1]] + 6*t[[-27]] - 6*t[[-28]] - t[[-54]] + t[[-55]]], {terms-55}]; t

a(n) = a(n-1) + 6*a(n-27) - 6*a(n-28) - a(n-54) + a(n-55), where the 55 initial terms can be computed using the Mathematica program.

G.f.: x^2*(73*x^53 +116*x^52 +100*x^51 +171*x^50 +104*x^49 +184*x^48 +57*x^47 +92*x^46 +55*x^45 +80*x^44 +88*x^43 +53*x^42 +139*x^41 +113*x^40 +139*x^39 +53*x^38 +88*x^37 +80*x^36 +55*x^35 +92*x^34 +57*x^33 +184*x^32 +104*x^31 +171*x^30 +100*x^29 +116*x^28 +73*x^27 -363*x^26 -568*x^25 -480*x^24 -797*x^23 -468*x^22 -792*x^21 -235*x^20 -368*x^19 -213*x^18 -300*x^17 -316*x^16 -183*x^15 -453*x^14 -339*x^13 -381*x^12 -135*x^11 -212*x^10 -180*x^9 -117*x^8 -184*x^7 -107*x^6 -312*x^5 -156*x^4 -229*x^3 -120*x^2 -128*x -75) / ((x -1)*(x^54 -6*x^27 +1)). - Colin Barker, May 18 2015

cp's OEIS Frontend

A129725 Nonnegative values x of solutions (x, y) to the Diophantine equation x^2+(x+521)^2 = y^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A160581 Decimal expansion of (601+276*sqrt(2))/457.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A160582 Decimal expansion of (213651 +31850*sqrt(2))/457^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

Formula

A201916 Nonnegative values x of solutions (x, y) to the Diophantine equation x^2+(x+2737)^2 = y^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula