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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
model:
9 z6 c8 s5 v$ R, L7 q2 gsets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
    轴上的距离、b表示region在y轴上的距离;            
; Q& @: ?8 I/ b, k/ L  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
/ o$ e- G( p# L% d6 G5 LT=
-1000 4 0
2 E' n- N- g3 T* ]0 -1000 4
0 0 -1000
;
1 z; S7 _' r: k0 F+ t9 b2 }5 \h=4;
8 ]" e; p7 h, j; I7 i& K. vg=1;
a=1,2,1;
, ]6 h; i! R6 Z; Y  A+ i$ Ab=1,1,1;
6 F; R) N5 W( C4 A! aenddata
  A: [5 v; ]/ ~* k7 K4 Z6 i" w! 关联度最大;
5 ^9 X0 ^& z2 q* j+ m5 A( Emax=@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));
  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
' f/ I$ t9 r/ g# T  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
* y( e! J. O1 x! p# O  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
9 ]! M: q- l+ m2 ]  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
& `- I7 t1 }) s8 Mend

zhangtao89524

loooog12

。。。。。