Euclid algorhitm

M=[[1,0],[0,1]]
a,b=2261,1275
A,B=a,b
def stR(x):
    if x<0: return '('+str(x)+')'
    else: return str(x)
while a<>0:
    factor=b/a
    a,b=b-a*factor,a
    M=[[M[1][0]-factor*M[0][0],M[1][1]-factor*M[0][1]],[M[0][0],M[0][1]]]
    print '(a,b)=',(a,b)
    print M
    print 'answer:', (M[1][0],M[1][1])
    print stR(A)+'*'+stR(M[1][0])+'+'+stR(B)+'*'+stR(M[1][1])+'=',A*M[1][0]+B*M[1][1]