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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
model:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
+ U* g* `4 L/ x    轴上的距离、b表示region在y轴上的距离;            
9 u2 C9 s" p7 l5 s& Z7 c  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
; m# u& C* K7 k9 j& P2 w3 r6 f: MT=
-1000 4 0
0 -1000 4
0 0 -1000
;
h=4;
g=1;
' g2 p5 @0 r. va=1,2,1;
& t2 ?5 E$ |2 Qb=1,1,1;
; v' j# v7 Z+ V1 D9 Jenddata
! 关联度最大;
& n  q5 x; D; z6 Z1 ?3 g. F; imax=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
- E, g- {' \  D0 e2 t  @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)))
       #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));
' }( T9 N5 `6 Y' |! P) Y6 ?, K  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
  s6 p& I/ q$ v* M* p  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
9 w) |# z9 V3 q( O9 a$ T: l  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
( s- p8 Q9 a- O. z" N4 r& Hend

loooog12

。。。。。

zhangtao89524