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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
) h6 g( g" s  Vmodel:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
5 M0 i# l9 l& g$ J. g    轴上的距离、b表示region在y轴上的距离;            
. y) _! j2 m$ k+ y0 o/ l  ?( K) f  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
" u/ ^  o6 n- z: J/ Ldata:  !关联度矩阵,要对称;
T=
% s; M# M! Y$ d-1000 4 0
0 -1000 4
0 0 -1000
;
h=4;
, E( f9 |# A% r  D* |; T0 y8 y+ Ng=1;
* ?( n) Y" L7 T: k# K6 za=1,2,1;
4 N* |4 t6 u0 F) m: {* O" L+ e/ S" u4 Pb=1,1,1;
enddata
! 关联度最大;
% j# a& j& k  o2 |. omax=@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)))
7 l' B9 }+ a9 W+ Q- h& 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))));
; w4 }' ^9 W  a  q: I# ?7 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)));  
6 i# L/ S2 f1 N% ~5 [  m4 fend

zhangtao89524

loooog12

。。。。。