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.

A339856 Primitive triples for integer-sided triangles whose sides a < b < c form a geometric progression.

This page as a plain text file.
%I A339856 #48 Jul 27 2023 17:30:13
%S A339856 4,6,9,9,12,16,16,20,25,25,30,36,25,35,49,25,40,64,36,42,49,49,56,64,
%T A339856 49,63,81,49,70,100,49,77,121,64,72,81,64,88,121,81,90,100,81,99,121,
%U A339856 81,117,169,81,126,196,100,110,121,100,130,169,121,132,144,121,143,169
%N A339856 Primitive triples for integer-sided triangles whose sides a < b < c form a geometric progression.
%C A339856 These triangles are called "geometric triangles" in Project Euler problem 370 (see link).
%C A339856 The triples are displayed in increasing lexicographic order (a, b, c).
%C A339856 Equivalently: triples of integer-sided triangles such that b^2 = a*c with a < c and gcd(a, c) = 1.
%C A339856 When a < b < c are in geometric progression with b = a*q, c = b*q, q is the constant, then 1 < q < (1+sqrt(5))/2 = phi = A001622 = 1.6180... (this bound is used in Maple code).
%C A339856 For each triple (a, b, c), there exists (r, s), 0 < r < s such that a = r^2, b = r*s, c = s^2, q = s/r.
%C A339856 Angle C < 90 degrees if 1 < q < sqrt(phi) and angle C > 90 degrees if sqrt(phi) < q < phi with sqrt(phi) = A139339 = 1.2720...
%C A339856 For k >= 2, each triple (a, b, c) of the form (k^2, k*(k+1), (k+1)^2) is (A008133(3k+1), A008133(3k+2), A008133(3k+3)).
%C A339856 Three geometrical properties about these triangles:
%C A339856   1) The sinus satisfy sin^2(B) = sin(A) * sin(C) with sin(A) < sin(B) < sin(C) that form a geometric progression.
%C A339856   2) The heights satisfy h_b^2 = h_a * h_c with h_c < h_b < h_a that form a geometric progression.
%C A339856   3) b^2 = 2 * R * h_b, with R = circumradius of the triangle ABC.
%H A339856 David A. Corneth, <a href="/A339856/b339856.txt">Table of n, a(n) for n = 1..10002</a>
%H A339856 Project Euler, <a href="https://projecteuler.net/problem=370">Problem 370: Geometric Triangles</a>.
%e A339856 The smallest such triangle is (4, 6, 9) with 4*9 = 6^2.
%e A339856 There exist four triangles with small side = 49 corresponding to triples (49, 56, 64), (49, 63, 81), (49, 70, 100) and (49, 77, 121).
%e A339856 The table begins:
%e A339856    4,  6,  9;
%e A339856    9, 12, 16;
%e A339856   16, 20, 25;
%e A339856   25, 30, 36;
%e A339856   25, 35, 49;
%e A339856   25, 40, 64;
%e A339856   36, 42, 49;
%e A339856   ...
%p A339856 for a from 1 to 300 do
%p A339856 for b from a+1 to floor((1+sqrt(5))/2 * a) do
%p A339856 for c from b+1 to floor((1+sqrt(5))/2 * b) do
%p A339856 k:=a*c;
%p A339856 if k=b^2 and igcd(a,b,c)=1 then print(a,b,c); end if;
%p A339856 end do;
%p A339856 end do;
%p A339856 end do;
%o A339856 (PARI) lista(nn) = {my(phi = (1+sqrt(5))/2); for (a=1, nn, for (b=a+1, floor(a*phi), for (c=b+1, floor(b*phi), if ((a*c == b^2) && (gcd([a,b,c])==1), print([a,b,c])););););} \\ _Michel Marcus_, Dec 25 2020
%o A339856 (PARI) upto(n) = my(res=List(), phi = (sqrt(5)+1) / 2); for(i = 2, sqrtint(n), for(j = i+1, (i*phi)\1, if(gcd(i, j)==1, listput(res, [i^2, i*j, j^2])))); concat(Vec(res)) \\ _David A. Corneth_, Dec 25 2020
%Y A339856 Cf. A339857 (smallest side), A339858 (middle side), A339859 (largest side), A339860 (perimeter).
%Y A339856 Cf. A336755 (similar for sides in arithmetic progression).
%Y A339856 Cf. A335893 (similar for angles in arithmetic progression).
%Y A339856 Cf. A001622 (phi), A139339 (sqrt(phi)), A008133.
%K A339856 nonn,easy,tabf
%O A339856 1,1
%A A339856 _Bernard Schott_, Dec 19 2020
%E A339856 Data corrected by _David A. Corneth_, Dec 25 2020