|  | 
 Explicit scriind un program mic:
 sage: N = 7220
 sage: N.divisors()
 [1, 2, 4, 5, 10, 19, 20, 38, 76, 95, 190, 361, 380, 722, 1444, 1805, 3610, 7220]
 
 sage:
 for a in N.divisors():
 b = ZZ(N/a)
 if gcd(a,b) == ZZ(19):    print "(a,b) = (%d,%d)" % (a,b)
 
 (a,b) = (19,380)
 (a,b) = (76,95)
 (a,b) = (95,76)
 (a,b) = (380,19)
 
 
 ---
 df (gauss)
 
 |