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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
% L. n4 J# K5 [  \7 A& z( _1 Pmodel:
* |: w* y0 C2 f: Csets:
: o" I$ O4 Q  F% r: Q  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
5 ^, q: S0 X, t    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
2 F# H9 [1 E1 F4 HT=
3 I* s0 E2 J6 L  {* z-1000 4 0
% @& Y& @& |% x" w0 -1000 4
0 0 -1000
0 K  d8 s. t! k+ `$ |/ w;
h=4;
0 Q* u5 o) c* h0 ?: Ag=1;
a=1,2,1;
5 a- u8 ~: ^8 F/ A$ tb=1,1,1;
enddata
  A3 ~5 v& Y4 [# X9 r: q2 p! 关联度最大;
% e! ?8 q. d  w$ [1 A- P) xmax=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
8 _% A9 B( _7 C* H& \6 j. M  Q  @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));
# Z, Q3 r  b* t: \8 |  @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))));
# M+ q# T! k$ N* l# e+ O: m9 T  @for(region(i):x(i)<h);
7 K5 `" u6 T* q" K8 ?+ n1 v8 i  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
# e: t; P8 ?8 I; P0 p1 Mend

zhangtao89524

loooog12

。。。。。