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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
4 C( w. @5 Z1 }model:
sets:
! T; J/ y1 f9 b# Z  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
    轴上的距离、b表示region在y轴上的距离;            
9 d6 \  m- O; `+ ^' Q! ^  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
+ w3 p$ A7 p# ]: _- n3 UT=
-1000 4 0
0 -1000 4
0 0 -1000
* P. Z) U9 b# L6 |' w" N7 r  X;
h=4;
g=1;
$ ]4 i, n! Q8 Za=1,2,1;
* s. M" S& G4 z/ i0 ub=1,1,1;
enddata
- q1 Q9 g$ z, }: S! 关联度最大;
- Q$ q9 S9 r) T& r5 Q+ _max=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
* D8 T& O. v/ C& g) c  @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)))
5 w) n. @8 E: C* Y' B       #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));
! H- T" {$ \6 q( I# b5 g: Y  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
, Z7 W' ~5 W3 F& y- O( G  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
2 P' N5 h4 P- n  @for(region(i):x(i)<h);
( E& f7 u' x% p5 E9 G9 @2 P  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
end

zhangtao89524

loooog12

。。。。。