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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
: P& L* F5 ]( S" G  z. Mmodel:
; U) C* k+ B6 K+ w2 Tsets:
: v  u: v5 r' h$ R+ ?9 I3 t  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
/ v4 I& Y* C* W. H- _    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
6 H1 h4 z! Y8 ?1 {2 Q. F. f' Vendsets
4 v- ?: \. Q! G: K' F: l: K4 sdata:  !关联度矩阵,要对称;
6 R2 [/ e$ W0 q% [0 v1 @T=
* ]$ a+ D$ p0 Q( ?% I6 F-1000 4 0
0 -1000 4
4 R$ ]0 u5 Z8 ]& o9 t9 g& n0 0 -1000
;
" H; d0 |* P# r7 |- bh=4;
- b& D  h" \! }8 D3 \# ~0 _+ ug=1;
a=1,2,1;
b=1,1,1;
enddata
! 关联度最大;
max=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
5 _. o. y% K1 U4 K* l: s* w  @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));
' G) k: ^2 d* w  t& w" N  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
$ _5 v2 Q( f. B9 Y% e5 F5 V  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
; o7 T, d8 v9 X6 m0 J( s: g  @for(region(i):x(i)<h);
+ P& G7 X: G( M" {# I+ H" H# Z1 @  @for(region(i):y(i)<g);
  E9 ]) B, y4 \5 Z% u7 q  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
/ I: K$ Z' ^, g0 k8 Q6 U; D7 Eend

zhangtao89524

loooog12

。。。。。