The method of Eratosthenes used to sieve out prime numbers is employed in this proof.
This sketch of a proof makes use of simple algebra only. This was the method by which Euler originally discovered the formula. There is a certain sieving property that we can use to our advantage:
Subtracting the second equation from the first we remove all elements that have a factor of 2:
Repeating for the next term:
Subtracting again we get:
where all elements having a factor of 3 or 2 (or both) are removed.
It can be seen that the right side is being sieved. Repeating infinitely for where is prime, we get:
Dividing both sides by everything but the ζ(s) we obtain:
This can be written more concisely as an infinite product over all primes p:
To make this proof rigorous, we need only to observe that when , the sieved right-hand side approaches 1, which follows immediately from the convergence of the Dirichlet series for .
An interesting result can be found for ζ(1), the harmonic series:
which can also be written as,
which is,
as,
thus,
While the series ratio test is inconclusive for the left-hand side it may be shown divergent by bounding logarithms. Similarly for the right-hand side the infinite coproduct of reals greater than one does not guarantee divergence, e.g.,
.
Instead, the denominator may be written in terms of the primorial numerator so that divergence is clear
Each factor (for a given prime p) in the product above can be expanded to a geometric series consisting of the reciprocal of p raised to multiples of s, as follows
When , this series converges absolutely. Hence we may take a finite number of factors, multiply them together, and rearrange terms. Taking all the primes p up to some prime number limit q, we have
where σ is the real part of s. By the fundamental theorem of arithmetic, the partial product when expanded out gives a sum consisting of those terms n−s where n is a product of primes less than or equal to q. The inequality results from the fact that therefore only integers larger than q can fail to appear in this expanded out partial product. Since the difference between the partial product and ζ(s) goes to zero when σ > 1, we have convergence in this region.