code: def dmp_ground_extract(f, g, u, K):
      """
      Extract common content from a pair of polynomials in ``K[X]``.
  
      Examples
      ========
  
      >>> from sympy.polys import ring, ZZ
      >>> R, x,y = ring("x,y", ZZ)
  
      >>> R.dmp_ground_extract(6*x*y + 12*x + 18, 4*x*y + 8*x + 12)
      (2, 3*x*y + 6*x + 9, 2*x*y + 4*x + 6)
  
      """
      fc = dmp_ground_content(f, u, K)
      gc = dmp_ground_content(g, u, K)
  
      gcd = K.gcd(fc, gc)
  
      if not K.is_one(gcd):
          f = dmp_quo_ground(f, gcd, u, K)
          g = dmp_quo_ground(g, gcd, u, K)
  
      return gcd, f, g
  
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 27 for calling a global function with example usage: dmp_ground_extract(f, g, u, K) and returns: gcd,  f,  g