PROGRAM DYNAMIC
      IMPLICIT INTEGER(A-Z)
      DOUBLE PRECISION TIMEA,TIMEB,TIMTOT,A1(30,30),U(3500000),SUM1,
     1                 A2(30,30),B(30),W1,W2,Z1,Z2,SUM2,TEMP,TEMP2,SUM,
     1                 SL,SR,DEFAYS,RDOM,EPS
      INTEGER V(3500000),S(30),SCOMP(30),IPERM(30)
C
C #####################################################################
C                    DYNAMIC PROGRAMMING ALGORITHM
C  OBJECTIVE: MAXIMIZE THE SUM OF THE PROXIMITIES IN THE UPPER TRIANGLE
C             OF THE MATRIX - DOMINANCE INDEX.  THIS VERSION OF THE
C             PROGRAM PROMPTS FOR INPUT OF THE OPTIMAL SINGLE-OBJECTIVE
C             FUNCTION VALUES, AND THEN PRODUCES 101 SOLUTIONS BY
C             LETTING THE WEIGHTS FOR EACH CRITERIONSPAN FROM 0 TO 1
C             IN INCREMENTS OF .01
C #####################################################################
C
      OPEN(1,FILE='ASYM.DAT')
      OPEN(3,FILE='SYM.DAT')
      OPEN(2,FILE='ASYM.OUT')
      EPS = 1.0D-06
      READ(1,*) N                              ! Read number of objects
      READ(1,*) ((A1(I,J),J=1,N),I=1,N)        ! Read proximity matrix
      READ(3,*) N                              ! Read number of objects
      READ(3,*) ((A2(I,J),J=1,N),I=1,N)        ! Read proximity matrix
      DO I = 1,N
        DO J = 1,N
          IF(I.NE.J) B(I) = B(I) + A2(I,J)
        END DO
      END DO
      WRITE(*,*) 'INPUT UNIOBJECTIVE OPTIMAL OBJECTIVE VALUE'
      WRITE(*,*) 'FOR THE SIGN INDEX'
      READ(*,*) Z1
      WRITE(*,*) 'INPUT UNIOBJECTIVE OPTIMAL OBJECTIVE VALUE'
      WRITE(*,*) 'FOR THE ABSOLUTE MAGNITUDE INDEX'
      READ(*,*) Z2
      CALL GETTIM (IHR, IMIN, ISEC, I100)
      CALL GETDAT (IYR, IMON, IDAY)
      TIMEA=DFLOAT(86400*IDAY+3600*IHR+60*IMIN+ISEC)+DFLOAT(I100)/100.
      DO 4000 IJK = 0,100
        W1 = IJK/100.
        W2 = 1.-W1
        IF(W1.LT.EPS) W1 = EPS
        IF(W2.LT.EPS) W2 = EPS
C
      NSUM = 2**N - 1
      DO I = 1,NSUM
        U(I) = -9.9D+12
      END DO
      DO I = 1,N
        INDEX = 2**(I-1)
        SUM = 0.0D0
        DO J = 1,N
          IF(I.NE.J) SUM = SUM + A1(I,J)
        END DO
        V(INDEX) = I
        U(INDEX) = (SUM/Z1)*W1 + (B(I)**2/Z2)*W2
      END DO
C
      MTC = 0
      DO 70 K = 1,N-1
        NFIRST = 0
        IF(MTC.EQ.0) THEN
          M2 = 0
          DO J = 1,K
            S(J) = M2 + J
          END DO
          IF(S(1).NE.N-K+1) MTC = 1
        END IF
 
 75     IF(MTC.EQ.0) GO TO 70
        IF(NFIRST.EQ.1) THEN
          IF(M2.LT.N-NH) NH = 0
          NH = NH + 1
          M2 = S(K+1-NH)
          DO J = 1,NH
            S(K+J-NH) = M2 + J
          END DO
          MTC = 1
          IF(S(1).EQ.N-K+1) MTC = 0
        END IF
        IF(NFIRST.EQ.0) NFIRST = 1
        INDEX = 0
        DO I = 1,K
          INDEX = INDEX + 2**(S(I)-1)
        END DO 
C      
        JJ = 1
        DO I = 1,N
          SCOMP(I) = 0
        END DO
        IDX = 0
        DO 76 I = 1,N
          DO J = 1,K
            IF(I.EQ.S(J)) GO TO 76
          END DO
          IDX = IDX + 1
          SCOMP(IDX) = I
 76     CONTINUE
C
        NK = N - K
        DO JJ = 1,NK
          J1 = SCOMP(JJ)
          SUM1 = 0.0D0
          DO I = 1,NK
            I1 = SCOMP(I)
            IF(I1.NE.J1) SUM1 = SUM1 + A1(J1,I1)
          END DO
          SUM2 = 0.0D0
          DO I = 1,K
            I1 = S(I)
            SUM2 = SUM2 + A2(J1,I1)
          END DO
          TEMP2 = (B(J1)-2*SUM2)**2
          TEMP = W1*(SUM1/Z1) + W2*(TEMP2/Z2) + U(INDEX)
          IDXTMP = INDEX + 2**(J1-1)
          IF(TEMP.GT.U(IDXTMP)) THEN
            U(IDXTMP) = TEMP
            V(IDXTMP) = J1
          END IF
        END DO
        GO TO 75
  70  CONTINUE
C
      IPERM(N) = V(NSUM)
      INDEX = NSUM
      LASTIN = IPERM(N)
      DO I = 1,N-1
        INDEX = INDEX - (2**(LASTIN-1))
        LASTIN = V(INDEX)
        IPERM(N-I) = LASTIN
      END DO
C
      RDOM = 0.
      DO I = 1,N-1
        II = IPERM(I)
        DO J = I+1,N
          JJ = IPERM(J)
          RDOM = RDOM + A1(II,JJ)
        END DO
      END DO
C
      DEFAYS = 0.
      DO I = 1,N
        II = IPERM(I)
        SL = 0.
        DO J = 1,I-1
          JJ = IPERM(J)
          SL = SL + A2(II,JJ)
        END DO
        SR = 0.
        DO J = I+1,N
          JJ = IPERM(J)
          SR = SR + A2(II,JJ)
        END DO
        DEFAYS = DEFAYS + (SL-SR)**2
      END DO
      WRITE(*,750) W1,W2,RDOM,DEFAYS,(IPERM(I),I=1,N)
      WRITE(2,750) W1,W2,RDOM,DEFAYS,(IPERM(I),I=1,N)
 4000 CONTINUE
      CALL GETTIM (IHR, IMIN, ISEC, I100)
      CALL GETDAT (IYR, IMON, IDAY)
      TIMEB=DFLOAT(86400*IDAY+3600*IHR+60*IMIN+ISEC)+DFLOAT(I100)/100.
      TIMTOT = TIMEB - TIMEA
C      WRITE(*,80) U(NSUM)
      WRITE(*,81) TIMTOT
C      WRITE(*,82) (IPERM(I),I=1,N)
C      WRITE(2,80) U(NSUM)
      WRITE(2,81) TIMTOT
C      WRITE(2,82) (IPERM(I),I=1,N)
 80   FORMAT(' OPTIMAL DOMINANCE INDEX  ',F15.5)
 81   FORMAT(' TOTAL CPU TIME (SECONDS) ',7X,F8.2)
 82   FORMAT(30I3)
 750  FORMAT(2F5.2,F5.0,F12.5,22I3)
      END