A001481 -id:A001481 - cp's OEIS frontend

A303656 Number of ways to write n as a^2 + b^2 + 3^c + 5^d, where a,b,c,d are nonnegative integers with a <= b.

0, 1, 1, 2, 1, 3, 2, 3, 2, 4, 3, 4, 2, 4, 4, 3, 2, 4, 4, 3, 2, 4, 3, 4, 1, 4, 5, 6, 4, 6, 5, 5, 6, 6, 5, 8, 4, 6, 6, 5, 4, 7, 5, 7, 5, 6, 4, 5, 3, 4, 7, 6, 7, 8, 5, 4, 7, 5, 5, 9, 3, 6, 5, 6, 4, 6, 5, 7, 7, 4, 5, 5, 5, 4, 6, 5, 6, 10, 5, 4, 5, 7, 4, 9, 2, 9, 8, 5, 6, 6

Offset: 1

Zhi-Wei Sun, Apr 27 2018

nonn

Conjecture: a(n) > 0 for all n > 1. In other words, any integer n > 1 can be written as the sum of two squares, a power of 3 and a power of 5.
It has been verified that a(n) > 0 for all n = 2..2*10^10.
It seems that any integer n > 1 also can be written as the sum of two squares, a power of 2 and a power of 3.
The author would like to offer 3500 US dollars as the prize for the first proof of his conjecture that a(n) > 0 for all n > 1. - Zhi-Wei Sun, Jun 05 2018
Jiao-Min Lin (a student at Nanjing University) has verified a(n) > 0 for all 1 < n <= 2.4*10^11. - Zhi-Wei Sun, Jul 30 2022

			a(2) = 1 with 2 = 0^2 + 0^2 + 3^0 + 5^0.
a(5) = 1 with 5 = 0^2 + 1^2 + 3^1 + 5^0.
a(25) = 1 with 25 = 1^2 + 4^2 + 3^1 + 5^1.

Zhi-Wei Sun, Table of n, a(n) for n = 1..100000
Zhi-Wei Sun, Refining Lagrange's four-square theorem, J. Number Theory 175(2017), 167-190.
Zhi-Wei Sun, New conjectures on representations of integers (I), Nanjing Univ. J. Math. Biquarterly 34(2017), no. 2, 97-120.
Zhi-Wei Sun, Restricted sums of four squares, arXiv:1701.05868 [math.NT], 2017-2018.

Cf. A000244, A000290, A000351, A001481, A273812, A302982, A302984, A303233, A303234, A303338, A303363, A303389, A303393, A303399, A303428, A303401, A303429, A303432, A303434, A303539, A303540, A303541, A303543, A303601, A303637, A303639, A303702, A303821.

SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]];
f[n_]:=f[n]=FactorInteger[n];
g[n_]:=g[n]=Sum[Boole[Mod[Part[Part[f[n],i],1],4]==3&&Mod[Part[Part[f[n],i],2],2]==1],{i,1,Length[f[n]]}]==0;
QQ[n_]:=QQ[n]=(n==0)||(n>0&&g[n]);
tab={};Do[r=0;Do[If[QQ[n-3^k-5^m],Do[If[SQ[n-3^k-5^m-x^2],r=r+1],{x,0,Sqrt[(n-3^k-5^m)/2]}]],{k,0,Log[3,n]},{m,0,If[n==3^k,-1,Log[5,n-3^k]]}];tab=Append[tab,r],{n,1,90}];Print[tab]

A303539 Number of ordered pairs (k, m) with 0 <= k <= m such that n - binomial(2k,k) - binomial(2m,m) can be written as the sum of two squares.

Original entry on oeis.org

0, 1, 2, 3, 2, 2, 3, 4, 3, 2, 3, 6, 4, 2, 2, 4, 4, 2, 2, 5, 5, 5, 4, 4, 4, 4, 4, 5, 4, 5, 4, 4, 3, 4, 3, 4, 4, 5, 6, 5, 5, 5, 4, 7, 3, 3, 4, 6, 4, 2, 3, 5, 6, 3, 4, 5, 6, 5, 2, 5, 4, 5, 3, 2, 4, 5, 4, 3, 3, 3, 6, 7, 5, 5, 6, 10, 6, 3, 4, 8

Offset: 1

Zhi-Wei Sun, Apr 25 2018

nonn

