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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
3 g& v( e9 j; V' ]: q) o) b! Smodel:
5 \' G4 x. v* ^7 ?) z( R  lsets:
( [: C- P5 v" W  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
. B& u1 t, f* f    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
7 I6 t1 u. R9 x. C" tdata:  !关联度矩阵,要对称;
7 f6 R5 }# h1 d; X9 V9 iT=
2 V- `% Q, k! m. K-1000 4 0
# f  {$ `0 R2 Q+ U' Q0 c7 Q2 ^0 -1000 4
0 0 -1000
9 C, P1 k& P* N( |;
# h1 j4 e. o& W8 g# qh=4;
* w) v# S) e% w# ]% y* C9 q6 jg=1;
a=1,2,1;
) a/ M% D8 m" ~6 y! Qb=1,1,1;
$ X/ G2 w( d4 ?& Cenddata
: H1 B. F  n2 j! 关联度最大;
$ V" B) s6 O1 {6 y2 Zmax=@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)))
4 e' @3 f. o* z: I+ z3 m9 V5 F" F       #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 c9 Z; Z. w$ R, P4 C  @for(region(i)： @for(region(j)|i #ne# j： @abs(x(i)-x(j))>=0.5*(a(i)+a(j))));
/ W4 j$ t5 F. y7 N& R  @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);
  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
1 u7 y& S4 d7 }) o7 send

zhangtao89524

loooog12

。。。。。