use sieved power sum in zeta for short power series; some refactoring