def extendEuclidean(a, b, s1=1, s2=0, t1=0, t2=1): if b: r=a%b return extendEuclidean(b, r, s2, s1-s2*(a//b), t2, t1-t2*(a//b)) return a, s1, t1