Conjecture: a(n) > 0 for all n > 1.
a(n) > 0 for all n = 2..10^10.
See also A303540 and A303541 for related sequences.

			a(2) = 1 with 2 - binomial(2*0,0) - binomial(2*0,0) = 0^2 + 0^2.
a(3) = 2 with 3 - binomial(2*0,0) - binomial(2*0,0) = 0^2 + 1^2 and 3 - binomial(2*0,0) - binomial(2*1,1) = 0^2 + 0^2.
a(5) = 2 with 5 - binomial(2*0,0) - binomial(2*1,1) = 1^2 + 1^2 and 5 - binomial(2*1,1) - binomial(2*1,1) = 0^2 + 1^2.

Zhi-Wei Sun, Table of n, a(n) for n = 1..10000
Zhi-Wei Sun, Refining Lagrange's four-square theorem, J. Number Theory 175(2017), 167-190.
Zhi-Wei Sun, New conjectures on representations of integers (I), Nanjing Univ. J. Math. Biquarterly 34(2017), no. 2, 97-120.
Zhi-Wei Sun, Restricted sums of four squares, arXiv:1701.05868 [math.NT], 2017-2018.

Cf. A000290, A000984, A001481, A303233, A303234, A303338, A303363, A303389, A303393, A303399, A303428, A303401, A303432, A303434, A303540, A303541, A303543.

c[n_]:=c[n]=Binomial[2n,n];
f[n_]:=f[n]=FactorInteger[n];
g[n_]:=g[n]=Sum[Boole[Mod[Part[Part[f[n],i],1],4]==3&&Mod[Part[Part[f[n],i],2],2]==1],{i,1,Length[f[n]]}]==0;
QQ[n_]:=QQ[n]=(n==0)||(n>0&&g[n]);
tab={};Do[r=0;k=0;Label[bb];If[c[k]>n,Goto[aa]];Do[If[QQ[n-c[k]-c[j]],r=r+1],{j,0,k}];k=k+1;Goto[bb];Label[aa];tab=Append[tab,r],{n,1,80}];Print[tab]

A303541 Numbers of the form k^2 + binomial(2*m,m) with k and m nonnegative integers.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 15, 17, 18, 20, 21, 22, 24, 26, 27, 29, 31, 36, 37, 38, 42, 45, 50, 51, 55, 56, 65, 66, 69, 70, 71, 74, 79, 82, 83, 84, 86, 87, 95, 101, 102, 106, 119, 120, 122, 123, 127, 134

Offset: 1

Zhi-Wei Sun, Apr 25 2018

nonn

The conjecture in A303540 has the following equivalent version: Each integer n > 1 can be written as the sum of two terms of the current sequence.

This has been verified for all n = 2..10^10.

			a(1) = 1 with 0^2 + binomial(2*0,0) = 1.
a(7) = 10 with 2^2 + binomial(2*2,2) = 10.
a(8) = 11 with 3^2 + binomial(2*1,1) = 11.

Zhi-Wei Sun, Table of n, a(n) for n = 1..10000
Zhi-Wei Sun, Refining Lagrange's four-square theorem, J. Number Theory 175(2017), 167-190.
Zhi-Wei Sun, New conjectures on representations of integers (I), Nanjing Univ. J. Math. Biquarterly 34(2017), no. 2, 97-120.
Zhi-Wei Sun, Restricted sums of four squares, arXiv:1701.05868 [math.NT], 2017-2018.

Cf. A000290, A000984, A001481, A303233, A303234, A303338, A303363, A303389, A303393, A303399, A303428, A303401, A303432, A303434, A303539, A303540, A303543.

c[n_]:=c[n]=Binomial[2n,n];
SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]];
tab={};n=0;Do[k=0;Label[bb];If[c[k]>m,Goto[aa]];If[SQ[m-c[k]],n=n+1;tab=Append[tab,m];Goto[aa],k=k+1;Goto[bb]];Label[aa],{m,1,134}];Print[tab]

A303543 Number of ways to write n as a^2 + b^2 + C(k) + C(m) with 0 <= a <= b and 0 < k <= m, where C(k) denotes the Catalan number binomial(2k,k)/(k+1).

Original entry on oeis.org

