A003325 -id:A003325 - cp's OEIS frontend

A000578 The cubes: a(n) = n^3.

0, 1, 8, 27, 64, 125, 216, 343, 512, 729, 1000, 1331, 1728, 2197, 2744, 3375, 4096, 4913, 5832, 6859, 8000, 9261, 10648, 12167, 13824, 15625, 17576, 19683, 21952, 24389, 27000, 29791, 32768, 35937, 39304, 42875, 46656, 50653, 54872, 59319, 64000, 68921, 74088, 79507

Offset: 0

N. J. A. Sloane

a(n) is the sum of the next n odd numbers; i.e., group the odd numbers so that the n-th group contains n elements like this: (1), (3, 5), (7, 9, 11), (13, 15, 17, 19), (21, 23, 25, 27, 29), ...; then each group sum = n^3 = a(n). Also the median of each group = n^2 = mean. As the sum of first n odd numbers is n^2 this gives another proof of the fact that the n-th partial sum = (n(n + 1)/2)^2. - Amarnath Murthy, Sep 14 2002
Total number of triangles resulting from criss-crossing cevians within a triangle so that two of its sides are each n-partitioned. - Lekraj Beedassy, Jun 02 2004. See Propp and Propp-Gubin for a proof.
Also structured triakis tetrahedral numbers (vertex structure 7) (cf. A100175 = alternate vertex); structured tetragonal prism numbers (vertex structure 7) (cf. A100177 = structured prisms); structured hexagonal diamond numbers (vertex structure 7) (cf. A100178 = alternate vertex; A000447 = structured diamonds); and structured trigonal anti-diamond numbers (vertex structure 7) (cf. A100188 = structured anti-diamonds). Cf. A100145 for more on structured polyhedral numbers. - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Schlaefli symbol for this polyhedron: {4, 3}.
Least multiple of n such that every partial sum is a square. - Amarnath Murthy, Sep 09 2005
Draw a regular hexagon. Construct points on each side of the hexagon such that these points divide each side into equally sized segments (i.e., a midpoint on each side or two points on each side placed to divide each side into three equally sized segments or so on), do the same construction for every side of the hexagon so that each side is equally divided in the same way. Connect all such points to each other with lines that are parallel to at least one side of the polygon. The result is a triangular tiling of the hexagon and the creation of a number of smaller regular hexagons. The equation gives the total number of regular hexagons found where n = the number of points drawn + 1. For example, if 1 point is drawn on each side then n = 1 + 1 = 2 and a(n) = 2^3 = 8 so there are 8 regular hexagons in total. If 2 points are drawn on each side then n = 2 + 1 = 3 and a(n) = 3^3 = 27 so there are 27 regular hexagons in total. - Noah Priluck (npriluck(AT)gmail.com), May 02 2007
The solutions of the Diophantine equation: (X/Y)^2 - X*Y = 0 are of the form: (n^3, n) with n >= 1. The solutions of the Diophantine equation: (m^2)*(X/Y)^2k - XY = 0 are of the form: (m*n^(2k + 1), m*n^(2k - 1)) with m >= 1, k >= 1 and n >= 1. The solutions of the Diophantine equation: (m^2)*(X/Y)^(2k + 1) - XY = 0 are of the form: (m*n^(k + 1), m*n^k) with m >= 1, k >= 1 and n >= 1. - Mohamed Bouhamida, Oct 04 2007
Except for the first two terms, the sequence corresponds to the Wiener indices of C_{2n} i.e., the cycle on 2n vertices (n > 1). - K.V.Iyer, Mar 16 2009
Totally multiplicative sequence with a(p) = p^3 for prime p. - Jaroslav Krizek, Nov 01 2009
Sums of rows of the triangle in A176271, n > 0. - Reinhard Zumkeller, Apr 13 2010
One of the 5 Platonic polyhedral (tetrahedral, cube, octahedral, dodecahedral and icosahedral) numbers (cf. A053012). - Daniel Forgues, May 14 2010
Numbers n for which order of torsion subgroup t of the elliptic curve y^2 = x^3 - n is t = 2. - Artur Jasinski, Jun 30 2010
The sequence with the lengths of the Pisano periods mod k is 1, 2, 3, 4, 5, 6, 7, 8, 3, 10, 11, 12, 13, 14, 15, 16, 17, 6, 19, 20, ... for k >= 1, apparently multiplicative and derived from A000027 by dividing every ninth term through 3. Cubic variant of A186646. - R. J. Mathar, Mar 10 2011
The number of atoms in a bcc (body-centered cubic) rhombic hexahedron with n atoms along one edge is n^3 (T. P. Martin, Shells of atoms, eq. (8)). - Brigitte Stepanov, Jul 02 2011
The inverse binomial transform yields the (finite) 0, 1, 6, 6 (third row in A019538 and A131689). - R. J. Mathar, Jan 16 2013
Twice the area of a triangle with vertices at (0, 0), (t(n - 1), t(n)), and (t(n), t(n - 1)), where t = A000217 are triangular numbers. - J. M. Bergot, Jun 25 2013
If n > 0 is not congruent to 5 (mod 6) then A010888(a(n)) divides a(n). - Ivan N. Ianakiev, Oct 16 2013
For n > 2, a(n) = twice the area of a triangle with vertices at points (binomial(n,3),binomial(n+2,3)), (binomial(n+1,3),binomial(n+1,3)), and (binomial(n+2,3),binomial(n,3)). - J. M. Bergot, Jun 14 2014
Determinants of the spiral knots S(4,k,(1,1,-1)). a(k) = det(S(4,k,(1,1,-1))). - Ryan Stees, Dec 14 2014
One of the oldest-known examples of this sequence is shown in the Senkereh tablet, BM 92698, which displays the first 32 terms in cuneiform. - Charles R Greathouse IV, Jan 21 2015
From Bui Quang Tuan, Mar 31 2015: (Start)
We construct a number triangle from the integers 1, 2, 3, ... 2*n-1 as follows. The first column contains all the integers 1, 2, 3, ... 2*n-1. Each succeeding column is the same as the previous column but without the first and last items. The last column contains only n. The sum of all the numbers in the triangle is n^3.
Here is the example for n = 4, where 1 + 2*2 + 3*3 + 4*4 + 3*5 + 2*6 + 7 = 64 = a(4):
  1
  2  2
  3  3  3
  4  4  4  4
  5  5  5
  6  6
  7
