code: def _gf_div(f, g, p):
      r"""
      Compute `\frac f g` modulo `p` for two univariate polynomials over
      `\mathbb Z_p`.
      """
      ring = f.ring
      densequo, denserem = gf_div(
          f.to_dense(), g.to_dense(), p, ring.domain
      )
      return ring.from_dense(densequo), ring.from_dense(denserem)
  
imports: from sympy.core.symbol import Dummy
  from sympy.ntheory import nextprime
  from sympy.ntheory.modular import crt
  from sympy.polys.domains import PolynomialRing
  from sympy.polys.galoistools import gf_gcd, gf_from_dict, gf_gcdex, gf_div, gf_lcm
  from sympy.polys.polyerrors import ModularGCDFailed
  from mpmath import sqrt
  import random
label_desc: Use image node_id 16 for calling a global function with example usage: _gf_div(f, g, p) and returns: ring,  ring