Cyclotomic Cosets calculation

The following code calculates the cyclotomic cosets that are need for
the calculation of idempotents modulo 1+x^n of some degree