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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
/ V6 V* N" m; Kmodel:
sets:
  region/1..3/:x,y,a,b; !x,y分别表示质心的横坐标、纵坐标,a表示region在x
3 }2 I" R2 ^8 v9 K: I# v    轴上的距离、b表示region在y轴上的距离;            
  link(region,region):T,z;  ! z(i,j)=1表示i和j紧邻，否则=0 ;
, p6 m% @* ^$ Q3 eendsets
- U8 M" X, w, U, ^$ Rdata:  !关联度矩阵,要对称;
" q1 k- I2 `: P2 ^T=
. h" C6 R; _) U7 S) z% a+ d3 l-1000 4 0
' V5 k/ w8 q) b* l: _( p0 -1000 4
0 0 -1000
; P6 N  q4 D0 _0 \# H" v;
" Z3 F: V$ L8 X+ r+ p4 @0 a) V6 [h=4;
g=1;
4 P! [; J, l& z* Y% o: o' pa=1,2,1;
b=1,1,1;
enddata
2 ?6 g- W4 T7 W$ `$ n6 N$ b! 关联度最大;
max=@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)))
       #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))));
. Z* Q$ A* e6 y  @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);
& c( I5 t2 L6 R8 K9 `  @for(region(i):y(i)<g);
- V, _, f+ c: w3 {5 Q) u  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
# p9 g* \0 `) V2 O8 kend

zhangtao89524

loooog12

。。。。。