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

hitman

运行几个小时，状态还是未知（unknown）,程序如下,模型附在最后，请大家帮我看看，谢谢！
0 t% w. f: c$ f# S! Bmodel:
  {: ~& j) q8 P: R9 ksets:
  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
* x# ~0 G4 c' }+ R! X2 A$ h- adata:  !关联度矩阵,要对称;
T=
0 k+ `2 U* n( }2 p2 R4 k) K-1000 4 0
3 ]; C; s) P8 p, f2 w0 -1000 4
0 0 -1000
/ g& S- W1 e% _9 B& j+ V  m;
h=4;
6 F2 o2 e/ D' ]; ?g=1;
, Z  R3 g/ y: s/ [a=1,2,1;
2 H: {4 ?0 U) L+ [# eb=1,1,1;
  b( {; F/ i# L% f9 e9 Eenddata
0 H/ e! y' d4 _! x/ Z  @! 关联度最大;
max=@sum(link(i,j)|i #ne# j: T(i,j)*z(i,j));
# S6 E! Y0 {- u. w( G4 C  @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))));
% Q+ v9 M4 M  ?7 f" b6 Z  @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);
  P, M$ A3 t5 n6 T0 l) [; o  @for(region(i):y(i)<g);
  @for(link(i,j)|i #ne# j： @bin(z(i,j)));  
  O. a7 E! |8 a: ?. L$ iend

zhangtao89524

loooog12

。。。。。