(End)
For n > 0, a(n) is the number of compositions of n+11 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016
Does not satisfy Benford's law [Ross, 2012]. - N. J. A. Sloane, Feb 08 2017
Number of inequivalent face colorings of the cube using at most n colors such that each color appears at least twice. - David Nacin, Feb 22 2017
Consider A = {a,b,c} a set with three distinct members. The number of subsets of A is 8, including {a,b,c} and the empty set. The number of subsets from each of those 8 subsets is 27. If the number of such iterations is n, then the total number of subsets is a(n-1). - Gregory L. Simay, Jul 27 2018
By Fermat's Last Theorem, these are the integers of the form x^k with the least possible value of k such that x^k = y^k + z^k never has a solution in positive integers x, y, z for that k. - Felix Fröhlich, Jul 27 2018

			For k=3, b(3) = 2 b(2) - b(1) = 4-1 = 3, so det(S(4,3,(1,1,-1))) = 3*3^2 = 27.
For n=3, a(3) = 3 + (3*0^2 + 3*0 + 3*1^2 + 3*1 + 3*2^2 + 3*2) = 27. - _Patrick J. McNab_, Mar 28 2016

Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See p. 191.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 43, 64, 81.
R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 255; 2nd. ed., p. 269. Worpitzky's identity (6.37).
Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 292.
T. Aaron Gulliver, "Sequences from cubes of integers", International Mathematical Journal, 4 (2003), no. 5, 439 - 445. See http://www.m-hikari.com/z2003.html for information about this journal. [I expanded the reference to make this easier to find. - N. J. A. Sloane, Feb 18 2019]
J. Propp and A. Propp-Gubin, "Counting Triangles in Triangles", Pi Mu Epsilon Journal (to appear).
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 6-7.
D. Wells, You Are A Mathematician, pp. 238-241, Penguin Books 1995.

