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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
+ ]& Y& q- ^2 H+ {  N& U+ cmodel:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
8 V& w+ {" m5 {    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
1 J7 @5 l, m2 H) hdata:  !关联度矩阵,要对称;
T=
/ i( G( k4 E, x# y9 o' I- o" |-1000 4 0
3 j, y; }7 |6 s7 r5 E# a: b0 -1000 4
0 0 -1000
9 Q5 B$ j" Q+ C;
h=4;
g=1;
a=1,2,1;
5 V9 p4 |- Q' S; Q" u; `' m  [b=1,1,1;
enddata
! 关联度最大;
/ {6 Y- v8 B( e: u/ dmax=@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 Y* X4 C% h       #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));
3 d+ Z3 A( o8 m" s5 e  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
1 `7 V. ^' k% I9 Z! K2 e  @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);
. M; j: _% `/ G6 l/ f& ?1 b. Q  @for(region(i):y(i)<g);
: T8 M; B3 V# J: ~+ F  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
, p" Y% x0 _+ @3 H* g4 U3 [  kend

loooog12

。。。。。

zhangtao89524