cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-9 of 9 results.

A063669 Hypotenuses of reciprocal Pythagorean triangles: number of solutions to 1/(12n)^2 = 1/b^2 + 1/c^2 [with b >= c > 0]; also number of values of A020885 (with repetitions) which divide n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 1, 1, 3, 1, 1, 1, 1, 2, 1, 1, 2, 1, 4, 1, 1, 1, 1, 4, 1, 1, 1, 1, 3, 1, 2, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 2, 1, 1, 1, 4, 1, 1, 1, 1, 2, 1, 1, 1, 1, 6, 1, 1, 1, 1, 2, 1, 1, 1, 1, 3, 1, 1, 1, 3, 2, 1, 1, 1, 1, 4, 2, 1, 1, 1, 2, 1, 1, 2, 1, 3, 1, 1, 1, 1, 5
Offset: 1

Views

Author

Henry Bottomley, Jul 28 2001

Keywords

Comments

Primitive reciprocal Pythagorean triangles 1/a^2 = 1/b^2 + 1/c^2 have a=fg, b=ef, c=eg where e^2 = f^2 + g^2; i.e., e,f,g represent the sides of primitive Pythagorean triangles. But the product of the two legs of primitive Pythagorean triangles are multiples of 12 and so the reciprocal of hypotenuses of reciprocal Pythagorean triangles are always multiples of 12 (A008594).

Examples

			a(1)=1 since 1/(12*1)^2 = 1/12^2 = 1/15^2 + 1/20^2;
a(70)=6 since 1/(12*70)^2 = 1/840^2 = 1/875^2 + 1/3000^2 = 1/888^2 + 1/2590^2 = 1/910^2 + 1/2184^2 = 1/952^2 + 1/1785^2 = 1/1050^2 + 1/1400^2 = 1/1160^2 + 1/1218^2.
Looking at A020885, 1 is divisible by 1, while 70 is divisible by 1, 5, 10, 14, 35 and again 35.
		

Crossrefs

A020883 Ordered long legs of primitive Pythagorean triangles.

Original entry on oeis.org

4, 12, 15, 21, 24, 35, 40, 45, 55, 56, 60, 63, 72, 77, 80, 84, 91, 99, 105, 112, 117, 120, 132, 140, 143, 144, 153, 156, 165, 168, 171, 176, 180, 187, 195, 208, 209, 220, 221, 224, 231, 240, 247, 252, 253, 255, 260, 264, 272, 273, 275, 285, 288, 299, 304, 308, 312, 323
Offset: 1

Views

Author

Keywords

Comments

Consider primitive Pythagorean triangles (A^2 + B^2 = C^2, (A, B) = 1, A < B); sequence gives values of B, sorted.
Any term in this sequence is given by f(m,n) = 2*m*n or g(m,n) = m^2 - n^2 where m and n are any two positive integers, m > 1, n < m, the greatest common divisor of m and n is 1, m and n are not both odd; e.g., f(m,n) = f(2,1) = 2*2*1 = 4. - Agola Kisira Odero, Apr 29 2016
All terms are composite. - Thomas Ordowski, Mar 12 2017
a(1) is the only power of 2. - Torlach Rush, Nov 08 2019
The first term appearing twice is 420 = a(75) = a(76) = A024410(1). - Giovanni Resta, Nov 11 2019
From Bernard Schott, May 05 2021: (Start)
Also, ordered sides a of primitive triples (a, b, c) for integer-sided triangles where side a is the harmonic mean of the 2 other sides b and c, i.e., 2/a = 1/b + 1/c with b < a < c.
Example: a(2) = 12, because the second triple is (12, 10, 15) with side a = 12, satisfying 2/12 = 1/10 + 1/15 and 15-12 < 10 < 15+12.
The first term appearing twice 420 corresponds to triples (420, 310, 651) and (420, 406, 435), the second one is 572 = a(101) = a(102) = A024410(2) and corresponds to triples (572, 407, 962) and (572, 455, 770). The terms that appear more than once in this sequence are in A024410.
For the corresponding primitive triples and miscellaneous properties and references, see A343891. (End)

