I am playing around with quotient rings in several variables and would like to get the congruence classes. How would I do that?
Let's say I have:
R.<x,y> = ZZ[]
S.<a,b> = R.quotient((x^2 + y^2, 17))
How do I get all congruence classes of R/[x^2 + y^2, 17]?