cast roots back to original ring

Author: c0rydoras

src/sage/rings/finite_rings/integer_mod_ring.py

@@ -1944,7 +1944,12 @@ def _roots_univariate_polynomial(self, f, ring=None, multiplicities=True, algori
 
         # Finite fields are a base case
         if self.is_field():
-            return f.change_ring(f.base_ring().field()).roots(multiplicities=False)
+            return list(
+                map(
+                    f.base_ring(),
+                    f.change_ring(f.base_ring().field()).roots(multiplicities=False),
+                )
+            )
 
         # Otherwise, find roots modulo each prime power
         fac = self.factored_order()