References

  • V. Lespinard & R. Pernet, Trigonométrie, Classe de Mathématiques élémentaires, programme 1962, problème B-337 p. 179, André Desvigne.

Crossrefs

Triangles with 2/a = 1/b + 1/c: A343891 (triples), A020883 (side a), A343892 (side b), A343893 (side c), A343894 (perimeter).

Programs

  • Maple
    for a from 4 to 325 do
    for b from floor(a/2)+1 to a-1 do
    c := a*b/(2*b-a);
    if c=floor(c) and igcd(a,b,c)=1 and c-bBernard Schott, May 05 2021

Extensions

Extended and corrected by David W. Wilson

A020884 Ordered short legs of primitive Pythagorean triangles.

Original entry on oeis.org

3, 5, 7, 8, 9, 11, 12, 13, 15, 16, 17, 19, 20, 20, 21, 23, 24, 25, 27, 28, 28, 29, 31, 32, 33, 33, 35, 36, 36, 37, 39, 39, 40, 41, 43, 44, 44, 45, 47, 48, 48, 49, 51, 51, 52, 52, 53, 55, 56, 57, 57, 59, 60, 60, 60, 61, 63, 64, 65, 65, 67, 68, 68, 69, 69, 71, 72, 73, 75, 75, 76, 76, 77
Offset: 1

Views

Author

Keywords

Comments

Consider primitive Pythagorean triangles (A^2 + B^2 = C^2, (A, B) = 1, A <= B); sequence gives values of A, sorted.
Union of A081874 and A081925. - Lekraj Beedassy, Jul 28 2006
Each term in this sequence is given by f(m,n) = m^2 - n^2 or g(m,n) = 2mn where m and n are relatively prime positive integers with m > n, m and n not both odd. For example, a(1) = f(2,1) = 2^2 - 1^2 = 3 and a(4) = g(4,1) = 2*4*1 = 8. - Agola Kisira Odero, Apr 29 2016
All powers of 2 greater than 4 (2^2) are terms, and are generated by the function g(m,n) = 2mn. - Torlach Rush, Nov 08 2019

Crossrefs

Cf. A009004, A020882, A020883, A020885, A020886. Different from A024352.
Cf. A024359 (gives the number of times n occurs).
Cf. A037213.

