

A262785


Number of ordered ways to write n as x^2 + y^2 + p*(p+d)/2, where 0 <= x <= y, d is 1 or 1, and p is prime.


7



1, 1, 3, 2, 3, 2, 3, 3, 1, 3, 5, 3, 2, 3, 3, 4, 2, 2, 6, 4, 4, 2, 6, 2, 2, 5, 2, 7, 4, 4, 4, 5, 3, 1, 7, 2, 5, 4, 4, 5, 4, 3, 3, 5, 1, 6, 4, 3, 1, 3, 5, 3, 8, 2, 7, 6, 3, 2, 4, 5, 3, 4, 2, 6, 5, 4, 5, 9, 2, 3, 9, 1, 5, 5, 7, 4, 3, 5, 5, 7, 3, 5, 7, 5, 3, 8, 4, 7, 4, 2, 9, 7, 6, 2, 9, 6, 1, 3, 3, 9
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,3


COMMENTS

Conjecture: (i) a(n) > 0 for all n > 0, and a(n) = 1 only for n = 1, 2, 9, 34, 45, 49, 72, 97, 241, 337, 538.
(ii) Any integer n > 9 can be written as x^2 + y^2 + z*(z+1), where x,y,z are nonnegative integers with z1 or z+1 prime.
In 2015, the author refined a result of Euler by proving that any positive integer can be written as the sum of two squares and a positive triangular number.


LINKS

ZhiWei Sun, Table of n, a(n) for n = 1..10000
ZhiWei Sun, Mixed sums of squares and triangular numbers, Acta Arith. 127(2007), 103113.
ZhiWei Sun, On universal sums ax^2+by^2+f(z), aT_x+bT_y+f(z) and aT_x+by^2+f(z), arXiv:1502.03056 [math.NT], 2015.


EXAMPLE

a(1) = 1 since 1 = 0^2 + 0^2 + 2*(21)/2 with 2 prime.
a(2) = 1 since 2 = 0^2 + 1^2 + 2*(21)/2 with 2 prime.
a(3) = 3 since 3 = 0^2 + 0^2 + 2*(2+1)/2 = 0^2 + 0^2 + 3*(31)/2 = 1^2 + 1^2 + 2*(21)/2 with 2 and 3 both prime.
a(9) = 1 since 9 = 2^2 + 2^2 + 2*(21)/2 with 2 prime.
a(34) = 1 since 34 = 2^2 + 3^2 + 7*(71)/2 with 7 prime.
a(45) = 1 since 45 = 1^2 + 4^2 + 7*(7+1)/2 with 7 prime.
a(49) = 1 since 49 = 3^2 + 5^2 + 5*(5+1)/2 with 5 prime.
a(72) = 1 since 72 = 1^2 + 4^2 + 11*(111)/2 with 11 prime.
a(97) = 1 since 97 = 1^2 + 9^2 + 5(5+1)/2 with 5 prime.
a(241) = 1 since 241 = 1^2 + 15^2 + 5*(5+1)/2 with 5 prime.
a(337) = 1 since 337 = 5^2 + 6^2 + 23*(23+1)/2 with 23 prime.
a(538) = 1 since 538 = 3^2 + 8^2 + 31*(311)/2 with 31 prime.


MATHEMATICA

SQ[n_]:=IntegerQ[Sqrt[n]]
f[d_, n_]:=Prime[n](Prime[n]+(1)^d)/2
Do[r=0; Do[If[SQ[nf[d, k]x^2], r=r+1], {d, 0, 1}, {k, 1, PrimePi[(Sqrt[8n+1](1)^d)/2]}, {x, 0, Sqrt[(nf[d, k])/2]}]; Print[n, " ", r]; Continue, {n, 1, 100}]


CROSSREFS

Cf. A000040, A000290, A000217, A001481, A254885, A262311, A262781.
Sequence in context: A049071 A168330 A176059 * A264843 A316290 A029211
Adjacent sequences: A262782 A262783 A262784 * A262786 A262787 A262788


KEYWORD

nonn


AUTHOR

ZhiWei Sun, Oct 01 2015


STATUS

approved