N. J. A. Sloane, Table of n, a(n) for n = 0..10000
H. Bottomley, Illustration of initial terms
British National Museum, Tablet 92698
N. Brothers, S. Evans, L. Taalman, L. Van Wyk, D. Witczak, and C. Yarnall, Spiral knots, Missouri J. of Math. Sci., 22 (2010).
M. DeLong, M. Russell, and J. Schrock, Colorability and determinants of T(m,n,r,s) twisted torus knots for n equiv. +/-1(mod m), Involve, Vol. 8 (2015), No. 3, 361-384.
Ralph Greenberg, Math For Poets
R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712. [Annotated scanned copy]
Milan Janjic, Enumerative Formulas for Some Functions on Finite Sets [Cached version at the Wayback Machine]
Hyun Kwang Kim, On Regular Polytope Numbers, Proc. Amer. Math. Soc., 131 (2002), 65-75. - fixed by _Felix Fröhlich_, Jun 16 2014
T. P. Martin, Shells of atoms, Phys. Reports, 273 (1996), 199-241, eq. (8).
Ed Pegg, Jr., Sequence Pictures, Math Games column, Dec 08 2003.
Ed Pegg, Jr., Sequence Pictures, Math Games column, Dec 08 2003 [Cached copy, with permission (pdf only)]
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
James Propp and Adam Propp-Gubin, Counting Triangles in Triangles, arXiv:2409.17117 [math.CO], 25 September 2024.
Kenneth A. Ross, First Digits of Squares and Cubes, Math. Mag. 85 (2012) 36-42. doi:10.4169/math.mag.85.1.36.
Eric Weisstein's World of Mathematics, Cubic Number, and Hex Pyramidal Number
Ronald Yannone, Hilbert Matrix Analyses
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
Index entries for "core" sequences
Index entries for sequences related to Benford's law

(1/12)*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, A006003, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
For sums of cubes, cf. A000537 (partial sums), A003072, A003325, A024166, A024670, A101102 (fifth partial sums).
Cf. A001158 (inverse Möbius transform), A007412 (complement), A030078(n) (cubes of primes), A048766, A058645 (binomial transform), A065876, A101094, A101097.
Subsequence of A145784.
Cf. A260260 (comment). - Bruno Berselli, Jul 22 2015
Cf. A000292 (tetrahedral numbers), A005900 (octahedral numbers), A006566 (dodecahedral numbers), A006564 (icosahedral numbers).
Cf. A098737 (main diagonal).

a000578 = (^ 3)
a000578_list = 0 : 1 : 8 : zipWith (+)
   (map (+ 6) a000578_list)
   (map (* 3) $ tail $ zipWith (-) (tail a000578_list) a000578_list)
-- Reinhard Zumkeller, Sep 05 2015, May 24 2012, Oct 22 2011

[ n^3 : n in [0..50] ]; // Wesley Ivan Hurt, Jun 14 2014

I:=[0,1,8,27]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..45]]; // Vincenzo Librandi, Jul 05 2014

A000578 := n->n^3;
seq(A000578(n), n=0..50);
isA000578 := proc(r)
    local p;
    if r = 0 or r =1 then
        true;
    else
        for p in ifactors(r)[2] do
            if op(2, p) mod 3 <> 0 then
                return false;
            end if;
        end do:
        true ;
    end if;
end proc: # R. J. Mathar, Oct 08 2013

