关于一道求解n-center的问题

heavenyhz

我刚刚开始接触linggo，下面代码是我写的一个求解n-center问题的代码。但是运行时，总说我end行错误，这是为什么呢？多谢指点！！

MODEL:
!2-center Problem
SETS:
   Marcket/ M1.. M12/:Wi, Hi, WDst, Dvy;
!Hi is the fixed operation cost at each marcket.
  Wi is Cij* #of load, where Cij is unit cost per distance;
   Center /C1.. C12/ ;
   LINKS (Marcket, Center): Dist, Xij, Sij,Connection;
   CloDist(Marcket,LINKS) : Dis21Cent    ;
ENDSETS
DATA:
    Hi = 200 210 220 300 310 260 150 180 220 305 160 270;
    Wi =    2   3    1    5    4   3    4   5   8    2   1    4;
  
   Dist =  0 15  7 54 45 28 18 37 44 38 27 50
          15  0  8 55 46 13 19 38 45 39 28 51
           7  8  0 47 38 21 11 30 37 31 20 43
          54 55 47  0  9 68 36 17 24 38 27 30
          45 46 38  9  0 59 27  8 15 29 18 21
          28 13 21 68 59  0 32 51 58 52 41 64
          18 19 11 36 27 32  0 19 26 20  9 33
          37 38 30 17  8 51 19  0  7 21 10 13
          44 45 37 24 15 58 26  7  0 28 17 20
          38 39 31 28 29 52 20 21 28  0 11 34
          27 28 20 27 18 41  9 10 17 11  0 23
          50 51 43 30 21 64 33 13 20 34 23  0;
   Connection = 0 0 1 0 0 0 0 0 0 0 0 0
                0 0 1 0 0 1 0 0 0 0 0 0
                1 1 0 0 0 0 1 0 0 0 0 0
                0 0 0 0 1 0 0 0 0 0 0 0
                0 0 0 1 0 0 0 1 0 0 0 0
                0 1 0 0 0 0 0 0 0 0 0 0
                0 0 1 0 0 0 0 0 0 0 1 0
                0 0 0 0 1 0 0 0 1 0 1 1
                0 0 0 0 0 0 0 1 0 0 0 0
                0 0 0 0 0 0 0 0 0 0 1 0
                0 0 0 0 0 0 1 1 0 1 0 0
                0 0 0 0 0 0 0 1 0 0 0 0;
         
ENDDATA     





!Sij<= Dist, Sij is the distance of the center from the lower vertix;
@FOR (LINKS (I,J):
           Sij(I,J) <= Connection (I,J)*Dist (I,J));

      !Xij is the arc that the center is located in. it can only happen in a  arc belonging to connection;
@for (LINKS (I,J): Xij(I,J) <= Connection(I,J));
!we need 2 center;
@FOR( LINKS (I,J):
      @SUM(LINKS(I,J): Xij )= 2);
@for (CloDist (K,I,J): Dis21Cent (K,I,J)=@if( Xij(I,J) #eq# 1,@smin(Dist(K,I)+Sij(I,J),Dist(K,J)+Dist(I,J)-Sij(I,J) ),99999999);
@for (Marcket (K): Dvy(K) = @min(CloDist(K,I,J): Dis21Cent(K,I,J)));
@for (Marcket(K): WDst(K)= Hi(K)+Wi(K)*Dvy(K));
MaxWDst = @Max (Marcket: WDst);
MaxWDst<=Z;
!Objective function;
MIN = Z;
!Xij,  is binery, it has value 1 only when a center is located between  i and j;
@FOR (LINKS: @bin (Xij));
END