login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Number of ways to write n as x + y + z with x, y, z positive integers such that 3*x^2*y^2 + 5*y^2*z^2 + 8*z^2*x^2 is a square.
5

%I #44 May 04 2021 07:37:56

%S 0,0,1,1,1,1,1,1,2,1,1,3,3,2,4,2,5,2,1,2,2,2,2,4,2,4,3,3,3,4,5,2,3,5,

%T 5,4,4,2,4,4,5,3,4,3,6,3,2,5,2,2,7,7,1,3,6,4,4,3,3,6,2,5,5,3,6,5,4,6,

%U 6,6,3,6,6,4,5,6,2,6,3,5,4,5,3,5,12,4,4,5,1,6,6,7,9,3,3,6,5,6,7,4

%N Number of ways to write n as x + y + z with x, y, z positive integers such that 3*x^2*y^2 + 5*y^2*z^2 + 8*z^2*x^2 is a square.

%C Conjecture 1: a(n) > 0 for all n > 2.

%C We have verified a(n) > 0 for all n = 3..10000. The conjecture holds if a(p) > 0 for every odd prime p. For any n > 0 we have a(3*n) > 0, since 3*n = n + n + n and 3 + 5 + 8 = 4^2.

%C It seems that a(n) = 1 only for n = 3..8, 10, 11, 19, 53, 89, 127, 178, 257, 461.

%C See also A343862 for similar conjectures.

%C Conjecture 1 holds for all n < 2^15. Note a(1823) = 1. - _Martin Ehrenstein_, May 03 2021

%H Martin Ehrenstein, <a href="/A340274/b340274.txt">Table of n, a(n) for n = 1..32767</a> (first 1500 terms from Zhi-Wei Sun)

%H Zhi-Wei Sun, <a href="https://doi.org/10.1016/j.jnt.2016.11.008">Refining Lagrange's four-square theorem</a>, J. Number Theory 175 (2017), 167-190. See also <a href="https://arxiv.org/abs/1604.06723">arXiv version</a>, arXiv:1604.06723 [math.NT], 2016-2017.

%e a(4) = 1 with 4 = 2 + 1 + 1 and 3*2^2*1^2 + 5*1^2*1^2 + 8*1^2*2^2 = 7^2.

%e a(19) = 1 with 19 = 9 + 9 + 1 and 3*9^2*9^2 + 5*9^2*1^2 + 8*1^2*9^2 = 144^2.

%e a(53) = 1 with 53 = 23 + 7 + 23 and 3*23^2*7^2 + 5*7^2*23^2 + 8*23^2*23^2 = 1564^2.

%e a(89) = 1 with 89 = 2 + 58 + 29 and 3*2^2*58^2 + 5*58^2*29^2 + 8*29^2*2^2 = 3770^2.

%e a(257) = 1 with 257 = 11 + 164 + 82 and 3*11^2*164^2 + 5*164^2*82^2 + 8*82^2*11^2 = 30340^2.

%e a(461) = 1 with 461 = 186 + 165 + 110 and 3*186^2*165^2 + 5*165^2*110^2 + 8*110^2*186^2 = 88440^2.

%t SQ[n_]:=IntegerQ[Sqrt[n]];

%t tab={};Do[r=0;Do[If[SQ[3x^2*y^2+(n-x-y)^2*(5*y^2+8*x^2)],r=r+1],{x,1,n-2},{y,1,n-1-x}];tab=Append[tab,r],{n,1,100}];Print[tab]

%Y Cf. A000041, A000290, A230121, A230747, A231168, A262357, A271719, A273278, A343862.

%K nonn

%O 1,9

%A _Zhi-Wei Sun_, Apr 24 2021