Move totient range to util.c, speed it up a little bit
Fix divisor sum on bigints
New data structure for phi(x,a): faster prime_count