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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
7 B$ `  n# e! C; ~model:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
    轴上的距离、b表示region在y轴上的距离;            
8 m+ I" H6 p, Z  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
4 |5 D1 ~1 ~4 q: [" t$ g/ |T=
) x- S7 j; E$ u( |-1000 4 0
4 Q* t* n, c$ Q% S, L' m0 -1000 4
0 0 -1000
: J- \  J7 h, J0 V) y% j0 |; q+ M8 M;
4 c3 j5 W9 l: \h=4;
# l4 `, s+ Q( F# o9 X* ]% R' Z) xg=1;
! O* {( @" H3 ja=1,2,1;
b=1,1,1;
' a5 B! y2 x: L: n% _enddata
! 关联度最大;
max=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
4 B5 T2 W' k- B+ N, v7 l% h; s$ |  @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));
( D7 J! Y9 x2 \+ ^% a" U/ o3 l. C  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
& A4 Y, m8 Y7 h8 J  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
9 n/ }$ q# G. L; S. L, x  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
end

zhangtao89524

loooog12

。。。。。