OFFSET
1,2
LINKS
Jean-François Alcover, Table of n, a(n) for n = 1..1000
MATHEMATICA
terms = 1000;
nmax = 12 terms;
okQ[n_] := AllTrue[FactorInteger[n][[All, 1]], MatchQ[Mod[#, 12], 1|11]&];
A072330 = Select[Range[nmax], okQ];
a[n_] := Module[{a, b, c, d, p}, d = If[n <= Length[A072330], A072330[[n]], Print["nmax = ", nmax, " insufficient"]; Exit[]]; If[n == 1, 1, For[b = 2 d, True, b++, a = b - d; c = b + d; p = (a + b + c)/2; If[IntegerQ[p] && IntegerQ[Sqrt[p (p - a) (p - b) (p - c)]] && GCD[a, b, c] == 1, Return[ Sqrt[b^2 - 4 d^2]/(2 Sqrt[3])]]]]];
a /@ Range[terms] (* Jean-François Alcover, Mar 07 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Nov 04 2003
EXTENSIONS
Extended by Ray Chandler, Jul 03 2004
STATUS
approved