lingo程序问题，算不出可行解，请教各位高手达人。谢谢！

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
model:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
! k% s: ?- k; }, F0 m' V    轴上的距离、b表示region在y轴上的距离;            
4 w% g$ d1 X0 j8 e  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
) {+ O; s" \8 _data:  !关联度矩阵,要对称;
T=
-1000 4 0
) _* j& n8 H# m9 `, I8 S2 A0 -1000 4
# |, x  B! @7 t4 ^& b0 0 -1000
1 j5 V+ ], Q* j;
% E3 D# N& C+ l; y6 g- O+ i  P1 \9 Ch=4;
4 W( y5 I. F$ f$ [g=1;
a=1,2,1;
" E- U' N3 x. v: o- h0 Mb=1,1,1;
( `) j4 N1 ?# ?  V. f" w  [: wenddata
! 关联度最大;
/ I& R+ ?, Q- B7 R, y0 Wmax=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
  @for(link(i,j)|i #ne# j:z(i,j)=@if((@abs(x(i)-x(j)) #eq# 0.5*(a(i)+a(j)) #and# @abs(y(i)-y(j))#lt# 0.5*(b(i)+b(j)))
8 m( K: y4 L9 @0 o: q       #or#(@abs(y(i)-y(j)) #eq# 0.5*(b(i)+b(j)) #and# @abs(x(i)-x(j)) #lt# 0.5*(a(i)+a(j))),1,0));
+ P+ q! K4 \; G1 X9 _( u  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
& j' g9 d- W  H/ J5 I2 D  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
( i) L0 i2 H- ^/ {  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
end

zhangtao89524

loooog12

。。。。。