problem with integer infeasible

Problems with modeling
Post Reply
ozge
User
User
Posts: 6
Joined: 2 months ago

problem with integer infeasible

Post by ozge » 2 months ago

Hi everyone, why do I get no solution? Can you help me again?

set i /0*7/;
alias (i,ip);
display i,ip;

set subset(i) /1*6/;
alias (subset,subsetp);
display subset, subsetp;

table d(i,ip) 'distances between cities'
$onDelim
,1,2,3,4,5,6
1,0.00,1.66,5.08,6.52,8.83,5.53
2,1.66,0.00,4.09,6.02,9.20,5.76
3,5.08,4.09,0.00,2.45,6.96,4.00
4,6.52,6.02,2.45,0.00,4.80,2.71
5,8.83,9.20,6.96,4.80,0.00,3.44
6,5.53,5.76,4.00,2.71,3.44,0.00



$offDelim
;

display d;
Scalar m 'total number of cities' / 6 /;

variables x(i,ip) 'if ip is visited after i 0-15 arası'
z 'total distance travelled' ;

binary variables x;

equations of
eq1
eq2
eq3(ip) 'city balance constraint'
eq4(subset)
eq5(subset);

of.. z=e= sum((subset,subsetp), x(subset,subsetp)*d(subset,subsetp));
eq1.. sum(subsetp, x('0',subsetp))=e= 1;
eq2.. sum(subset, x(subset,'7'))=e=1;
eq3(subsetp).. (sum( subset $(not sameas(subset,subsetp)), x(subset,subsetp)) - sum(subset $(not sameas(subset,subsetp)), x(subsetp,subset)) ) =e= 0;
eq4(subset).. sum(subsetp $(not sameas(subset,subsetp)), x(subset,subsetp) ) =g= m*x('0', subset);
eq5(subset).. sum(subsetp $(not sameas(subset,subsetp)), x(subsetp,subset) )=g= m*x(subset,'7');

model tsp2test / all /;
solve tsp2test minimizing z using mip;

display x.l, x.m, z.l;
Attachments
Ekran Resmi 2018-10-02 15.08.30.png

Post Reply