# Thread: Sums of Powers of Positive Integers

1. ## Sums of Powers of Positive Integers

Let's consider sums of powers of positive integers.
$\sum_{k=1}^n 1 = n \\ \sum_{k=1}^n k = \frac12 n (n+1) \\ \sum_{k=1}^n k^2 = \frac16 n (n+1) (2n+1) \\ \sum_{k=1}^n k^3 = \frac14 n^2 (n+1)^2 \\ \sum_{k=1}^n k^4 = \frac{1}{30} n (n+1) (2n+1) (3n^2 + 3n - 1)$

Not very much of a pattern. Let's look at something similar.

$\sum_{k=1}^n k (k+1) \cdots (k+m) = \frac{1}{m+2} n (n+1) \cdots (n+m+1)$

Easy to prove. But Abramowitz and Stegun have some general formulas for positive-integer powers, though they use something called Bernoulli polynomials, Bn(x). These are defined by

$\frac{t e^{x t}}{e^t - 1} = \sum_{n=0}^{\infty} B_n(x) \frac{t^n} {n!}$

Bernoulli numbers Bn are Bn(0). Some related polynomials are the Euler polynomials:

$\frac{2e^{x t}}{e^t + 1} = \sum_{n=0}^{\infty} E_n(x) \frac{t^n} {n!}$

However, the corresponding Euler numbers are En = 2n En(1/2).

The sum-of-powers formula:

$\sum_{k=0}^n k^m = \frac{ B_{m+1}(n+1) - B_{m+1} }{m+1}$

2. Bernoulli and Euler numbers fill out the series expressions for the trigonometric functions. The best-known series are

$\sin x = \sum_{k=0}^\infty (-1)^k \frac{x^{2k+1}}{(2k+1)!} \\ \cos x = \sum_{k=0}^\infty (-1)^k \frac{x^{2k}}{(2k)!}$

For the remaining four functions, the series are
$\tan x = \sum_{k=0}^\infty (-1)^{k+1} 2^{2k} (2^{2k}-1) B_{2k} \frac{x^{2k-1}}{(2k)!} \\ \cot x = \sum_{k=0}^\infty (-1)^{k} 2^{2k} B_{2k} \frac{x^{2k-1}}{(2k)!} \\ \sec x = \sum_{k=0}^\infty (-1)^k E_{2k} \frac{x^{2k}}{(2k)!} \\ \csc x = \sum_{k=0}^\infty (-1)^{k+1} 2 (2^{2k-1}-1) B_{2k} \frac{x^{2k-1}}{(2k)!} \\$

For products and powers of trigonometric functions, one can use trigonometric identities and differentiation, like: cos(x)2 = (1/2)*(1 + cos(2x)) and sec(x)2 = (d/dx) tan(x).

3. Negative powers of positive integers have some interesting features. Let's consider reciprocals first.

$\sum_{k=1}^n \frac{1}{k} \to \log n + \gamma$

as n -> oo, where the last term is the Euler-Mascheroni constant, 0.577216...

It can be generalized to

$\sum_{k=1}^n \frac{(\log k)^m}{k} \to \frac{(\log n)^{m+1}}{m+1} + \gamma_m$

For negative powers less than -1, the series will converge, and we get something called the "Riemann zeta function" in the power:

$\zeta(s) = \sum_{n=1}^\infty \frac{1}{n^s} = \prod_{p} \frac{1}{1 - p^{-s}} = \frac{1}{\Gamma(s)} \int_{x=0}^\infty \frac{x^{s-1}}{e^x - 1} dx$

where the p's are all the prime numbers. There are some interrelationships, like

$\sum_{n=0}^\infty \frac{1}{(2n+1)^s} = (1 - 2^{-s}) \zeta(s) \\ \sum_{n=1}^\infty \frac{(-1)^{n+1}}{n^s} = (1 - 2^{1-s}) \zeta(s) = \frac{1}{\Gamma(s)} \int_{x=0}^\infty \frac{x^{s-1}}{e^x + 1} dx$

4. The Riemann zeta function has the series expansion

$\zeta(s) = \frac{1}{s-1} + \sum_{n=0}^\infty (-1)^n \gamma_n \frac{(s-1)^n}{n!}$

It diverges for s = 1, as one might expect, but for s < 1, the series will still converge, despite its defining sum failing to converge. This is from something called "analytic continuation", something which is roughly going around poles or divergent spots like s = 1 here, going around them in the complex plane.

One finds a remarkable sort of reflection theorem:

$\zeta(s) = 2^s \pi^{s-1} \sin ((1/2) \pi s) \Gamma(1-s) \zeta(1-s)$

It gives values of the function for s < 1, values where the function's series diverges. Some special values, for positive integers n:

