|
%I #11 Dec 14 2015 12:29:26
%S 3,7,8,5,7,8,5,19,21,7,13,15,7,37,40,8,13,15,9,61,65,11,31,35,11,91,
%T 96,13,43,48,13,127,133,15,169,176,16,19,21,16,49,55,17,73,80,17,217,
%U 225,19,91,99,19,271,280,21,331,341,23,133,143,23,397,408
%N Primitive Eisenstein triples: (a,b,c) in lexicographic order such that a^2 + b^2 - a*b - c^2 = 0, a < b < c, and gcd(a, b) = 1.
%C The sides of a primitive 60-degree integer triangle.
%H Colin Barker, <a href="/A264826/b264826.txt">Table of n, a(n) for n = 1..9999</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PythagoreanTriple.html">Pythagorean Triple</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Eisenstein_triple">Eisenstein triple</a>
%o (PARI)
%o pt60(a) = {
%o my(L=List(), n=-3*a^2, f, g, b, c);
%o fordiv(n, f,
%o g=n\f;
%o if(f>g && (g+f)%2==0 && (f-g)%4==0,
%o b=(f-g)\4; c=((f+g)\2+a)\2;
%o if(c>0 && a<b && gcd(a, c)==1, listput(L, [a,b,c]))
%o )
%o );
%o Vec(L)
%o }
%o concat(concat(vector(30, a, pt60(a))))
%Y Cf. A089025, A121992, A201223, A263728, A264827.
%K nonn,tabf
%O 1,1
%A _Colin Barker_, Nov 26 2015
|
