def test_ideal_log_relation():
K = QuadraticField(-41)
q = (K * 3).factor()[0][0]
C_K = K.class_group()
t = ideal_log_relation(q)
e = q.ideal_class_log()
assert q == t * prod(gi**ei.sage() for gi, ei in zip(C_K.gens_ideals(), e))
