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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
model:
, Q0 U% M% R* T$ j; ~9 V4 d1 `+ V9 ssets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
. S' K; L" g* S; F3 H    轴上的距离、b表示region在y轴上的距离;            
# w  G1 u- _& [5 W% q; t! v  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
8 o. v3 l: @+ [T=
-1000 4 0
0 -1000 4
0 0 -1000
;
0 w( \! T6 `5 Q) Y9 qh=4;
- Q* M) |0 L% H- i, K& M+ Yg=1;
a=1,2,1;
7 D& C- R8 k5 f( P: T7 Rb=1,1,1;
4 t' P# E& Q+ c5 P% d% e- N0 O1 genddata
% ~/ m$ u5 o/ ?' U2 y6 }0 {! 关联度最大;
4 z  E* K7 Y# }0 \" S- P, h3 ymax=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
2 j9 s9 |* a, |7 P  @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)))
6 e- ?! }( N* u; s9 s       #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));
  @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)));  
. k+ g$ Q6 A. O! d0 v+ oend

zhangtao89524

loooog12

。。。。。