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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
1 h" a/ |" f8 {8 J5 Z7 B* b2 _! H% }9 wmodel:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
) I' P+ {5 x$ I. Ldata:  !关联度矩阵,要对称;
! b7 p8 i7 }+ _' w) r1 BT=
-1000 4 0
: E  o# g' D. y' p+ B6 C& h0 -1000 4
' K- Y. S, }# h# @* B0 0 -1000
& l# ?1 U+ @, S7 ^2 e4 T2 m4 F;
h=4;
g=1;
a=1,2,1;
b=1,1,1;
enddata
! 关联度最大;
+ |6 h8 h( J0 Smax=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
3 `. x% Z/ Q9 x! n# Z  @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)))
3 c' z1 g" |' _5 M& G, 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));
8 D8 |: x& q& n: L' _  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
7 V8 n* W$ s/ k' F9 b& b  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
: |# q( R( L+ j: ?" {  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
) m, L* {8 E/ ^4 U- `$ {end

zhangtao89524

loooog12

。。。。。