Programs

  • Haskell
    a020884 n = a020884_list !! (n-1)
    a020884_list = f 1 1 where
       f u v | v > uu `div` 2        = f (u + 1) (u + 2)
             | gcd u v > 1 || w == 0 = f u (v + 2)
             | otherwise             = u : f u (v + 2)
             where uu = u ^ 2; w = a037213 (uu + v ^ 2)
    -- Reinhard Zumkeller, Nov 09 2012
  • Mathematica
    shortLegs = {}; amx = 99; Do[For[b = a + 1, b < (a^2/2), c = (a^2 + b^2)^(1/2); If[c == IntegerPart[c] && GCD[a, b, c] == 1, AppendTo[shortLegs, a]]; b = b + 2], {a, 3, amx}]; shortLegs (* Vladimir Joseph Stephan Orlovsky, Aug 07 2008 *)
    Take[Union[Sort/@({Times@@#,(Last[#]^2-First[#]^2)/2}&/@(Select[Subsets[Range[1,101,2],{2}],GCD@@#==1&]))][[;;,1]],80] (* Harvey P. Dale, Feb 06 2025 *)

Extensions

Extended and corrected by David W. Wilson

A020886 Ordered semiperimeters of primitive Pythagorean triangles.

Original entry on oeis.org

6, 15, 20, 28, 35, 42, 45, 63, 66, 72, 77, 88, 91, 99, 104, 110, 117, 120, 130, 143, 153, 156, 165, 170, 187, 190, 195, 204, 209, 210, 221, 228, 231, 238, 247, 255, 266, 272, 273, 276, 285, 299, 304, 322, 323, 325, 336, 342, 345, 350, 357, 368, 378, 391, 399
Offset: 1

Views

Author

Keywords

Comments

k is in this sequence iff A078926(k) > 0.
Also, ordered sides c of primitive triples (a, b, c) for integer-sided triangles where side a is the harmonic mean of the 2 other sides b and c, i.e., 2/a = 1/b + 1/c with b < a < c (A343893). - Bernard Schott, May 06 2021
a(n) are the ordered radii of inscribed circles in squares, from which the tangents to the circles are cut off by primitive Pythagorean triangles. - Alexander M. Domashenko, Oct 17 2024

Crossrefs

Subsequence of A005279.
Triangles with 2/a = 1/b + 1/c: A343891 (triples), A020883 (side a), A343892 (side b), A343893 (side c), A343894 (perimeter).

Programs

  • Maple
    isA020886 := proc(an) local r::integer,s::integer ; for r from floor((an/2)^(1/2)) to floor(an^(1/2)) do for s from r-1 to 1 by -2 do if r*(r+s) = an and gcd(r,s) < 2 then RETURN(true) ; fi ; if r*(r+s) < an then break ; fi ; od ; od : RETURN(false) ; end : for n from 2 to 400 do if isA020886(n) then printf("%d,",n) ; fi ; od ; # R. J. Mathar, Jun 08 2006
  • Mathematica
    A078926[n_] := Sum[Boole[n < d^2 < 2n && CoprimeQ[d, n/d]], {d, Divisors[ n/2^IntegerExponent[n, 2]]}];
    Select[Range[1000], A078926[#]>0&] (* Jean-François Alcover, Mar 23 2020 *)
  • PARI
    is(n,f=factor(n))=my(P=apply(i->f[i,1]^f[i,2],[2-n%2..#f~]),nn=2*n); forvec(v=vector(#P,i,[0,1]), my(d=prod(i=1,#v,P[i]^v[i]),d2=d^2); if(d2n, return(1))); 0
    list(lim)=my(v=List()); forfactored(n=6,lim\1, if(is(n[1],n[2]), listput(v,n[1]))); Vec(v) \\ Charles R Greathouse IV, Feb 03 2023

Formula

a(n) = A024364(n)/2.

A024406 Ordered areas of primitive Pythagorean triangles.

Original entry on oeis.org

6, 30, 60, 84, 180, 210, 210, 330, 504, 546, 630, 840, 924, 990, 1224, 1320, 1386, 1560, 1710, 1716, 2310, 2340, 2574, 2730, 2730, 3036, 3570, 3900, 4080, 4290, 4620, 4914, 5016, 5610, 5814, 6090, 6630, 7140, 7440, 7854, 7956, 7980, 7980, 8970, 8976, 9690
Offset: 1

Views

Author

Keywords

Comments

This sequence also gives Fibonacci's congruous numbers (or congrua) divided by 4 with multiplicities, not regarding leg exchange in the underlying primitive Pythagorean triangle. See A258150 and the example. - Wolfdieter Lang, Jun 14 2015
The squarefree part of an entry which is not squarefree is a primitive congruent number from A006991 belonging to a Pythagorean triangle with rational (not all integer) side lengths (and its companion obtained by exchanging the legs). See the W. Lang link. - Wolfdieter Lang, Oct 25 2016

Examples

			a(6) = a(7) = 210 corresponds to the area (in some squared length unit) of the primitive Pythagorean triangles (21, 20, 29) and (35, 12, 37). Fibonacci's congruum C = 840 = 210*4 belongs to the two triples [x, y, z] = [29, 41, 1] and [37, 47, 23], solving x^2 + C = y^2 and x^2 - C = z^2. - _Wolfdieter Lang_, Jun 14 2015
a(5) = 180 = 6^2*5 lead to the primitive congruent number A006991(1) = 5 from the primitive Pythagorean triangle [9, 40, 41] after division by 6: [3/2, 20/3, 41/6]. See the link for the other nonsquarefree a(n) numbers. - _Wolfdieter Lang_, Oct 25 2016
		

Crossrefs

Formula

a(n) = 6*A020885(n). - Lekraj Beedassy, Apr 30 2004
a(n) = A121728(n)*A121729(n)/2. - M. F. Hasler, Apr 16 2020

A063011 Ordered products of the sides of primitive Pythagorean triangles.

Original entry on oeis.org

60, 780, 2040, 4200, 12180, 14760, 15540, 40260, 65520, 66780, 92820, 120120, 189840, 192720, 199980, 235620, 277680, 354960, 453960, 497640, 595140, 619020, 643500, 1021020, 1063860, 1075620, 1265880, 1484340, 1609080, 1761540
Offset: 1

Views

Author

Henry Bottomley, Jul 26 2001

Keywords

Comments

It is an open question whether any two distinct Pythagorean triples can have the same product of their sides.

Examples

			a(1)=3*4*5=60; a(2)=5*12*13=780 (rather than 6*8*10=480, which would not be primitive).
		

Crossrefs

Programs

  • Mathematica
    k=17000000;lst={};Do[Do[If[IntegerQ[a=Sqrt[c^2-b^2]]&&GCD[a,b,c]==1,If[a>=b,Break[]];x=a*b*c;If[x<=k,AppendTo[lst,x]]],{b,c-1,4,-1}],{c,5,700,1}];Union@lst (* Vladimir Joseph Stephan Orlovsky, Sep 05 2009 *)
    With[{nn=50},Take[(Times@@#)Sqrt[#[[1]]^2+#[[2]]^2]&/@Union[Sort/@ ({Times@@#, (Last[#]^2-First[#]^2)/2}&/@(Select[Subsets[Range[1,nn+1,2],{2}],GCD@@#==1&]))]//Union,nn]] (* Harvey P. Dale, Jun 08 2018 *)

A258150 Triangle of Fibonacci's congruum (congruous) numbers divided by 24 based on primitive Pythagorean triangles. Areas divided by 6 of these triangles.

Original entry on oeis.org

1, 0, 5, 10, 0, 14, 0, 35, 0, 30, 35, 0, 0, 0, 55, 0, 105, 0, 154, 0, 91, 84, 0, 220, 0, 260, 0, 140, 0, 231, 0, 390, 0, 0, 0, 204, 165, 0, 455, 0, 0, 0, 595, 0, 285, 0, 429, 0, 770, 0, 935, 0, 836, 0, 385, 286, 0, 0, 0, 1190, 0, 1330, 0, 0, 0, 506
Offset: 2

Views

Author

Wolfdieter Lang, Jun 11 2015

Keywords

Comments

The problem is: given a square, find a positive integer that, whether added to or subtracted from that square, yields a square. That is, both x^2 + C = y^2 and x^2 - C = z^2. Equivalently: z^2 + C = x^2 and x^2 + C = y^2 (squares in arithmetic progression). This is treated in Fibonacci's 'The book of squares' (Liber Quadratorum (1225) but for rational x,y,z). See the Sigler reference, Proposition 14, pp. 53-74 (note that the formulation of this problem on p. 53 is not correct, 'from a square' should read 'from the same square'). See also van der Waerden, pp. 40-42, and A. Weil, pp. 13-14. The desired number C is called a congruum by Fibonacci (a congruous number in Sigler's translation).
For the history of this problem, see Dickson, pp. 459-472 (he uses the (misleading) term congruent number).
The following solution is based on primitive Pythagorean triangles. (Fibonacci's solution is based on sums of odd squares.) The triangle T(n, m) = 24*C(n, m) will have 0's for those (n, m) not leading to primitive Pythagorean triples.
Addition of the two equations, substitution of y = u + v > 0 and z = |u - v| and division by 2 leads to x^2 = u^2 + v^2. Consider primitive Pythagorean triples (u, v, x) with even v which are pairwise relatively prime. Then also GCD(u,v,x) = 1. A common factor f for u, v and x would lead to a multiplication by f^2 on both sides of the two equations. For primitive Pythagorean triples see A249866. One has u = n^2 - m^2, v = 2*n*m and x = n^2 + m^2 with GCD(n, m) = 1 , n > m >= 1, n + m odd. Then C = C(n, m) = 4*n*m*(n^2 - m^2) = 2*v(n, m)*u(n, m). This is four times the area of the Pythagorean triangle. C is divisible by 4! = 24 (see A020885). Define T(n, m) = C(n, m)/4!, for 2 <= m + 1 <= n. This is the area of the corresponding primitive Pythagorean triangle divided by 6.
The corresponding x = x(n, m), y = y(n, m) and z(n, m) number triangles are given in A222946, A225949 and A258149 respectively.
T(n, m) = n*m*(n^2 - m^2)/6, for m = 1, 2, ..., n-1, has for n >= 2 the minimum value at m = 1, and the next largest value appears for n >= 3 at m = n-1. Note all (n, m) pairs are considered here. The proof of the first part is easy. The proof of T(n, m) - T(n, n-1) > 0, for m = 2, 3, ..., n-2, and n >= 3, is equivalent to n^2*(m-2) + 3*n > m^3 +1 and this is easy to prove with n >= m+2 and m >= 2. Therefore the triangle T(n, m) with 0's attains for even n the smallest nonzero row entry at m = 1, and for odd n the smallest nonzero row entry appears at m = n-1 (last entry).
This allows us to find (after solution of two cubic equations for even and odd n, named ne = ne(N) and no = no(N)) a row number nmin(N) = max(ne(n), no(N)) such that N will not appear in any row n > nmin(N).
The original problem posed to Fibonacci by Giovanni di Palermo (Master John of Palermo) was to find a [rational] square that when increased or decreased by 5 gives a square. Fibonacci gave the solution in his Liber Quadratorum in Proposition 17 (see Sigler, pp. 77-81) as x^2 = (41/12)^2 = 1681/144, y^2 = (49/12)^2 = 2401/144 and z^2 = (31/12)^2 = 961/144. This corresponds to the integer quartet (C; x, y, z) = (720; 41, 49, 31) corresponding to the primitive Pythagorean triple [9, 40, 41]. See the examples for (n, m) = (5, 4).
The numbers without zeros, in nondecreasing order, are given in A020885 = A024406/6.
Comments from Eric Snyder, Feb 07 2023: (Start)
If m+n > 3 and not divisible by 3, then m+n | T(n,m).
Additionally, if 2n-1 > 3 and not divisible by 3, then 2n-1 = 6k+-1, and T(n,n-1) = (2n-1)*P(-+k), where P(-+k) is a generalized pentagonal number (A001318). For example, T(6,5) = 11*P(-2) = 11*5.
T(n,n-1) = A000330(n-1) for n>=2. (End)

Examples

			The triangle T(n, m) begins:
n\m   1   2   3   4    5   6    7   8   9  10    11
2:    1
3:    0   5
4:   10   0  14
5:    0  35   0  30
6:   35   0   0   0   55
7:    0 105   0 154    0  91
8:   84   0 220   0  260   0  140
9:    0 231   0 390    0   0    0 204
10: 165   0 455   0    0   0  595   0 285
11:   0 429   0 770    0 935    0 836   0 385
12: 286   0   0   0 1190   0 1330   0   0   0   506
...
The smallest nonzero number for each row with even n is T(n, 1), and for odd n it is T(n, n-1).
The above mentioned nmin(N) will for N = 300 be 12.
Therefore, no number > 300 will appear for rows with n > 12.
-----------------------------------------------------
The corresponding quartets (C; x, y, z) are:
n=2:  (24; 5, 7, 1),
n=3:  (120; 13, 17, 7),
n=4:  (240; 17, 23, 7), (336; 25, 31, 17),
n=5:  (840; 29, 41, 1), (720; 41, 49, 31),
n=6:  (840; 37, 47, 23), (1320; 61, 71, 49),
n=7:  (2520; 53, 73, 17), (3696; 65, 89, 23),
      (2184; 85, 97, 71),
n=8:  (2016; 65, 79, 47), (5280; 73, 103, 7),
      (6240; 89, 119, 41), (3360; 113, 127, 97),
n=9:  (5544; 85, 113, 41), (9360; 97, 137, 7),
      (4896; 145, 161, 127),
n=10: (3960; 101, 119, 79), (10920; 109, 151, 31),
      (14280; 149, 191, 89), (6840; 181, 199, 161),
n=11: (10296; 125, 161, 73), (18480; 137, 193, 17),
      (22440; 157, 217, 47), (20064; 185, 233, 119),
      (9240; 221, 241, 199),
n=12: (6864; 145, 167, 119), (28560; 169, 239, 1),
      (31920; 193, 263, 73), (12144; 265, 287, 241),
...
-----------------------------------------------------
The corresponding primitive Pythagorean triples
(u, v, x) are:
n=2:  [3, 4, 5],
n=3:  [5, 12, 13],
n=4:  [15, 8, 17], [7, 24, 25],
n=5:  [21, 20, 29],[9, 40, 41],
n=6:  [35, 12, 37], [11, 60, 61],
n=7:  [45, 28, 53], [33, 56, 65],
      [13, 84, 85],
n=8:  [63, 16, 65], [55, 48, 73],
      [39, 80, 89], [15, 112, 113],
n=9:  [77, 36, 85], [65, 72, 97],
      [17, 144, 145],
n=10: [99, 20, 101], [91, 60, 109],
      [51, 140, 149], [19, 180, 181],
n=11: [117, 44, 125], [105, 88, 137],
      [85, 132, 157], [57, 176, 185],
      [21, 220, 221],
n=12: [143, 24, 145], [119, 120, 169],
      [95, 168, 193], [23, 264, 265],
...
		

References

  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 2, 1920, pp. 459-472.
  • L. E. Sigler, Leonardo Pisano, Fibonacci, The book of squares, Academic Press, 1987.
  • B. L. van der Waerden, A History of Algebra, Springer, 1985, pp. 40-42.
  • André Weil, Number Theory, An approach through history, From Hammurapi to Legendre, Birkhäuser, 1984, pp. 13-14.

Crossrefs

Programs

  • Mathematica
    T[n_, m_] /; 2 <= m+1 <= n && OddQ[n+m] && CoprimeQ[n, m] := n*m*(n^2 - m^2)/6; T[, ] = 0; Table[T[n, m], {n, 2, 12}, {m, 1, n-1}] // Flatten (* Jean-François Alcover, Jun 16 2015, after given formula *)

Formula

T(n, m) = n*m*(n^2 - m^2)/6 if 2 <= m+1 <= n, n+m odd, GCD(n, m) = 1 and 0 otherwise.

A258151 Areas of primitive Pythagorean triangles divided by 6, in increasing order without multiple entries.

Original entry on oeis.org

1, 5, 10, 14, 30, 35, 55, 84, 91, 105, 140, 154, 165, 204, 220, 231, 260, 285, 286, 385, 390, 429, 455, 506, 595, 650, 680, 715, 770, 819, 836, 935, 969, 1015, 1105, 1190, 1240, 1309, 1326, 1330, 1495, 1496, 1615, 1729, 1771, 1785, 1820, 1925
Offset: 1

Views

Author

Wolfdieter Lang, Jun 14 2015

Keywords

Comments

See A020885 for this sequence with multiplicities. See A024365 for the areas with multiplicities.
This sequence gives also Fibonacci's congruous numbers divided by 24 without multiple entries. See A258150.

Examples

			See A020885.
		

Crossrefs

A263573 Intersection of A024365 and A129912.

Original entry on oeis.org

6, 30, 60, 180, 210, 2310, 4620, 60060, 510510, 10810800, 116396280, 200560490130, 401120980260
Offset: 1

Views

Author

Bill McEachen, Oct 21 2015

Keywords

Comments

The two sequences involve areas of primitive Pythagorean triples and primorial products. Intersections are only considered once (no repeats). Conjecture: the sequence is infinite.
Conjecture: The next two entries are a(12) = 200560490130, a(13) = 401120980260.
From G. C. Greubel, Dec 29 2015: (Start)
6|a(n) for n>=1,
30|a(n) for n>=2,
a(n)/6 = {1, 5, 10, 30, 35, 385, 770, 10010, ...} is a subset of values found in A008706.
(End)
a(12) and a(13) confirmed. a(14) > 2*10^31, if it exists. - Giovanni Resta, Mar 31 2017

Examples

			A024365 begins {6, 30, 60, 84, 180, 210, 210, 330, 504, 546, 630, 840, 924, 990, 1224, 1320, 1386, 1560, 1710, 1716, 2310, ...}.
A129912 begins {1, 2, 6, 12, 30, 60, 180, 210, 360, 420, 1260, 2310, 2520, ...}.
So, common entries encountered are {6, 30, 60, 180, 210, 2310, ...}.
Specifically, we see that A024365(1) = A129912(3), A024365(2) = A129912(5), A024365(3) = A129912(6), A024365(5) = A129912(7).
These are then the first four entries of the sequence (6, 30, 60, 180).
		

Crossrefs

Programs

  • Mathematica
    s = 6 Take[Sort[(Times @@ #)/12 & /@ ({Times @@ #, (Last[#]^2 - First[#]^2)/2} & /@ Select[Subsets[Range[1, 3600, 2], {2}], GCD @@ # == 1 &])], 1800]; f[m_] := f[m] = Union[Times @@@ Subsets[FoldList[Times, 1, Prime[Range[m]]]]][[1 ;; 100]]; f[10]; f[m = 11]; While[f[m] != f[m - 1], m++]; t = f[m]; Intersection[s, t] (* Michael De Vlieger, Oct 22 2015, after Harvey P. Dale at A020885 and Jean-François Alcover at A129912 *) (* or *)
    ok[n_] := Block[{a, f = Power @@@ FactorInteger[2 n]}, SelectFirst[ Subsets[f, {1, Floor[ Length[f]/2]}], (a = Times @@ #; IntegerQ@ Sqrt[a^2 + (2 n/a)^2]) &, {}] != {}]; pr[n_] := Product[ Prime[n+1-i]^i, {i, n}];  upto[mx_] := Block[{ric, j=1}, ric[n_, ip_, ex_] := If[n < mx, Block[{p = Prime[ip + 1]}, If[ex == 1 && ok[n], Sow@ n]; ric[n p^ex, ip + 1, ex]; If[ex > 1, ric[n p^(ex - 1), ip+1, ex-1]]]]; Sort@ Reap[ While[pr[j] < mx, ric[2^j, 1, j]; j++]][[2, 1]]]; upto[10^12] (* much faster, Giovanni Resta, Mar 31 2017 *)
  • PARI
    \\note: code does not generate the sequence, just checks for a matching PPT entry
    genit(area)={myMax=floor(sqrt(2*area));i5=myMax;endless=0;soln=List();
    while(i5>=2,dun=0;j=2.*myVal/i5; k=floor(j); if(j>k, dun=1 );if(dun<1,
    c=sqrt(i5^2 + k^2);w=floor(c);if(c>w,dun=1); if(dun<1,if(gcd(k,i5)>1,dun=1 ));
    if(dun<1,listput(soln,k); listput(soln,i5);listput(soln,w);listsort(soln);
    print("soln a,b,c = ", soln[1],"  ",soln[2],"  ",soln[3] );dun=2;break ));
    i5--;endless++);if(i5<=2&&dun<1,print("no solution ") );if(i5>2&&dun<2,
    print("max iteration limit was hit ",endless) );print (endless);}
    (C++)
    #include 
    #include 
    using namespace std;
    int main(){ifstream fin1,fin2;
    int myValue,myValue2,ptr,fptr,i5,j5;
    unsigned long list1[9999]={0};
    unsigned long list2[999]={0};
    unsigned long final[31]={0};
    fin1.open("A024365.txt"); fin2.open("A129912.txt");
    ptr=1;
    while(ptr<9999)
    {fin1>> myValue;fin1.get();list1[ptr]=myValue;
        if(ptr<999)
           {fin2>> myValue2;fin2.get();list2[ptr]=myValue2;}
        ptr++;}
    fin1.close();fin2.close();fptr=1;
    for(i5=1;i5<9990;i5++)
    {for(j5=1;j5<999;j5++){
    if(list1[i5]==list2[j5] )
    {
        fptr++;
        if(fptr>30){break;}
        final[fptr]=list1[i5];
        cout << final[fptr] << ",";
        break;
    }}if(fptr>30){break;}}}

Extensions

a(12)-a(13) from Giovanni Resta, Mar 31 2017
Showing 1-9 of 9 results.