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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
% L/ n/ X! {, Q+ @+ _( o: D9 fmodel:
sets:
. `: Q/ U& O( {- z, w  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
$ N8 h! t$ a% r9 x7 Y4 S: @9 I    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
endsets
data:  !关联度矩阵,要对称;
$ A- B; o0 q2 h: eT=
8 ^& h7 ~. ~' ~1 y7 j-1000 4 0
0 -1000 4
% U2 K7 i3 n% ?8 ~0 G0 0 -1000
% E8 E$ A. m3 Y: F2 |7 g/ [;
/ y4 J+ Y9 }- I+ t, Jh=4;
# ^) f) `# E# Z( R5 r- o9 Pg=1;
a=1,2,1;
. y" P  ]6 a. g8 w( }" ~7 _b=1,1,1;
; `0 e1 v6 J% t& G0 C) u% s' V3 Fenddata
, {' r% e/ }2 y$ }! e+ S1 f( i! 关联度最大;
) [5 G' Y8 U7 C2 ^5 ^  lmax=@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 E. A" H* P/ O       #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)));  
end

zhangtao89524

loooog12

。。。。。