levels.pi

import cp.

main =>
     % matrix of values X*Y=Z
     T={
        {2,50,100},
        {1,6,6},
        {1,7,7},
        {3,1,3}
     },
     [X1,X2] :: [T[I,1] : I in 1..len(T)],
     %
     X1 #!= X2,
     Error1 #= sum([abs(X1*T[I,2]-T[I,1]*T[I,2]) : I in 1..len(T)]),
     Error2 #= sum([abs(X2*T[I,2]-T[I,1]*T[I,2]) : I in 1..len(T)]),
     
     solve($[min(Error1+Error2)],[X1,X2]),
     writeln([X1,X2]).