$\zeta(0) = \frac12 ,\ \zeta(1) = \infty ,\ \zeta(-2n) = 0 ,\ \zeta(1-2n) = - \frac{B_{2n}}{2n} ,\\ \zeta(2n) = \frac{(2\pi)^{2n}}{2 (2n)!} (-1)^{n+1} B_{2n}$

Another special value:
$\sum_{n=1}^\infty \frac{(-1)^{n+1}}{n} = \log 2$

5. It is evident that the Riemann zeta function has zeros -2, -4, -6, ..., its "trivial" ones. There is a theorem about its remaining zeros, its "nontrivial" ones. It is that these zeros all have real part 1/2 -- the This theorem has never been proved, despite a lot of effort and despite failure to find counterexamples. It has implications in, among other things, the distribution of prime numbers.

6. Thus, one gets

$\sum_{n=1}^\infty \frac{1}{n^2} = \frac{\pi^2}{6} \\ \sum_{n=1}^\infty \frac{1}{n^4} = \frac{\pi^4}{90}$

$\xi(s) = \frac12 \pi^{-s/2} s (s-1) \Gamma(s/2) \zeta(s)$

It satisfies $\xi(s) = \xi(1 - s)$

There are ways of constructing series expansions that converge over smaller values of the real part of the arg. The defining series converges only for Re(s) > 1, and this series converges for Re(s) > 0:

$\zeta(s) = \frac{1}{1 - 2^{1-s}} \sum_{n=1}^\infty \frac{(-1)^{n+1}}{n^s}$

and these "Dirichlet series", converging for Re(s) > 0 and Re(s) > -1:

$\zeta(s) = \frac{1}{s-1} \sum_{n=1}^\infty \left( \frac{n}{(n+1)^s} - \frac{n - s}{n^s} \right)$

$\zeta(s) = \frac{1}{s-1} \sum_{n=1}^\infty \frac{n(n+1)}{2} \left( \frac{2n+3+s}{(n+1)^{s+2}} - \frac{2n - 1 - s}{n^{s+2}} \right)$

8. But there is an interesting series for the reciprocal of the Riemann zeta function:

$\frac{1}{\zeta(s)} = \sum_{n=1}^\infty \frac{\mu(n)}{n^s}$

where μ(n) is the Möbius (Moebius) mu function. Defined for positive-integer n, it has values
• n = 1: 1
• n has maximum power 1 of all its primes (it's square-free): (-1)^(number of primes)
• n's factors have at least one prime power greater than 1: 0

Thus, mu(1) = 1, mu(2) = -1, mu(3) = -1, mu(4) = 0, mu(5) = -1, mu(6) = 1, ...

-

This function appears in an interesting context. An algebraic field is a combination of additionlike and multiplicationlike operations on some set where both are commutative (abelian) groups, the first over the entire set, and the second over all the set but the additive identity (0).

All the finite fields are known, and they were discovered by Evariste Galois. They all have numbers of elements that are powers of primes: GF(p^n) for prime p. The fields GF(p) are integers 0, 1, ..., p-1 with addition and multiplication both modulo p. The fields GF(p^n) for n > 1 are more complicated.

Their elements are polynomials in some variable with coefficients in GF(p). Thus, GF(4) has elements 0, 1, x, 1+x for variable x. Addition and multiplication are defined in the usual ways for polynomials, though polynomial multiplication involves taking the remainder from division by some degree-n irreducible monic polynomial. "Irreducible" meaning that it cannot be factored, and monic meaning that the highest term has coefficient 1. Let us see what degree-2 irreducible monic polynomials are available for GF(4).

x^2 = x*x (no), x^2 + 1 = (x+1)^2 (no), x^2 + x = x*(x+1) (no), and x^2 + x + 1 (yes)

For the second one, remember that the coefficients are in GF(2), and 1 + 1 = 0 in GF(2).

For larger fields, one has several such polynomials available, and how many there are is given by this formula:

$\frac1n \sum_{d | n} \mu(d) p^{n/d}$

for all d evenly dividing n.

9. For GF(8), one has
x^3 (no), x^3 + 1 = (x^2+x+1)*(x+1) (no), x^3 + x = x*(x+1)^2 (no), x^3 + x + 1 (yes), x^3 + x^2 = x^2*(x+1) (no), x^3 + x^2 + 1 (yes), x^3 + x^2 + x = x*(x^2+x+1) (no), x^3+x^2+x+1 = (x+1)^3 (no)

Two irreducible ones.

For GF(9), one has

x^2 (no), x^2 + 1 (yes), x^2 + 2 = (x+1)*(x+2) (no), x^2 + x = x*(x+1) (no), x^2 + x + 1 = (x+2)^2 (no), x^2 + x + 2 (yes), x^2 + 2x = x*(x+2) (no), x^2 + 2x + 1 = (x+1)^2 (no), x^2 + 2x + 2 (yes)

Three irreducible ones.

#### Posting Permissions

• You may not post new threads
• You may not post replies
• You may not post attachments
• You may not edit your posts
•