divisor function sigma

Calculate the divisor function sigma_k(n) for positive integer n

sigma(n,0) is equal to divisor_count(n)
sigma(n,1) is equal to sum(divisors(n))
sigma(n,k) is equal to sum([x**k for x in divisors(n)])

Default for k is 1.

References
==========

- http://en.wikipedia.org/wiki/Divisor_function
While the other two lines are true, this formula is true in every case (they reduce to this when `k=1` or `k=0`).  It might be clearer to just write this formula, and then note in a paragraph below that it reduces to the other things for these special values of `k`.