A118337 -id:A118337 - cp's OEIS frontend

A156572 Squares of the form k^2+(k+23)^2 with integer k.

289, 529, 1369, 4225, 13225, 42025, 139129, 444889, 1423249, 4721929, 15108769, 48344209, 160402225, 513249025, 1642275625, 5448949489, 17435353849, 55789022809, 185103876169, 592288777609, 1895184495649, 6288082836025

Offset: 1

Klaus Brockhaus, Feb 11 2009

Square roots of k^2+(k+17)^2 are in A156567, values k are in A118337.

			4225 = 65^2 is of the form k^2+(k+23)^2 with k = 33: 33^2+56^2 = 4225. Hence 4225 is in the sequence.

G. C. Greubel, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,34,-34,0,-1,1).

Cf. A156567, A156575 (first trisection), A156573 (second trisection), A156574 (third trisection).

Cf. A118337, A156035 (decimal expansion of 3+2*sqrt(2)), A156164 (decimal expansion of 17+12*sqrt(2)), A156571 (decimal expansion of (27+10*sqrt(2))/23), A157472 (decimal expansion of (627+238*sqrt(2))/23^2).

LinearRecurrence[{1,0,34,-34,0,-1,1}, {289,529,1369,4225,13225,42025,139129}, 30] (* Harvey P. Dale, Mar 21 2020 *)

{forstep(n=-8, 1800000, [1, 3], if(issquare(a=2*n*(n+23)+529), print1(a, ",")))}

def f(n,p,q): return p*chebyshev_U(n,17) - q*chebyshev_U(n-1,17)
def a(n):
    if (n%3==0): return -289*bool(n==0) + (1/4)*(529 + 3*f(n/3, 209, 5457))
    elif (n%3==1): return (1/4)*(529 + 3*f((n-1)/3, 209, 1649))
    else: return (1/4)*(529 + 3*f((n-2)/3, 529, 529))
[a(n) for n in (1..30)] # G. C. Greubel, Jan 04 2022

a(n) = 34*a(n-3) - a(n-6) - 4232 for n > 6; a(1)=289, a(2)=529, a(3)=1369, a(4)=4225, a(5)=13225, a(6)=42025.
a(n) = A156567(n)^2.
G.f.: x*(289 +240*x +840*x^2 -6970*x^3 +840*x^4 +240*x^5 +289*x^6)/((1-x)*(1 -34*x^3 +x^6)).
Limit_{n -> infinity} a(n)/a(n-3) = 17 + 12*sqrt(2).
Limit_{n -> infinity} a(n)/a(n-1) = ((627 + 238*sqrt(2))/23^2)^2 for n mod 3 = 1.
Limit_{n -> infinity} a(n)/a(n-1) = ((27 + 10*sqrt(2))/23)^2 for n mod 3 = {0, 2}.
a(n) = -289*[n=0] + (529/4) + (3/4)*( f(n/3, 209, 5457)*(n mod 3 = 1) + f((n-1)/3, 209, 1649)*(n mod 3 = 1) + f((n-2)/2, 529, 529)*(n mod 3 = 2) ), where f(n, p, q) = p*ChebyshevU(n, 17) - q*ChebyshevU(n-1, 17). - G. C. Greubel, Jan 04 2022

Revised by Klaus Brockhaus, Feb 16 2009

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

A157472 Decimal expansion of (627 + 238*sqrt(2))/23^2.

Original entry on oeis.org

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

Offset: 1

Klaus Brockhaus, Mar 12 2009

			(627 + 238*sqrt(2))/23^2 = 1.82151763297693123178...

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

Cf. A118337, A156567, A002193 (decimal expansion of sqrt(2)), A156035 (decimal expansion of 3+2*sqrt(2)), A156571 (decimal expansion of (27+10*sqrt(2))/23).

SetDefaultRealField(RealField(100)); (627+238*sqrt(2))/23^2; // G. C. Greubel, Aug 17 2018

RealDigits[(627+238Sqrt[2])/23^2,10,120][[1]] (* Harvey P. Dale, Feb 03 2015 *)

(627 + 238*sqrt(2))/23^2 \\ G. C. Greubel, Aug 17 2018

Equals (34 + 7*sqrt(2))/(34 - 7*sqrt(2)) = (3+2*sqrt(2))*(10- 2*sqrt(2) )^2/(10+2*sqrt(2))^2.
Equals lim_{n -> infinity} b(n)/b(n-1) for n mod 3 = 0, where b is A118337.
Equals lim_{n -> infinity} b(n)/b(n-1) for n mod 3 = 1, where b is A156567.

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

A156572 Squares of the form k^2+(k+23)^2 with integer k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

Extensions

A157472 Decimal expansion of (627 + 238*sqrt(2))/23^2.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Magma

Mathematica

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