In number theory, the sum of the first n cubes is the square of the nth triangular number. That is,

The same equation may be written more compactly using the mathematical notation for summation:

This identity is sometimes called **Nicomachus's theorem**, after Nicomachus of Gerasa (c. 60 – c. 120 CE).

Nicomachus, at the end of Chapter 20 of his *Introduction to Arithmetic*, pointed out that if one writes a list of the odd numbers, the first is the cube of 1, the sum of the next two is the cube of 2, the sum of the next three is the cube of 3, and so on. He does not go further than this, but from this it follows that the sum of the first n cubes equals the sum of the first odd numbers, that is, the odd numbers from 1 to . The average of these numbers is obviously , and there are of them, so their sum is

Many early mathematicians have studied and provided proofs of Nicomachus's theorem. Stroeker (1995) claims that "every student of number theory surely must have marveled at this miraculous fact". Pengelley (2002) finds references to the identity not only in the works of Nicomachus in what is now Jordan in the first century CE, but also in those of Aryabhata in India in the fifth century, and in those of Al-Karaji circa 1000 in Persia. Bressoud (2004) mentions several additional early mathematical works on this formula, by Al-Qabisi (tenth century Arabia), Gersonides (circa 1300 France), and Nilakantha Somayaji (circa 1500 India); he reproduces Nilakantha's visual proof.

The sequence of squared triangular numbers is^{[1]}

These numbers can be viewed as figurate numbers, a four-dimensional hyperpyramidal generalization of the triangular numbers and square pyramidal numbers.

As Stein (1971) observes, these numbers also count the number of rectangles with horizontal and vertical sides formed in an *n* × *n* grid. For instance, the points of a 4 × 4 grid (or a square made up of three smaller squares on a side) can form 36 different rectangles. The number of squares in a square grid is similarly counted by the square pyramidal numbers.

The identity also admits a natural probabilistic interpretation as follows. Let *X*, *Y*, *Z*, *W* be four integer numbers independently and uniformly chosen at random between 1 and n. Then, the probability that W is the largest of the four numbers equals the probability that Y is at least as large as X and that W is at least as large as Z. That is, . For any particular value of W, the combinations of X, Y, and Z that make W largest form a cube 1 ≤ *X*, *Y*, *Z* ≤ *n* so (adding the size of this cube over all choices of W) the number of combinations of *X*, *Y*, *Z*, *W* for which W is largest is a sum of cubes, the left hand side of the Nichomachus identity. The sets of pairs (*X*,*Y*) with *X* ≤ *Y* and of pairs (*Z*,*W*) with *Z* ≤ *W* form isosceles right triangles, and the set counted by the right hand side of the equation of probabilities is the Cartesian product of these two triangles, so its size is the square of a triangular number on the right hand side of the Nichomachus identity. The probabilities themselves are respectively the left and right sides of the Nichomachus identity, normalized to make probabilities by dividing both sides by *n*^{4}.^{[citation needed]}

Charles Wheatstone (1854) gives a particularly simple derivation, by expanding each cube in the sum into a set of consecutive odd numbers. He begins by giving the identity

That identity is related to triangular numbers in the following way:

and thus the summands forming start off just after those forming all previous values up to .
Applying this property, along with another well-known identity:

produces the following derivation:

Row (1893) obtains another proof by summing the numbers in a square multiplication table in two different ways. The sum of the th row is times a triangular number, from which it follows that the sum of all the rows is the square of a triangular number. Alternatively, one can decompose the table into a sequence of nested gnomons, each consisting of the products in which the larger of the two terms is some fixed value. The sum within each gmonon is a cube, so the sum of the whole table is a sum of cubes.

In the more recent mathematical literature, Edmonds (1957) provides a proof using summation by parts. Stein (1971) uses the rectangle-counting interpretation of these numbers to form a geometric proof of the identity (see also Benjamin, Quinn & Wurtz 2006); he observes that it may also be proved easily (but uninformatively) by induction, and states that Toeplitz (1963) provides "an interesting old Arabic proof". Kanim (2004) provides a purely visual proof, Benjamin & Orrison (2002) provide two additional proofs, and Nelsen (1993) gives seven geometric proofs.

A similar result to Nicomachus's theorem holds for all power sums, namely that odd power sums (sums of odd powers) are a polynomial in triangular numbers.
These are called Faulhaber polynomials, of which the sum of cubes is the simplest and most elegant example.
However, in no other case is one power sum a square of another.^{[2]}

Stroeker (1995) studies more general conditions under which the sum of a consecutive sequence of cubes forms a square. Garrett & Hummel (2004) and Warnaar (2004) study polynomial analogues of the square triangular number formula, in which series of polynomials add to the square of another polynomial.