A232437 Numbers whose square is expressible in only one way as x^2+xy+y^2, with x and y > 0.

7, 13, 14, 19, 21, 26, 28, 31, 35, 37, 38, 39, 42, 43, 52, 56, 57, 61, 62, 63, 65, 67, 70, 73, 74, 76, 77, 78, 79, 84, 86, 93, 95, 97, 103, 104, 105, 109, 111, 112, 114, 117, 119, 122, 124, 126, 127, 129, 130, 134, 139, 140, 143, 146, 148, 151, 152, 154, 155, 156, 157, 158, 161

Offset

1,1

Comments

Analog of A084645 for 120-degree angle triangles with integer sides.

Numbers with exactly one prime divisor of the form 6k+1 with multiplicity one.

Primitive elements of A050931.

Links

Ray Chandler, Table of n, a(n) for n = 1..10000

Index entries for sequences related to A2 = hexagonal = triangular lattice

Formula

Terms are obtained by the products A230780(k)*A002476(p) for k, p > 0, ordered by increasing values.

Example

a(1) = 7 as 7^2 = 3^2 + 3*5 + 5^2.

Mathematica

r[k_] := Reduce[x>0 && y>0 && k^2 == x^2 + x y + y^2, {x, y}, Integers];
selQ[k_] := Which[rk = r[k]; rk === False, False, rk[[0]] === And && Length[rk] == 2, False, rk[[0]] === Or && Length[rk] == 2, True, True, False];
Select[Range[1000], selQ] (* Jean-François Alcover, Feb 20 2020 *)

Crossrefs

Cf. A002476, A050931, A230780, A232436 (subsequence).

Cf. A084645, A232437, A248599, A254063, A254064.

Sequence in context: A013651 A050931 A072864 * A275201 A274558 A120100

Adjacent sequences: A232434 A232435 A232436 * A232438 A232439 A232440

Keyword

nonn

Author

Jean-Christophe Hervé, Nov 24 2013

Status

approved