%I #21 Apr 09 2017 03:32:56
%S 168,2688,13608,43008,105000,108000,217728,403368,688128,1102248,
%T 1680000,1728000
%N Integer areas A of the integer-sided triangles such that the length of the inradius and the circumradius are both a perfect square.
%C Subset of A208984.
%C The areas of the primitive triangles of sides (a, b, c) and inradius, circumradius equals respectively to r and R are 672, 108000, ... The sides of the nonprimitive triangles are of the form (a*k^2, b*k^2, c*k^2) with r' = r*k^2 and R' = R*k^2 where r', R' are respectively the inradius and the circumradius of the nonprimitive triangles. The areas A' of the nonprimitive triangles are A' = A*k^4.
%C The area A of a triangle whose sides have lengths a, b, and c is given by Heron's formula: A = sqrt(s*(s-a)*(s-b)*(s-c)), where s = (a+b+c)/2. The inradius r is given by r = A/s and the circumradius is given by R = abc/4A.
%C The following table gives the first values (A, a, b, c, r, R).
%C +---------+------+------+------+-----+------+
%C | A | a | b | c | r | R |
%C +---------+------+------+------+-----+------+
%C | 168 | 14 | 30 | 40 | 4 | 25 |
%C | 2688 | 56 | 120 | 160 | 16 | 100 |
%C | 13608 | 126 | 270 | 360 | 36 | 225 |
%C | 43008 | 224 | 480 | 640 | 64 | 400 |
%C | 105000 | 350 | 750 | 1000 | 100 | 625 |
%C | 108000 | 480 | 510 | 510 | 144 | 289 |
%C | 217728 | 504 | 1080 | 1440 | 144 | 900 |
%C | 403368 | 686 | 1470 | 1960 | 196 | 1225 |
%C | 688128 | 896 | 1920 | 2560 | 256 | 1600 |
%C | 1102248 | 1134 | 2430 | 3240 | 324 | 2025 |
%C | 1680000 | 1400 | 3000 | 4000 | 400 | 2500 |
%C | 1728000 | 1920 | 2040 | 2040 | 576 | 1156 |
%C +---------+------+------+------+-----+------+
%H Mohammad K. Azarian, <a href="http://www.jstor.org/stable/25678790">Solution of problem 125: Circumradius and Inradius</a>, Math Horizons, Vol. 16, No. 2 (Nov. 2008), p. 32.
%H Eric W. Weisstein, <a href="http://mathworld.wolfram.com/Circumradius.html">MathWorld: Circumradius</a>
%H Eric W. Weisstein, <a href="http://mathworld.wolfram.com/Inradius.html">MathWorld: Inradius</a>
%e a(1) = 168 because, for (a,b,c) = (14, 30, 40) => s= (14 + 30 + 40)/2 = 42, and
%e A = sqrt(42*(42-14)*(42-30)*(42-40)) = sqrt(28224) = 168;
%e R = abc/4A = 14*30*40/(4*168) = 25;
%e r = A/s = 168/42 = 4.
%t nn=2000;Do[s=(a+b+c)/2;If[IntegerQ[s],area2=s (s-a)(s-b)(s-c);If[0<area2&&IntegerQ[Sqrt[area2]]&&IntegerQ[Sqrt[Sqrt[area2]/s]]&&IntegerQ[Sqrt[a*b*c/(4*Sqrt[area2])]],Print[Sqrt[area2]," ",a " ",b," ",c," ", Sqrt[area2]/s," ", a*b*c/(4*Sqrt[area2])]]],{a,nn},{b,a},{c,b}]
%Y Cf. A188158, A208984.
%K nonn,hard
%O 1,1
%A _Michel Lagneau_, Nov 22 2013
|