code: def test_issue_25793():
      R, x = ring("x", ZZ)
      f = x - 4851  # failure starts for values more than 4850
      g = f * (2 * x + 1)
      H, cff, cfg = R.dup_zz_heu_gcd(f, g)
      assert H == f
  
imports: from sympy.polys.rings import ring
  from sympy.polys.domains import ZZ
  from sympy.polys.heuristicgcd import heugcd
label_desc: Use image node_id 4 for calling a global function with example usage: test_issue_25793() without return types