b(n)=my(x=factor(polcyclo(n,2)/gcd(polcyclo(n,2),n)));for(k=1,omega(polcyclo(n,2)/gcd(polcyclo(n,2),n))-1,print1("M( ",n," )C: ",x[k,1]);print());print1("M( ",n," )D");print()
c(n)=gcd(polcyclo(4*n,2),2^n-2^((2*n+2)/4)+1)
d(n)=gcd(polcyclo(4*n,2),2^n+2^((2*n+2)/4)+1)
e(n)=my(x=factor(c(n)/gcd(c(n),4*n)));for(k=1,omega(c(n)/gcd(c(n),4*n))-1,print1("M( ",4*n,"A )C: ",x[k,1]);print());print1("M( ",4*n,"A )D");print()
f(n)=my(x=factor(d(n)/gcd(d(n),4*n)));for(k=1,omega(d(n)/gcd(d(n),4*n))-1,print1("M( ",4*n,"B )C: ",x[k,1]);print());print1("M( ",4*n,"B )D");print()
for(k=1,1024,if(k%8==4,e(k/4);f(k/4),b(k)))