0, 1, 2, 3, 2, 3, 4, 4, 2, 3, 5, 5, 2, 3, 5, 5, 4, 3, 6, 8, 4, 3, 6, 6, 3, 3, 5, 7, 6, 3, 4, 8, 5, 2, 6, 7, 3, 4, 5, 5, 6, 4, 5, 10, 6, 4, 7, 8, 4, 2, 7, 9, 9, 5, 7, 11, 8, 2, 5, 11, 5, 4, 4, 8, 8, 4, 6, 11, 10, 3, 6, 8, 5, 5, 6, 7, 6, 6, 5, 9

Offset: 1

Zhi-Wei Sun, Apr 25 2018

nonn

Conjecture: a(n) > 0 for all n > 1. In other words, any integer n > 1 can be written as the sum of two squares and two Catalan numbers.

This is similar to the author's conjecture in A303540. It has been verified that a(n) > 0 for all n = 2..10^9.

			a(2) = 1 with 2 = 0^2 + 0^2 + C(1) + C(1).
a(3) = 2 with 3 = 0^2 + 1^2 + C(1) + C(1) = 0^2 + 0^2 + C(1) + C(2).

Zhi-Wei Sun, Table of n, a(n) for n = 1..10000
Zhi-Wei Sun, Refining Lagrange's four-square theorem, J. Number Theory 175(2017), 167-190.
Zhi-Wei Sun, New conjectures on representations of integers (I), Nanjing Univ. J. Math. Biquarterly 34(2017), no. 2, 97-120.
Zhi-Wei Sun, Restricted sums of four squares, arXiv:1701.05868 [math.NT], 2017-2018.

Cf. A000108, A000290, A001481, A303233, A303234, A303338, A303363, A303389, A303393, A303399, A303428, A303401, A303432, A303434, A303539, A303540, A303541, A303601.

SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]];
c[n_]:=c[n]=Binomial[2n,n]/(n+1);
f[n_]:=f[n]=FactorInteger[n];
g[n_]:=g[n]=Sum[Boole[Mod[Part[Part[f[n],i],1],4]==3&&Mod[Part[Part[f[n],i],2],2]==1],{i,1,Length[f[n]]}]==0;
QQ[n_]:=QQ[n]=(n==0)||(n>0&&g[n]);
tab={};Do[r=0;k=1;Label[bb];If[c[k]>n,Goto[aa]];Do[If[QQ[n-c[k]-c[j]],Do[If[SQ[n-c[k]-c[j]-x^2],r=r+1],{x,0,Sqrt[(n-c[k]-c[j])/2]}]],{j,1,k}];k=k+1;Goto[bb];Label[aa];tab=Append[tab,r],{n,1,80}];Print[tab]

A007692 Numbers that are the sum of 2 nonzero squares in 2 or more ways.

Original entry on oeis.org

50, 65, 85, 125, 130, 145, 170, 185, 200, 205, 221, 250, 260, 265, 290, 305, 325, 338, 340, 365, 370, 377, 410, 425, 442, 445, 450, 481, 485, 493, 500, 505, 520, 530, 533, 545, 565, 578, 580, 585, 610, 625, 629, 650, 680, 685, 689, 697, 725, 730, 740, 745, 754, 765

Offset: 1

N. J. A. Sloane

A025426(a(n)) > 1. - Reinhard Zumkeller, Aug 16 2011
For the question that is in the link AskNRICH Archive: It is easy to show that (a^2 + b^2)*(c^2 + d^2) = (a*c + b*d)^2 + (a*d - b*c)^2 = (a*d + b*c)^2 + (a*c - b*d)^2. So terms of this sequence can be generated easily. 5 is the least number of the form a^2 + b^2 where a and b distinct positive integers and this is a list sequence. This is the why we observe that there are many terms which are divisible by 5. - Altug Alkan, May 16 2016
Square roots of square terms: {25, 50, 65, 75, 85, 100, 125, 130, 145, 150, 169, 170, 175, 185, 195, 200, 205, 221, 225, 250, 255, 260, 265, 275, 289, 290, 300, 305, ...}. They are also listed by A009177. - Michael De Vlieger, May 16 2016

			50 is a term since 1^2 + 7^2 and 5^2 + 5^2 equal 50.
25 is not a term since though 3^2 + 4^2 = 25, 25 is square, i.e., 0^2 + 5^2 = 25, leaving it with only one possible sum of 2 nonzero squares.
625 is a term since it is the sum of squares of {0,25}, {7,24}, and {15,20}.

