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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
  U- `2 _/ {6 l& O# l2 S" Y$ ymodel:
sets:
  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 ;
endsets
data:  !关联度矩阵,要对称;
# D/ A/ c; `1 V1 ^7 \5 a5 w) }  jT=
4 q' X2 r$ A; }/ `-1000 4 0
  F3 e$ s# d! v) d. A7 Z0 -1000 4
0 0 -1000
/ N" w/ ~( U. \3 P" o. \3 x;
h=4;
, N3 W! U9 x5 y8 u# a6 Jg=1;
a=1,2,1;
b=1,1,1;
enddata
6 h1 B/ T2 S# [. L" x) @4 z! 关联度最大;
  A% L- C- L: H  Q% n1 E7 m3 ?max=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
$ |2 q5 p9 G1 N  @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)))
) f5 L# r- r* M& n6 D       #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));
+ N  g6 [7 \: E2 m: D  @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))));
  @for(region(i):x(i)<h);
  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
  e' S4 x3 I3 I- p$ C. Mend

zhangtao89524

loooog12

。。。。。