Function solve_linear_diophantine 

pub fn solve_linear_diophantine<T>(
    a: T,
    b: T,
    c: T,
) -> Option<LinearDiophantineSolution<T>>
where
    T: Signed,

Solve ax + by = c