Table[n^3, {n, 0, 30}] (* Stefan Steinerberger, Apr 01 2006 *)
CoefficientList[Series[x (1 + 4 x + x^2)/(1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Jul 05 2014 *)
Accumulate[Table[3n^2+3n+1,{n,0,20}]] (* or *) LinearRecurrence[{4,-6,4,-1},{1,8,27,64},20](* Harvey P. Dale, Aug 18 2018 *)

A000578(n):=n^3$
makelist(A000578(n),n,0,30); /* Martin Ettl, Nov 03 2012 */

A000578(n)=n^3 \\ M. F. Hasler, Apr 12 2008

is(n)=ispower(n,3) \\ Charles R Greathouse IV, Feb 20 2012

A000578_list, m = [], [6, -6, 1, 0]
for _ in range(10**2):
    A000578_list.append(m[-1])
    for i in range(3):
        m[i+1] += m[i] # Chai Wah Wu, Dec 15 2015

(define (A000578 n) (* n n n)) ;; Antti Karttunen, Oct 06 2017

a(n) = Sum_{i=0..n-1} A003215(i).
Multiplicative with a(p^e) = p^(3e). - David W. Wilson, Aug 01 2001
G.f.: x*(1+4*x+x^2)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
Dirichlet generating function: zeta(s-3). - Franklin T. Adams-Watters, Sep 11 2005, Amarnath Murthy, Sep 09 2005
E.g.f.: (1+3*x+x^2)*x*exp(x). - Franklin T. Adams-Watters, Sep 11 2005 - Amarnath Murthy, Sep 09 2005
a(n) = Sum_{i=1..n} (Sum_{j=i..n+i-1} A002024(j,i)). - Reinhard Zumkeller, Jun 24 2007
a(n) = lcm(n, (n - 1)^2) - (n - 1)^2. E.g.: lcm(1, (1 - 1)^2) - (1 - 1)^2 = 0, lcm(2, (2 - 1)^2) - (2 - 1)^2 = 1, lcm(3, (3 - 1)^2) - (3 - 1)^2 = 8, ... - Mats Granvik, Sep 24 2007
Starting (1, 8, 27, 64, 125, ...), = binomial transform of [1, 7, 12, 6, 0, 0, 0, ...]. - Gary W. Adamson, Nov 21 2007
a(n) = A007531(n) + A000567(n). - Reinhard Zumkeller, Sep 18 2009
a(n) = binomial(n+2,3) + 4*binomial(n+1,3) + binomial(n,3). [Worpitzky's identity for cubes. See. e.g., Graham et al., eq. (6.37). - Wolfdieter Lang, Jul 17 2019]
a(n) = n + 6*binomial(n+1,3) = binomial(n,1)+6*binomial(n+1,3). - Ron Knott, Jun 10 2019
A010057(a(n)) = 1. - Reinhard Zumkeller, Oct 22 2011
a(n) = A000537(n) - A000537(n-1), difference between 2 squares of consecutive triangular numbers. - Pierre CAMI, Feb 20 2012
a(n) = A048395(n) - 2*A006002(n). - J. M. Bergot, Nov 25 2012
a(n) = 1 + 7*(n-1) + 6*(n-1)*(n-2) + (n-1)*(n-2)*(n-3). - Antonio Alberto Olivares, Apr 03 2013
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6. - Ant King Apr 29 2013
a(n) = A000330(n) + Sum_{i=1..n-1} A014105(i), n >= 1. - Ivan N. Ianakiev, Sep 20 2013
a(k) = det(S(4,k,(1,1,-1))) = k*b(k)^2, where b(1)=1, b(2)=2, b(k) = 2*b(k-1) - b(k-2) = b(2)*b(k-1) - b(k-2). - Ryan Stees, Dec 14 2014
For n >= 1, a(n) = A152618(n-1) + A033996(n-1). - Bui Quang Tuan, Apr 01 2015
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Jon Tavasanis, Feb 21 2016
a(n) = n + Sum_{j=0..n-1} Sum_{k=1..2} binomial(3,k)*j^(3-k). - Patrick J. McNab, Mar 28 2016
a(n) = A000292(n-1) * 6 + n. - Zhandos Mambetaliyev, Nov 24 2016
a(n) = n*binomial(n+1, 2) + 2*binomial(n+1, 3) + binomial(n,3). - Tony Foster III, Nov 14 2017
From Amiram Eldar, Jul 02 2020: (Start)
Sum_{n>=1} 1/a(n) = zeta(3) (A002117).
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*zeta(3)/4 (A197070). (End)
From Amiram Eldar, Jan 20 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)/Pi.
Product_{n>=2} (1 - 1/a(n)) = cosh(sqrt(3)*Pi/2)/(3*Pi). (End)
a(n) = Sum_{d|n} sigma_3(d)*mu(n/d) = Sum_{d|n} A001158(d)*A008683(n/d). Moebius transform of sigma_3(n). - Ridouane Oudra, Apr 15 2021

A000404 Numbers that are the sum of 2 nonzero squares.

Original entry on oeis.org

2, 5, 8, 10, 13, 17, 18, 20, 25, 26, 29, 32, 34, 37, 40, 41, 45, 50, 52, 53, 58, 61, 65, 68, 72, 73, 74, 80, 82, 85, 89, 90, 97, 98, 100, 101, 104, 106, 109, 113, 116, 117, 122, 125, 128, 130, 136, 137, 145, 146, 148, 149, 153, 157, 160, 162, 164, 169, 170, 173, 178

Offset: 1

N. J. A. Sloane and J. H. Conway

From the formula it is easy to see that if k is in this sequence, then so are all odd powers of k. - T. D. Noe, Jan 13 2009
Also numbers whose cubes are the sum of two nonzero squares. - Joe Namnath and Lawrence Sze
A line perpendicular to y=mx has its first integral y-intercept at a^2+b^2. The remaining ones for that slope are multiples of that primitive value. - Larry J Zimmermann, Aug 19 2010
The primes in this sequence are sequence A002313.
Complement of A018825; A025426(a(n)) > 0; A063725(a(n)) > 0. - Reinhard Zumkeller, Aug 16 2011
If the two squares are not equal, then any power is still in the sequence: if k = x^2 + y^2 with x != y, then k^2 = (x^2-y^2)^2 + (2xy)^2 and k^3 = (x(x^2-3y^2))^2 + (y(3x^2-y^2))^2, etc. - Carmine Suriano, Jul 13 2012
There are never more than 3 consecutive terms that differ by 1. Triples of consecutive terms that differ by 1 occur infinitely many times, for example, 2(k^2 + k)^2, (k^2 - 1)^2 + (k^2 + 2 k)^2, and (k^2 + k - 1)^2 + (k^2 + k + 1)^2 for any integer k > 1. - Ivan Neretin, Mar 16 2017 [Corrected by Jerzy R Borysowicz, Apr 14 2017]
Number of terms less than 10^k, k=1,2,3,...: 3, 34, 308, 2690, 23873, 215907, 1984228, ... - Muniru A Asiru, Feb 01 2018
The squares in this sequence are the squares of the so-called hypotenuse numbers A009003. - M. F. Hasler, Jun 20 2025

			25 = 3^2 + 4^2, therefore 25 is a term. Note that also 25^3 = 15625 = 44^2 + 117^2, therefore 15625 is a term.

David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989.
GCHQ, The GCHQ Puzzle Book, Penguin, 2016. See page 103.
E. Grosswald, Representations of Integers as Sums of Squares. Springer-Verlag, NY, 1985, p. 75, Theorem 4, with Theorem 2, p. 15.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, p. 219, th. 251, 252.
Ian Stewart, "Game, Set and Math", Chapter 8, 'Close Encounters of the Fermat Kind', Penguin Books, Ed. 1991, pp. 107-124.

T. D. Noe, Table of n, a(n) for n = 1..10000
J. M. De Koninck and V. Ouellet, On the n-th element of a set of positive integers, Annales Univ. Sci. Budapest Sect. Comput. 44 (2015), 153-164. See 2. on p. 162.
Étienne Fouvry, Claude Levesque, and Michel Waldschmidt, Representation of integers by cyclotomic binary forms, arXiv:1712.09019 [math.NT], 2017.
Joshua Harrington, Lenny Jones, and Alicia Lamarche, Representing integers as the sum of two squares in the ring Z_n, arXiv:1404.0187 [math.NT], 2014.
David Rabahy, Google Sheets
G. Xiao, Two squares.
Reinhard Zumkeller, Illustration for A084888 and A000404.
Index entries for sequences related to sums of squares

A001481 gives another version (allowing for zero squares).
Cf. A004431 (2 distinct squares), A063725 (number of representations), A024509 (numbers with multiplicity), A025284, A018825. Also A050803, A050801, A001105, A033431, A084888, A000578, A000290, A057961, A232499, A007692.
Cf. A003325 (analog for cubes), A003336 (analog for 4th powers).
Cf. A009003 (square roots of the squares in this sequence).
Column k=2 of A336725.
Cf. A355237, A355238.

P:=List([1..10^4],i->i^2);;
A000404 := Set(Flat(List(P, i->List(P, j -> i+j)))); # Muniru A Asiru, Feb 01 2018

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

lst:=[]; for n in [1..178] do f:=Factorization(n); if IsSquare(n) then for m in [1..#f] do d:=f[m]; if d[1] mod 4 eq 1 then Append(~lst, n); break; end if; end for; else t:=0; for m in [1..#f] do d:=f[m]; if d[1] mod 4 eq 3 and d[2] mod 2 eq 1 then t:=1; break; end if; end for; if t eq 0 then Append(~lst, n); end if; end if; end for; lst; // Arkadiusz Wesolowski, Feb 16 2017

nMax:=178: A:={}: for i to floor(sqrt(nMax)) do for j to floor(sqrt(nMax)) do if i^2+j^2 <= nMax then A := `union`(A, {i^2+j^2}) else  end if end do end do: A; # Emeric Deutsch, Jan 02 2017

nMax=1000; n2=Floor[Sqrt[nMax-1]]; Union[Flatten[Table[a^2+b^2, {a,n2}, {b,a,Floor[Sqrt[nMax-a^2]]}]]]
Select[Range@ 200, Length[PowersRepresentations[#, 2, 2] /. {0, } -> Nothing] > 0 &] (* _Michael De Vlieger, Mar 24 2016 *)
Module[{upto=200},Select[Union[Total/@Tuples[Range[Sqrt[upto]]^2,2]],#<= upto&]] (* Harvey P. Dale, Sep 18 2021 *)

is_A000404(n)= for( i=1,#n=factor(n)~%4, n[1,i]==3 && n[2,i]%2 && return); n && ( vecmin(n[1,])==1 || (n[1,1]==2 && n[2,1]%2)) \\ M. F. Hasler, Feb 07 2009

list(lim)=my(v=List(),x2); lim\=1; for(x=1,sqrtint(lim-1), x2=x^2; for(y=1,sqrtint(lim-x2), listput(v,x2+y^2))); Set(v) \\ Charles R Greathouse IV, Apr 30 2016

from itertools import count, islice
from sympy import factorint
def A000404_gen(startvalue=1): # generator of terms >= startvalue
    for n in count(max(startvalue,1)):
        c = False
        for p in (f:=factorint(n)):
            if (q:= p & 3)==3 and f[p]&1:
                break
            elif q == 1:
                c = True
        else:
            if c or f.get(2,0)&1:
                yield n
A000404_list = list(islice(A000404_gen(),30)) # Chai Wah Wu, Jul 01 2022

Let k = 2^t * p_1^a_1 * p_2^a_2 * ... * p_r^a_r * q_1^b_1 * q_2^b_2 * ... * q_s^b_s with t >= 0, a_i >= 0 for i=1..r, where p_i == 1 (mod 4) for i=1..r and q_j == -1 (mod 4) for j=1..s. Then k is a term iff 1) b_j == 0 (mod 2) for j=1..s and 2) r > 0 or t == 1 (mod 2) (or both).
From Charles R Greathouse IV, Nov 18 2022: (Start)
a(n) ~ k*n*sqrt(log n), where k = 1.3085... = 1/A064533.
There are B(x) = (x/sqrt(log x)) * (K + B2/log x + O(1/log^2 x)) terms of this sequence up to x, where K = A064533 and B2 = A227158. (End)

Edited by Ralf Stephan, Nov 15 2004
Typo in formula corrected by M. F. Hasler, Feb 07 2009
Erroneous Mathematica program fixed by T. D. Noe, Aug 07 2009
PARI code fixed for versions > 2.5 by M. F. Hasler, Jan 01 2013

A001235 Taxi-cab numbers: sums of 2 cubes in more than 1 way.

Original entry on oeis.org

1729, 4104, 13832, 20683, 32832, 39312, 40033, 46683, 64232, 65728, 110656, 110808, 134379, 149389, 165464, 171288, 195841, 216027, 216125, 262656, 314496, 320264, 327763, 373464, 402597, 439101, 443889, 513000, 513856, 515375, 525824, 558441, 593047, 684019, 704977

Offset: 1

N. J. A. Sloane

From Wikipedia: "1729 is known as the Hardy-Ramanujan number after a famous anecdote of the British mathematician G. H. Hardy regarding a hospital visit to the Indian mathematician Srinivasa Ramanujan. In Hardy's words: 'I remember once going to see him when he was ill at Putney. I had ridden in taxi cab number 1729 and remarked that the number seemed to me rather a dull one, and that I hoped it was not an unfavorable omen. "No," he replied, "it is a very interesting number; it is the smallest number expressible as the sum of two cubes in two different ways."'"
A011541 gives another version of "taxicab numbers".
If n is in this sequence, then n*k^3 is also in this sequence for all k > 0. So this sequence is obviously infinite. - Altug Alkan, May 09 2016

			4104 belongs to the sequence as 4104 = 2^3 + 16^3 = 9^3 + 15^3.

R. K. Guy, Unsolved Problems in Number Theory, Section D1.
G. H. Hardy, Ramanujan, Cambridge Univ. Press, 1940, p. 12.
Ya. I. Perelman, Algebra can be fun, pp. 142-143.
H. W. Richmond, On integers which satisfy the equation t^3 +- x^3 +- y^3 +- z^3, Trans. Camb. Phil. Soc., 22 (1920), 389-403, see p. 402.
D. Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, 165.

Shahar Amitai, Table of n, a(n) for n = 1..30000 (terms a(1)-a(4724) from T. D. Noe, terms a(4725)-a(10000) from Zak Seidov).
Shahar Amitai, Python code to generate all taxicab numbers up to N.
J. Charles-É, Recreomath, Ramanujan's Number.
A. Grinstein, Ramanujan and 1729, University of Melbourne Dept. of Math and Statistics Newsletter: Issue 3, 1998.
Henk Koppelaar, Peyman Nasehpour, and Maarten Looijen, Symmetry between Series if Entangled by Sums, Preprints.org, 2024.
Istanbul Bilgi University, Ramanujan and Hardy's Taxi
Christopher Lane, The First ten Ta(2) and their double distinct cubic sums representations, Find Ramanujan's Taxi Number using JavaScript. [WayBack Machine copy]
J. Leech, Some solutions of Diophantine equations, Proc. Camb. Phil. Soc., 53 (1957), 778-780.
J. Loy, The Hardy-Ramanujan Number.
Mia Muessig, Julia code for finding general taxicab numbers
Ken Ono and Sarah Trebat-Leder, The 1729 K3 surface, arXiv:1510.00735 [math.NT], 2015.
Eric Weisstein's World of Mathematics, Cubic Number
Eric Weisstein's World of Mathematics, Diophantine Equation 3rd Powers
Eric Weisstein's World of Mathematics, Taxicab Number
D. W. Wilson, The Fifth Taxicab Number is 48988659276962496, J. Integer Sequences, Vol. 2, 1999, #99.1.9.

Subsequence of A003325.
Cf. A007692, A008917, A011541, A018786, A018850 (primitive solutions), A051347 (allows negatives), A343708, A360619.
Solutions in greater numbers of ways:
(>2): A018787 (A003825 for primitive, A023050 for coprime),
(>3): A023051 (A003826 for primitive),
(>4): A051167 (A155057 for primitive).

Select[Range[750000],Length[PowersRepresentations[#,2,3]]>1&] (* Harvey P. Dale, Nov 25 2014, with correction by Zak Seidov, Jul 13 2015 *)

is(n)=my(t);for(k=ceil((n/2)^(1/3)),(n-.4)^(1/3),if(ispower(n-k^3,3),if(t,return(1),t=1)));0 \\ Charles R Greathouse IV, Jul 15 2011

T=thueinit(x^3+1,1);
is(n)=my(v=thue(T,n)); sum(i=1,#v,v[i][1]>=0 && v[i][2]>=v[i][1])>1 \\ Charles R Greathouse IV, May 09 2016

A003072 Numbers that are the sum of 3 positive cubes.

Original entry on oeis.org

3, 10, 17, 24, 29, 36, 43, 55, 62, 66, 73, 80, 81, 92, 99, 118, 127, 129, 134, 136, 141, 153, 155, 160, 179, 190, 192, 197, 216, 218, 225, 232, 244, 251, 253, 258, 270, 277, 281, 288, 307, 314, 342, 344, 345, 349, 352, 359, 368, 371, 375, 378, 397, 405, 408, 415, 433, 434

Offset: 1

N. J. A. Sloane, David W. Wilson

A119977 is a subsequence; if m is a term then there exists at least one k>0 such that m-k^3 is a term of A003325. - Reinhard Zumkeller, Jun 03 2006
A025456(a(n)) > 0. - Reinhard Zumkeller, Apr 23 2009
Davenport proved that a(n) << n^(54/47 + e) for every e > 0. - Charles R Greathouse IV, Mar 26 2012

			a(11) = 73 = 1^3 + 2^3 + 4^3, which is sum of three cubes.
a(15) = 99 = 2^3 + 3^3 + 4^3, which is sum of three cubes.

K. D. Bajpai, Table of n, a(n) for n = 1..12955 (first 1000 terms from T. D. Noe)
H. Davenport, Sums of three positive cubes, J. London Math. Soc., 25 (1950), 339-343. Coll. Works III p. 999.
Eric Weisstein's World of Mathematics, Cubic Number
Index entries for sequences related to sums of cubes

Subsequence of A004825.
Cf. A003325, A024981, A057904 (complement), A010057, A000578, A023042 (subsequence of cubes).
Cf. A###### (x, y) = Numbers that are the sum of x nonzero y-th powers:
- squares: A000404 (2, 2), A000408 (3, 2), A000414 (4, 2), A047700 (5, 2);
- cubes: A003325 (2, 3), A003072 (3, 3), A003327 (4, 3), A003328 (5, 3), A003329 (6, 3), A003330 (7, 3), A003331 (8, 3), A003332 (9, 3), A003333 (10, 3), A003334 (11, 3), A003335 (12, 3);
- fourth powers: A003336 (2, 4), A003337 (3, 4), A003338 (4, 4), A003339 (5, 4), A003340 (6, 4), A003341 (7, 4), A003342 (8, 4), A003343 (9, 4), A003344 (10, 4), A003345 (11, 4), A003346 (12, 4);
- fifth powers: A003347 (2, 5), A003348 (3, 5), A003349 (4, 5), A003350 (5, 5), A003351 (6, 5), A003352 (7, 5), A003353 (8, 5), A003354 (9, 5), A003355 (10, 5), A003356 (11, 5), A003357 (12, 5);
- sixth powers: A003358 (2, 6), A003359 (3, 6), A003360 (4, 6), A003361 (5, 6), A003362 (6, 6), A003363 (7, 6), A003364 (8, 6), A003365 (9, 6), A003366 (10, 6), A003367 (11, 6), A003368 (12, 6);
- seventh powers: A003369 (2, 7), A003370 (3, 7), A003371 (4, 7), A003372 (5, 7), A003373 (6, 7), A003374 (7, 7), A003375 (8, 7), A003376 (9, 7), A003377 (10, 7), A003378 (11, 7), A003379 (12, 7);
- eighth powers: A003380 (2, 8), A003381 (3, 8), A003382 (4, 8), A003383 (5, 8), A003384 (6, 8), A003385 (7, 8), A003386 (8, 8), A003387 (9, 8), A003388 (10, 8), A003389 (11, 8), A003390 (12, 8);
- ninth powers: A003391 (2, 9), A003392 (3, 9), A003393 (4, 9), A003394 (5, 9), A003395 (6, 9), A003396 (7, 9), A003397 (8, 9), A003398 (9, 9), A003399 (10, 9), A004800 (11, 9), A004801 (12, 9);
- tenth powers: A004802 (2, 10), A004803 (3, 10), A004804 (4, 10), A004805 (5, 10), A004806 (6, 10), A004807 (7, 10), A004808 (8, 10), A004809 (9, 10), A004810 (10, 10), A004811 (11, 10), A004812 (12, 10);
- eleventh powers: A004813 (2, 11), A004814 (3, 11), A004815 (4, 11), A004816 (5, 11), A004817 (6, 11), A004818 (7, 11), A004819 (8, 11), A004820 (9, 11), A004821 (10, 11), A004822 (11, 11), A004823 (12, 11).

a003072 n = a003072_list !! (n-1)
a003072_list = filter c3 [1..] where
   c3 x = any (== 1) $ map (a010057 . fromInteger) $
                       takeWhile (> 0) $ map (x -) $ a003325_list
-- Reinhard Zumkeller, Mar 24 2012

isA003072 := proc(n)
    local x,y,z;
    for x from 1 do
        if 3*x^3 > n then
            return false;
        end if;
        for y from x do
            if x^3+2*y^3 > n then
                break;
            end if;
            if isA000578(n-x^3-y^3) then
                return true;
            end if;
        end do:
    end do:
end proc:
for n from 1 to 1000 do
    if isA003072(n) then
        printf("%d,",n) ;
    end if;
end do: # R. J. Mathar, Jan 23 2016

Select[Range[435], (p = PowersRepresentations[#, 3, 3]; (Select[p, #[[1]] > 0 && #[[2]] > 0 && #[[3]] > 0 &] != {})) &] (* Jean-François Alcover, Apr 29 2011 *)
With[{upto=500},Select[Union[Total/@Tuples[Range[Floor[Surd[upto-2,3]]]^3,3]],#<=upto&]] (* Harvey P. Dale, Oct 25 2021 *)

sum(n=1,11,x^(n^3),O(x^1400))^3 /* Then [i|i<-[1..#%],polcoef(%,i)] gives the list of powers with nonzero coefficient. - M. F. Hasler, Aug 02 2020 */

list(lim)=my(v=List(),k,t); lim\=1; for(x=1,sqrtnint(lim-2,3), for(y=1, min(sqrtnint(lim-x^3-1,3),x), k=x^3+y^3; for(z=1,min(sqrtnint(lim-k,3), y), listput(v, k+z^3)))); Set(v) \\ Charles R Greathouse IV, Sep 14 2015

{n: A025456(n) >0}. - R. J. Mathar, Jun 15 2018

Incorrect program removed by David A. Corneth, Aug 01 2020

cp's OEIS Frontend

A024666 Positions of odd numbers in A003325.

Views

Author

Keywords

A024668 Position of 1 + n^3 in A003325.

Views

Author

Keywords

A024669 Position of n^3 + (n+1)^3 in A003325.

Views

Author

Keywords

A139211 Partial sums of A003325.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Extensions

A143186 A triangle sequence of holes in a triangle sequence numbering of n^3+m^3 (A003325) by comparing the n level to the 2*n level.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A273186 Taxi-cab numbers (A001235) n such that n+1 is the sum of two positive cubes (A003325).

Views

Author

Keywords

Comments

Examples

Crossrefs

A000578 The cubes: a(n) = n^3.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Scheme

Formula

A000404 Numbers that are the sum of 2 nonzero squares.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions