A218794 Numbers that can be written as p^2 + 3pq + q^2 with primes p < q.

31, 59, 79, 95, 121, 179, 191, 229, 251, 295, 311, 389, 395, 401, 451, 479, 491, 541, 569, 671, 695, 719, 745, 809, 899, 971, 1019, 1061, 1109, 1111, 1121, 1151, 1249, 1271, 1301, 1409, 1451, 1499, 1595, 1619, 1661, 1711, 1919, 1931, 1949, 1991, 2059, 2105, 2111, 2141, 2195, 2201, 2245

Offset

1,1

Comments

This is a subsequence of A218793, with the restriction that p < q, excluding terms of the form 5p^2 unless they would have another decomposition of the given form.

Sequence A218771 is the subsequence of primes in this sequence.

Links

Table of n, a(n) for n=1..53.

Example

a(1) = 31 = p^2+3pq+q^2 for p=2, q=3.
a(20) = 671 = p^2+3pq+q^2 for (p,q)=(2,23) and (5,19) is the least term to allow more than 1 decomposition. See A218795 for more such terms.

Mathematica

With[{nn=60}, Take[Union[#[[1]]^2+3Times@@#+#[[2]]^2&/@Subsets[Prime[ Range[ Floor[nn/3]]], {2}]], nn]] (* Harvey P. Dale, Apr 08 2013 *)

Prog

(PARI) is_A218794(n, v=0)={ /* set v=1 to count number of decompositions, and v=2 to print them */ my(r, c=0); forprime( q=1, sqrtint((n-1)\5), issquare(4*n+5*q^2, &r) || next; isprime((r-3*q)/2) || next; v || return(1); v>1 && print1([q, (r-3*q)/2]", "); c++); c}

Crossrefs

Sequence in context: A193362 A293345 A096453 * A218771 A109356 A061251

Adjacent sequences: A218791 A218792 A218793 * A218795 A218796 A218797

Keyword

nonn

Author

M. F. Hasler, Nov 05 2012

Status

approved