def dmp_content(f, u, K):
 (cont, v) = (dmp_LC(f, K), (u - 1))
 if dmp_zero_p(f, u):
  return cont
 for c in f[1:]:
  cont = dmp_gcd(cont, c, v, K)
  if dmp_one_p(cont, v, K):
   break
 if K.is_negative(dmp_ground_LC(cont, v, K)):
  return dmp_neg(cont, v, K)
 else:
  return cont

def dmp_content ( f , u , K ) : ( cont , v ) = ( dmp_LC ( f , K ) , ( u - 1 ) ) if dmp_zero_p ( f , u ) : return cont for c in f [ 1 : ] : cont = dmp_gcd ( cont , c , v , K ) if dmp_one_p ( cont , v , K ) : break if K . is_negative ( dmp_ground_LC ( cont , v , K ) ) : return dmp_neg ( cont , v , K ) else : return cont
returns gcd of multivariate coefficients .