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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
4 I( G3 g% s- r8 E, ~: R/ v6 amodel:
, N; V: F1 m% p) ], K. {sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
. O; ~# h$ F. T8 E8 X6 B7 A1 x    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
0 E) z6 U$ \7 R1 dendsets
data:  !关联度矩阵,要对称;
9 v. _% R  G+ U- \& y$ rT=
# F6 B; B/ T3 c3 _-1000 4 0
0 -1000 4
$ T, Z' V8 ~# z  T2 x0 W' i0 0 -1000
;
h=4;
* Y* U- M" o: ^: _* H* G7 Z* }g=1;
; m( D8 Q: C; w5 xa=1,2,1;
b=1,1,1;
enddata
! 关联度最大;
max=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
+ x7 h) w# {+ q2 b  @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));
' ]/ F7 `$ s' T' F% Z5 ^" B5 `  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
5 D; @; G+ C+ U  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
/ G, }# R5 g* @6 X4 l  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
$ s* j* E5 e4 }/ M3 lend

zhangtao89524

loooog12

。。。。。