sigma divisor function

use closed form for geometric partial sum