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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
, w% J" q1 |! W3 C7 X& `0 H! ?5 qmodel:
sets:
. q. |( s- o" E! Y% R2 z+ b1 R* K  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
' b" Y1 E) g9 N1 z% i0 D, z    轴上的距离、b表示region在y轴上的距离;            
3 Q; {# c8 I' T2 `  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
T=
-1000 4 0
0 -1000 4
  S3 {  i& {/ N/ |+ I7 m/ G) O0 0 -1000
;
h=4;
g=1;
( u% i- J4 Z) N- G% O; P+ va=1,2,1;
  m* w. M( _+ [- [$ n% O  H- Eb=1,1,1;
+ u: R, a& [% a/ y: T. D* renddata
5 Y5 S( L5 X6 Q9 k! 关联度最大;
' P$ w# F' R- Z! v# \& `" `max=@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)))
       #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));
, p+ X7 k" K# ~  m. T  @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))));
  @for(region(i):x(i)<h);
4 C2 N" `# G* Y+ G  @for(region(i):y(i)<g);
2 v0 D1 W3 o/ J4 D, W  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
end

zhangtao89524

loooog12

。。。。。