def equation_diophantienne(a,b,c):
    p=pgcd(a,b)
    x0=c*bezout(a,b)[0]/p
    y0=c*bezout(a,b)[1]/p
    return [(x0,y0),(b/p,"k +",x0,";"-a/p,"k +",y0)]