cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A022549 Sum of a square and a nonnegative cube.

Original entry on oeis.org

0, 1, 2, 4, 5, 8, 9, 10, 12, 16, 17, 24, 25, 26, 27, 28, 31, 33, 36, 37, 43, 44, 49, 50, 52, 57, 63, 64, 65, 68, 72, 73, 76, 80, 81, 82, 89, 91, 100, 101, 108, 113, 121, 122, 125, 126, 127, 128, 129, 134, 141, 144, 145
Offset: 1

Views

Author

Keywords

Comments

It appears that there are no modular constraints on this sequence; i.e., every residue class of every integer has representatives here. - Franklin T. Adams-Watters, Dec 03 2009
A045634(a(n)) > 0. - Reinhard Zumkeller, Jul 17 2010

Crossrefs

Complement of A022550; A002760 and A179509 are subsequences.

Programs

  • Mathematica
    q=30; imax=q^2; Select[Union[Flatten[Table[x^2+y^3, {y,0,q^(2/3)}, {x,0,q}]]], #<=imax&] (* Vladimir Joseph Stephan Orlovsky, Apr 20 2011 *)
  • PARI
    is(n)=for(k=0,sqrtnint(n,3), if(issquare(n-k^3), return(1))); 0 \\ Charles R Greathouse IV, Aug 24 2020
    
  • PARI
    list(lim)=my(v=List(),t); for(k=0,sqrtnint(lim\=1,3), t=k^3; for(n=0,sqrtint(lim-t), listput(v,t+n^2))); Set(v) \\ Charles R Greathouse IV, Aug 24 2020

A045634 Number of ways in which n can be partitioned as a sum of a square and cube.

Original entry on oeis.org

1, 2, 1, 0, 1, 1, 0, 0, 1, 2, 1, 0, 1, 0, 0, 0, 1, 2, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 2, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 2, 2, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 2, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 2
Offset: 0

Views

Author

Keywords

Comments

a(A022550(n))=0; a(A179509(n))=1; a(A022549(n))>0; a(A060861(n))=n. [From Reinhard Zumkeller, Jul 17 2010]

Examples

			a(9)=2 because 9=2^3+1^2 and 9=3^2+0^3.
		

Crossrefs

Programs

  • Maple
    M:=100; M2:=M^2; t0:=array(0..M2);
    for i from 0 to M2 do t0[i]:=0; od:
    for a from 0 to M do for b from 0 to M do
    i:=a^2+b^3; if i <= M2 then t0[i]:=t0[i]+1; fi; od: od:
    [seq(t0[i],i=0..M2)];
  • Mathematica
    max = 100; Clear[a]; a[_] = 0;
    Do[n = i^2 + j^3; a[n] += 1, {i, 0, Sqrt[max]}, {j, 0, max^(1/3)}];
    Table[a[n], {n, 0, max}] (* Jean-François Alcover, Aug 02 2018 *)

Extensions

More terms from Erich Friedman
Showing 1-2 of 2 results.