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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
# ?+ v, [5 |* }9 W* ]5 R9 amodel:
sets:
# B+ |5 y) z9 n: u2 F* u, B  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 ;
  q% j; R1 B: s( ]7 Hendsets
( ?+ q2 w4 C) V/ Rdata:  !关联度矩阵,要对称;
! R9 Y2 l  H+ pT=
-1000 4 0
0 -1000 4
0 0 -1000
* g* w) n! Y% |! \;
  g/ w# j, X+ B, e' Vh=4;
" z2 {0 X1 c  ng=1;
a=1,2,1;
b=1,1,1;
/ a8 S( _" S, ienddata
% c" y9 u, w+ Y, n! 关联度最大;
7 ?, V7 U! d0 lmax=@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)))
" r, ]! p3 N- M! ^. ?7 z8 p+ @2 @5 }       #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));
1 |) s8 J, B# n! X& B6 L  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
" G  T& E+ ?$ A3 Q7 q  @for(region(i)： @for(region(j)|i #ne# j： @abs(y(i)-y(j))>=0.5*(b(i)+b(j))));
( ~# C, x+ R7 x$ ?3 m  @for(region(i):x(i)<h);
4 v, a' b' Z, e  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
, b3 `$ g5 F4 U- H7 D8 O4 P# Mend

zhangtao89524

loooog12

。。。。。