code: def dmp_ground_content(f, u, K):
      """
      Compute the GCD of coefficients of ``f`` in ``K[X]``.
  
      Examples
      ========
  
      >>> from sympy.polys import ring, ZZ, QQ
  
      >>> R, x,y = ring("x,y", ZZ)
      >>> f = 2*x*y + 6*x + 4*y + 12
  
      >>> R.dmp_ground_content(f)
      2
  
      >>> R, x,y = ring("x,y", QQ)
      >>> f = 2*x*y + 6*x + 4*y + 12
  
      >>> R.dmp_ground_content(f)
      2
  
      """
      from sympy.polys.domains import QQ
  
      if not u:
          return dup_content(f, K)
  
      if dmp_zero_p(f, u):
          return K.zero
  
      cont, v = K.zero, u - 1
  
      if K == QQ:
          for c in f:
              cont = K.gcd(cont, dmp_ground_content(c, v, K))
      else:
          for c in f:
              cont = K.gcd(cont, dmp_ground_content(c, v, K))
  
              if K.is_one(cont):
                  break
  
      return cont
  
imports: from sympy.polys.densearith import dup_add_term, dmp_add_term, dup_lshift, dup_add, dmp_add, dup_sub, dmp_sub, dup_mul, dmp_mul, dup_sqr, dup_div, dup_rem, dmp_rem, dmp_expand, dup_mul_ground, dmp_mul_ground, dup_quo_ground, dmp_quo_ground, dup_exquo_ground, dmp_exquo_ground
  from sympy.polys.densebasic import dup_strip, dmp_strip, dup_convert, dmp_convert, dup_degree, dmp_degree, dmp_to_dict, dmp_from_dict, dup_LC, dmp_LC, dmp_ground_LC, dup_TC, dmp_TC, dmp_zero, dmp_ground, dmp_zero_p, dup_to_raw_dict, dup_from_raw_dict, dmp_zeros
  from sympy.polys.polyerrors import MultivariatePolynomialError, DomainError
  from sympy.utilities import variations
  from math import ceil, log
label_desc: Use image node_id 23 for calling a global function with example usage: dmp_ground_content(f, u, K) and returns: cont, dup_content, K