@vv thanks for the answer and the code. Actually I should have mentioned that I'm looking for a nonnegative integer solution. So I tried to modify your code by adding assume=nonnegint into minimize, but this function does not accept the option.