Ming-Sun Li, Kathryn Robertson, Thomas J. Osler, Abdul Hassen, Christopher S. Simons and Marcus Wright, "On numbers equal to the sum of two squares in more than one way", Mathematics and Computer Education, 43 (2009), 102 - 108.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
D. Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, 125.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
AskNRICH Archive, Numbers expressible as the sum of 2 squares in more than one way
D. J. C. Mackay and S. Mahajan, Numbers that are Sums of Squares in Several Ways
G. Xiao, Two squares
Index entries for sequences related to sums of squares

Subsequence of A001481. A subsequence is A025285 (2 ways).

Cf. A004431, A118882, A000404, A018825, A025284 (one way).

import Data.List (findIndices)
a007692 n = a007692_list !! (n-1)
a007692_list = findIndices (> 1) a025426_list
-- Reinhard Zumkeller, Aug 16 2011

Select[Range@ 800, Length@ Select[PowersRepresentations[#, 2, 2], First@ # != 0 &] > 1 &] (* Michael De Vlieger, May 16 2016 *)

isA007692(n)=nb = 0; lim = sqrtint(n); for (x=1, lim, if ((n-x^2 >= x^2) && issquare(n-x^2), nb++); ); nb >= 2; \\ Altug Alkan, May 16 2016

is(n)=my(t); if(n<9, return(0)); for(k=sqrtint(n\2-1)+1,sqrtint(n-1), if(issquare(n-k^2)&&t++>1, return(1))); 0 \\ Charles R Greathouse IV, Jun 08 2016

A303637 Number of ways to write n as x^2 + y^2 + 2^z + 5*2^w, where x,y,z,w are nonnegative integers with x <= y.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 2, 2, 2, 2, 4, 3, 4, 5, 5, 5, 4, 4, 5, 4, 5, 9, 8, 6, 6, 9, 7, 6, 8, 8, 10, 8, 4, 8, 5, 7, 9, 12, 9, 6, 10, 9, 11, 10, 8, 16, 12, 8, 9, 12, 9, 11, 12, 11, 9, 10, 12, 14, 10, 10

Offset: 1

Zhi-Wei Sun, Apr 27 2018

nonn

Conjecture: a(n) > 0 for all n > 5.
This has been verified for n up to 5*10^9. Note that 321256731 cannot be written as x^2 + (2*y)^2 + 2^z + 5*2^w with x,y,z,w nonnegative integers.
In contrast, Crocker proved in 2008 that there are infinitely many positive integers not representable as the sum of two squares and at most two powers of 2.
570143 cannot be written as x^2 + y^2 + 2^z + 3*2^w with x,y,z,w nonnegative integers, and 2284095 cannot be written as x^2 + y^2 + 2^z + 7*2^w with x,y,z,w nonnegative integers.
Jiao-Min Lin (a student at Nanjing University) found a counterexample to the conjecture: a(18836421387) = 0. - Zhi-Wei Sun, Jul 21 2022

			a(6) = 1 with 6 = 0^2 + 0^2 + 2^0 + 5*2^0.
a(8) = 2 with 8 = 1^2 + 1^2 + 2^0 + 5*2^0 = 0^2 + 1^2 + 2^1 + 5*2^0.
a(9) = 2 with 9 = 1^2 + 1^2 + 2^1 + 5*2^0 = 0^2 + 0^2 + 2^2 + 5*2^0.
a(10) = 2 with 10 = 0^2 + 2^2 + 2^0 + 5*2^0 = 0^2 + 1^2 + 2^2 + 5*2^0.

R. C. Crocker, On the sum of two squares and two powers of k, Colloq. Math. 112(2008), 235-267.

Zhi-Wei Sun, Table of n, a(n) for n = 1..10000
Zhi-Wei Sun, Refining Lagrange's four-square theorem, J. Number Theory 175(2017), 167-190.
Zhi-Wei Sun, New conjectures on representations of integers (I), Nanjing Univ. J. Math. Biquarterly 34(2017), no. 2, 97-120.
Zhi-Wei Sun, Restricted sums of four squares, arXiv:1701.05868 [math.NT], 2017-2018.

Cf. A000079, A000290, A001481, A273812, A302982, A302984, A303233, A303234, A303338, A303363, A303389, A303393, A303399, A303428, A303401, A303432, A303434, A303539, A303540, A303541, A303543, A303601, A303639, A303656.

SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]];
f[n_]:=f[n]=FactorInteger[n];
g[n_]:=g[n]=Sum[Boole[Mod[Part[Part[f[n],i],1],4]==3&&Mod[Part[Part[f[n],i],2],2]==1],{i,1,Length[f[n]]}]==0;
QQ[n_]:=QQ[n]=(n==0)||(n>0&&g[n]);
tab={};Do[r=0;Do[If[QQ[n-5*2^k-2^m],Do[If[SQ[n-5*2^k-2^m-x^2],r=r+1],{x,0,Sqrt[(n-5*2^k-2^m)/2]}]],{k,0,Log[2,n/5]},{m,0,If[n/5==2^k,-1,Log[2,n-5*2^k]]}];tab=Append[tab,r],{n,1,60}];Print[tab]

A020756 Numbers that are the sum of two triangular numbers.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 9, 10, 11, 12, 13, 15, 16, 18, 20, 21, 22, 24, 25, 27, 28, 29, 30, 31, 34, 36, 37, 38, 39, 42, 43, 45, 46, 48, 49, 51, 55, 56, 57, 58, 60, 61, 64, 65, 66, 67, 69, 70, 72, 73, 76, 78, 79, 81, 83, 84, 87, 88, 90, 91, 92, 93, 94, 97, 99, 100, 101, 102, 105, 106, 108

Offset: 1

David W. Wilson

The possible sums of a square and a promic, i.e., x^2+n(n+1), e.g., 3^2 + 2*3 = 9 + 6 = 15 is present. - Jon Perry, May 28 2003
A052343(a(n)) > 0; union of A118139 and A119345. - Reinhard Zumkeller, May 15 2006
Also union of A051533 and A000217. - Ant King, Nov 29 2010

T. D. Noe, Table of n, a(n) for n = 1..10000
John A. Ewell, On Sums of Triangular Numbers and Sums of Squares, The American Mathematical Monthly, 99:8 (October 1992), pp. 752-757.
T. Khovanova, K. Knop, and A. Radul, Baron Munchhausen's Sequence, J. Int. Seq. 13 (2010) # 10.8.7.
L. K. Mork, Keith Sullivan, Trenton Vogt, and Darin J. Ulness, A group theoretical approach to the partitioning of integers: Application to triangular numbers, squares, and centered polygonal numbers, Australasian J. Comb. (2021) Vol. 80, No. 3, 305-321.

Complement of A020757.
Cf. A051533 (sums of two positive triangular numbers), A001481 (sums of two squares), A002378, A000217.
Cf. A052343.

a020756 n = a020756_list !! (n-1)
a020756_list = filter ((> 0) . a052343) [0..]
-- Reinhard Zumkeller, Jul 25 2014

q[k_] := If[! Head[Reduce[m (m + 1) + n (n + 1) == 2 k && 0 <= m && 0 <= n, {m, n}, Integers]] === Symbol, k, {}]; DeleteCases[Table[q[i], {i, 0, 108}], {}] (* Ant King, Nov 29 2010 *)
Take[Union[Total/@Tuples[Accumulate[Range[0,20]],2]],80] (* Harvey P. Dale, May 02 2012 *)

v=vector(200); vc=0; for (x=0,10, for (y=0,10,v[vc++ ]=x^2+y*(y+1))); v=vecsort(v); v

is(n)=my(f=factor(4*n+1));for(i=1,#f~,if(f[i,1]%4==3 && f[i,2]%2, return(0))); 1 \\ Charles R Greathouse IV, Jul 05 2013

Numbers n such that 4n+1 is the sum of two squares, i.e. such that 4n+1 is in A001481. Hence n is a member if and only if 4n+1 = odd square * product of distinct primes of form 4k+1. (Fred Helenius and others, Dec 18 2004)
Equivalently, we may say that a positive integer n can be partitioned into a sum of two triangular numbers if and only if every 4 k + 3 prime factor in the canonical form of 4 n + 1 occurs with an even exponent. - Ant King, Nov 29 2010
Also, the values of n for which 8n+2 can be partitioned into a sum of two squares of natural numbers. - Ant King, Nov 29 2010
Closed under the operation f(x, y) = 4*x*y + x + y.

Entry revised by N. J. A. Sloane, Dec 20 2004

A000925 Number of ordered ways of writing n as a sum of 2 squares of nonnegative integers.

Original entry on oeis.org

1, 2, 1, 0, 2, 2, 0, 0, 1, 2, 2, 0, 0, 2, 0, 0, 2, 2, 1, 0, 2, 0, 0, 0, 0, 4, 2, 0, 0, 2, 0, 0, 1, 0, 2, 0, 2, 2, 0, 0, 2, 2, 0, 0, 0, 2, 0, 0, 0, 2, 3, 0, 2, 2, 0, 0, 0, 0, 2, 0, 0, 2, 0, 0, 2, 4, 0, 0, 2, 0, 0, 0, 1, 2, 2, 0, 0, 0, 0, 0, 2, 2, 2, 0, 0, 4, 0, 0, 0, 2, 2, 0, 0, 0, 0, 0, 0, 2, 1, 0, 4

Offset: 0

Jacques Haubrich (jhaubrich(AT)freeler.nl)

A. Das and A. C. Melissinos, Quantum Mechanics: A Modern Introduction, Gordon and Breach, 1986, p. 47.
E. Grosswald, Representations of Integers as Sums of Squares. Springer-Verlag, NY, 1985.

Cf. A000290, A010052, A000161, A247367.

a000925 n = sum $ map (a010052 . (n -)) $ takeWhile (<= n) a000290_list
-- Reinhard Zumkeller, Sep 14 2014

a[n_] := (pr = PowersRepresentations[n, 2, 2]; Count[Union[Join[pr, Reverse /@ pr]], {j_ /; j >= 0, k_ /; k >= 0}]); a /@ Range[0, 100] (* Jean-François Alcover, Apr 05 2011 *)
nn = 100; t = CoefficientList[Series[Sum[x^k^2, {k, 0, Sqrt[nn]}]^2, {x, 0, nn}], x] (* T. D. Noe, Apr 05 2011 *)
SquareQ[n_] := IntegerQ[Sqrt[n]]; Table[Count[FrobeniusSolve[{1, 1}, n], {?SquareQ}], {n, 0, 100}] (* Robert G. Wilson v, Apr 15 2017 *)

a(n)=sum(i=0,n,sum(j=0,n,if(i^2+j^2-n,0,1)))

Coefficient of q^k in (1/4)*(1 + theta_3(0, q))^2.

a(A001481(n))>0; a(A022544(n))=0. - Benoit Cloitre, Apr 20 2003

A025284 Numbers that are the sum of 2 nonzero squares in exactly 1 way.

Original entry on oeis.org

2, 5, 8, 10, 13, 17, 18, 20, 25, 26, 29, 32, 34, 37, 40, 41, 45, 52, 53, 58, 61, 68, 72, 73, 74, 80, 82, 89, 90, 97, 98, 100, 101, 104, 106, 109, 113, 116, 117, 122, 128, 136, 137, 146, 148, 149, 153, 157, 160, 162, 164, 169, 173, 178, 180, 181, 193, 194, 197, 202, 208, 212

Offset: 1

David W. Wilson

nonn

A025426(a(n)) = 1. - Reinhard Zumkeller, Aug 16 2011

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Square Number
G. Xiao, Two squares
Index entries for sequences related to sums of squares

Cf. A000404, A025426, A081324, A125022, A018825, A001481, A007692, A243148.

import Data.List (elemIndices)
a025284 n = a025284_list !! (n-1)
a025284_list = elemIndices 1 a025426_list
-- Reinhard Zumkeller, Aug 16 2011

selQ[n_] := Length[ Select[ PowersRepresentations[n, 2, 2], Times @@ # != 0 &]] == 1; Select[Range[300], selQ] (* Jean-François Alcover, Oct 03 2013 *)
b[n_, i_, k_, t_] := b[n, i, k, t] = If[n == 0, If[t == 0, 1, 0], If[i<1 || t<1, 0, b[n, i - 1, k, t] + If[i^2 > n, 0, b[n - i^2, i, k, t - 1]]]];
T[n_, k_] := b[n, Sqrt[n] // Floor, k, k];
Position[Table[T[n, 2], {n, 0, 300}], 1] - 1 // Flatten (* Jean-François Alcover, Nov 06 2020, after Alois P. Heinz in A243148 *)

A243148(a(n),2) = 1. - Alois P. Heinz, Feb 25 2019

A133388 Largest integer m such that n-m^2 is a square, or 0 if no such m exists.

Original entry on oeis.org

1, 1, 0, 2, 2, 0, 0, 2, 3, 3, 0, 0, 3, 0, 0, 4, 4, 3, 0, 4, 0, 0, 0, 0, 5, 5, 0, 0, 5, 0, 0, 4, 0, 5, 0, 6, 6, 0, 0, 6, 5, 0, 0, 0, 6, 0, 0, 0, 7, 7, 0, 6, 7, 0, 0, 0, 0, 7, 0, 0, 6, 0, 0, 8, 8, 0, 0, 8, 0, 0, 0, 6, 8, 7, 0, 0, 0, 0, 0, 8, 9, 9, 0, 0, 9, 0, 0, 0, 8, 9, 0, 0, 0, 0, 0, 0, 9, 7, 0, 10

Offset: 1

M. F. Hasler, Nov 23 2007

nonn

The sequence could be extended to a(0) = 0.

We could have defined a(n) = -1 instead of 0 if n is not sum of two squares, and then include unambiguously a(0) = 0. At present, a(n) = 0 <=> A000161(n) = 0.

			a(3) = 0 since 3 cannot be written as sum of 2 perfect squares;
a(5) = 2 since 5 = 2^2 + 1^2.

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Cf. A000161, A001481.

a:= proc(n) local t, d;
      for t from 0 do d:= n-t^2;
        if d<0 then break elif issqr(d) then return isqrt(d) fi
      od; 0
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, May 14 2015

a[n_] := Module[{m, d, s}, For[m = 0, True, m++, d = n - m^2; If[d < 0, Break[], s = Sqrt[d]; If[IntegerQ[s], Return[s]]]]; 0];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, May 18 2018, after Alois P. Heinz *)

sum2sqr(n)={ if(n>1, my(L=List(), f, p=1); for(i=1, matsize(f=factor(n))[1], if(f[i,1]%4==1, listput(L, [qfbsolve(Qfb(1,0,1), f[i,1])*[1, I]~, f[i,2]] ),/*elseif*/ f[i,1]==2, p = (1+I)^f[i,2],/*elseif*/ bittest(f[i,2],0), return([]),/*else*/ p *= f[i,1]^(f[i,2]\2))); L=apply(s->vector(s[2]+1, j, s[1]^(s[2]+1-j)*conj(s[1])^(j-1)), L); my(S=List()); forvec(T=vector(#L, i, [1,#L[i]]), listput(S, prod( j=1, #T, L[j][T[j]] ))); Set(apply(f->vecsort(abs([real(f), imag(f)])), Set(S)*p)), if(n<0, [], [[0, n]]))} \\ updated by M. F. Hasler, May 12 2018. (If PARI version 2.12.x returns an error, append [1] to qfbsolve(...) above. - M. F. Hasler, Dec 12 2019)
apply( A133388=n->if(n=sum2sqr(n),vecmax(Mat(n~))), [1..50]) \\ This sequence: maximum

from sympy.solvers.diophantine.diophantine import diop_DN
def A133388(n): return max((a for a, b in diop_DN(-1,n)),default=0) # Chai Wah Wu, Sep 08 2022

a(n) = max( sup { max(a,b) | a^2+b^2 = n ; a,b in Z }, 0 )
a(A022544(j))=0, j>0. - R. J. Mathar, Jun 17 2009
a(n^2) = a(n^2 + 1) = n, for all n. Conversely, whenever a(n) = a(n+1), then n = k^2. - M. F. Hasler, Sep 02 2018

cp's OEIS Frontend

A303656 Number of ways to write n as a^2 + b^2 + 3^c + 5^d, where a,b,c,d are nonnegative integers with a <= b.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A303539 Number of ordered pairs (k, m) with 0 <= k <= m such that n - binomial(2*k,k) - binomial(2*m,m) can be written as the sum of two squares.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A303541 Numbers of the form k^2 + binomial(2*m,m) with k and m nonnegative integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A303543 Number of ways to write n as a^2 + b^2 + C(k) + C(m) with 0 <= a <= b and 0 < k <= m, where C(k) denotes the Catalan number binomial(2k,k)/(k+1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A007692 Numbers that are the sum of 2 nonzero squares in 2 or more ways.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

A303637 Number of ways to write n as x^2 + y^2 + 2^z + 5*2^w, where x,y,z,w are nonnegative integers with x <= y.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

A020756 Numbers that are the sum of two triangular numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Formula

Extensions

A303539 Number of ordered pairs (k, m) with 0 <= k <= m such that n - binomial(2k,k) - binomial(2m,m) can be written as the sum of two squares.