A Heuristic Nonlinear Integer Programming Solver Applied to an Eight-Department/Facility Quadratic Assignment Problem

Jsun Yui Wong

The computer program listed below tries to solve the 8-department problem of Nugent, Vollmann, and Ruml (1968).

0 DEFINT A-Z
3 DEFSNG H,S,P,M
5 DIM X(66),A(66)
8 DIM HS(33,33),HR(22,22)
11 FOR IZ=1 TO 8
14 FOR JZ=1 TO 8
16 READ HR(IZ,JZ)
18 NEXT JZ
19 NEXT IZ
24 DATA 999,1,2,3,1,2,3,4
25 DATA 1,999,1,2,2,1,2,3
26 DATA 2,1,999,1,3,2,1,2
27 DATA 3,2,1,999,4,3,2,1
28 DATA 1,2,3,4,999,1,2,3
29 DATA 2,1,2,3,1,999,1,2
30 DATA 3,2,1,2,2,1,999,1
31 DATA 4,3,2,1,3,2,1,999
49 FOR IZ=1 TO 7
50 FOR JZ=IZ+1 TO 8
51 READ HS(IZ,JZ)
52 NEXT JZ
53 NEXT IZ
54 DATA 5,2,4,1,0,0,6
55 DATA 3,0,2,2,2,0
56 DATA 0,0,0,0,5
57 DATA 5,2,2,10
58 DATA 10,0,0
59 DATA 5,1
60 DATA 10
74 FOR JJJJ=-32000 TO 32000
75 RANDOMIZE JJJJ
76 M=-999999999999#
81 FOR ZI=1 TO 8
82 A(ZI)=ZI
83 NEXT ZI
126 IMAR=10+RND*50
128 FOR I=1 TO IMAR
129 FOR KK=1 TO 8
131 X(KK)=A(KK)
132 NEXT KK
223 IJL=1+FIX(RND*8)
224 IJM=1+FIX(RND*8)
236 X(IJL)=A(IJM):X(IJM)=A(IJL)
1631 SUMUL=0
1632 FOR IS1=1 TO 7
1634 FOR IS2=IS1+1 TO 8
1637 SUMUL=SUMUL+HS(IS1,IS2)*HR(X(IS1),X(IS2))
1638 NEXT IS2
1639 NEXT IS1
1646 P=-SUMUL
1656 IF P<=M THEN 1670
1657 FOR KEW=1 TO 8
1658 A(KEW)=X(KEW)
1659 NEXT KEW
1661 M=P
1666 GOTO 128
1670 NEXT I
1890 IF M>-110 THEN 1912 ELSE 1999
1912 PRINT A(1),A(2),A(3),A(4),A(5)
1929 PRINT A(6),A(7),A(8),M,JJJJ
1999 NEXT JJJJ

This BASIC computer program was run with Microsoft's GW BASIC 3.11 interpreter, which is not a compiler. Its output through JJJJ=-31938 is presented below. What follows is a manual copy from the computer monitor.

6 5 1 3 7
8 4 2 -109 -31997

6 5 1 7 8
4 3 2 -107 -31995

7 8 4 2 6
5 1 3 -109 -31963

7 8 4 2 6
5 1 3 -109 -31961

3 4 8 2 1
5 6 7 -107 -31960

2 1 5 3 4
8 7 6 -107 -31957

7 8 4 6 5
1 2 3 -107 -31952

6 5 1 7 8
4 3 2 -107 -31938

107 is optimal, Nugent et al. (1968). Interpreted in accordance with line 1912 and line 1929, the output above was obtained in three seconds of running on a personal computer with an Intel 2.66 GHz. chip and with the IBM basica/D interpreter.

References

Hillier, F. S. (1963), "Quantitative tools for plant layout analysis," J. Indust. Eng. 14, 33-40.

Loiola, E. M., N. M. M. de Abreu, P. O. Boaventura-Netto, P. Hahn, and T. Querido
(2007), "A survey for the quadratic assignment problem," European Journal of Operational Research 176, 657-690.

Microsoft Corp. BASIC. Second Edition (May 1982), Version 1.10. Boca Raton, Florida: IBM Corp., Personal Computer, P. O. Box 1328-C, Boca Raton, Florida 33432, 1981.

Nugent, C. E., T. E. Vollmann, and J. Ruml (1968), "An experimental comparisons of techniques for the assignment of facilities to locations," Operations Research 16, 150-173.

Assigning Forty-Two Activities to Conflicting Facilities Subject to Less-Than-or-Equal-to Constraints

Jsun Yui Wong

The scheduling problems of Carlson and Nemhauser [2], Davis, Devine, and Lutz [3], and Arnold, Beckwith, and Jones [1] are related. The difference between the original Carlson and Nemhauser problem [2] and the problem of Davis, Devine, and Lutz [3] is "that their problem did not consider conflicting facilities," Davis, Devine, and Lutz [3, p. 225]. For the newer problem, "a cost is incurred if two activities are assigned to the same facility, or to two facilities whose operation is in conflict," Davis et al. [3, p. 224]. Line 32 and line 33 of the computer program below say the operation of facility 1 and facility 2 is in conflict; line 33 and line 34 say the operation of facility 2 and facility 3 is in conflict.

There is the constraint that period one can handle only five courses or less, and there is also the constraint that period six can handle only five courses or less. These two restrictions exemplify the less-than-or-equal-to constraints on pages 1096 of Arnold, Beckwith, and Jones [1] and are handled mainly by line 366 through line 400 of the program listed below. Carlson and Nemhausser [2] "seek to find a conflict-minimizing schedule, but do not require that the schedule be balanced," [1, p. 1096].

Like Lutz, Devine, Kumin, and Smith [4], the present paper has another type of restrictions. Courses 1 through 14 may be assigned to time blocks 1 through 4, courses 15 through 23 may be assigned in time blocks 2 through 4, and courses 24 through 42 may be assigned to time blocks 4 through 6. These restrictions are handled by line 341 through line 363.

The conflict costs used in line 401 through line 5981 of the computer program below come from Nugent, Vollmann, and Ruml [5] and Skorin-Kapov [6].

In this paper, A(1)=2, A(2)=2, A(3)=2,..., and A(42)=2 stand for course 1 in time block 2, course 2 in time block 2, course 3 in time block 2,..., and course 42 in time block 2, respectively.

The computer program below is modelled after the computer program of the post "Assigning Activities to Conflicting Facilities To Minimize Conflict Cost" posted in the present blog on March 7, 2008.

4 DEFINT A-Z
5 DEFDBL M,T,P
7 DIM A(56),X(56),T(1333)
9 DIM TBM(6,6)
32 TBM(1,1)=1:TBM(1,2)=1
33 TBM(2,1)=1:TBM(2,2)=1:TBM(2,3)=1
34 TBM(3,2)=1:TBM(3,3)=1
35 TBM(4,4)=1
36 TBM(5,5)=1
37 TBM(6,6)=1
42 FOR JJJJ=-32000 TO 32000
44 RANDOMIZE JJJJ
46 M=-1D+17
61 FOR KX=1 TO 14
63 A(KX)=2+FIX(RND*3)
65 NEXT KX
71 FOR KY=15 TO 23
73 A(KY)=2+FIX(RND*3)
75 NEXT KY
81 FOR KZ=24 TO 42
83 A(KZ)=4+FIX(RND*2)
85 NEXT KZ
111 IMAR=10+FIX(RND*1500)
128 FOR I=1 TO IMAR
129 FOR K=1 TO 42
301 X(K)=A(K)
302 NEXT K
311 IF RND<.33 THEN 341 ELSE IF RND<.32 THEN 351 ELSE 361
341 IAP1=1+FIX(RND*14)
343 X(IAP1)=1+FIX(RND*4)
346 GOTO 366
351 IAP2=15+FIX(RND*9)
353 X(IAP2)=2+FIX(RND*3)
356 GOTO 366
361 IAP3=24+FIX(RND*19)
363 X(IAP3)=4+FIX(RND*3)
366 SUM1=0:SUM2=0
371 FOR I11=1 TO 42
381 IF X(I11)=1 THEN 391 ELSE IF X(I11)=6 THEN 397 ELSE 398
391 SUM1=SUM1+1
393 GOTO 398
397 SUM2=SUM2+6
398 NEXT I11
399 IF SUM1>5 THEN 6670
400 IF SUM2>30 THEN 6670
401 T(1)=3*10000*(1/(10000*(1-TBM(X(1),X(2)))+1))
402 T(2)=2*10000*(1/(10000*(1-TBM(X(1),X(3)))+1))
403 T(3)=0*10000*(1/(10000*(1-TBM(X(1),X(4)))+1))
404 T(4)=0*10000*(1/(10000*(1-TBM(X(1),X(5)))+1))
405 T(5)=2*10000*(1/(10000*(1-TBM(X(1),X(6)))+1))
406 T(6)=10*10000*(1/(10000*(1-TBM(X(1),X(7)))+1))
407 T(7)=5*10000*(1/(10000*(1-TBM(X(1),X(8)))+1))
408 T(8)=0*10000*(1/(10000*(1-TBM(X(1),X(9)))+1))
409 T(9)=5*10000*(1/(10000*(1-TBM(X(1),X(10)))+1))
410 T(10)=2*10000*(1/(10000*(1-TBM(X(1),X(11)))+1))
411 T(11)=5*10000*(1/(10000*(1-TBM(X(1),X(12)))+1))
412 T(12)=0*10000*(1/(10000*(1-TBM(X(1),X(13)))+1))
413 T(13)=0*10000*(1/(10000*(1-TBM(X(1),X(14)))+1))
414 T(14)=2*10000*(1/(10000*(1-TBM(X(1),X(15)))+1))
415 T(15)=0*10000*(1/(10000*(1-TBM(X(1),X(16)))+1))
416 T(16)=5*10000*(1/(10000*(1-TBM(X(1),X(17)))+1))
417 T(17)=6*10000*(1/(10000*(1-TBM(X(1),X(18)))+1))
418 T(18)=3*10000*(1/(10000*(1-TBM(X(1),X(19)))+1))
419 T(19)=0*10000*(1/(10000*(1-TBM(X(1),X(20)))+1))
420 T(20)=1*10000*(1/(10000*(1-TBM(X(1),X(21)))+1))
421 T(21)=10*10000*(1/(10000*(1-TBM(X(1),X(22)))+1))
422 T(22)=0*10000*(1/(10000*(1-TBM(X(1),X(23)))+1))
423 T(23)=10*10000*(1/(10000*(1-TBM(X(1),X(24)))+1))
424 T(24)=2*10000*(1/(10000*(1-TBM(X(1),X(25)))+1))
425 T(25)=1*10000*(1/(10000*(1-TBM(X(1),X(26)))+1))
426 T(26)=1*10000*(1/(10000*(1-TBM(X(1),X(27)))+1))
427 T(27)=1*10000*(1/(10000*(1-TBM(X(1),X(28)))+1))
428 T(28)=0*10000*(1/(10000*(1-TBM(X(1),X(29)))+1))
429 T(29)=1*10000*(1/(10000*(1-TBM(X(1),X(30)))+1))
430 T(30)=4*10000*(1/(10000*(1-TBM(X(2),X(3)))+1))
431 T(31)=0*10000*(1/(10000*(1-TBM(X(2),X(4)))+1))
432 T(32)=10*10000*(1/(10000*(1-TBM(X(2),X(5)))+1))
433 T(33)=4*10000*(1/(10000*(1-TBM(X(2),X(6)))+1))
434 T(34)=0*10000*(1/(10000*(1-TBM(X(2),X(7)))+1))
435 T(35)=0*10000*(1/(10000*(1-TBM(X(2),X(8)))+1))
436 T(36)=2*10000*(1/(10000*(1-TBM(X(2),X(9)))+1))
437 T(37)=2*10000*(1/(10000*(1-TBM(X(2),X(10)))+1))
438 T(38)=1*10000*(1/(10000*(1-TBM(X(2),X(11)))+1))
439 T(39)=0*10000*(1/(10000*(1-TBM(X(2),X(12)))+1))
440 T(40)=5*10000*(1/(10000*(1-TBM(X(2),X(13)))+1))
441 T(41)=0*10000*(1/(10000*(1-TBM(X(2),X(14)))+1))
442 T(42)=0*10000*(1/(10000*(1-TBM(X(2),X(15)))+1))
443 T(43)=0*10000*(1/(10000*(1-TBM(X(2),X(16)))+1))
444 T(44)=0*10000*(1/(10000*(1-TBM(X(2),X(17)))+1))
445 T(45)=2*10000*(1/(10000*(1-TBM(X(2),X(18)))+1))
446 T(46)=0*10000*(1/(10000*(1-TBM(X(2),X(19)))+1))
447 T(47)=1*10000*(1/(10000*(1-TBM(X(2),X(20)))+1))
448 T(48)=6*10000*(1/(10000*(1-TBM(X(2),X(21)))+1))
449 T(49)=1*10000*(1/(10000*(1-TBM(X(2),X(22)))+1))
450 T(50)=0*10000*(1/(10000*(1-TBM(X(2),X(23)))+1))
451 T(51)=1*10000*(1/(10000*(1-TBM(X(2),X(24)))+1))
452 T(52)=2*10000*(1/(10000*(1-TBM(X(2),X(25)))+1))
453 T(53)=2*10000*(1/(10000*(1-TBM(X(2),X(26)))+1))
454 T(54)=5*10000*(1/(10000*(1-TBM(X(2),X(27)))+1))
455 T(55)=1*10000*(1/(10000*(1-TBM(X(2),X(28)))+1))
456 T(56)=10*10000*(1/(10000*(1-TBM(X(2),X(29)))+1))
457 T(57)=5*10000*(1/(10000*(1-TBM(X(2),X(30)))+1))
458 T(58)=3*10000*(1/(10000*(1-TBM(X(3),X(4)))+1))
459 T(59)=4*10000*(1/(10000*(1-TBM(X(3),X(5)))+1))
460 T(60)=0*10000*(1/(10000*(1-TBM(X(3),X(6)))+1))
461 T(61)=5*10000*(1/(10000*(1-TBM(X(3),X(7)))+1))
462 T(62)=5*10000*(1/(10000*(1-TBM(X(3),X(8)))+1))
463 T(63)=5*10000*(1/(10000*(1-TBM(X(3),X(9)))+1))
464 T(64)=1*10000*(1/(10000*(1-TBM(X(3),X(10)))+1))
465 T(65)=4*10000*(1/(10000*(1-TBM(X(3),X(11)))+1))
466 T(66)=1*10000*(1/(10000*(1-TBM(X(3),X(12)))+1))
467 T(67)=0*10000*(1/(10000*(1-TBM(X(3),X(13)))+1))
468 T(68)=4*10000*(1/(10000*(1-TBM(X(3),X(14)))+1))
469 T(69)=0*10000*(1/(10000*(1-TBM(X(3),X(15)))+1))
470 T(70)=4*10000*(1/(10000*(1-TBM(X(3),X(16)))+1))
471 T(71)=0*10000*(1/(10000*(1-TBM(X(3),X(17)))+1))
472 T(72)=6*10000*(1/(10000*(1-TBM(X(3),X(18)))+1))
473 T(73)=3*10000*(1/(10000*(1-TBM(X(3),X(19)))+1))
474 T(74)=2*10000*(1/(10000*(1-TBM(X(3),X(20)))+1))
475 T(75)=5*10000*(1/(10000*(1-TBM(X(3),X(21)))+1))
476 T(76)=5*10000*(1/(10000*(1-TBM(X(3),X(22)))+1))
477 T(77)=2*10000*(1/(10000*(1-TBM(X(3),X(23)))+1))
478 T(78)=1*10000*(1/(10000*(1-TBM(X(3),X(24)))+1))
479 T(79)=0*10000*(1/(10000*(1-TBM(X(3),X(25)))+1))
480 T(80)=0*10000*(1/(10000*(1-TBM(X(3),X(26)))+1))
481 T(81)=3*10000*(1/(10000*(1-TBM(X(3),X(27)))+1))
482 T(82)=1*10000*(1/(10000*(1-TBM(X(3),X(28)))+1))
483 T(83)=0*10000*(1/(10000*(1-TBM(X(3),X(29)))+1))
484 T(84)=2*10000*(1/(10000*(1-TBM(X(3),X(30)))+1))
485 T(85)=0*10000*(1/(10000*(1-TBM(X(4),X(5)))+1))
486 T(86)=0*10000*(1/(10000*(1-TBM(X(4),X(6)))+1))
487 T(87)=0*10000*(1/(10000*(1-TBM(X(4),X(7)))+1))
488 T(88)=2*10000*(1/(10000*(1-TBM(X(4),X(8)))+1))
489 T(89)=2*10000*(1/(10000*(1-TBM(X(4),X(9)))+1))
490 T(90)=0*10000*(1/(10000*(1-TBM(X(4),X(10)))+1))
491 T(91)=6*10000*(1/(10000*(1-TBM(X(4),X(11)))+1))
492 T(92)=0*10000*(1/(10000*(1-TBM(X(4),X(12)))+1))
493 T(93)=2*10000*(1/(10000*(1-TBM(X(4),X(13)))+1))
494 T(94)=5*10000*(1/(10000*(1-TBM(X(4),X(14)))+1))
495 T(95)=2*10000*(1/(10000*(1-TBM(X(4),X(15)))+1))
496 T(96)=5*10000*(1/(10000*(1-TBM(X(4),X(16)))+1))
497 T(97)=1*10000*(1/(10000*(1-TBM(X(4),X(17)))+1))
498 T(98)=1*10000*(1/(10000*(1-TBM(X(4),X(18)))+1))
499 T(99)=1*10000*(1/(10000*(1-TBM(X(4),X(19)))+1))
500 T(100)=1*10000*(1/(10000*(1-TBM(X(4),X(20)))+1))
501 T(101)=2*10000*(1/(10000*(1-TBM(X(4),X(21)))+1))
502 T(102)=2*10000*(1/(10000*(1-TBM(X(4),X(22)))+1))
503 T(103)=4*10000*(1/(10000*(1-TBM(X(4),X(23)))+1))
504 T(104)=0*10000*(1/(10000*(1-TBM(X(4),X(24)))+1))
505 T(105)=2*10000*(1/(10000*(1-TBM(X(4),X(25)))+1))
506 T(106)=0*10000*(1/(10000*(1-TBM(X(4),X(26)))+1))
507 T(107)=2*10000*(1/(10000*(1-TBM(X(4),X(27)))+1))
508 T(108)=2*10000*(1/(10000*(1-TBM(X(4),X(28)))+1))
509 T(109)=5*10000*(1/(10000*(1-TBM(X(4),X(29)))+1))
510 T(110)=5*10000*(1/(10000*(1-TBM(X(4),X(30)))+1))
511 T(111)=5*10000*(1/(10000*(1-TBM(X(5),X(6)))+1))
512 T(112)=2*10000*(1/(10000*(1-TBM(X(5),X(7)))+1))
513 T(113)=0*10000*(1/(10000*(1-TBM(X(5),X(8)))+1))
514 T(114)=0*10000*(1/(10000*(1-TBM(X(5),X(9)))+1))
515 T(115)=0*10000*(1/(10000*(1-TBM(X(5),X(10)))+1))
516 T(116)=0*10000*(1/(10000*(1-TBM(X(5),X(11)))+1))
517 T(117)=2*10000*(1/(10000*(1-TBM(X(5),X(12)))+1))
518 T(118)=0*10000*(1/(10000*(1-TBM(X(5),X(13)))+1))
519 T(119)=0*10000*(1/(10000*(1-TBM(X(5),X(14)))+1))
520 T(120)=0*10000*(1/(10000*(1-TBM(X(5),X(15)))+1))
521 T(121)=0*10000*(1/(10000*(1-TBM(X(5),X(16)))+1))
522 T(122)=2*10000*(1/(10000*(1-TBM(X(5),X(17)))+1))
523 T(123)=1*10000*(1/(10000*(1-TBM(X(5),X(18)))+1))
524 T(124)=0*10000*(1/(10000*(1-TBM(X(5),X(19)))+1))
525 T(125)=0*10000*(1/(10000*(1-TBM(X(5),X(20)))+1))
526 T(126)=2*10000*(1/(10000*(1-TBM(X(5),X(21)))+1))
527 T(127)=0*10000*(1/(10000*(1-TBM(X(5),X(22)))+1))
528 T(128)=5*10000*(1/(10000*(1-TBM(X(5),X(23)))+1))
529 T(129)=1*10000*(1/(10000*(1-TBM(X(5),X(24)))+1))
530 T(130)=0*10000*(1/(10000*(1-TBM(X(5),X(25)))+1))
531 T(131)=2*10000*(1/(10000*(1-TBM(X(5),X(26)))+1))
532 T(132)=1*10000*(1/(10000*(1-TBM(X(5),X(27)))+1))
533 T(133)=0*10000*(1/(10000*(1-TBM(X(5),X(28)))+1))
534 T(134)=2*10000*(1/(10000*(1-TBM(X(5),X(29)))+1))
535 T(135)=1*10000*(1/(10000*(1-TBM(X(5),X(30)))+1))
536 T(136)=1*10000*(1/(10000*(1-TBM(X(6),X(7)))+1))
537 T(137)=2*10000*(1/(10000*(1-TBM(X(6),X(8)))+1))
538 T(138)=2*10000*(1/(10000*(1-TBM(X(6),X(9)))+1))
539 T(139)=1*10000*(1/(10000*(1-TBM(X(6),X(10)))+1))
540 T(140)=4*10000*(1/(10000*(1-TBM(X(6),X(11)))+1))
541 T(141)=10*10000*(1/(10000*(1-TBM(X(6),X(12)))+1))
542 T(142)=10*10000*(1/(10000*(1-TBM(X(6),X(13)))+1))
543 T(143)=2*10000*(1/(10000*(1-TBM(X(6),X(14)))+1))
544 T(144)=5*10000*(1/(10000*(1-TBM(X(6),X(15)))+1))
545 T(145)=5*10000*(1/(10000*(1-TBM(X(6),X(16)))+1))
546 T(146)=0*10000*(1/(10000*(1-TBM(X(6),X(17)))+1))
547 T(147)=5*10000*(1/(10000*(1-TBM(X(6),X(18)))+1))
548 T(148)=0*10000*(1/(10000*(1-TBM(X(6),X(19)))+1))
549 T(149)=0*10000*(1/(10000*(1-TBM(X(6),X(20)))+1))
550 T(150)=0*10000*(1/(10000*(1-TBM(X(6),X(21)))+1))
551 T(151)=10*10000*(1/(10000*(1-TBM(X(6),X(22)))+1))
552 T(152)=0*10000*(1/(10000*(1-TBM(X(6),X(23)))+1))
553 T(153)=0*10000*(1/(10000*(1-TBM(X(6),X(24)))+1))
554 T(154)=0*10000*(1/(10000*(1-TBM(X(6),X(25)))+1))
555 T(155)=4*10000*(1/(10000*(1-TBM(X(6),X(26)))+1))
556 T(156)=0*10000*(1/(10000*(1-TBM(X(6),X(27)))+1))
557 T(157)=10*10000*(1/(10000*(1-TBM(X(6),X(28)))+1))
558 T(158)=1*10000*(1/(10000*(1-TBM(X(6),X(29)))+1))
559 T(159)=1*10000*(1/(10000*(1-TBM(X(6),X(30)))+1))
560 T(160)=10*10000*(1/(10000*(1-TBM(X(7),X(8)))+1))
561 T(161)=10*10000*(1/(10000*(1-TBM(X(7),X(9)))+1))
562 T(162)=5*10000*(1/(10000*(1-TBM(X(7),X(10)))+1))
563 T(163)=10*10000*(1/(10000*(1-TBM(X(7),X(11)))+1))
564 T(164)=10*10000*(1/(10000*(1-TBM(X(7),X(12)))+1))
565 T(165)=6*10000*(1/(10000*(1-TBM(X(7),X(13)))+1))
566 T(166)=0*10000*(1/(10000*(1-TBM(X(7),X(14)))+1))
567 T(167)=0*10000*(1/(10000*(1-TBM(X(7),X(15)))+1))
568 T(168)=10*10000*(1/(10000*(1-TBM(X(7),X(16)))+1))
569 T(169)=2*10000*(1/(10000*(1-TBM(X(7),X(17)))+1))
570 T(170)=1*10000*(1/(10000*(1-TBM(X(7),X(18)))+1))
571 T(171)=10*10000*(1/(10000*(1-TBM(X(7),X(19)))+1))
572 T(172)=1*10000*(1/(10000*(1-TBM(X(7),X(20)))+1))
573 T(173)=5*10000*(1/(10000*(1-TBM(X(7),X(21)))+1))
574 T(174)=5*10000*(1/(10000*(1-TBM(X(7),X(22)))+1))
575 T(175)=2*10000*(1/(10000*(1-TBM(X(7),X(23)))+1))
576 T(176)=3*10000*(1/(10000*(1-TBM(X(7),X(24)))+1))
577 T(177)=5*10000*(1/(10000*(1-TBM(X(7),X(25)))+1))
578 T(178)=0*10000*(1/(10000*(1-TBM(X(7),X(26)))+1))
579 T(179)=2*10000*(1/(10000*(1-TBM(X(7),X(27)))+1))
580 T(180)=0*10000*(1/(10000*(1-TBM(X(7),X(28)))+1))
581 T(181)=1*10000*(1/(10000*(1-TBM(X(7),X(29)))+1))
582 T(182)=3*10000*(1/(10000*(1-TBM(X(7),X(30)))+1))
583 T(183)=1*10000*(1/(10000*(1-TBM(X(8),X(9)))+1))
584 T(184)=3*10000*(1/(10000*(1-TBM(X(8),X(10)))+1))
585 T(185)=5*10000*(1/(10000*(1-TBM(X(8),X(11)))+1))
586 T(186)=0*10000*(1/(10000*(1-TBM(X(8),X(12)))+1))
587 T(187)=0*10000*(1/(10000*(1-TBM(X(8),X(13)))+1))
588 T(188)=0*10000*(1/(10000*(1-TBM(X(8),X(14)))+1))
589 T(189)=2*10000*(1/(10000*(1-TBM(X(8),X(15)))+1))
590 T(190)=4*10000*(1/(10000*(1-TBM(X(8),X(16)))+1))
591 T(191)=5*10000*(1/(10000*(1-TBM(X(8),X(17)))+1))
592 T(192)=2*10000*(1/(10000*(1-TBM(X(8),X(18)))+1))
593 T(193)=10*10000*(1/(10000*(1-TBM(X(8),X(19)))+1))
594 T(194)=6*10000*(1/(10000*(1-TBM(X(8),X(20)))+1))
595 T(195)=0*10000*(1/(10000*(1-TBM(X(8),X(21)))+1))
596 T(196)=5*10000*(1/(10000*(1-TBM(X(8),X(22)))+1))
597 T(197)=5*10000*(1/(10000*(1-TBM(X(8),X(23)))+1))
598 T(198)=2*10000*(1/(10000*(1-TBM(X(8),X(24)))+1))
599 T(199)=5*10000*(1/(10000*(1-TBM(X(8),X(25)))+1))
600 T(200)=0*10000*(1/(10000*(1-TBM(X(8),X(26)))+1))
601 T(201)=5*10000*(1/(10000*(1-TBM(X(8),X(27)))+1))
602 T(202)=5*10000*(1/(10000*(1-TBM(X(8),X(28)))+1))
603 T(203)=0*10000*(1/(10000*(1-TBM(X(8),X(29)))+1))
604 T(204)=2*10000*(1/(10000*(1-TBM(X(8),X(30)))+1))
605 T(205)=10*10000*(1/(10000*(1-TBM(X(9),X(10)))+1))
606 T(206)=2*10000*(1/(10000*(1-TBM(X(9),X(11)))+1))
607 T(207)=1*10000*(1/(10000*(1-TBM(X(9),X(12)))+1))
608 T(208)=5*10000*(1/(10000*(1-TBM(X(9),X(13)))+1))
609 T(209)=2*10000*(1/(10000*(1-TBM(X(9),X(14)))+1))
610 T(210)=0*10000*(1/(10000*(1-TBM(X(9),X(15)))+1))
611 T(211)=3*10000*(1/(10000*(1-TBM(X(9),X(16)))+1))
612 T(212)=0*10000*(1/(10000*(1-TBM(X(9),X(17)))+1))
613 T(213)=2*10000*(1/(10000*(1-TBM(X(9),X(18)))+1))
614 T(214)=0*10000*(1/(10000*(1-TBM(X(9),X(19)))+1))
615 T(215)=0*10000*(1/(10000*(1-TBM(X(9),X(20)))+1))
616 T(216)=4*10000*(1/(10000*(1-TBM(X(9),X(21)))+1))
617 T(217)=0*10000*(1/(10000*(1-TBM(X(9),X(22)))+1))
618 T(218)=5*10000*(1/(10000*(1-TBM(X(9),X(23)))+1))
619 T(219)=2*10000*(1/(10000*(1-TBM(X(9),X(24)))+1))
620 T(220)=0*10000*(1/(10000*(1-TBM(X(9),X(25)))+1))
621 T(221)=5*10000*(1/(10000*(1-TBM(X(9),X(26)))+1))
622 T(222)=2*10000*(1/(10000*(1-TBM(X(9),X(27)))+1))
623 T(223)=2*10000*(1/(10000*(1-TBM(X(9),X(28)))+1))
624 T(224)=5*10000*(1/(10000*(1-TBM(X(9),X(29)))+1))
625 T(225)=2*10000*(1/(10000*(1-TBM(X(9),X(30)))+1))
626 T(226)=5*10000*(1/(10000*(1-TBM(X(10),X(11)))+1))
627 T(227)=5*10000*(1/(10000*(1-TBM(X(10),X(12)))+1))
628 T(228)=6*10000*(1/(10000*(1-TBM(X(10),X(13)))+1))
629 T(229)=0*10000*(1/(10000*(1-TBM(X(10),X(14)))+1))
630 T(230)=1*10000*(1/(10000*(1-TBM(X(10),X(15)))+1))
631 T(231)=5*10000*(1/(10000*(1-TBM(X(10),X(16)))+1))
632 T(232)=5*10000*(1/(10000*(1-TBM(X(10),X(17)))+1))
633 T(233)=0*10000*(1/(10000*(1-TBM(X(10),X(18)))+1))
634 T(234)=5*10000*(1/(10000*(1-TBM(X(10),X(19)))+1))
635 T(235)=2*10000*(1/(10000*(1-TBM(X(10),X(20)))+1))
636 T(236)=3*10000*(1/(10000*(1-TBM(X(10),X(21)))+1))
637 T(237)=5*10000*(1/(10000*(1-TBM(X(10),X(22)))+1))
638 T(238)=0*10000*(1/(10000*(1-TBM(X(10),X(23)))+1))
639 T(239)=5*10000*(1/(10000*(1-TBM(X(10),X(24)))+1))
640 T(240)=2*10000*(1/(10000*(1-TBM(X(10),X(25)))+1))
641 T(241)=10*10000*(1/(10000*(1-TBM(X(10),X(26)))+1))
642 T(242)=10*10000*(1/(10000*(1-TBM(X(10),X(27)))+1))
643 T(243)=1*10000*(1/(10000*(1-TBM(X(10),X(28)))+1))
644 T(244)=5*10000*(1/(10000*(1-TBM(X(10),X(29)))+1))
645 T(245)=2*10000*(1/(10000*(1-TBM(X(10),X(30)))+1))
646 T(246)=0*10000*(1/(10000*(1-TBM(X(11),X(12)))+1))
647 T(247)=0*10000*(1/(10000*(1-TBM(X(11),X(13)))+1))
648 T(248)=1*10000*(1/(10000*(1-TBM(X(11),X(14)))+1))
649 T(249)=2*10000*(1/(10000*(1-TBM(X(11),X(15)))+1))
650 T(250)=1*10000*(1/(10000*(1-TBM(X(11),X(16)))+1))
651 T(251)=0*10000*(1/(10000*(1-TBM(X(11),X(17)))+1))
652 T(252)=2*10000*(1/(10000*(1-TBM(X(11),X(18)))+1))
653 T(253)=0*10000*(1/(10000*(1-TBM(X(11),X(19)))+1))
654 T(254)=0*10000*(1/(10000*(1-TBM(X(11),X(20)))+1))
655 T(255)=0*10000*(1/(10000*(1-TBM(X(11),X(21)))+1))
656 T(256)=6*10000*(1/(10000*(1-TBM(X(11),X(22)))+1))
657 T(257)=6*10000*(1/(10000*(1-TBM(X(11),X(23)))+1))
658 T(258)=0*10000*(1/(10000*(1-TBM(X(11),X(24)))+1))
659 T(259)=4*10000*(1/(10000*(1-TBM(X(11),X(25)))+1))
660 T(260)=5*10000*(1/(10000*(1-TBM(X(11),X(26)))+1))
661 T(261)=3*10000*(1/(10000*(1-TBM(X(11),X(27)))+1))
662 T(262)=2*10000*(1/(10000*(1-TBM(X(11),X(28)))+1))
663 T(263)=2*10000*(1/(10000*(1-TBM(X(11),X(29)))+1))
664 T(264)=10*10000*(1/(10000*(1-TBM(X(11),X(30)))+1))
665 T(265)=5*10000*(1/(10000*(1-TBM(X(12),X(13)))+1))
666 T(266)=5*10000*(1/(10000*(1-TBM(X(12),X(14)))+1))
667 T(267)=2*10000*(1/(10000*(1-TBM(X(12),X(15)))+1))
668 T(268)=0*10000*(1/(10000*(1-TBM(X(12),X(16)))+1))
669 T(269)=0*10000*(1/(10000*(1-TBM(X(12),X(17)))+1))
670 T(270)=0*10000*(1/(10000*(1-TBM(X(12),X(18)))+1))
671 T(271)=0*10000*(1/(10000*(1-TBM(X(12),X(19)))+1))
672 T(272)=2*10000*(1/(10000*(1-TBM(X(12),X(20)))+1))
673 T(273)=0*10000*(1/(10000*(1-TBM(X(12),X(21)))+1))
674 T(274)=4*10000*(1/(10000*(1-TBM(X(12),X(22)))+1))
675 T(275)=5*10000*(1/(10000*(1-TBM(X(12),X(23)))+1))
676 T(276)=10*10000*(1/(10000*(1-TBM(X(12),X(24)))+1))
677 T(277)=1*10000*(1/(10000*(1-TBM(X(12),X(25)))+1))
678 T(278)=0*10000*(1/(10000*(1-TBM(X(12),X(26)))+1))
679 T(279)=0*10000*(1/(10000*(1-TBM(X(12),X(27)))+1))
680 T(280)=0*10000*(1/(10000*(1-TBM(X(12),X(28)))+1))
681 T(281)=0*10000*(1/(10000*(1-TBM(X(12),X(29)))+1))
682 T(282)=1*10000*(1/(10000*(1-TBM(X(12),X(30)))+1))
683 T(283)=2*10000*(1/(10000*(1-TBM(X(13),X(14)))+1))
684 T(284)=0*10000*(1/(10000*(1-TBM(X(13),X(15)))+1))
685 T(285)=4*10000*(1/(10000*(1-TBM(X(13),X(16)))+1))
686 T(286)=2*10000*(1/(10000*(1-TBM(X(13),X(17)))+1))
687 T(287)=2*10000*(1/(10000*(1-TBM(X(13),X(18)))+1))
688 T(288)=1*10000*(1/(10000*(1-TBM(X(13),X(19)))+1))
689 T(289)=0*10000*(1/(10000*(1-TBM(X(13),X(20)))+1))
690 T(290)=6*10000*(1/(10000*(1-TBM(X(13),X(21)))+1))
691 T(291)=2*10000*(1/(10000*(1-TBM(X(13),X(22)))+1))
692 T(292)=1*10000*(1/(10000*(1-TBM(X(13),X(23)))+1))
693 T(293)=5*10000*(1/(10000*(1-TBM(X(13),X(24)))+1))
694 T(294)=5*10000*(1/(10000*(1-TBM(X(13),X(25)))+1))
695 T(295)=0*10000*(1/(10000*(1-TBM(X(13),X(26)))+1))
696 T(296)=0*10000*(1/(10000*(1-TBM(X(13),X(27)))+1))
697 T(297)=1*10000*(1/(10000*(1-TBM(X(13),X(28)))+1))
698 T(298)=5*10000*(1/(10000*(1-TBM(X(13),X(29)))+1))
699 T(299)=5*10000*(1/(10000*(1-TBM(X(13),X(30)))+1))
700 T(300)=2*10000*(1/(10000*(1-TBM(X(14),X(15)))+1))
701 T(301)=1*10000*(1/(10000*(1-TBM(X(14),X(16)))+1))
702 T(302)=0*10000*(1/(10000*(1-TBM(X(14),X(17)))+1))
703 T(303)=5*10000*(1/(10000*(1-TBM(X(14),X(18)))+1))
704 T(304)=3*10000*(1/(10000*(1-TBM(X(14),X(19)))+1))
705 T(305)=10*10000*(1/(10000*(1-TBM(X(14),X(20)))+1))
706 T(306)=0*10000*(1/(10000*(1-TBM(X(14),X(21)))+1))
707 T(307)=0*10000*(1/(10000*(1-TBM(X(14),X(22)))+1))
708 T(308)=4*10000*(1/(10000*(1-TBM(X(14),X(23)))+1))
709 T(309)=2*10000*(1/(10000*(1-TBM(X(14),X(24)))+1))
710 T(310)=0*10000*(1/(10000*(1-TBM(X(14),X(25)))+1))
711 T(311)=0*10000*(1/(10000*(1-TBM(X(14),X(26)))+1))
712 T(312)=4*10000*(1/(10000*(1-TBM(X(14),X(27)))+1))
713 T(313)=2*10000*(1/(10000*(1-TBM(X(14),X(28)))+1))
714 T(314)=5*10000*(1/(10000*(1-TBM(X(14),X(29)))+1))
715 T(315)=5*10000*(1/(10000*(1-TBM(X(14),X(30)))+1))
716 T(316)=4*10000*(1/(10000*(1-TBM(X(15),X(16)))+1))
717 T(317)=5*10000*(1/(10000*(1-TBM(X(15),X(17)))+1))
718 T(318)=1*10000*(1/(10000*(1-TBM(X(15),X(18)))+1))
719 T(319)=0*10000*(1/(10000*(1-TBM(X(15),X(19)))+1))
720 T(320)=1*10000*(1/(10000*(1-TBM(X(15),X(20)))+1))
721 T(321)=0*10000*(1/(10000*(1-TBM(X(15),X(21)))+1))
722 T(322)=5*10000*(1/(10000*(1-TBM(X(15),X(22)))+1))
723 T(323)=0*10000*(1/(10000*(1-TBM(X(15),X(23)))+1))
724 T(324)=2*10000*(1/(10000*(1-TBM(X(15),X(24)))+1))
725 T(325)=0*10000*(1/(10000*(1-TBM(X(15),X(25)))+1))
726 T(326)=0*10000*(1/(10000*(1-TBM(X(15),X(26)))+1))
727 T(327)=5*10000*(1/(10000*(1-TBM(X(15),X(27)))+1))
728 T(328)=1*10000*(1/(10000*(1-TBM(X(15),X(28)))+1))
729 T(329)=1*10000*(1/(10000*(1-TBM(X(15),X(29)))+1))
730 T(330)=0*10000*(1/(10000*(1-TBM(X(15),X(30)))+1))
731 T(331)=0*10000*(1/(10000*(1-TBM(X(16),X(17)))+1))
732 T(332)=3*10000*(1/(10000*(1-TBM(X(16),X(18)))+1))
733 T(333)=0*10000*(1/(10000*(1-TBM(X(16),X(19)))+1))
734 T(334)=2*10000*(1/(10000*(1-TBM(X(16),X(20)))+1))
735 T(335)=2*10000*(1/(10000*(1-TBM(X(16),X(21)))+1))
736 T(336)=0*10000*(1/(10000*(1-TBM(X(16),X(22)))+1))
737 T(337)=2*10000*(1/(10000*(1-TBM(X(16),X(23)))+1))
738 T(338)=0*10000*(1/(10000*(1-TBM(X(16),X(24)))+1))
739 T(339)=5*10000*(1/(10000*(1-TBM(X(16),X(25)))+1))
740 T(340)=0*10000*(1/(10000*(1-TBM(X(16),X(26)))+1))
741 T(341)=5*10000*(1/(10000*(1-TBM(X(16),X(27)))+1))
742 T(342)=2*10000*(1/(10000*(1-TBM(X(16),X(28)))+1))
743 T(343)=5*10000*(1/(10000*(1-TBM(X(16),X(29)))+1))
744 T(344)=10*10000*(1/(10000*(1-TBM(X(16),X(30)))+1))
745 T(345)=2*10000*(1/(10000*(1-TBM(X(17),X(18)))+1))
746 T(346)=2*10000*(1/(10000*(1-TBM(X(17),X(19)))+1))
747 T(347)=0*10000*(1/(10000*(1-TBM(X(17),X(20)))+1))
748 T(348)=0*10000*(1/(10000*(1-TBM(X(17),X(21)))+1))
749 T(349)=0*10000*(1/(10000*(1-TBM(X(17),X(22)))+1))
750 T(350)=6*10000*(1/(10000*(1-TBM(X(17),X(23)))+1))
751 T(351)=5*10000*(1/(10000*(1-TBM(X(17),X(24)))+1))
752 T(352)=3*10000*(1/(10000*(1-TBM(X(17),X(25)))+1))
753 T(353)=5*10000*(1/(10000*(1-TBM(X(17),X(26)))+1))
754 T(354)=0*10000*(1/(10000*(1-TBM(X(17),X(27)))+1))
755 T(355)=0*10000*(1/(10000*(1-TBM(X(17),X(28)))+1))
756 T(356)=5*10000*(1/(10000*(1-TBM(X(17),X(29)))+1))
757 T(357)=1*10000*(1/(10000*(1-TBM(X(17),X(30)))+1))
758 T(358)=5*10000*(1/(10000*(1-TBM(X(18),X(19)))+1))
759 T(359)=1*10000*(1/(10000*(1-TBM(X(18),X(20)))+1))
760 T(360)=2*10000*(1/(10000*(1-TBM(X(18),X(21)))+1))
761 T(361)=10*10000*(1/(10000*(1-TBM(X(18),X(22)))+1))
762 T(362)=10*10000*(1/(10000*(1-TBM(X(18),X(23)))+1))
763 T(363)=4*10000*(1/(10000*(1-TBM(X(18),X(24)))+1))
764 T(364)=0*10000*(1/(10000*(1-TBM(X(18),X(25)))+1))
765 T(365)=0*10000*(1/(10000*(1-TBM(X(18),X(26)))+1))
766 T(366)=5*10000*(1/(10000*(1-TBM(X(18),X(27)))+1))
767 T(367)=0*10000*(1/(10000*(1-TBM(X(18),X(28)))+1))
768 T(368)=0*10000*(1/(10000*(1-TBM(X(18),X(29)))+1))
769 T(369)=0*10000*(1/(10000*(1-TBM(X(18),X(30)))+1))
770 T(370)=0*10000*(1/(10000*(1-TBM(X(19),X(20)))+1))
771 T(371)=5*10000*(1/(10000*(1-TBM(X(19),X(21)))+1))
772 T(372)=5*10000*(1/(10000*(1-TBM(X(19),X(22)))+1))
773 T(373)=1*10000*(1/(10000*(1-TBM(X(19),X(23)))+1))
774 T(374)=0*10000*(1/(10000*(1-TBM(X(19),X(24)))+1))
775 T(375)=5*10000*(1/(10000*(1-TBM(X(19),X(25)))+1))
776 T(376)=2*10000*(1/(10000*(1-TBM(X(19),X(26)))+1))
777 T(377)=1*10000*(1/(10000*(1-TBM(X(19),X(27)))+1))
778 T(378)=2*10000*(1/(10000*(1-TBM(X(19),X(28)))+1))
779 T(379)=10*10000*(1/(10000*(1-TBM(X(19),X(29)))+1))
780 T(380)=10*10000*(1/(10000*(1-TBM(X(19),X(30)))+1))
781 T(381)=5*10000*(1/(10000*(1-TBM(X(20),X(21)))+1))
782 T(382)=2*10000*(1/(10000*(1-TBM(X(20),X(22)))+1))
783 T(383)=1*10000*(1/(10000*(1-TBM(X(20),X(23)))+1))
784 T(384)=3*10000*(1/(10000*(1-TBM(X(20),X(24)))+1))
785 T(385)=1*10000*(1/(10000*(1-TBM(X(20),X(25)))+1))
786 T(386)=5*10000*(1/(10000*(1-TBM(X(20),X(26)))+1))
787 T(387)=6*10000*(1/(10000*(1-TBM(X(20),X(27)))+1))
788 T(388)=5*10000*(1/(10000*(1-TBM(X(20),X(28)))+1))
789 T(389)=5*10000*(1/(10000*(1-TBM(X(20),X(29)))+1))
790 T(390)=3*10000*(1/(10000*(1-TBM(X(20),X(30)))+1))
791 T(391)=4*10000*(1/(10000*(1-TBM(X(21),X(22)))+1))
792 T(392)=0*10000*(1/(10000*(1-TBM(X(21),X(23)))+1))
793 T(393)=1*10000*(1/(10000*(1-TBM(X(21),X(24)))+1))
794 T(394)=0*10000*(1/(10000*(1-TBM(X(21),X(25)))+1))
795 T(395)=0*10000*(1/(10000*(1-TBM(X(21),X(26)))+1))
796 T(396)=0*10000*(1/(10000*(1-TBM(X(21),X(27)))+1))
797 T(397)=5*10000*(1/(10000*(1-TBM(X(21),X(28)))+1))
798 T(398)=0*10000*(1/(10000*(1-TBM(X(21),X(29)))+1))
799 T(399)=0*10000*(1/(10000*(1-TBM(X(21),X(30)))+1))
800 T(400)=5*10000*(1/(10000*(1-TBM(X(22),X(23)))+1))
801 T(401)=0*10000*(1/(10000*(1-TBM(X(22),X(24)))+1))
802 T(402)=4*10000*(1/(10000*(1-TBM(X(22),X(25)))+1))
803 T(403)=4*10000*(1/(10000*(1-TBM(X(22),X(26)))+1))
804 T(404)=5*10000*(1/(10000*(1-TBM(X(22),X(27)))+1))
805 T(405)=0*10000*(1/(10000*(1-TBM(X(22),X(28)))+1))
806 T(406)=2*10000*(1/(10000*(1-TBM(X(22),X(29)))+1))
807 T(407)=5*10000*(1/(10000*(1-TBM(X(22),X(30)))+1))
808 T(408)=0*10000*(1/(10000*(1-TBM(X(23),X(24)))+1))
809 T(409)=4*10000*(1/(10000*(1-TBM(X(23),X(25)))+1))
810 T(410)=4*10000*(1/(10000*(1-TBM(X(23),X(26)))+1))
811 T(411)=1*10000*(1/(10000*(1-TBM(X(23),X(27)))+1))
812 T(412)=0*10000*(1/(10000*(1-TBM(X(23),X(28)))+1))
813 T(413)=2*10000*(1/(10000*(1-TBM(X(23),X(29)))+1))
814 T(414)=2*10000*(1/(10000*(1-TBM(X(23),X(30)))+1))
815 T(415)=5*10000*(1/(10000*(1-TBM(X(24),X(25)))+1))
816 T(416)=5*10000*(1/(10000*(1-TBM(X(24),X(26)))+1))
817 T(417)=0*10000*(1/(10000*(1-TBM(X(24),X(27)))+1))
818 T(418)=1*10000*(1/(10000*(1-TBM(X(24),X(28)))+1))
819 T(419)=0*10000*(1/(10000*(1-TBM(X(24),X(29)))+1))
820 T(420)=0*10000*(1/(10000*(1-TBM(X(24),X(30)))+1))
821 T(421)=1*10000*(1/(10000*(1-TBM(X(25),X(26)))+1))
822 T(422)=0*10000*(1/(10000*(1-TBM(X(25),X(27)))+1))
823 T(423)=10*10000*(1/(10000*(1-TBM(X(25),X(28)))+1))
824 T(424)=1*10000*(1/(10000*(1-TBM(X(25),X(29)))+1))
825 T(425)=0*10000*(1/(10000*(1-TBM(X(25),X(30)))+1))
826 T(426)=0*10000*(1/(10000*(1-TBM(X(26),X(27)))+1))
827 T(427)=0*10000*(1/(10000*(1-TBM(X(26),X(28)))+1))
828 T(428)=0*10000*(1/(10000*(1-TBM(X(26),X(29)))+1))
829 T(429)=0*10000*(1/(10000*(1-TBM(X(26),X(30)))+1))
830 T(430)=0*10000*(1/(10000*(1-TBM(X(27),X(28)))+1))
831 T(431)=0*10000*(1/(10000*(1-TBM(X(27),X(29)))+1))
832 T(432)=10*10000*(1/(10000*(1-TBM(X(27),X(30)))+1))
833 T(433)=2*10000*(1/(10000*(1-TBM(X(28),X(29)))+1))
834 T(434)=2*10000*(1/(10000*(1-TBM(X(28),X(30)))+1))
835 T(435)=2*10000*(1/(10000*(1-TBM(X(29),X(30)))+1))
1429 T(436)=1*10000*(1/(10000*(1-TBM(X(1),X(31)))+1))
1430 T(437)=0*10000*(1/(10000*(1-TBM(X(1),X(32)))+1))
1431 T(438)=4*10000*(1/(10000*(1-TBM(X(1),X(33)))+1))
1432 T(439)=0*10000*(1/(10000*(1-TBM(X(1),X(34)))+1))
1433 T(440)=3*10000*(1/(10000*(1-TBM(X(1),X(35)))+1))
1434 T(441)=5*10000*(1/(10000*(1-TBM(X(1),X(36)))+1))
1435 T(442)=2*10000*(1/(10000*(1-TBM(X(1),X(37)))+1))
1436 T(443)=4*10000*(1/(10000*(1-TBM(X(1),X(38)))+1))
1437 T(444)=0*10000*(1/(10000*(1-TBM(X(1),X(39)))+1))
1438 T(445)=0*10000*(1/(10000*(1-TBM(X(1),X(40)))+1))
2557 T(446)=1*10000*(1/(10000*(1-TBM(X(2),X(31)))+1))
2558 T(447)=0*10000*(1/(10000*(1-TBM(X(2),X(32)))+1))
2559 T(448)=0*10000*(1/(10000*(1-TBM(X(2),X(33)))+1))
2560 T(449)=1*10000*(1/(10000*(1-TBM(X(2),X(34)))+1))
2561 T(450)=5*10000*(1/(10000*(1-TBM(X(2),X(35)))+1))
2562 T(451)=2*10000*(1/(10000*(1-TBM(X(2),X(36)))+1))
2563 T(452)=2*10000*(1/(10000*(1-TBM(X(2),X(37)))+1))
2564 T(453)=0*10000*(1/(10000*(1-TBM(X(2),X(38)))+1))
2565 T(454)=5*10000*(1/(10000*(1-TBM(X(2),X(39)))+1))
2566 T(455)=4*10000*(1/(10000*(1-TBM(X(2),X(40)))+1))
2684 T(456)=0*10000*(1/(10000*(1-TBM(X(3),X(31)))+1))
2685 T(457)=0*10000*(1/(10000*(1-TBM(X(3),X(32)))+1))
2686 T(458)=0*10000*(1/(10000*(1-TBM(X(3),X(33)))+1))
2687 T(459)=0*10000*(1/(10000*(1-TBM(X(3),X(34)))+1))
2688 T(460)=1*10000*(1/(10000*(1-TBM(X(3),X(35)))+1))
2689 T(461)=5*10000*(1/(10000*(1-TBM(X(3),X(36)))+1))
2690 T(462)=10*10000*(1/(10000*(1-TBM(X(3),X(37)))+1))
2691 T(463)=5*10000*(1/(10000*(1-TBM(X(3),X(38)))+1))
2692 T(464)=1*10000*(1/(10000*(1-TBM(X(3),X(39)))+1))
2693 T(465)=0*10000*(1/(10000*(1-TBM(X(3),X(40)))+1))
2740 T(467)=5*10000*(1/(10000*(1-TBM(X(4),X(31)))+1))
2741 T(468)=2*10000*(1/(10000*(1-TBM(X(4),X(32)))+1))
2742 T(469)=0*10000*(1/(10000*(1-TBM(X(4),X(33)))+1))
2743 T(470)=0*10000*(1/(10000*(1-TBM(X(4),X(34)))+1))
2744 T(471)=5*10000*(1/(10000*(1-TBM(X(4),X(35)))+1))
2745 T(472)=10*10000*(1/(10000*(1-TBM(X(4),X(36)))+1))
2746 T(473)=0*10000*(1/(10000*(1-TBM(X(4),X(37)))+1))
2747 T(474)=4*10000*(1/(10000*(1-TBM(X(4),X(38)))+1))
2748 T(475)=5*10000*(1/(10000*(1-TBM(X(4),X(39)))+1))
2749 T(476)=2*10000*(1/(10000*(1-TBM(X(4),X(40)))+1))
2835 T(487)=1*10000*(1/(10000*(1-TBM(X(5),X(31)))+1))
2836 T(488)=2*10000*(1/(10000*(1-TBM(X(5),X(32)))+1))
2837 T(489)=1*10000*(1/(10000*(1-TBM(X(5),X(33)))+1))
2838 T(490)=5*10000*(1/(10000*(1-TBM(X(5),X(34)))+1))
2839 T(491)=2*10000*(1/(10000*(1-TBM(X(5),X(35)))+1))
2840 T(492)=5*10000*(1/(10000*(1-TBM(X(5),X(36)))+1))
2841 T(493)=10*10000*(1/(10000*(1-TBM(X(5),X(37)))+1))
2842 T(494)=5*10000*(1/(10000*(1-TBM(X(5),X(38)))+1))
2843 T(495)=10*10000*(1/(10000*(1-TBM(X(5),X(39)))+1))
2844 T(496)=0*10000*(1/(10000*(1-TBM(X(5),X(40)))+1))
2959 T(497)=2*10000*(1/(10000*(1-TBM(X(6),X(31)))+1))
2960 T(498)=1*10000*(1/(10000*(1-TBM(X(6),X(32)))+1))
2961 T(499)=0*10000*(1/(10000*(1-TBM(X(6),X(33)))+1))
2962 T(500)=1*10000*(1/(10000*(1-TBM(X(6),X(34)))+1))
2963 T(501)=1*10000*(1/(10000*(1-TBM(X(6),X(35)))+1))
2964 T(502)=0*10000*(1/(10000*(1-TBM(X(6),X(36)))+1))
2965 T(503)=5*10000*(1/(10000*(1-TBM(X(6),X(37)))+1))
2966 T(504)=0*10000*(1/(10000*(1-TBM(X(6),X(38)))+1))
2967 T(505)=10*10000*(1/(10000*(1-TBM(X(6),X(39)))+1))
2968 T(506)=0*10000*(1/(10000*(1-TBM(X(6),X(40)))+1))
3052 T(507)=5*10000*(1/(10000*(1-TBM(X(7),X(31)))+1))
3053 T(508)=10*10000*(1/(10000*(1-TBM(X(7),X(32)))+1))
3054 T(509)=5*10000*(1/(10000*(1-TBM(X(7),X(33)))+1))
3055 T(510)=5*10000*(1/(10000*(1-TBM(X(7),X(34)))+1))
3056 T(511)=10*10000*(1/(10000*(1-TBM(X(7),X(35)))+1))
3057 T(512)=0*10000*(1/(10000*(1-TBM(X(7),X(36)))+1))
3058 T(513)=0*10000*(1/(10000*(1-TBM(X(7),X(37)))+1))
3059 T(514)=1*10000*(1/(10000*(1-TBM(X(7),X(38)))+1))
3060 T(515)=2*10000*(1/(10000*(1-TBM(X(7),X(39)))+1))
3061 T(516)=0*10000*(1/(10000*(1-TBM(X(7),X(40)))+1))
3101 T(517)=5*10000*(1/(10000*(1-TBM(X(8),X(31)))+1))
3102 T(518)=0*10000*(1/(10000*(1-TBM(X(8),X(32)))+1))
3103 T(519)=0*10000*(1/(10000*(1-TBM(X(8),X(33)))+1))
3104 T(520)=2*10000*(1/(10000*(1-TBM(X(8),X(34)))+1))
3105 T(521)=0*10000*(1/(10000*(1-TBM(X(8),X(35)))+1))
3106 T(522)=2*10000*(1/(10000*(1-TBM(X(8),X(36)))+1))
3107 T(523)=1*10000*(1/(10000*(1-TBM(X(8),X(37)))+1))
3108 T(524)=1*10000*(1/(10000*(1-TBM(X(8),X(38)))+1))
3109 T(525)=2*10000*(1/(10000*(1-TBM(X(8),X(39)))+1))
3110 T(526)=10*10000*(1/(10000*(1-TBM(X(8),X(40)))+1))
3201 T(527)=5*10000*(1/(10000*(1-TBM(X(9),X(31)))+1))
3202 T(528)=0*10000*(1/(10000*(1-TBM(X(9),X(32)))+1))
3203 T(529)=10*10000*(1/(10000*(1-TBM(X(9),X(33)))+1))
3204 T(530)=1*10000*(1/(10000*(1-TBM(X(9),X(34)))+1))
3205 T(531)=0*10000*(1/(10000*(1-TBM(X(9),X(35)))+1))
3206 T(532)=3*10000*(1/(10000*(1-TBM(X(9),X(36)))+1))
3207 T(533)=0*10000*(1/(10000*(1-TBM(X(9),X(37)))+1))
3208 T(534)=0*10000*(1/(10000*(1-TBM(X(9),X(38)))+1))
3209 T(535)=1*10000*(1/(10000*(1-TBM(X(9),X(39)))+1))
3210 T(536)=2*10000*(1/(10000*(1-TBM(X(9),X(40)))+1))
3301 T(537)=5*10000*(1/(10000*(1-TBM(X(10),X(31)))+1))
3302 T(538)=2*10000*(1/(10000*(1-TBM(X(10),X(32)))+1))
3303 T(539)=1*10000*(1/(10000*(1-TBM(X(10),X(33)))+1))
3304 T(540)=0*10000*(1/(10000*(1-TBM(X(10),X(34)))+1))
3305 T(541)=1*10000*(1/(10000*(1-TBM(X(10),X(35)))+1))
3306 T(542)=1*10000*(1/(10000*(1-TBM(X(10),X(36)))+1))
3307 T(543)=1*10000*(1/(10000*(1-TBM(X(10),X(37)))+1))
3308 T(544)=5*10000*(1/(10000*(1-TBM(X(10),X(38)))+1))
3309 T(545)=6*10000*(1/(10000*(1-TBM(X(10),X(39)))+1))
3310 T(546)=3*10000*(1/(10000*(1-TBM(X(10),X(40)))+1))
3401 T(547)=2*10000*(1/(10000*(1-TBM(X(11),X(31)))+1))
3402 T(548)=0*10000*(1/(10000*(1-TBM(X(11),X(32)))+1))
3403 T(549)=0*10000*(1/(10000*(1-TBM(X(11),X(33)))+1))
3404 T(550)=0*10000*(1/(10000*(1-TBM(X(11),X(34)))+1))
3405 T(551)=0*10000*(1/(10000*(1-TBM(X(11),X(35)))+1))
3406 T(552)=0*10000*(1/(10000*(1-TBM(X(11),X(36)))+1))
3407 T(553)=0*10000*(1/(10000*(1-TBM(X(11),X(37)))+1))
3408 T(554)=0*10000*(1/(10000*(1-TBM(X(11),X(38)))+1))
3409 T(555)=5*10000*(1/(10000*(1-TBM(X(11),X(39)))+1))
3410 T(556)=1*10000*(1/(10000*(1-TBM(X(11),X(40)))+1))
3501 T(557)=5*10000*(1/(10000*(1-TBM(X(12),X(31)))+1))
3503 T(558)=0*10000*(1/(10000*(1-TBM(X(12),X(33)))+1))
3504 T(559)=0*10000*(1/(10000*(1-TBM(X(12),X(34)))+1))
3505 T(560)=0*10000*(1/(10000*(1-TBM(X(12),X(35)))+1))
3506 T(561)=2*10000*(1/(10000*(1-TBM(X(12),X(36)))+1))
3507 T(562)=0*10000*(1/(10000*(1-TBM(X(12),X(37)))+1))
3508 T(563)=0*10000*(1/(10000*(1-TBM(X(12),X(38)))+1))
3509 T(564)=6*10000*(1/(10000*(1-TBM(X(12),X(39)))+1))
3510 T(565)=1*10000*(1/(10000*(1-TBM(X(12),X(40)))+1))
3601 T(566)=0*10000*(1/(10000*(1-TBM(X(13),X(31)))+1))
3602 T(567)=1*10000*(1/(10000*(1-TBM(X(13),X(32)))+1))
3603 T(568)=1*10000*(1/(10000*(1-TBM(X(13),X(33)))+1))
3604 T(569)=1*10000*(1/(10000*(1-TBM(X(13),X(34)))+1))
3605 T(570)=0*10000*(1/(10000*(1-TBM(X(13),X(35)))+1))
3606 T(571)=2*10000*(1/(10000*(1-TBM(X(13),X(36)))+1))
3607 T(572)=0*10000*(1/(10000*(1-TBM(X(13),X(37)))+1))
3608 T(573)=2*10000*(1/(10000*(1-TBM(X(13),X(38)))+1))
3609 T(574)=1*10000*(1/(10000*(1-TBM(X(13),X(39)))+1))
3610 T(575)=1*10000*(1/(10000*(1-TBM(X(13),X(40)))+1))
3701 T(576)=2*10000*(1/(10000*(1-TBM(X(14),X(31)))+1))
3702 T(577)=1*10000*(1/(10000*(1-TBM(X(14),X(32)))+1))
3703 T(578)=3*10000*(1/(10000*(1-TBM(X(14),X(33)))+1))
3704 T(579)=2*10000*(1/(10000*(1-TBM(X(14),X(34)))+1))
3705 T(580)=5*10000*(1/(10000*(1-TBM(X(14),X(35)))+1))
3706 T(581)=0*10000*(1/(10000*(1-TBM(X(14),X(36)))+1))
3707 T(582)=10*10000*(1/(10000*(1-TBM(X(14),X(37)))+1))
3708 T(583)=2*10000*(1/(10000*(1-TBM(X(14),X(38)))+1))
3709 T(584)=10*10000*(1/(10000*(1-TBM(X(14),X(39)))+1))
3710 T(585)=0*10000*(1/(10000*(1-TBM(X(14),X(40)))+1))
3801 T(586)=0*10000*(1/(10000*(1-TBM(X(15),X(31)))+1))
3802 T(587)=0*10000*(1/(10000*(1-TBM(X(15),X(32)))+1))
3803 T(588)=0*10000*(1/(10000*(1-TBM(X(15),X(33)))+1))
3804 T(589)=0*10000*(1/(10000*(1-TBM(X(15),X(34)))+1))
3805 T(590)=5*10000*(1/(10000*(1-TBM(X(15),X(35)))+1))
3806 T(591)=1*10000*(1/(10000*(1-TBM(X(15),X(36)))+1))
3807 T(592)=0*10000*(1/(10000*(1-TBM(X(15),X(37)))+1))
3808 T(593)=3*10000*(1/(10000*(1-TBM(X(15),X(38)))+1))
3809 T(594)=5*10000*(1/(10000*(1-TBM(X(15),X(39)))+1))
3810 T(595)=0*10000*(1/(10000*(1-TBM(X(15),X(40)))+1))
3901 T(596)=5*10000*(1/(10000*(1-TBM(X(16),X(31)))+1))
3902 T(597)=0*10000*(1/(10000*(1-TBM(X(16),X(32)))+1))
3903 T(598)=5*10000*(1/(10000*(1-TBM(X(16),X(33)))+1))
3904 T(599)=5*10000*(1/(10000*(1-TBM(X(16),X(34)))+1))
3905 T(600)=0*10000*(1/(10000*(1-TBM(X(16),X(35)))+1))
3906 T(601)=10*10000*(1/(10000*(1-TBM(X(16),X(36)))+1))
3907 T(602)=0*10000*(1/(10000*(1-TBM(X(16),X(37)))+1))
3908 T(603)=2*10000*(1/(10000*(1-TBM(X(16),X(38)))+1))
3909 T(604)=5*10000*(1/(10000*(1-TBM(X(16),X(39)))+1))
3910 T(605)=0*10000*(1/(10000*(1-TBM(X(16),X(40)))+1))
4001 T(606)=0*10000*(1/(10000*(1-TBM(X(17),X(31)))+1))
4002 T(607)=10*10000*(1/(10000*(1-TBM(X(17),X(32)))+1))
4003 T(608)=2*10000*(1/(10000*(1-TBM(X(17),X(33)))+1))
4004 T(609)=5*10000*(1/(10000*(1-TBM(X(17),X(34)))+1))
4005 T(610)=4*10000*(1/(10000*(1-TBM(X(17),X(35)))+1))
4006 T(611)=10*10000*(1/(10000*(1-TBM(X(17),X(36)))+1))
4007 T(612)=4*10000*(1/(10000*(1-TBM(X(17),X(37)))+1))
4008 T(613)=3*10000*(1/(10000*(1-TBM(X(17),X(38)))+1))
4009 T(614)=5*10000*(1/(10000*(1-TBM(X(17),X(39)))+1))
4010 T(615)=0*10000*(1/(10000*(1-TBM(X(17),X(40)))+1))
4169 T(616)=5*10000*(1/(10000*(1-TBM(X(18),X(31)))+1))
4170 T(617)=0*10000*(1/(10000*(1-TBM(X(18),X(32)))+1))
4171 T(618)=0*10000*(1/(10000*(1-TBM(X(18),X(33)))+1))
4172 T(619)=10*10000*(1/(10000*(1-TBM(X(18),X(34)))+1))
4173 T(620)=5*10000*(1/(10000*(1-TBM(X(18),X(35)))+1))
4174 T(621)=0*10000*(1/(10000*(1-TBM(X(18),X(36)))+1))
4175 T(622)=2*10000*(1/(10000*(1-TBM(X(18),X(37)))+1))
4176 T(623)=2*10000*(1/(10000*(1-TBM(X(18),X(38)))+1))
4177 T(624)=10*10000*(1/(10000*(1-TBM(X(18),X(39)))+1))
4178 T(625)=2*10000*(1/(10000*(1-TBM(X(18),X(40)))+1))
4280 T(626)=1*10000*(1/(10000*(1-TBM(X(19),X(31)))+1))
4281 T(627)=2*10000*(1/(10000*(1-TBM(X(19),X(32)))+1))
4282 T(628)=0*10000*(1/(10000*(1-TBM(X(19),X(33)))+1))
4283 T(629)=10*10000*(1/(10000*(1-TBM(X(19),X(34)))+1))
4284 T(630)=0*10000*(1/(10000*(1-TBM(X(19),X(35)))+1))
4285 T(631)=1*10000*(1/(10000*(1-TBM(X(19),X(36)))+1))
4286 T(632)=0*10000*(1/(10000*(1-TBM(X(19),X(37)))+1))
4287 T(633)=6*10000*(1/(10000*(1-TBM(X(19),X(38)))+1))
4288 T(634)=2*10000*(1/(10000*(1-TBM(X(19),X(39)))+1))
4289 T(635)=0*10000*(1/(10000*(1-TBM(X(19),X(40)))+1))
4390 T(636)=5*10000*(1/(10000*(1-TBM(X(20),X(31)))+1))
4391 T(637)=1*10000*(1/(10000*(1-TBM(X(20),X(32)))+1))
4392 T(638)=5*10000*(1/(10000*(1-TBM(X(20),X(33)))+1))
4393 T(639)=2*10000*(1/(10000*(1-TBM(X(20),X(34)))+1))
4394 T(640)=5*10000*(1/(10000*(1-TBM(X(20),X(35)))+1))
4395 T(641)=1*10000*(1/(10000*(1-TBM(X(20),X(36)))+1))
4396 T(642)=0*10000*(1/(10000*(1-TBM(X(20),X(37)))+1))
4397 T(643)=1*10000*(1/(10000*(1-TBM(X(20),X(38)))+1))
4398 T(644)=0*10000*(1/(10000*(1-TBM(X(20),X(39)))+1))
4399 T(645)=2*10000*(1/(10000*(1-TBM(X(20),X(40)))+1))
4499 T(646)=2*10000*(1/(10000*(1-TBM(X(21),X(31)))+1))
4500 T(647)=0*10000*(1/(10000*(1-TBM(X(21),X(32)))+1))
4501 T(648)=0*10000*(1/(10000*(1-TBM(X(21),X(33)))+1))
4502 T(649)=4*10000*(1/(10000*(1-TBM(X(21),X(34)))+1))
4503 T(650)=10*10000*(1/(10000*(1-TBM(X(21),X(35)))+1))
4504 T(651)=0*10000*(1/(10000*(1-TBM(X(21),X(36)))+1))
4505 T(652)=5*10000*(1/(10000*(1-TBM(X(21),X(37)))+1))
4506 T(653)=0*10000*(1/(10000*(1-TBM(X(21),X(38)))+1))
4507 T(654)=0*10000*(1/(10000*(1-TBM(X(21),X(39)))+1))
4508 T(655)=0*10000*(1/(10000*(1-TBM(X(21),X(40)))+1))
4537 T(656)=2*10000*(1/(10000*(1-TBM(X(22),X(31)))+1))
4538 T(657)=10*10000*(1/(10000*(1-TBM(X(22),X(32)))+1))
4539 T(658)=1*10000*(1/(10000*(1-TBM(X(22),X(33)))+1))
4540 T(659)=10*10000*(1/(10000*(1-TBM(X(22),X(34)))+1))
4541 T(660)=5*10000*(1/(10000*(1-TBM(X(22),X(35)))+1))
4542 T(661)=1*10000*(1/(10000*(1-TBM(X(22),X(36)))+1))
4543 T(662)=1*10000*(1/(10000*(1-TBM(X(22),X(37)))+1))
4544 T(663)=1*10000*(1/(10000*(1-TBM(X(22),X(38)))+1))
4545 T(664)=2*10000*(1/(10000*(1-TBM(X(22),X(39)))+1))
4546 T(665)=0*10000*(1/(10000*(1-TBM(X(22),X(40)))+1))
4614 T(666)=10*10000*(1/(10000*(1-TBM(X(23),X(31)))+1))
4615 T(667)=5*10000*(1/(10000*(1-TBM(X(23),X(32)))+1))
4616 T(668)=2*10000*(1/(10000*(1-TBM(X(23),X(33)))+1))
4617 T(669)=5*10000*(1/(10000*(1-TBM(X(23),X(34)))+1))
4618 T(670)=5*10000*(1/(10000*(1-TBM(X(23),X(35)))+1))
4619 T(671)=2*10000*(1/(10000*(1-TBM(X(23),X(36)))+1))
4620 T(672)=0*10000*(1/(10000*(1-TBM(X(23),X(37)))+1))
4621 T(673)=5*10000*(1/(10000*(1-TBM(X(23),X(38)))+1))
4622 T(674)=5*10000*(1/(10000*(1-TBM(X(23),X(39)))+1))
4623 T(675)=10*10000*(1/(10000*(1-TBM(X(23),X(40)))+1))
4720 T(676)=0*10000*(1/(10000*(1-TBM(X(24),X(31)))+1))
4721 T(677)=1*10000*(1/(10000*(1-TBM(X(24),X(32)))+1))
4722 T(678)=2*10000*(1/(10000*(1-TBM(X(24),X(33)))+1))
4723 T(679)=1*10000*(1/(10000*(1-TBM(X(24),X(34)))+1))
4724 T(680)=2*10000*(1/(10000*(1-TBM(X(24),X(35)))+1))
4725 T(681)=2*10000*(1/(10000*(1-TBM(X(24),X(36)))+1))
4726 T(682)=5*10000*(1/(10000*(1-TBM(X(24),X(37)))+1))
4727 T(683)=6*10000*(1/(10000*(1-TBM(X(24),X(38)))+1))
4728 T(684)=10*10000*(1/(10000*(1-TBM(X(24),X(39)))+1))
4729 T(685)=2*10000*(1/(10000*(1-TBM(X(24),X(40)))+1))
4825 T(686)=5*10000*(1/(10000*(1-TBM(X(25),X(31)))+1))
4826 T(687)=10*10000*(1/(10000*(1-TBM(X(25),X(32)))+1))
4827 T(688)=10*10000*(1/(10000*(1-TBM(X(25),X(33)))+1))
4828 T(689)=10*10000*(1/(10000*(1-TBM(X(25),X(34)))+1))
4829 T(690)=1*10000*(1/(10000*(1-TBM(X(25),X(35)))+1))
4830 T(691)=2*10000*(1/(10000*(1-TBM(X(25),X(36)))+1))
4831 T(692)=2*10000*(1/(10000*(1-TBM(X(25),X(37)))+1))
4832 T(693)=2*10000*(1/(10000*(1-TBM(X(25),X(38)))+1))
4833 T(694)=2*10000*(1/(10000*(1-TBM(X(25),X(39)))+1))
4834 T(695)=5*10000*(1/(10000*(1-TBM(X(25),X(40)))+1))
4929 T(696)=4*10000*(1/(10000*(1-TBM(X(26),X(31)))+1))
4930 T(697)=0*10000*(1/(10000*(1-TBM(X(26),X(32)))+1))
4931 T(698)=0*10000*(1/(10000*(1-TBM(X(26),X(33)))+1))
4932 T(699)=5*10000*(1/(10000*(1-TBM(X(26),X(34)))+1))
4933 T(700)=5*10000*(1/(10000*(1-TBM(X(26),X(35)))+1))
4934 T(701)=2*10000*(1/(10000*(1-TBM(X(26),X(36)))+1))
4935 T(702)=5*10000*(1/(10000*(1-TBM(X(26),X(37)))+1))
4936 T(703)=10*10000*(1/(10000*(1-TBM(X(26),X(38)))+1))
4937 T(704)=0*10000*(1/(10000*(1-TBM(X(26),X(39)))+1))
4938 T(705)=5*10000*(1/(10000*(1-TBM(X(26),X(40)))+1))
5032 T(706)=5*10000*(1/(10000*(1-TBM(X(27),X(31)))+1))
5033 T(707)=5*10000*(1/(10000*(1-TBM(X(27),X(32)))+1))
5034 T(708)=0*10000*(1/(10000*(1-TBM(X(27),X(33)))+1))
5035 T(709)=2*10000*(1/(10000*(1-TBM(X(27),X(34)))+1))
5036 T(710)=3*10000*(1/(10000*(1-TBM(X(27),X(35)))+1))
5037 T(711)=10*10000*(1/(10000*(1-TBM(X(27),X(36)))+1))
5038 T(712)=1*10000*(1/(10000*(1-TBM(X(27),X(37)))+1))
5039 T(713)=1*10000*(1/(10000*(1-TBM(X(27),X(38)))+1))
5040 T(714)=0*10000*(1/(10000*(1-TBM(X(27),X(39)))+1))
5041 T(715)=0*10000*(1/(10000*(1-TBM(X(27),X(40)))+1))
5134 T(716)=2*10000*(1/(10000*(1-TBM(X(28),X(31)))+1))
5135 T(717)=2*10000*(1/(10000*(1-TBM(X(28),X(32)))+1))
5136 T(718)=0*10000*(1/(10000*(1-TBM(X(28),X(33)))+1))
5137 T(719)=0*10000*(1/(10000*(1-TBM(X(28),X(34)))+1))
5138 T(720)=0*10000*(1/(10000*(1-TBM(X(28),X(35)))+1))
5139 T(721)=0*10000*(1/(10000*(1-TBM(X(28),X(36)))+1))
5140 T(722)=1*10000*(1/(10000*(1-TBM(X(28),X(37)))+1))
5141 T(723)=0*10000*(1/(10000*(1-TBM(X(28),X(38)))+1))
5142 T(724)=5*10000*(1/(10000*(1-TBM(X(28),X(39)))+1))
5143 T(725)=10*10000*(1/(10000*(1-TBM(X(28),X(40)))+1))
5235 REM
5236 T(726)=5*10000*(1/(10000*(1-TBM(X(29),X(31)))+1))
5237 T(727)=10*10000*(1/(10000*(1-TBM(X(29),X(32)))+1))
5238 T(728)=1*10000*(1/(10000*(1-TBM(X(29),X(33)))+1))
5239 T(729)=10*10000*(1/(10000*(1-TBM(X(29),X(34)))+1))
5240 T(730)=2*10000*(1/(10000*(1-TBM(X(29),X(35)))+1))
5241 T(731)=10*10000*(1/(10000*(1-TBM(X(29),X(36)))+1))
5242 T(732)=2*10000*(1/(10000*(1-TBM(X(29),X(37)))+1))
5243 T(733)=5*10000*(1/(10000*(1-TBM(X(29),X(38)))+1))
5244 T(734)=5*10000*(1/(10000*(1-TBM(X(29),X(39)))+1))
5245 T(735)=0*10000*(1/(10000*(1-TBM(X(29),X(40)))+1))
5301 T(801)=2*10000*(1/(10000*(1-TBM(X(30),X(31)))+1))
5302 T(802)=2*10000*(1/(10000*(1-TBM(X(30),X(32)))+1))
5303 T(803)=1*10000*(1/(10000*(1-TBM(X(30),X(33)))+1))
5304 T(804)=0*10000*(1/(10000*(1-TBM(X(30),X(34)))+1))
5305 T(805)=1*10000*(1/(10000*(1-TBM(X(30),X(35)))+1))
5306 T(806)=5*10000*(1/(10000*(1-TBM(X(30),X(36)))+1))
5307 T(807)=5*10000*(1/(10000*(1-TBM(X(30),X(37)))+1))
5308 T(808)=10*10000*(1/(10000*(1-TBM(X(30),X(38)))+1))
5309 T(809)=0*10000*(1/(10000*(1-TBM(X(30),X(39)))+1))
5310 T(810)=2*10000*(1/(10000*(1-TBM(X(30),X(40)))+1))
5331 T(811)=1*10000*(1/(10000*(1-TBM(X(31),X(32)))+1))
5332 T(812)=2*10000*(1/(10000*(1-TBM(X(31),X(33)))+1))
5333 T(813)=1*10000*(1/(10000*(1-TBM(X(31),X(34)))+1))
5334 T(814)=1*10000*(1/(10000*(1-TBM(X(31),X(35)))+1))
5335 T(815)=10*10000*(1/(10000*(1-TBM(X(31),X(36)))+1))
5336 T(816)=10*10000*(1/(10000*(1-TBM(X(31),X(37)))+1))
5337 T(817)=0*10000*(1/(10000*(1-TBM(X(31),X(38)))+1))
5338 T(818)=0*10000*(1/(10000*(1-TBM(X(31),X(39)))+1))
5339 T(819)=5*10000*(1/(10000*(1-TBM(X(31),X(40)))+1))
5401 T(821)=0*10000*(1/(10000*(1-TBM(X(32),X(33)))+1))
5402 T(822)=0*10000*(1/(10000*(1-TBM(X(32),X(34)))+1))
5403 T(823)=1*10000*(1/(10000*(1-TBM(X(32),X(35)))+1))
5404 T(824)=2*10000*(1/(10000*(1-TBM(X(32),X(36)))+1))
5405 T(825)=5*10000*(1/(10000*(1-TBM(X(32),X(37)))+1))
5406 T(826)=10*10000*(1/(10000*(1-TBM(X(32),X(38)))+1))
5407 T(827)=1*10000*(1/(10000*(1-TBM(X(32),X(39)))+1))
5408 T(828)=6*10000*(1/(10000*(1-TBM(X(32),X(40)))+1))
5511 T(831)=2*10000*(1/(10000*(1-TBM(X(33),X(34)))+1))
5512 T(832)=2*10000*(1/(10000*(1-TBM(X(33),X(35)))+1))
5513 T(833)=10*10000*(1/(10000*(1-TBM(X(33),X(36)))+1))
5514 T(834)=2*10000*(1/(10000*(1-TBM(X(33),X(37)))+1))
5515 T(835)=10*10000*(1/(10000*(1-TBM(X(33),X(38)))+1))
5516 T(836)=2*10000*(1/(10000*(1-TBM(X(33),X(39)))+1))
5517 T(837)=4*10000*(1/(10000*(1-TBM(X(33),X(40)))+1))
5601 T(841)=6*10000*(1/(10000*(1-TBM(X(34),X(35)))+1))
5602 T(842)=0*10000*(1/(10000*(1-TBM(X(34),X(36)))+1))
5603 T(843)=0*10000*(1/(10000*(1-TBM(X(34),X(37)))+1))
5604 T(844)=1*10000*(1/(10000*(1-TBM(X(34),X(38)))+1))
5605 T(845)=4*10000*(1/(10000*(1-TBM(X(34),X(39)))+1))
5606 T(846)=4*10000*(1/(10000*(1-TBM(X(34),X(40)))+1))
5701 T(851)=2*10000*(1/(10000*(1-TBM(X(35),X(36)))+1))
5702 T(852)=5*10000*(1/(10000*(1-TBM(X(35),X(37)))+1))
5703 T(853)=5*10000*(1/(10000*(1-TBM(X(35),X(38)))+1))
5704 T(854)=2*10000*(1/(10000*(1-TBM(X(35),X(39)))+1))
5705 T(855)=1*10000*(1/(10000*(1-TBM(X(35),X(40)))+1))
5751 T(861)=0*10000*(1/(10000*(1-TBM(X(36),X(37)))+1))
5752 T(862)=0*10000*(1/(10000*(1-TBM(X(36),X(38)))+1))
5753 T(863)=10*10000*(1/(10000*(1-TBM(X(36),X(39)))+1))
5754 T(864)=2*10000*(1/(10000*(1-TBM(X(36),X(40)))+1))
5781 T(865)=6*10000*(1/(10000*(1-TBM(X(37),X(38)))+1))
5782 T(866)=2*10000*(1/(10000*(1-TBM(X(37),X(39)))+1))
5783 T(867)=2*10000*(1/(10000*(1-TBM(X(37),X(40)))+1))
5801 T(871)=6*10000*(1/(10000*(1-TBM(X(38),X(39)))+1))
5802 T(872)=0*10000*(1/(10000*(1-TBM(X(38),X(40)))+1))
5851 T(881)=5*10000*(1/(10000*(1-TBM(X(39),X(40)))+1))
5901 T(901)=2*10000*(1/(10000*(1-TBM(X(1),X(41)))+1))
5902 T(902)=10*10000*(1/(10000*(1-TBM(X(1),X(42)))+1))
5903 T(903)=5*10000*(1/(10000*(1-TBM(X(2),X(41)))+1))
5904 T(904)=4*10000*(1/(10000*(1-TBM(X(2),X(42)))+1))
5905 T(905)=1*10000*(1/(10000*(1-TBM(X(3),X(41)))+1))
5906 T(906)=5*10000*(1/(10000*(1-TBM(X(3),X(42)))+1))
5907 T(907)=6*10000*(1/(10000*(1-TBM(X(4),X(41)))+1))
5908 T(908)=5*10000*(1/(10000*(1-TBM(X(4),X(42)))+1))
5909 T(909)=0*10000*(1/(10000*(1-TBM(X(5),X(41)))+1))
5910 T(910)=2*10000*(1/(10000*(1-TBM(X(5),X(42)))+1))
5911 T(911)=0*10000*(1/(10000*(1-TBM(X(6),X(41)))+1))
5912 T(912)=0*10000*(1/(10000*(1-TBM(X(6),X(42)))+1))
5913 T(913)=1*10000*(1/(10000*(1-TBM(X(7),X(41)))+1))
5914 T(914)=2*10000*(1/(10000*(1-TBM(X(7),X(42)))+1))
5915 T(915)=5*10000*(1/(10000*(1-TBM(X(8),X(41)))+1))
5916 T(916)=0*10000*(1/(10000*(1-TBM(X(8),X(42)))+1))
5917 T(917)=5*10000*(1/(10000*(1-TBM(X(9),X(41)))+1))
5918 T(918)=0*10000*(1/(10000*(1-TBM(X(9),X(42)))+1))
5919 T(919)=2*10000*(1/(10000*(1-TBM(X(10),X(41)))+1))
5920 T(920)=5*10000*(1/(10000*(1-TBM(X(10),X(42)))+1))
5921 T(921)=5*10000*(1/(10000*(1-TBM(X(11),X(41)))+1))
5922 T(922)=0*10000*(1/(10000*(1-TBM(X(11),X(42)))+1))
5923 T(923)=10*10000*(1/(10000*(1-TBM(X(12),X(41)))+1))
5924 T(924)=0*10000*(1/(10000*(1-TBM(X(12),X(42)))+1))
5925 T(925)=1*10000*(1/(10000*(1-TBM(X(13),X(41)))+1))
5926 T(926)=4*10000*(1/(10000*(1-TBM(X(13),X(42)))+1))
5927 T(927)=2*10000*(1/(10000*(1-TBM(X(14),X(41)))+1))
5928 T(928)=0*10000*(1/(10000*(1-TBM(X(14),X(42)))+1))
5929 T(929)=0*10000*(1/(10000*(1-TBM(X(15),X(41)))+1))
5930 T(930)=0*10000*(1/(10000*(1-TBM(X(15),X(42)))+1))
5931 T(931)=2*10000*(1/(10000*(1-TBM(X(16),X(41)))+1))
5932 T(932)=0*10000*(1/(10000*(1-TBM(X(16),X(42)))+1))
5933 T(933)=5*10000*(1/(10000*(1-TBM(X(17),X(41)))+1))
5934 T(934)=5*10000*(1/(10000*(1-TBM(X(17),X(42)))+1))
5935 T(935)=0*10000*(1/(10000*(1-TBM(X(18),X(41)))+1))
5936 T(936)=10*10000*(1/(10000*(1-TBM(X(18),X(42)))+1))
5937 T(937)=10*10000*(1/(10000*(1-TBM(X(19),X(41)))+1))
5938 T(938)=0*10000*(1/(10000*(1-TBM(X(19),X(42)))+1))
5939 T(939)=0*10000*(1/(10000*(1-TBM(X(20),X(41)))+1))
5940 T(940)=0*10000*(1/(10000*(1-TBM(X(20),X(42)))+1))
5941 T(941)=2*10000*(1/(10000*(1-TBM(X(21),X(41)))+1))
5942 T(942)=1*10000*(1/(10000*(1-TBM(X(21),X(42)))+1))
5943 T(943)=0*10000*(1/(10000*(1-TBM(X(22),X(41)))+1))
5944 T(944)=5*10000*(1/(10000*(1-TBM(X(22),X(42)))+1))
5945 T(945)=1*10000*(1/(10000*(1-TBM(X(23),X(41)))+1))
5946 T(946)=0*10000*(1/(10000*(1-TBM(X(23),X(42)))+1))
5947 T(947)=0*10000*(1/(10000*(1-TBM(X(24),X(41)))+1))
5948 T(948)=0*10000*(1/(10000*(1-TBM(X(24),X(42)))+1))
5949 T(949)=2*10000*(1/(10000*(1-TBM(X(25),X(41)))+1))
5950 T(950)=1*10000*(1/(10000*(1-TBM(X(25),X(42)))+1))
5951 T(951)=0*10000*(1/(10000*(1-TBM(X(26),X(41)))+1))
5952 T(952)=5*10000*(1/(10000*(1-TBM(X(26),X(42)))+1))
5953 T(953)=10*10000*(1/(10000*(1-TBM(X(27),X(41)))+1))
5954 T(954)=0*10000*(1/(10000*(1-TBM(X(27),X(42)))+1))
5955 T(955)=1*10000*(1/(10000*(1-TBM(X(28),X(41)))+1))
5956 T(956)=0*10000*(1/(10000*(1-TBM(X(28),X(42)))+1))
5957 T(957)=10*10000*(1/(10000*(1-TBM(X(29),X(41)))+1))
5958 T(958)=10*10000*(1/(10000*(1-TBM(X(29),X(42)))+1))
5959 T(959)=0*10000*(1/(10000*(1-TBM(X(30),X(41)))+1))
5960 T(960)=2*10000*(1/(10000*(1-TBM(X(30),X(42)))+1))
5961 T(961)=1*10000*(1/(10000*(1-TBM(X(31),X(41)))+1))
5962 T(962)=1*10000*(1/(10000*(1-TBM(X(31),X(42)))+1))
5963 T(963)=0*10000*(1/(10000*(1-TBM(X(32),X(41)))+1))
5964 T(964)=0*10000*(1/(10000*(1-TBM(X(32),X(42)))+1))
5965 T(965)=3*10000*(1/(10000*(1-TBM(X(33),X(41)))+1))
5966 T(966)=5*10000*(1/(10000*(1-TBM(X(33),X(42)))+1))
5967 T(967)=10*10000*(1/(10000*(1-TBM(X(34),X(41)))+1))
5968 T(968)=0*10000*(1/(10000*(1-TBM(X(34),X(42)))+1))
5969 T(969)=5*10000*(1/(10000*(1-TBM(X(35),X(41)))+1))
5970 T(970)=0*10000*(1/(10000*(1-TBM(X(35),X(42)))+1))
5971 T(971)=2*10000*(1/(10000*(1-TBM(X(36),X(41)))+1))
5972 T(972)=1*10000*(1/(10000*(1-TBM(X(36),X(42)))+1))
5973 T(973)=1*10000*(1/(10000*(1-TBM(X(37),X(41)))+1))
5974 T(974)=3*10000*(1/(10000*(1-TBM(X(37),X(42)))+1))
5975 T(975)=10*10000*(1/(10000*(1-TBM(X(38),X(41)))+1))
5976 T(976)=5*10000*(1/(10000*(1-TBM(X(38),X(42)))+1))
5977 T(977)=0*10000*(1/(10000*(1-TBM(X(39),X(41)))+1))
5978 T(978)=10*10000*(1/(10000*(1-TBM(X(39),X(42)))+1))
5979 T(979)=0*10000*(1/(10000*(1-TBM(X(40),X(41)))+1))
5980 T(980)=1*10000*(1/(10000*(1-TBM(X(40),X(42)))+1))
5981 T(981)=4*10000*(1/(10000*(1-TBM(X(41),X(42)))+1))
6151 P1NEW=0
6152 FOR KAU7=1 TO 981
6153 P1NEW=P1NEW+T(KAU7)
6154 NEXT KAU7
6450 P=-P1NEW
6451 IF P<=M THEN 6670
6657 FOR KEW=1 TO 42
6658 A(KEW)=X(KEW)
6659 NEXT KEW
6661 M=P
6662 PP1=P1NEW
6669 GOTO 128
6670 NEXT I
6890 IF M>-2700000! THEN 6941 ELSE 6999
6941 PRINT A(1),A(2),A(3),A(4),A(5),A(6),A(7),A(8),A(9),A(10)
6942 PRINT A(11),A(12),A(13),A(14),A(15),A(16),A(17),A(18),A(19),A(20)
6943 PRINT A(21),A(22),A(23),A(24),A(25),A(26),A(27),A(28),A(29),A(30)
6944 PRINT A(31),A(32),A(33),A(34),A(35),A(36),A(37),A(38),A(39),A(40)
6945 PRINT A(41),A(42)
6951 PRINT JJJJ,M
6999 NEXT JJJJ

This BASIC computer program was run with the IBM basica/D interpreter, and its output through JJJJ=-21195 is presented below.

3 3 1 1 4
1 1 4 3 1
3 4 3 3 4
3 3 4 3 3
4 3 3 5 6
4 5 5 4 4
5 4 4 5 5
5 6 5 6 6
6 5
-30974 -2672059.794020599

1 3 1 3 3
4 4 1 3 4
3 3 1 1 4
3 3 4 3 3
4 3 4 5 4
5 6 5 6 5
5 5 5 5 6
4 4 4 6 6
5 5
-28948 -2692057.794220578

4 3 4 3 3
1 1 1 3 1
3 3 4 1 4
3 3 3 3 3
4 3 4 5 6
4 4 5 4 5
5 5 4 5 5
6 6 6 4 6
5 5
-24437 -2602066.793320668

1 3 1 3 3
4 4 1 3 4
3 3 1 1 4
3 3 3 3 3
4 3 4 6 4
5 6 5 6 5
5 5 5 5 5
4 4 6 5 6
4 4
-22969 -2652061.793820618

3 1 4 4 3
1 1 4 3 1
3 4 3 1 4
3 3 4 3 3
4 3 3 5 6
4 5 5 4 4
5 4 4 5 5
6 6 5 6 5
6 5
-21897 -2652061.793820619

4 3 1 3 3
1 1 4 3 1
3 3 4 1 4
3 3 3 3 3
4 3 4 5 6
4 6 5 6 4
5 4 4 5 5
5 6 5 4 6
4 5
-21195 -2582068.793120688

Copied manually from the computer monitor and interpreted in accordance with line 6941 through 6951, the output above was produced in 22 hours on a personal computer with an Intel 2.66 GHz. chip and the IBM interpreter, which is slower than the corresponding compiler.

References

[1] L. R. Arnold, R. E. Beckwith, C. M. Jones, Scheduling the 41st-ORSA-meeting sessions: the visiting-fireman problem, II, Operations Research 21 (1973) 1095-1103.

[2] R. C. Carlson, G. L. Nemhauser, Scheduling to minimize interaction cost, Operations Research 14 (1966) 52-58.

[3] F. E. Davis, M. D. Devine, R. P. Lutz, Scheduling activities among conflicting facilities to minimize conflict, Mathematical Programming 6 (1974) 224-228.

[4] R. P. Lutz, M. D. Devine, H. J. Kumin, W. C. Smith, An application of operations research to school desegregation, Management Science 19 (1972) P100-P109.

[5] C. E. Nugent, T. E. Vollmann, J. Ruml, An experimental comparison of techniques for the assignment of facilities to locations, Operations Research 16 (1968) 150-173.

[6] J. Skorin-Kapov, Tabu Search Applied to the Quadratic Assignment Problem, Journal on Computing 2 (1990) 33-45.





Assigning Forty-Two Activities to Conflicting Facilities Subject to Less-Than-or-Equal-to Constraints

Jsun Yui Wong

The scheduling problems of Carlson and Nemhauser [2], Davis, Devine, and Lutz [3], and Arnold, Beckwith, and Jones [1] are related. The difference between the original Carlson and Nemhauser problem [2] and the problem of Davis, Devine, and Lutz [3] is "that their problem did not consider conflicting facilities," Davis, Devine, and Lutz [3, p. 225]. For the newer problem, "a cost is incurred if two activities are assigned to the same facility, or to two facilities whose operation is in conflict," Davis et al. [3, p. 224]. Line 32 and line 33 of the computer program below say the operation of facility 1 and facility 2 is in conflict; line 33 and line 34 say the operation of facility 2 and facility 3 is in conflict.

There is the constraint that period one can handle only five courses or less, and there is also the constraint that period six can handle only five courses or less. These two restrictions exemplify the less-than-or-equal-to constraints on pages 1096 of Arnold, Beckwith, and Jones [1] and are handled mainly by line 366 through line 400 of the program listed below. Carlson and Nemhausser [2] "seek to find a conflict-minimizing schedule, but do not require that the schedule be balanced," [1, p. 1096].

Like Lutz, Devine, Kumin, and Smith [4], the present paper has another type of restrictions. Courses 1 through 14 may be assigned to time blocks 1 through 4, courses 15 through 23 may be assigned in time blocks 2 through 4, and courses 24 through 42 may be assigned to time blocks 4 through 6. These restrictions are handled by line 341 through line 363.

The conflict costs used in line 401 through line 5981 of the computer program below come from Nugent, Vollmann, and Ruml [5] and Skorin-Kapov [6].

In this paper, A(1)=2, A(2)=2, A(3)=2,..., and A(42)=2 stand for course 1 in time block 2, course 2 in time block 2, course 3 in time block 2,..., and course 42 in time block 2, respectively.

The computer program below is modelled after the computer program of the post "Assigning Activities to Conflicting Facilities To Minimize Conflict Cost" posted in the present blog on March 7, 2008.

4 DEFINT A-Z
5 DEFDBL M,T,P
7 DIM A(56),X(56),T(1333)
9 DIM TBM(6,6)
32 TBM(1,1)=1:TBM(1,2)=1
33 TBM(2,1)=1:TBM(2,2)=1:TBM(2,3)=1
34 TBM(3,2)=1:TBM(3,3)=1
35 TBM(4,4)=1
36 TBM(5,5)=1
37 TBM(6,6)=1
42 FOR JJJJ=-32000 TO 32000
44 RANDOMIZE JJJJ
46 M=-1D+17
61 FOR KX=1 TO 14
63 A(KX)=2+FIX(RND*3)
65 NEXT KX
71 FOR KY=15 TO 23
73 A(KY)=2+FIX(RND*3)
75 NEXT KY
81 FOR KZ=24 TO 42
83 A(KZ)=4+FIX(RND*2)
85 NEXT KZ
111 IMAR=10+FIX(RND*1500)
128 FOR I=1 TO IMAR
129 FOR K=1 TO 42
301 X(K)=A(K)
302 NEXT K
311 IF RND<.33 THEN 341 ELSE IF RND<.32 THEN 351 ELSE 361
341 IAP1=1+FIX(RND*14)
343 X(IAP1)=1+FIX(RND*4)
346 GOTO 366
351 IAP2=15+FIX(RND*9)
353 X(IAP2)=2+FIX(RND*3)
356 GOTO 366
361 IAP3=24+FIX(RND*19)
363 X(IAP3)=4+FIX(RND*3)
366 SUM1=0:SUM2=0
371 FOR I11=1 TO 42
381 IF X(I11)=1 THEN 391 ELSE IF X(I11)=6 THEN 397 ELSE 398
391 SUM1=SUM1+1
393 GOTO 398
397 SUM2=SUM2+6
398 NEXT I11
399 IF SUM1>5 THEN 6670
400 IF SUM2>30 THEN 6670
401 T(1)=3*10000*(1/(10000*(1-TBM(X(1),X(2)))+1))
402 T(2)=2*10000*(1/(10000*(1-TBM(X(1),X(3)))+1))
403 T(3)=0*10000*(1/(10000*(1-TBM(X(1),X(4)))+1))
404 T(4)=0*10000*(1/(10000*(1-TBM(X(1),X(5)))+1))
405 T(5)=2*10000*(1/(10000*(1-TBM(X(1),X(6)))+1))
406 T(6)=10*10000*(1/(10000*(1-TBM(X(1),X(7)))+1))
407 T(7)=5*10000*(1/(10000*(1-TBM(X(1),X(8)))+1))
408 T(8)=0*10000*(1/(10000*(1-TBM(X(1),X(9)))+1))
409 T(9)=5*10000*(1/(10000*(1-TBM(X(1),X(10)))+1))
410 T(10)=2*10000*(1/(10000*(1-TBM(X(1),X(11)))+1))
411 T(11)=5*10000*(1/(10000*(1-TBM(X(1),X(12)))+1))
412 T(12)=0*10000*(1/(10000*(1-TBM(X(1),X(13)))+1))
413 T(13)=0*10000*(1/(10000*(1-TBM(X(1),X(14)))+1))
414 T(14)=2*10000*(1/(10000*(1-TBM(X(1),X(15)))+1))
415 T(15)=0*10000*(1/(10000*(1-TBM(X(1),X(16)))+1))
416 T(16)=5*10000*(1/(10000*(1-TBM(X(1),X(17)))+1))
417 T(17)=6*10000*(1/(10000*(1-TBM(X(1),X(18)))+1))
418 T(18)=3*10000*(1/(10000*(1-TBM(X(1),X(19)))+1))
419 T(19)=0*10000*(1/(10000*(1-TBM(X(1),X(20)))+1))
420 T(20)=1*10000*(1/(10000*(1-TBM(X(1),X(21)))+1))
421 T(21)=10*10000*(1/(10000*(1-TBM(X(1),X(22)))+1))
422 T(22)=0*10000*(1/(10000*(1-TBM(X(1),X(23)))+1))
423 T(23)=10*10000*(1/(10000*(1-TBM(X(1),X(24)))+1))
424 T(24)=2*10000*(1/(10000*(1-TBM(X(1),X(25)))+1))
425 T(25)=1*10000*(1/(10000*(1-TBM(X(1),X(26)))+1))
426 T(26)=1*10000*(1/(10000*(1-TBM(X(1),X(27)))+1))
427 T(27)=1*10000*(1/(10000*(1-TBM(X(1),X(28)))+1))
428 T(28)=0*10000*(1/(10000*(1-TBM(X(1),X(29)))+1))
429 T(29)=1*10000*(1/(10000*(1-TBM(X(1),X(30)))+1))
430 T(30)=4*10000*(1/(10000*(1-TBM(X(2),X(3)))+1))
431 T(31)=0*10000*(1/(10000*(1-TBM(X(2),X(4)))+1))
432 T(32)=10*10000*(1/(10000*(1-TBM(X(2),X(5)))+1))
433 T(33)=4*10000*(1/(10000*(1-TBM(X(2),X(6)))+1))
434 T(34)=0*10000*(1/(10000*(1-TBM(X(2),X(7)))+1))
435 T(35)=0*10000*(1/(10000*(1-TBM(X(2),X(8)))+1))
436 T(36)=2*10000*(1/(10000*(1-TBM(X(2),X(9)))+1))
437 T(37)=2*10000*(1/(10000*(1-TBM(X(2),X(10)))+1))
438 T(38)=1*10000*(1/(10000*(1-TBM(X(2),X(11)))+1))
439 T(39)=0*10000*(1/(10000*(1-TBM(X(2),X(12)))+1))
440 T(40)=5*10000*(1/(10000*(1-TBM(X(2),X(13)))+1))
441 T(41)=0*10000*(1/(10000*(1-TBM(X(2),X(14)))+1))
442 T(42)=0*10000*(1/(10000*(1-TBM(X(2),X(15)))+1))
443 T(43)=0*10000*(1/(10000*(1-TBM(X(2),X(16)))+1))
444 T(44)=0*10000*(1/(10000*(1-TBM(X(2),X(17)))+1))
445 T(45)=2*10000*(1/(10000*(1-TBM(X(2),X(18)))+1))
446 T(46)=0*10000*(1/(10000*(1-TBM(X(2),X(19)))+1))
447 T(47)=1*10000*(1/(10000*(1-TBM(X(2),X(20)))+1))
448 T(48)=6*10000*(1/(10000*(1-TBM(X(2),X(21)))+1))
449 T(49)=1*10000*(1/(10000*(1-TBM(X(2),X(22)))+1))
450 T(50)=0*10000*(1/(10000*(1-TBM(X(2),X(23)))+1))
451 T(51)=1*10000*(1/(10000*(1-TBM(X(2),X(24)))+1))
452 T(52)=2*10000*(1/(10000*(1-TBM(X(2),X(25)))+1))
453 T(53)=2*10000*(1/(10000*(1-TBM(X(2),X(26)))+1))
454 T(54)=5*10000*(1/(10000*(1-TBM(X(2),X(27)))+1))
455 T(55)=1*10000*(1/(10000*(1-TBM(X(2),X(28)))+1))
456 T(56)=10*10000*(1/(10000*(1-TBM(X(2),X(29)))+1))
457 T(57)=5*10000*(1/(10000*(1-TBM(X(2),X(30)))+1))
458 T(58)=3*10000*(1/(10000*(1-TBM(X(3),X(4)))+1))
459 T(59)=4*10000*(1/(10000*(1-TBM(X(3),X(5)))+1))
460 T(60)=0*10000*(1/(10000*(1-TBM(X(3),X(6)))+1))
461 T(61)=5*10000*(1/(10000*(1-TBM(X(3),X(7)))+1))
462 T(62)=5*10000*(1/(10000*(1-TBM(X(3),X(8)))+1))
463 T(63)=5*10000*(1/(10000*(1-TBM(X(3),X(9)))+1))
464 T(64)=1*10000*(1/(10000*(1-TBM(X(3),X(10)))+1))
465 T(65)=4*10000*(1/(10000*(1-TBM(X(3),X(11)))+1))
466 T(66)=1*10000*(1/(10000*(1-TBM(X(3),X(12)))+1))
467 T(67)=0*10000*(1/(10000*(1-TBM(X(3),X(13)))+1))
468 T(68)=4*10000*(1/(10000*(1-TBM(X(3),X(14)))+1))
469 T(69)=0*10000*(1/(10000*(1-TBM(X(3),X(15)))+1))
470 T(70)=4*10000*(1/(10000*(1-TBM(X(3),X(16)))+1))
471 T(71)=0*10000*(1/(10000*(1-TBM(X(3),X(17)))+1))
472 T(72)=6*10000*(1/(10000*(1-TBM(X(3),X(18)))+1))
473 T(73)=3*10000*(1/(10000*(1-TBM(X(3),X(19)))+1))
474 T(74)=2*10000*(1/(10000*(1-TBM(X(3),X(20)))+1))
475 T(75)=5*10000*(1/(10000*(1-TBM(X(3),X(21)))+1))
476 T(76)=5*10000*(1/(10000*(1-TBM(X(3),X(22)))+1))
477 T(77)=2*10000*(1/(10000*(1-TBM(X(3),X(23)))+1))
478 T(78)=1*10000*(1/(10000*(1-TBM(X(3),X(24)))+1))
479 T(79)=0*10000*(1/(10000*(1-TBM(X(3),X(25)))+1))
480 T(80)=0*10000*(1/(10000*(1-TBM(X(3),X(26)))+1))
481 T(81)=3*10000*(1/(10000*(1-TBM(X(3),X(27)))+1))
482 T(82)=1*10000*(1/(10000*(1-TBM(X(3),X(28)))+1))
483 T(83)=0*10000*(1/(10000*(1-TBM(X(3),X(29)))+1))
484 T(84)=2*10000*(1/(10000*(1-TBM(X(3),X(30)))+1))
485 T(85)=0*10000*(1/(10000*(1-TBM(X(4),X(5)))+1))
486 T(86)=0*10000*(1/(10000*(1-TBM(X(4),X(6)))+1))
487 T(87)=0*10000*(1/(10000*(1-TBM(X(4),X(7)))+1))
488 T(88)=2*10000*(1/(10000*(1-TBM(X(4),X(8)))+1))
489 T(89)=2*10000*(1/(10000*(1-TBM(X(4),X(9)))+1))
490 T(90)=0*10000*(1/(10000*(1-TBM(X(4),X(10)))+1))
491 T(91)=6*10000*(1/(10000*(1-TBM(X(4),X(11)))+1))
492 T(92)=0*10000*(1/(10000*(1-TBM(X(4),X(12)))+1))
493 T(93)=2*10000*(1/(10000*(1-TBM(X(4),X(13)))+1))
494 T(94)=5*10000*(1/(10000*(1-TBM(X(4),X(14)))+1))
495 T(95)=2*10000*(1/(10000*(1-TBM(X(4),X(15)))+1))
496 T(96)=5*10000*(1/(10000*(1-TBM(X(4),X(16)))+1))
497 T(97)=1*10000*(1/(10000*(1-TBM(X(4),X(17)))+1))
498 T(98)=1*10000*(1/(10000*(1-TBM(X(4),X(18)))+1))
499 T(99)=1*10000*(1/(10000*(1-TBM(X(4),X(19)))+1))
500 T(100)=1*10000*(1/(10000*(1-TBM(X(4),X(20)))+1))
501 T(101)=2*10000*(1/(10000*(1-TBM(X(4),X(21)))+1))
502 T(102)=2*10000*(1/(10000*(1-TBM(X(4),X(22)))+1))
503 T(103)=4*10000*(1/(10000*(1-TBM(X(4),X(23)))+1))
504 T(104)=0*10000*(1/(10000*(1-TBM(X(4),X(24)))+1))
505 T(105)=2*10000*(1/(10000*(1-TBM(X(4),X(25)))+1))
506 T(106)=0*10000*(1/(10000*(1-TBM(X(4),X(26)))+1))
507 T(107)=2*10000*(1/(10000*(1-TBM(X(4),X(27)))+1))
508 T(108)=2*10000*(1/(10000*(1-TBM(X(4),X(28)))+1))
509 T(109)=5*10000*(1/(10000*(1-TBM(X(4),X(29)))+1))
510 T(110)=5*10000*(1/(10000*(1-TBM(X(4),X(30)))+1))
511 T(111)=5*10000*(1/(10000*(1-TBM(X(5),X(6)))+1))
512 T(112)=2*10000*(1/(10000*(1-TBM(X(5),X(7)))+1))
513 T(113)=0*10000*(1/(10000*(1-TBM(X(5),X(8)))+1))
514 T(114)=0*10000*(1/(10000*(1-TBM(X(5),X(9)))+1))
515 T(115)=0*10000*(1/(10000*(1-TBM(X(5),X(10)))+1))
516 T(116)=0*10000*(1/(10000*(1-TBM(X(5),X(11)))+1))
517 T(117)=2*10000*(1/(10000*(1-TBM(X(5),X(12)))+1))
518 T(118)=0*10000*(1/(10000*(1-TBM(X(5),X(13)))+1))
519 T(119)=0*10000*(1/(10000*(1-TBM(X(5),X(14)))+1))
520 T(120)=0*10000*(1/(10000*(1-TBM(X(5),X(15)))+1))
521 T(121)=0*10000*(1/(10000*(1-TBM(X(5),X(16)))+1))
522 T(122)=2*10000*(1/(10000*(1-TBM(X(5),X(17)))+1))
523 T(123)=1*10000*(1/(10000*(1-TBM(X(5),X(18)))+1))
524 T(124)=0*10000*(1/(10000*(1-TBM(X(5),X(19)))+1))
525 T(125)=0*10000*(1/(10000*(1-TBM(X(5),X(20)))+1))
526 T(126)=2*10000*(1/(10000*(1-TBM(X(5),X(21)))+1))
527 T(127)=0*10000*(1/(10000*(1-TBM(X(5),X(22)))+1))
528 T(128)=5*10000*(1/(10000*(1-TBM(X(5),X(23)))+1))
529 T(129)=1*10000*(1/(10000*(1-TBM(X(5),X(24)))+1))
530 T(130)=0*10000*(1/(10000*(1-TBM(X(5),X(25)))+1))
531 T(131)=2*10000*(1/(10000*(1-TBM(X(5),X(26)))+1))
532 T(132)=1*10000*(1/(10000*(1-TBM(X(5),X(27)))+1))
533 T(133)=0*10000*(1/(10000*(1-TBM(X(5),X(28)))+1))
534 T(134)=2*10000*(1/(10000*(1-TBM(X(5),X(29)))+1))
535 T(135)=1*10000*(1/(10000*(1-TBM(X(5),X(30)))+1))
536 T(136)=1*10000*(1/(10000*(1-TBM(X(6),X(7)))+1))
537 T(137)=2*10000*(1/(10000*(1-TBM(X(6),X(8)))+1))
538 T(138)=2*10000*(1/(10000*(1-TBM(X(6),X(9)))+1))
539 T(139)=1*10000*(1/(10000*(1-TBM(X(6),X(10)))+1))
540 T(140)=4*10000*(1/(10000*(1-TBM(X(6),X(11)))+1))
541 T(141)=10*10000*(1/(10000*(1-TBM(X(6),X(12)))+1))
542 T(142)=10*10000*(1/(10000*(1-TBM(X(6),X(13)))+1))
543 T(143)=2*10000*(1/(10000*(1-TBM(X(6),X(14)))+1))
544 T(144)=5*10000*(1/(10000*(1-TBM(X(6),X(15)))+1))
545 T(145)=5*10000*(1/(10000*(1-TBM(X(6),X(16)))+1))
546 T(146)=0*10000*(1/(10000*(1-TBM(X(6),X(17)))+1))
547 T(147)=5*10000*(1/(10000*(1-TBM(X(6),X(18)))+1))
548 T(148)=0*10000*(1/(10000*(1-TBM(X(6),X(19)))+1))
549 T(149)=0*10000*(1/(10000*(1-TBM(X(6),X(20)))+1))
550 T(150)=0*10000*(1/(10000*(1-TBM(X(6),X(21)))+1))
551 T(151)=10*10000*(1/(10000*(1-TBM(X(6),X(22)))+1))
552 T(152)=0*10000*(1/(10000*(1-TBM(X(6),X(23)))+1))
553 T(153)=0*10000*(1/(10000*(1-TBM(X(6),X(24)))+1))
554 T(154)=0*10000*(1/(10000*(1-TBM(X(6),X(25)))+1))
555 T(155)=4*10000*(1/(10000*(1-TBM(X(6),X(26)))+1))
556 T(156)=0*10000*(1/(10000*(1-TBM(X(6),X(27)))+1))
557 T(157)=10*10000*(1/(10000*(1-TBM(X(6),X(28)))+1))
558 T(158)=1*10000*(1/(10000*(1-TBM(X(6),X(29)))+1))
559 T(159)=1*10000*(1/(10000*(1-TBM(X(6),X(30)))+1))
560 T(160)=10*10000*(1/(10000*(1-TBM(X(7),X(8)))+1))
561 T(161)=10*10000*(1/(10000*(1-TBM(X(7),X(9)))+1))
562 T(162)=5*10000*(1/(10000*(1-TBM(X(7),X(10)))+1))
563 T(163)=10*10000*(1/(10000*(1-TBM(X(7),X(11)))+1))
564 T(164)=10*10000*(1/(10000*(1-TBM(X(7),X(12)))+1))
565 T(165)=6*10000*(1/(10000*(1-TBM(X(7),X(13)))+1))
566 T(166)=0*10000*(1/(10000*(1-TBM(X(7),X(14)))+1))
567 T(167)=0*10000*(1/(10000*(1-TBM(X(7),X(15)))+1))
568 T(168)=10*10000*(1/(10000*(1-TBM(X(7),X(16)))+1))
569 T(169)=2*10000*(1/(10000*(1-TBM(X(7),X(17)))+1))
570 T(170)=1*10000*(1/(10000*(1-TBM(X(7),X(18)))+1))
571 T(171)=10*10000*(1/(10000*(1-TBM(X(7),X(19)))+1))
572 T(172)=1*10000*(1/(10000*(1-TBM(X(7),X(20)))+1))
573 T(173)=5*10000*(1/(10000*(1-TBM(X(7),X(21)))+1))
574 T(174)=5*10000*(1/(10000*(1-TBM(X(7),X(22)))+1))
575 T(175)=2*10000*(1/(10000*(1-TBM(X(7),X(23)))+1))
576 T(176)=3*10000*(1/(10000*(1-TBM(X(7),X(24)))+1))
577 T(177)=5*10000*(1/(10000*(1-TBM(X(7),X(25)))+1))
578 T(178)=0*10000*(1/(10000*(1-TBM(X(7),X(26)))+1))
579 T(179)=2*10000*(1/(10000*(1-TBM(X(7),X(27)))+1))
580 T(180)=0*10000*(1/(10000*(1-TBM(X(7),X(28)))+1))
581 T(181)=1*10000*(1/(10000*(1-TBM(X(7),X(29)))+1))
582 T(182)=3*10000*(1/(10000*(1-TBM(X(7),X(30)))+1))
583 T(183)=1*10000*(1/(10000*(1-TBM(X(8),X(9)))+1))
584 T(184)=3*10000*(1/(10000*(1-TBM(X(8),X(10)))+1))
585 T(185)=5*10000*(1/(10000*(1-TBM(X(8),X(11)))+1))
586 T(186)=0*10000*(1/(10000*(1-TBM(X(8),X(12)))+1))
587 T(187)=0*10000*(1/(10000*(1-TBM(X(8),X(13)))+1))
588 T(188)=0*10000*(1/(10000*(1-TBM(X(8),X(14)))+1))
589 T(189)=2*10000*(1/(10000*(1-TBM(X(8),X(15)))+1))
590 T(190)=4*10000*(1/(10000*(1-TBM(X(8),X(16)))+1))
591 T(191)=5*10000*(1/(10000*(1-TBM(X(8),X(17)))+1))
592 T(192)=2*10000*(1/(10000*(1-TBM(X(8),X(18)))+1))
593 T(193)=10*10000*(1/(10000*(1-TBM(X(8),X(19)))+1))
594 T(194)=6*10000*(1/(10000*(1-TBM(X(8),X(20)))+1))
595 T(195)=0*10000*(1/(10000*(1-TBM(X(8),X(21)))+1))
596 T(196)=5*10000*(1/(10000*(1-TBM(X(8),X(22)))+1))
597 T(197)=5*10000*(1/(10000*(1-TBM(X(8),X(23)))+1))
598 T(198)=2*10000*(1/(10000*(1-TBM(X(8),X(24)))+1))
599 T(199)=5*10000*(1/(10000*(1-TBM(X(8),X(25)))+1))
600 T(200)=0*10000*(1/(10000*(1-TBM(X(8),X(26)))+1))
601 T(201)=5*10000*(1/(10000*(1-TBM(X(8),X(27)))+1))
602 T(202)=5*10000*(1/(10000*(1-TBM(X(8),X(28)))+1))
603 T(203)=0*10000*(1/(10000*(1-TBM(X(8),X(29)))+1))
604 T(204)=2*10000*(1/(10000*(1-TBM(X(8),X(30)))+1))
605 T(205)=10*10000*(1/(10000*(1-TBM(X(9),X(10)))+1))
606 T(206)=2*10000*(1/(10000*(1-TBM(X(9),X(11)))+1))
607 T(207)=1*10000*(1/(10000*(1-TBM(X(9),X(12)))+1))
608 T(208)=5*10000*(1/(10000*(1-TBM(X(9),X(13)))+1))
609 T(209)=2*10000*(1/(10000*(1-TBM(X(9),X(14)))+1))
610 T(210)=0*10000*(1/(10000*(1-TBM(X(9),X(15)))+1))
611 T(211)=3*10000*(1/(10000*(1-TBM(X(9),X(16)))+1))
612 T(212)=0*10000*(1/(10000*(1-TBM(X(9),X(17)))+1))
613 T(213)=2*10000*(1/(10000*(1-TBM(X(9),X(18)))+1))
614 T(214)=0*10000*(1/(10000*(1-TBM(X(9),X(19)))+1))
615 T(215)=0*10000*(1/(10000*(1-TBM(X(9),X(20)))+1))
616 T(216)=4*10000*(1/(10000*(1-TBM(X(9),X(21)))+1))
617 T(217)=0*10000*(1/(10000*(1-TBM(X(9),X(22)))+1))
618 T(218)=5*10000*(1/(10000*(1-TBM(X(9),X(23)))+1))
619 T(219)=2*10000*(1/(10000*(1-TBM(X(9),X(24)))+1))
620 T(220)=0*10000*(1/(10000*(1-TBM(X(9),X(25)))+1))
621 T(221)=5*10000*(1/(10000*(1-TBM(X(9),X(26)))+1))
622 T(222)=2*10000*(1/(10000*(1-TBM(X(9),X(27)))+1))
623 T(223)=2*10000*(1/(10000*(1-TBM(X(9),X(28)))+1))
624 T(224)=5*10000*(1/(10000*(1-TBM(X(9),X(29)))+1))
625 T(225)=2*10000*(1/(10000*(1-TBM(X(9),X(30)))+1))
626 T(226)=5*10000*(1/(10000*(1-TBM(X(10),X(11)))+1))
627 T(227)=5*10000*(1/(10000*(1-TBM(X(10),X(12)))+1))
628 T(228)=6*10000*(1/(10000*(1-TBM(X(10),X(13)))+1))
629 T(229)=0*10000*(1/(10000*(1-TBM(X(10),X(14)))+1))
630 T(230)=1*10000*(1/(10000*(1-TBM(X(10),X(15)))+1))
631 T(231)=5*10000*(1/(10000*(1-TBM(X(10),X(16)))+1))
632 T(232)=5*10000*(1/(10000*(1-TBM(X(10),X(17)))+1))
633 T(233)=0*10000*(1/(10000*(1-TBM(X(10),X(18)))+1))
634 T(234)=5*10000*(1/(10000*(1-TBM(X(10),X(19)))+1))
635 T(235)=2*10000*(1/(10000*(1-TBM(X(10),X(20)))+1))
636 T(236)=3*10000*(1/(10000*(1-TBM(X(10),X(21)))+1))
637 T(237)=5*10000*(1/(10000*(1-TBM(X(10),X(22)))+1))
638 T(238)=0*10000*(1/(10000*(1-TBM(X(10),X(23)))+1))
639 T(239)=5*10000*(1/(10000*(1-TBM(X(10),X(24)))+1))
640 T(240)=2*10000*(1/(10000*(1-TBM(X(10),X(25)))+1))
641 T(241)=10*10000*(1/(10000*(1-TBM(X(10),X(26)))+1))
642 T(242)=10*10000*(1/(10000*(1-TBM(X(10),X(27)))+1))
643 T(243)=1*10000*(1/(10000*(1-TBM(X(10),X(28)))+1))
644 T(244)=5*10000*(1/(10000*(1-TBM(X(10),X(29)))+1))
645 T(245)=2*10000*(1/(10000*(1-TBM(X(10),X(30)))+1))
646 T(246)=0*10000*(1/(10000*(1-TBM(X(11),X(12)))+1))
647 T(247)=0*10000*(1/(10000*(1-TBM(X(11),X(13)))+1))
648 T(248)=1*10000*(1/(10000*(1-TBM(X(11),X(14)))+1))
649 T(249)=2*10000*(1/(10000*(1-TBM(X(11),X(15)))+1))
650 T(250)=1*10000*(1/(10000*(1-TBM(X(11),X(16)))+1))
651 T(251)=0*10000*(1/(10000*(1-TBM(X(11),X(17)))+1))
652 T(252)=2*10000*(1/(10000*(1-TBM(X(11),X(18)))+1))
653 T(253)=0*10000*(1/(10000*(1-TBM(X(11),X(19)))+1))
654 T(254)=0*10000*(1/(10000*(1-TBM(X(11),X(20)))+1))
655 T(255)=0*10000*(1/(10000*(1-TBM(X(11),X(21)))+1))
656 T(256)=6*10000*(1/(10000*(1-TBM(X(11),X(22)))+1))
657 T(257)=6*10000*(1/(10000*(1-TBM(X(11),X(23)))+1))
658 T(258)=0*10000*(1/(10000*(1-TBM(X(11),X(24)))+1))
659 T(259)=4*10000*(1/(10000*(1-TBM(X(11),X(25)))+1))
660 T(260)=5*10000*(1/(10000*(1-TBM(X(11),X(26)))+1))
661 T(261)=3*10000*(1/(10000*(1-TBM(X(11),X(27)))+1))
662 T(262)=2*10000*(1/(10000*(1-TBM(X(11),X(28)))+1))
663 T(263)=2*10000*(1/(10000*(1-TBM(X(11),X(29)))+1))
664 T(264)=10*10000*(1/(10000*(1-TBM(X(11),X(30)))+1))
665 T(265)=5*10000*(1/(10000*(1-TBM(X(12),X(13)))+1))
666 T(266)=5*10000*(1/(10000*(1-TBM(X(12),X(14)))+1))
667 T(267)=2*10000*(1/(10000*(1-TBM(X(12),X(15)))+1))
668 T(268)=0*10000*(1/(10000*(1-TBM(X(12),X(16)))+1))
669 T(269)=0*10000*(1/(10000*(1-TBM(X(12),X(17)))+1))
670 T(270)=0*10000*(1/(10000*(1-TBM(X(12),X(18)))+1))
671 T(271)=0*10000*(1/(10000*(1-TBM(X(12),X(19)))+1))
672 T(272)=2*10000*(1/(10000*(1-TBM(X(12),X(20)))+1))
673 T(273)=0*10000*(1/(10000*(1-TBM(X(12),X(21)))+1))
674 T(274)=4*10000*(1/(10000*(1-TBM(X(12),X(22)))+1))
675 T(275)=5*10000*(1/(10000*(1-TBM(X(12),X(23)))+1))
676 T(276)=10*10000*(1/(10000*(1-TBM(X(12),X(24)))+1))
677 T(277)=1*10000*(1/(10000*(1-TBM(X(12),X(25)))+1))
678 T(278)=0*10000*(1/(10000*(1-TBM(X(12),X(26)))+1))
679 T(279)=0*10000*(1/(10000*(1-TBM(X(12),X(27)))+1))
680 T(280)=0*10000*(1/(10000*(1-TBM(X(12),X(28)))+1))
681 T(281)=0*10000*(1/(10000*(1-TBM(X(12),X(29)))+1))
682 T(282)=1*10000*(1/(10000*(1-TBM(X(12),X(30)))+1))
683 T(283)=2*10000*(1/(10000*(1-TBM(X(13),X(14)))+1))
684 T(284)=0*10000*(1/(10000*(1-TBM(X(13),X(15)))+1))
685 T(285)=4*10000*(1/(10000*(1-TBM(X(13),X(16)))+1))
686 T(286)=2*10000*(1/(10000*(1-TBM(X(13),X(17)))+1))
687 T(287)=2*10000*(1/(10000*(1-TBM(X(13),X(18)))+1))
688 T(288)=1*10000*(1/(10000*(1-TBM(X(13),X(19)))+1))
689 T(289)=0*10000*(1/(10000*(1-TBM(X(13),X(20)))+1))
690 T(290)=6*10000*(1/(10000*(1-TBM(X(13),X(21)))+1))
691 T(291)=2*10000*(1/(10000*(1-TBM(X(13),X(22)))+1))
692 T(292)=1*10000*(1/(10000*(1-TBM(X(13),X(23)))+1))
693 T(293)=5*10000*(1/(10000*(1-TBM(X(13),X(24)))+1))
694 T(294)=5*10000*(1/(10000*(1-TBM(X(13),X(25)))+1))
695 T(295)=0*10000*(1/(10000*(1-TBM(X(13),X(26)))+1))
696 T(296)=0*10000*(1/(10000*(1-TBM(X(13),X(27)))+1))
697 T(297)=1*10000*(1/(10000*(1-TBM(X(13),X(28)))+1))
698 T(298)=5*10000*(1/(10000*(1-TBM(X(13),X(29)))+1))
699 T(299)=5*10000*(1/(10000*(1-TBM(X(13),X(30)))+1))
700 T(300)=2*10000*(1/(10000*(1-TBM(X(14),X(15)))+1))
701 T(301)=1*10000*(1/(10000*(1-TBM(X(14),X(16)))+1))
702 T(302)=0*10000*(1/(10000*(1-TBM(X(14),X(17)))+1))
703 T(303)=5*10000*(1/(10000*(1-TBM(X(14),X(18)))+1))
704 T(304)=3*10000*(1/(10000*(1-TBM(X(14),X(19)))+1))
705 T(305)=10*10000*(1/(10000*(1-TBM(X(14),X(20)))+1))
706 T(306)=0*10000*(1/(10000*(1-TBM(X(14),X(21)))+1))
707 T(307)=0*10000*(1/(10000*(1-TBM(X(14),X(22)))+1))
708 T(308)=4*10000*(1/(10000*(1-TBM(X(14),X(23)))+1))
709 T(309)=2*10000*(1/(10000*(1-TBM(X(14),X(24)))+1))
710 T(310)=0*10000*(1/(10000*(1-TBM(X(14),X(25)))+1))
711 T(311)=0*10000*(1/(10000*(1-TBM(X(14),X(26)))+1))
712 T(312)=4*10000*(1/(10000*(1-TBM(X(14),X(27)))+1))
713 T(313)=2*10000*(1/(10000*(1-TBM(X(14),X(28)))+1))
714 T(314)=5*10000*(1/(10000*(1-TBM(X(14),X(29)))+1))
715 T(315)=5*10000*(1/(10000*(1-TBM(X(14),X(30)))+1))
716 T(316)=4*10000*(1/(10000*(1-TBM(X(15),X(16)))+1))
717 T(317)=5*10000*(1/(10000*(1-TBM(X(15),X(17)))+1))
718 T(318)=1*10000*(1/(10000*(1-TBM(X(15),X(18)))+1))
719 T(319)=0*10000*(1/(10000*(1-TBM(X(15),X(19)))+1))
720 T(320)=1*10000*(1/(10000*(1-TBM(X(15),X(20)))+1))
721 T(321)=0*10000*(1/(10000*(1-TBM(X(15),X(21)))+1))
722 T(322)=5*10000*(1/(10000*(1-TBM(X(15),X(22)))+1))
723 T(323)=0*10000*(1/(10000*(1-TBM(X(15),X(23)))+1))
724 T(324)=2*10000*(1/(10000*(1-TBM(X(15),X(24)))+1))
725 T(325)=0*10000*(1/(10000*(1-TBM(X(15),X(25)))+1))
726 T(326)=0*10000*(1/(10000*(1-TBM(X(15),X(26)))+1))
727 T(327)=5*10000*(1/(10000*(1-TBM(X(15),X(27)))+1))
728 T(328)=1*10000*(1/(10000*(1-TBM(X(15),X(28)))+1))
729 T(329)=1*10000*(1/(10000*(1-TBM(X(15),X(29)))+1))
730 T(330)=0*10000*(1/(10000*(1-TBM(X(15),X(30)))+1))
731 T(331)=0*10000*(1/(10000*(1-TBM(X(16),X(17)))+1))
732 T(332)=3*10000*(1/(10000*(1-TBM(X(16),X(18)))+1))
733 T(333)=0*10000*(1/(10000*(1-TBM(X(16),X(19)))+1))
734 T(334)=2*10000*(1/(10000*(1-TBM(X(16),X(20)))+1))
735 T(335)=2*10000*(1/(10000*(1-TBM(X(16),X(21)))+1))
736 T(336)=0*10000*(1/(10000*(1-TBM(X(16),X(22)))+1))
737 T(337)=2*10000*(1/(10000*(1-TBM(X(16),X(23)))+1))
738 T(338)=0*10000*(1/(10000*(1-TBM(X(16),X(24)))+1))
739 T(339)=5*10000*(1/(10000*(1-TBM(X(16),X(25)))+1))
740 T(340)=0*10000*(1/(10000*(1-TBM(X(16),X(26)))+1))
741 T(341)=5*10000*(1/(10000*(1-TBM(X(16),X(27)))+1))
742 T(342)=2*10000*(1/(10000*(1-TBM(X(16),X(28)))+1))
743 T(343)=5*10000*(1/(10000*(1-TBM(X(16),X(29)))+1))
744 T(344)=10*10000*(1/(10000*(1-TBM(X(16),X(30)))+1))
745 T(345)=2*10000*(1/(10000*(1-TBM(X(17),X(18)))+1))
746 T(346)=2*10000*(1/(10000*(1-TBM(X(17),X(19)))+1))
747 T(347)=0*10000*(1/(10000*(1-TBM(X(17),X(20)))+1))
748 T(348)=0*10000*(1/(10000*(1-TBM(X(17),X(21)))+1))
749 T(349)=0*10000*(1/(10000*(1-TBM(X(17),X(22)))+1))
750 T(350)=6*10000*(1/(10000*(1-TBM(X(17),X(23)))+1))
751 T(351)=5*10000*(1/(10000*(1-TBM(X(17),X(24)))+1))
752 T(352)=3*10000*(1/(10000*(1-TBM(X(17),X(25)))+1))
753 T(353)=5*10000*(1/(10000*(1-TBM(X(17),X(26)))+1))
754 T(354)=0*10000*(1/(10000*(1-TBM(X(17),X(27)))+1))
755 T(355)=0*10000*(1/(10000*(1-TBM(X(17),X(28)))+1))
756 T(356)=5*10000*(1/(10000*(1-TBM(X(17),X(29)))+1))
757 T(357)=1*10000*(1/(10000*(1-TBM(X(17),X(30)))+1))
758 T(358)=5*10000*(1/(10000*(1-TBM(X(18),X(19)))+1))
759 T(359)=1*10000*(1/(10000*(1-TBM(X(18),X(20)))+1))
760 T(360)=2*10000*(1/(10000*(1-TBM(X(18),X(21)))+1))
761 T(361)=10*10000*(1/(10000*(1-TBM(X(18),X(22)))+1))
762 T(362)=10*10000*(1/(10000*(1-TBM(X(18),X(23)))+1))
763 T(363)=4*10000*(1/(10000*(1-TBM(X(18),X(24)))+1))
764 T(364)=0*10000*(1/(10000*(1-TBM(X(18),X(25)))+1))
765 T(365)=0*10000*(1/(10000*(1-TBM(X(18),X(26)))+1))
766 T(366)=5*10000*(1/(10000*(1-TBM(X(18),X(27)))+1))
767 T(367)=0*10000*(1/(10000*(1-TBM(X(18),X(28)))+1))
768 T(368)=0*10000*(1/(10000*(1-TBM(X(18),X(29)))+1))
769 T(369)=0*10000*(1/(10000*(1-TBM(X(18),X(30)))+1))
770 T(370)=0*10000*(1/(10000*(1-TBM(X(19),X(20)))+1))
771 T(371)=5*10000*(1/(10000*(1-TBM(X(19),X(21)))+1))
772 T(372)=5*10000*(1/(10000*(1-TBM(X(19),X(22)))+1))
773 T(373)=1*10000*(1/(10000*(1-TBM(X(19),X(23)))+1))
774 T(374)=0*10000*(1/(10000*(1-TBM(X(19),X(24)))+1))
775 T(375)=5*10000*(1/(10000*(1-TBM(X(19),X(25)))+1))
776 T(376)=2*10000*(1/(10000*(1-TBM(X(19),X(26)))+1))
777 T(377)=1*10000*(1/(10000*(1-TBM(X(19),X(27)))+1))
778 T(378)=2*10000*(1/(10000*(1-TBM(X(19),X(28)))+1))
779 T(379)=10*10000*(1/(10000*(1-TBM(X(19),X(29)))+1))
780 T(380)=10*10000*(1/(10000*(1-TBM(X(19),X(30)))+1))
781 T(381)=5*10000*(1/(10000*(1-TBM(X(20),X(21)))+1))
782 T(382)=2*10000*(1/(10000*(1-TBM(X(20),X(22)))+1))
783 T(383)=1*10000*(1/(10000*(1-TBM(X(20),X(23)))+1))
784 T(384)=3*10000*(1/(10000*(1-TBM(X(20),X(24)))+1))
785 T(385)=1*10000*(1/(10000*(1-TBM(X(20),X(25)))+1))
786 T(386)=5*10000*(1/(10000*(1-TBM(X(20),X(26)))+1))
787 T(387)=6*10000*(1/(10000*(1-TBM(X(20),X(27)))+1))
788 T(388)=5*10000*(1/(10000*(1-TBM(X(20),X(28)))+1))
789 T(389)=5*10000*(1/(10000*(1-TBM(X(20),X(29)))+1))
790 T(390)=3*10000*(1/(10000*(1-TBM(X(20),X(30)))+1))
791 T(391)=4*10000*(1/(10000*(1-TBM(X(21),X(22)))+1))
792 T(392)=0*10000*(1/(10000*(1-TBM(X(21),X(23)))+1))
793 T(393)=1*10000*(1/(10000*(1-TBM(X(21),X(24)))+1))
794 T(394)=0*10000*(1/(10000*(1-TBM(X(21),X(25)))+1))
795 T(395)=0*10000*(1/(10000*(1-TBM(X(21),X(26)))+1))
796 T(396)=0*10000*(1/(10000*(1-TBM(X(21),X(27)))+1))
797 T(397)=5*10000*(1/(10000*(1-TBM(X(21),X(28)))+1))
798 T(398)=0*10000*(1/(10000*(1-TBM(X(21),X(29)))+1))
799 T(399)=0*10000*(1/(10000*(1-TBM(X(21),X(30)))+1))
800 T(400)=5*10000*(1/(10000*(1-TBM(X(22),X(23)))+1))
801 T(401)=0*10000*(1/(10000*(1-TBM(X(22),X(24)))+1))
802 T(402)=4*10000*(1/(10000*(1-TBM(X(22),X(25)))+1))
803 T(403)=4*10000*(1/(10000*(1-TBM(X(22),X(26)))+1))
804 T(404)=5*10000*(1/(10000*(1-TBM(X(22),X(27)))+1))
805 T(405)=0*10000*(1/(10000*(1-TBM(X(22),X(28)))+1))
806 T(406)=2*10000*(1/(10000*(1-TBM(X(22),X(29)))+1))
807 T(407)=5*10000*(1/(10000*(1-TBM(X(22),X(30)))+1))
808 T(408)=0*10000*(1/(10000*(1-TBM(X(23),X(24)))+1))
809 T(409)=4*10000*(1/(10000*(1-TBM(X(23),X(25)))+1))
810 T(410)=4*10000*(1/(10000*(1-TBM(X(23),X(26)))+1))
811 T(411)=1*10000*(1/(10000*(1-TBM(X(23),X(27)))+1))
812 T(412)=0*10000*(1/(10000*(1-TBM(X(23),X(28)))+1))
813 T(413)=2*10000*(1/(10000*(1-TBM(X(23),X(29)))+1))
814 T(414)=2*10000*(1/(10000*(1-TBM(X(23),X(30)))+1))
815 T(415)=5*10000*(1/(10000*(1-TBM(X(24),X(25)))+1))
816 T(416)=5*10000*(1/(10000*(1-TBM(X(24),X(26)))+1))
817 T(417)=0*10000*(1/(10000*(1-TBM(X(24),X(27)))+1))
818 T(418)=1*10000*(1/(10000*(1-TBM(X(24),X(28)))+1))
819 T(419)=0*10000*(1/(10000*(1-TBM(X(24),X(29)))+1))
820 T(420)=0*10000*(1/(10000*(1-TBM(X(24),X(30)))+1))
821 T(421)=1*10000*(1/(10000*(1-TBM(X(25),X(26)))+1))
822 T(422)=0*10000*(1/(10000*(1-TBM(X(25),X(27)))+1))
823 T(423)=10*10000*(1/(10000*(1-TBM(X(25),X(28)))+1))
824 T(424)=1*10000*(1/(10000*(1-TBM(X(25),X(29)))+1))
825 T(425)=0*10000*(1/(10000*(1-TBM(X(25),X(30)))+1))
826 T(426)=0*10000*(1/(10000*(1-TBM(X(26),X(27)))+1))
827 T(427)=0*10000*(1/(10000*(1-TBM(X(26),X(28)))+1))
828 T(428)=0*10000*(1/(10000*(1-TBM(X(26),X(29)))+1))
829 T(429)=0*10000*(1/(10000*(1-TBM(X(26),X(30)))+1))
830 T(430)=0*10000*(1/(10000*(1-TBM(X(27),X(28)))+1))
831 T(431)=0*10000*(1/(10000*(1-TBM(X(27),X(29)))+1))
832 T(432)=10*10000*(1/(10000*(1-TBM(X(27),X(30)))+1))
833 T(433)=2*10000*(1/(10000*(1-TBM(X(28),X(29)))+1))
834 T(434)=2*10000*(1/(10000*(1-TBM(X(28),X(30)))+1))
835 T(435)=2*10000*(1/(10000*(1-TBM(X(29),X(30)))+1))
1429 T(436)=1*10000*(1/(10000*(1-TBM(X(1),X(31)))+1))
1430 T(437)=0*10000*(1/(10000*(1-TBM(X(1),X(32)))+1))
1431 T(438)=4*10000*(1/(10000*(1-TBM(X(1),X(33)))+1))
1432 T(439)=0*10000*(1/(10000*(1-TBM(X(1),X(34)))+1))
1433 T(440)=3*10000*(1/(10000*(1-TBM(X(1),X(35)))+1))
1434 T(441)=5*10000*(1/(10000*(1-TBM(X(1),X(36)))+1))
1435 T(442)=2*10000*(1/(10000*(1-TBM(X(1),X(37)))+1))
1436 T(443)=4*10000*(1/(10000*(1-TBM(X(1),X(38)))+1))
1437 T(444)=0*10000*(1/(10000*(1-TBM(X(1),X(39)))+1))
1438 T(445)=0*10000*(1/(10000*(1-TBM(X(1),X(40)))+1))
2557 T(446)=1*10000*(1/(10000*(1-TBM(X(2),X(31)))+1))
2558 T(447)=0*10000*(1/(10000*(1-TBM(X(2),X(32)))+1))
2559 T(448)=0*10000*(1/(10000*(1-TBM(X(2),X(33)))+1))
2560 T(449)=1*10000*(1/(10000*(1-TBM(X(2),X(34)))+1))
2561 T(450)=5*10000*(1/(10000*(1-TBM(X(2),X(35)))+1))
2562 T(451)=2*10000*(1/(10000*(1-TBM(X(2),X(36)))+1))
2563 T(452)=2*10000*(1/(10000*(1-TBM(X(2),X(37)))+1))
2564 T(453)=0*10000*(1/(10000*(1-TBM(X(2),X(38)))+1))
2565 T(454)=5*10000*(1/(10000*(1-TBM(X(2),X(39)))+1))
2566 T(455)=4*10000*(1/(10000*(1-TBM(X(2),X(40)))+1))
2684 T(456)=0*10000*(1/(10000*(1-TBM(X(3),X(31)))+1))
2685 T(457)=0*10000*(1/(10000*(1-TBM(X(3),X(32)))+1))
2686 T(458)=0*10000*(1/(10000*(1-TBM(X(3),X(33)))+1))
2687 T(459)=0*10000*(1/(10000*(1-TBM(X(3),X(34)))+1))
2688 T(460)=1*10000*(1/(10000*(1-TBM(X(3),X(35)))+1))
2689 T(461)=5*10000*(1/(10000*(1-TBM(X(3),X(36)))+1))
2690 T(462)=10*10000*(1/(10000*(1-TBM(X(3),X(37)))+1))
2691 T(463)=5*10000*(1/(10000*(1-TBM(X(3),X(38)))+1))
2692 T(464)=1*10000*(1/(10000*(1-TBM(X(3),X(39)))+1))
2693 T(465)=0*10000*(1/(10000*(1-TBM(X(3),X(40)))+1))
2740 T(467)=5*10000*(1/(10000*(1-TBM(X(4),X(31)))+1))
2741 T(468)=2*10000*(1/(10000*(1-TBM(X(4),X(32)))+1))
2742 T(469)=0*10000*(1/(10000*(1-TBM(X(4),X(33)))+1))
2743 T(470)=0*10000*(1/(10000*(1-TBM(X(4),X(34)))+1))
2744 T(471)=5*10000*(1/(10000*(1-TBM(X(4),X(35)))+1))
2745 T(472)=10*10000*(1/(10000*(1-TBM(X(4),X(36)))+1))
2746 T(473)=0*10000*(1/(10000*(1-TBM(X(4),X(37)))+1))
2747 T(474)=4*10000*(1/(10000*(1-TBM(X(4),X(38)))+1))
2748 T(475)=5*10000*(1/(10000*(1-TBM(X(4),X(39)))+1))
2749 T(476)=2*10000*(1/(10000*(1-TBM(X(4),X(40)))+1))
2835 T(487)=1*10000*(1/(10000*(1-TBM(X(5),X(31)))+1))
2836 T(488)=2*10000*(1/(10000*(1-TBM(X(5),X(32)))+1))
2837 T(489)=1*10000*(1/(10000*(1-TBM(X(5),X(33)))+1))
2838 T(490)=5*10000*(1/(10000*(1-TBM(X(5),X(34)))+1))
2839 T(491)=2*10000*(1/(10000*(1-TBM(X(5),X(35)))+1))
2840 T(492)=5*10000*(1/(10000*(1-TBM(X(5),X(36)))+1))
2841 T(493)=10*10000*(1/(10000*(1-TBM(X(5),X(37)))+1))
2842 T(494)=5*10000*(1/(10000*(1-TBM(X(5),X(38)))+1))
2843 T(495)=10*10000*(1/(10000*(1-TBM(X(5),X(39)))+1))
2844 T(496)=0*10000*(1/(10000*(1-TBM(X(5),X(40)))+1))
2959 T(497)=2*10000*(1/(10000*(1-TBM(X(6),X(31)))+1))
2960 T(498)=1*10000*(1/(10000*(1-TBM(X(6),X(32)))+1))
2961 T(499)=0*10000*(1/(10000*(1-TBM(X(6),X(33)))+1))
2962 T(500)=1*10000*(1/(10000*(1-TBM(X(6),X(34)))+1))
2963 T(501)=1*10000*(1/(10000*(1-TBM(X(6),X(35)))+1))
2964 T(502)=0*10000*(1/(10000*(1-TBM(X(6),X(36)))+1))
2965 T(503)=5*10000*(1/(10000*(1-TBM(X(6),X(37)))+1))
2966 T(504)=0*10000*(1/(10000*(1-TBM(X(6),X(38)))+1))
2967 T(505)=10*10000*(1/(10000*(1-TBM(X(6),X(39)))+1))
2968 T(506)=0*10000*(1/(10000*(1-TBM(X(6),X(40)))+1))
3052 T(507)=5*10000*(1/(10000*(1-TBM(X(7),X(31)))+1))
3053 T(508)=10*10000*(1/(10000*(1-TBM(X(7),X(32)))+1))
3054 T(509)=5*10000*(1/(10000*(1-TBM(X(7),X(33)))+1))
3055 T(510)=5*10000*(1/(10000*(1-TBM(X(7),X(34)))+1))
3056 T(511)=10*10000*(1/(10000*(1-TBM(X(7),X(35)))+1))
3057 T(512)=0*10000*(1/(10000*(1-TBM(X(7),X(36)))+1))
3058 T(513)=0*10000*(1/(10000*(1-TBM(X(7),X(37)))+1))
3059 T(514)=1*10000*(1/(10000*(1-TBM(X(7),X(38)))+1))
3060 T(515)=2*10000*(1/(10000*(1-TBM(X(7),X(39)))+1))
3061 T(516)=0*10000*(1/(10000*(1-TBM(X(7),X(40)))+1))
3101 T(517)=5*10000*(1/(10000*(1-TBM(X(8),X(31)))+1))
3102 T(518)=0*10000*(1/(10000*(1-TBM(X(8),X(32)))+1))
3103 T(519)=0*10000*(1/(10000*(1-TBM(X(8),X(33)))+1))
3104 T(520)=2*10000*(1/(10000*(1-TBM(X(8),X(34)))+1))
3105 T(521)=0*10000*(1/(10000*(1-TBM(X(8),X(35)))+1))
3106 T(522)=2*10000*(1/(10000*(1-TBM(X(8),X(36)))+1))
3107 T(523)=1*10000*(1/(10000*(1-TBM(X(8),X(37)))+1))
3108 T(524)=1*10000*(1/(10000*(1-TBM(X(8),X(38)))+1))
3109 T(525)=2*10000*(1/(10000*(1-TBM(X(8),X(39)))+1))
3110 T(526)=10*10000*(1/(10000*(1-TBM(X(8),X(40)))+1))
3201 T(527)=5*10000*(1/(10000*(1-TBM(X(9),X(31)))+1))
3202 T(528)=0*10000*(1/(10000*(1-TBM(X(9),X(32)))+1))
3203 T(529)=10*10000*(1/(10000*(1-TBM(X(9),X(33)))+1))
3204 T(530)=1*10000*(1/(10000*(1-TBM(X(9),X(34)))+1))
3205 T(531)=0*10000*(1/(10000*(1-TBM(X(9),X(35)))+1))
3206 T(532)=3*10000*(1/(10000*(1-TBM(X(9),X(36)))+1))
3207 T(533)=0*10000*(1/(10000*(1-TBM(X(9),X(37)))+1))
3208 T(534)=0*10000*(1/(10000*(1-TBM(X(9),X(38)))+1))
3209 T(535)=1*10000*(1/(10000*(1-TBM(X(9),X(39)))+1))
3210 T(536)=2*10000*(1/(10000*(1-TBM(X(9),X(40)))+1))
3301 T(537)=5*10000*(1/(10000*(1-TBM(X(10),X(31)))+1))
3302 T(538)=2*10000*(1/(10000*(1-TBM(X(10),X(32)))+1))
3303 T(539)=1*10000*(1/(10000*(1-TBM(X(10),X(33)))+1))
3304 T(540)=0*10000*(1/(10000*(1-TBM(X(10),X(34)))+1))
3305 T(541)=1*10000*(1/(10000*(1-TBM(X(10),X(35)))+1))
3306 T(542)=1*10000*(1/(10000*(1-TBM(X(10),X(36)))+1))
3307 T(543)=1*10000*(1/(10000*(1-TBM(X(10),X(37)))+1))
3308 T(544)=5*10000*(1/(10000*(1-TBM(X(10),X(38)))+1))
3309 T(545)=6*10000*(1/(10000*(1-TBM(X(10),X(39)))+1))
3310 T(546)=3*10000*(1/(10000*(1-TBM(X(10),X(40)))+1))
3401 T(547)=2*10000*(1/(10000*(1-TBM(X(11),X(31)))+1))
3402 T(548)=0*10000*(1/(10000*(1-TBM(X(11),X(32)))+1))
3403 T(549)=0*10000*(1/(10000*(1-TBM(X(11),X(33)))+1))
3404 T(550)=0*10000*(1/(10000*(1-TBM(X(11),X(34)))+1))
3405 T(551)=0*10000*(1/(10000*(1-TBM(X(11),X(35)))+1))
3406 T(552)=0*10000*(1/(10000*(1-TBM(X(11),X(36)))+1))
3407 T(553)=0*10000*(1/(10000*(1-TBM(X(11),X(37)))+1))
3408 T(554)=0*10000*(1/(10000*(1-TBM(X(11),X(38)))+1))
3409 T(555)=5*10000*(1/(10000*(1-TBM(X(11),X(39)))+1))
3410 T(556)=1*10000*(1/(10000*(1-TBM(X(11),X(40)))+1))
3501 T(557)=5*10000*(1/(10000*(1-TBM(X(12),X(31)))+1))
3503 T(558)=0*10000*(1/(10000*(1-TBM(X(12),X(33)))+1))
3504 T(559)=0*10000*(1/(10000*(1-TBM(X(12),X(34)))+1))
3505 T(560)=0*10000*(1/(10000*(1-TBM(X(12),X(35)))+1))
3506 T(561)=2*10000*(1/(10000*(1-TBM(X(12),X(36)))+1))
3507 T(562)=0*10000*(1/(10000*(1-TBM(X(12),X(37)))+1))
3508 T(563)=0*10000*(1/(10000*(1-TBM(X(12),X(38)))+1))
3509 T(564)=6*10000*(1/(10000*(1-TBM(X(12),X(39)))+1))
3510 T(565)=1*10000*(1/(10000*(1-TBM(X(12),X(40)))+1))
3601 T(566)=0*10000*(1/(10000*(1-TBM(X(13),X(31)))+1))
3602 T(567)=1*10000*(1/(10000*(1-TBM(X(13),X(32)))+1))
3603 T(568)=1*10000*(1/(10000*(1-TBM(X(13),X(33)))+1))
3604 T(569)=1*10000*(1/(10000*(1-TBM(X(13),X(34)))+1))
3605 T(570)=0*10000*(1/(10000*(1-TBM(X(13),X(35)))+1))
3606 T(571)=2*10000*(1/(10000*(1-TBM(X(13),X(36)))+1))
3607 T(572)=0*10000*(1/(10000*(1-TBM(X(13),X(37)))+1))
3608 T(573)=2*10000*(1/(10000*(1-TBM(X(13),X(38)))+1))
3609 T(574)=1*10000*(1/(10000*(1-TBM(X(13),X(39)))+1))
3610 T(575)=1*10000*(1/(10000*(1-TBM(X(13),X(40)))+1))
3701 T(576)=2*10000*(1/(10000*(1-TBM(X(14),X(31)))+1))
3702 T(577)=1*10000*(1/(10000*(1-TBM(X(14),X(32)))+1))
3703 T(578)=3*10000*(1/(10000*(1-TBM(X(14),X(33)))+1))
3704 T(579)=2*10000*(1/(10000*(1-TBM(X(14),X(34)))+1))
3705 T(580)=5*10000*(1/(10000*(1-TBM(X(14),X(35)))+1))
3706 T(581)=0*10000*(1/(10000*(1-TBM(X(14),X(36)))+1))
3707 T(582)=10*10000*(1/(10000*(1-TBM(X(14),X(37)))+1))
3708 T(583)=2*10000*(1/(10000*(1-TBM(X(14),X(38)))+1))
3709 T(584)=10*10000*(1/(10000*(1-TBM(X(14),X(39)))+1))
3710 T(585)=0*10000*(1/(10000*(1-TBM(X(14),X(40)))+1))
3801 T(586)=0*10000*(1/(10000*(1-TBM(X(15),X(31)))+1))
3802 T(587)=0*10000*(1/(10000*(1-TBM(X(15),X(32)))+1))
3803 T(588)=0*10000*(1/(10000*(1-TBM(X(15),X(33)))+1))
3804 T(589)=0*10000*(1/(10000*(1-TBM(X(15),X(34)))+1))
3805 T(590)=5*10000*(1/(10000*(1-TBM(X(15),X(35)))+1))
3806 T(591)=1*10000*(1/(10000*(1-TBM(X(15),X(36)))+1))
3807 T(592)=0*10000*(1/(10000*(1-TBM(X(15),X(37)))+1))
3808 T(593)=3*10000*(1/(10000*(1-TBM(X(15),X(38)))+1))
3809 T(594)=5*10000*(1/(10000*(1-TBM(X(15),X(39)))+1))
3810 T(595)=0*10000*(1/(10000*(1-TBM(X(15),X(40)))+1))
3901 T(596)=5*10000*(1/(10000*(1-TBM(X(16),X(31)))+1))
3902 T(597)=0*10000*(1/(10000*(1-TBM(X(16),X(32)))+1))
3903 T(598)=5*10000*(1/(10000*(1-TBM(X(16),X(33)))+1))
3904 T(599)=5*10000*(1/(10000*(1-TBM(X(16),X(34)))+1))
3905 T(600)=0*10000*(1/(10000*(1-TBM(X(16),X(35)))+1))
3906 T(601)=10*10000*(1/(10000*(1-TBM(X(16),X(36)))+1))
3907 T(602)=0*10000*(1/(10000*(1-TBM(X(16),X(37)))+1))
3908 T(603)=2*10000*(1/(10000*(1-TBM(X(16),X(38)))+1))
3909 T(604)=5*10000*(1/(10000*(1-TBM(X(16),X(39)))+1))
3910 T(605)=0*10000*(1/(10000*(1-TBM(X(16),X(40)))+1))
4001 T(606)=0*10000*(1/(10000*(1-TBM(X(17),X(31)))+1))
4002 T(607)=10*10000*(1/(10000*(1-TBM(X(17),X(32)))+1))
4003 T(608)=2*10000*(1/(10000*(1-TBM(X(17),X(33)))+1))
4004 T(609)=5*10000*(1/(10000*(1-TBM(X(17),X(34)))+1))
4005 T(610)=4*10000*(1/(10000*(1-TBM(X(17),X(35)))+1))
4006 T(611)=10*10000*(1/(10000*(1-TBM(X(17),X(36)))+1))
4007 T(612)=4*10000*(1/(10000*(1-TBM(X(17),X(37)))+1))
4008 T(613)=3*10000*(1/(10000*(1-TBM(X(17),X(38)))+1))
4009 T(614)=5*10000*(1/(10000*(1-TBM(X(17),X(39)))+1))
4010 T(615)=0*10000*(1/(10000*(1-TBM(X(17),X(40)))+1))
4169 T(616)=5*10000*(1/(10000*(1-TBM(X(18),X(31)))+1))
4170 T(617)=0*10000*(1/(10000*(1-TBM(X(18),X(32)))+1))
4171 T(618)=0*10000*(1/(10000*(1-TBM(X(18),X(33)))+1))
4172 T(619)=10*10000*(1/(10000*(1-TBM(X(18),X(34)))+1))
4173 T(620)=5*10000*(1/(10000*(1-TBM(X(18),X(35)))+1))
4174 T(621)=0*10000*(1/(10000*(1-TBM(X(18),X(36)))+1))
4175 T(622)=2*10000*(1/(10000*(1-TBM(X(18),X(37)))+1))
4176 T(623)=2*10000*(1/(10000*(1-TBM(X(18),X(38)))+1))
4177 T(624)=10*10000*(1/(10000*(1-TBM(X(18),X(39)))+1))
4178 T(625)=2*10000*(1/(10000*(1-TBM(X(18),X(40)))+1))
4280 T(626)=1*10000*(1/(10000*(1-TBM(X(19),X(31)))+1))
4281 T(627)=2*10000*(1/(10000*(1-TBM(X(19),X(32)))+1))
4282 T(628)=0*10000*(1/(10000*(1-TBM(X(19),X(33)))+1))
4283 T(629)=10*10000*(1/(10000*(1-TBM(X(19),X(34)))+1))
4284 T(630)=0*10000*(1/(10000*(1-TBM(X(19),X(35)))+1))
4285 T(631)=1*10000*(1/(10000*(1-TBM(X(19),X(36)))+1))
4286 T(632)=0*10000*(1/(10000*(1-TBM(X(19),X(37)))+1))
4287 T(633)=6*10000*(1/(10000*(1-TBM(X(19),X(38)))+1))
4288 T(634)=2*10000*(1/(10000*(1-TBM(X(19),X(39)))+1))
4289 T(635)=0*10000*(1/(10000*(1-TBM(X(19),X(40)))+1))
4390 T(636)=5*10000*(1/(10000*(1-TBM(X(20),X(31)))+1))
4391 T(637)=1*10000*(1/(10000*(1-TBM(X(20),X(32)))+1))
4392 T(638)=5*10000*(1/(10000*(1-TBM(X(20),X(33)))+1))
4393 T(639)=2*10000*(1/(10000*(1-TBM(X(20),X(34)))+1))
4394 T(640)=5*10000*(1/(10000*(1-TBM(X(20),X(35)))+1))
4395 T(641)=1*10000*(1/(10000*(1-TBM(X(20),X(36)))+1))
4396 T(642)=0*10000*(1/(10000*(1-TBM(X(20),X(37)))+1))
4397 T(643)=1*10000*(1/(10000*(1-TBM(X(20),X(38)))+1))
4398 T(644)=0*10000*(1/(10000*(1-TBM(X(20),X(39)))+1))
4399 T(645)=2*10000*(1/(10000*(1-TBM(X(20),X(40)))+1))
4499 T(646)=2*10000*(1/(10000*(1-TBM(X(21),X(31)))+1))
4500 T(647)=0*10000*(1/(10000*(1-TBM(X(21),X(32)))+1))
4501 T(648)=0*10000*(1/(10000*(1-TBM(X(21),X(33)))+1))
4502 T(649)=4*10000*(1/(10000*(1-TBM(X(21),X(34)))+1))
4503 T(650)=10*10000*(1/(10000*(1-TBM(X(21),X(35)))+1))
4504 T(651)=0*10000*(1/(10000*(1-TBM(X(21),X(36)))+1))
4505 T(652)=5*10000*(1/(10000*(1-TBM(X(21),X(37)))+1))
4506 T(653)=0*10000*(1/(10000*(1-TBM(X(21),X(38)))+1))
4507 T(654)=0*10000*(1/(10000*(1-TBM(X(21),X(39)))+1))
4508 T(655)=0*10000*(1/(10000*(1-TBM(X(21),X(40)))+1))
4537 T(656)=2*10000*(1/(10000*(1-TBM(X(22),X(31)))+1))
4538 T(657)=10*10000*(1/(10000*(1-TBM(X(22),X(32)))+1))
4539 T(658)=1*10000*(1/(10000*(1-TBM(X(22),X(33)))+1))
4540 T(659)=10*10000*(1/(10000*(1-TBM(X(22),X(34)))+1))
4541 T(660)=5*10000*(1/(10000*(1-TBM(X(22),X(35)))+1))
4542 T(661)=1*10000*(1/(10000*(1-TBM(X(22),X(36)))+1))
4543 T(662)=1*10000*(1/(10000*(1-TBM(X(22),X(37)))+1))
4544 T(663)=1*10000*(1/(10000*(1-TBM(X(22),X(38)))+1))
4545 T(664)=2*10000*(1/(10000*(1-TBM(X(22),X(39)))+1))
4546 T(665)=0*10000*(1/(10000*(1-TBM(X(22),X(40)))+1))
4614 T(666)=10*10000*(1/(10000*(1-TBM(X(23),X(31)))+1))
4615 T(667)=5*10000*(1/(10000*(1-TBM(X(23),X(32)))+1))
4616 T(668)=2*10000*(1/(10000*(1-TBM(X(23),X(33)))+1))
4617 T(669)=5*10000*(1/(10000*(1-TBM(X(23),X(34)))+1))
4618 T(670)=5*10000*(1/(10000*(1-TBM(X(23),X(35)))+1))
4619 T(671)=2*10000*(1/(10000*(1-TBM(X(23),X(36)))+1))
4620 T(672)=0*10000*(1/(10000*(1-TBM(X(23),X(37)))+1))
4621 T(673)=5*10000*(1/(10000*(1-TBM(X(23),X(38)))+1))
4622 T(674)=5*10000*(1/(10000*(1-TBM(X(23),X(39)))+1))
4623 T(675)=10*10000*(1/(10000*(1-TBM(X(23),X(40)))+1))
4720 T(676)=0*10000*(1/(10000*(1-TBM(X(24),X(31)))+1))
4721 T(677)=1*10000*(1/(10000*(1-TBM(X(24),X(32)))+1))
4722 T(678)=2*10000*(1/(10000*(1-TBM(X(24),X(33)))+1))
4723 T(679)=1*10000*(1/(10000*(1-TBM(X(24),X(34)))+1))
4724 T(680)=2*10000*(1/(10000*(1-TBM(X(24),X(35)))+1))
4725 T(681)=2*10000*(1/(10000*(1-TBM(X(24),X(36)))+1))
4726 T(682)=5*10000*(1/(10000*(1-TBM(X(24),X(37)))+1))
4727 T(683)=6*10000*(1/(10000*(1-TBM(X(24),X(38)))+1))
4728 T(684)=10*10000*(1/(10000*(1-TBM(X(24),X(39)))+1))
4729 T(685)=2*10000*(1/(10000*(1-TBM(X(24),X(40)))+1))
4825 T(686)=5*10000*(1/(10000*(1-TBM(X(25),X(31)))+1))
4826 T(687)=10*10000*(1/(10000*(1-TBM(X(25),X(32)))+1))
4827 T(688)=10*10000*(1/(10000*(1-TBM(X(25),X(33)))+1))
4828 T(689)=10*10000*(1/(10000*(1-TBM(X(25),X(34)))+1))
4829 T(690)=1*10000*(1/(10000*(1-TBM(X(25),X(35)))+1))
4830 T(691)=2*10000*(1/(10000*(1-TBM(X(25),X(36)))+1))
4831 T(692)=2*10000*(1/(10000*(1-TBM(X(25),X(37)))+1))
4832 T(693)=2*10000*(1/(10000*(1-TBM(X(25),X(38)))+1))
4833 T(694)=2*10000*(1/(10000*(1-TBM(X(25),X(39)))+1))
4834 T(695)=5*10000*(1/(10000*(1-TBM(X(25),X(40)))+1))
4929 T(696)=4*10000*(1/(10000*(1-TBM(X(26),X(31)))+1))
4930 T(697)=0*10000*(1/(10000*(1-TBM(X(26),X(32)))+1))
4931 T(698)=0*10000*(1/(10000*(1-TBM(X(26),X(33)))+1))
4932 T(699)=5*10000*(1/(10000*(1-TBM(X(26),X(34)))+1))
4933 T(700)=5*10000*(1/(10000*(1-TBM(X(26),X(35)))+1))
4934 T(701)=2*10000*(1/(10000*(1-TBM(X(26),X(36)))+1))
4935 T(702)=5*10000*(1/(10000*(1-TBM(X(26),X(37)))+1))
4936 T(703)=10*10000*(1/(10000*(1-TBM(X(26),X(38)))+1))
4937 T(704)=0*10000*(1/(10000*(1-TBM(X(26),X(39)))+1))
4938 T(705)=5*10000*(1/(10000*(1-TBM(X(26),X(40)))+1))
5032 T(706)=5*10000*(1/(10000*(1-TBM(X(27),X(31)))+1))
5033 T(707)=5*10000*(1/(10000*(1-TBM(X(27),X(32)))+1))
5034 T(708)=0*10000*(1/(10000*(1-TBM(X(27),X(33)))+1))
5035 T(709)=2*10000*(1/(10000*(1-TBM(X(27),X(34)))+1))
5036 T(710)=3*10000*(1/(10000*(1-TBM(X(27),X(35)))+1))
5037 T(711)=10*10000*(1/(10000*(1-TBM(X(27),X(36)))+1))
5038 T(712)=1*10000*(1/(10000*(1-TBM(X(27),X(37)))+1))
5039 T(713)=1*10000*(1/(10000*(1-TBM(X(27),X(38)))+1))
5040 T(714)=0*10000*(1/(10000*(1-TBM(X(27),X(39)))+1))
5041 T(715)=0*10000*(1/(10000*(1-TBM(X(27),X(40)))+1))
5134 T(716)=2*10000*(1/(10000*(1-TBM(X(28),X(31)))+1))
5135 T(717)=2*10000*(1/(10000*(1-TBM(X(28),X(32)))+1))
5136 T(718)=0*10000*(1/(10000*(1-TBM(X(28),X(33)))+1))
5137 T(719)=0*10000*(1/(10000*(1-TBM(X(28),X(34)))+1))
5138 T(720)=0*10000*(1/(10000*(1-TBM(X(28),X(35)))+1))
5139 T(721)=0*10000*(1/(10000*(1-TBM(X(28),X(36)))+1))
5140 T(722)=1*10000*(1/(10000*(1-TBM(X(28),X(37)))+1))
5141 T(723)=0*10000*(1/(10000*(1-TBM(X(28),X(38)))+1))
5142 T(724)=5*10000*(1/(10000*(1-TBM(X(28),X(39)))+1))
5143 T(725)=10*10000*(1/(10000*(1-TBM(X(28),X(40)))+1))
5235 REM
5236 T(726)=5*10000*(1/(10000*(1-TBM(X(29),X(31)))+1))
5237 T(727)=10*10000*(1/(10000*(1-TBM(X(29),X(32)))+1))
5238 T(728)=1*10000*(1/(10000*(1-TBM(X(29),X(33)))+1))
5239 T(729)=10*10000*(1/(10000*(1-TBM(X(29),X(34)))+1))
5240 T(730)=2*10000*(1/(10000*(1-TBM(X(29),X(35)))+1))
5241 T(731)=10*10000*(1/(10000*(1-TBM(X(29),X(36)))+1))
5242 T(732)=2*10000*(1/(10000*(1-TBM(X(29),X(37)))+1))
5243 T(733)=5*10000*(1/(10000*(1-TBM(X(29),X(38)))+1))
5244 T(734)=5*10000*(1/(10000*(1-TBM(X(29),X(39)))+1))
5245 T(735)=0*10000*(1/(10000*(1-TBM(X(29),X(40)))+1))
5301 T(801)=2*10000*(1/(10000*(1-TBM(X(30),X(31)))+1))
5302 T(802)=2*10000*(1/(10000*(1-TBM(X(30),X(32)))+1))
5303 T(803)=1*10000*(1/(10000*(1-TBM(X(30),X(33)))+1))
5304 T(804)=0*10000*(1/(10000*(1-TBM(X(30),X(34)))+1))
5305 T(805)=1*10000*(1/(10000*(1-TBM(X(30),X(35)))+1))
5306 T(806)=5*10000*(1/(10000*(1-TBM(X(30),X(36)))+1))
5307 T(807)=5*10000*(1/(10000*(1-TBM(X(30),X(37)))+1))
5308 T(808)=10*10000*(1/(10000*(1-TBM(X(30),X(38)))+1))
5309 T(809)=0*10000*(1/(10000*(1-TBM(X(30),X(39)))+1))
5310 T(810)=2*10000*(1/(10000*(1-TBM(X(30),X(40)))+1))
5331 T(811)=1*10000*(1/(10000*(1-TBM(X(31),X(32)))+1))
5332 T(812)=2*10000*(1/(10000*(1-TBM(X(31),X(33)))+1))
5333 T(813)=1*10000*(1/(10000*(1-TBM(X(31),X(34)))+1))
5334 T(814)=1*10000*(1/(10000*(1-TBM(X(31),X(35)))+1))
5335 T(815)=10*10000*(1/(10000*(1-TBM(X(31),X(36)))+1))
5336 T(816)=10*10000*(1/(10000*(1-TBM(X(31),X(37)))+1))
5337 T(817)=0*10000*(1/(10000*(1-TBM(X(31),X(38)))+1))
5338 T(818)=0*10000*(1/(10000*(1-TBM(X(31),X(39)))+1))
5339 T(819)=5*10000*(1/(10000*(1-TBM(X(31),X(40)))+1))
5401 T(821)=0*10000*(1/(10000*(1-TBM(X(32),X(33)))+1))
5402 T(822)=0*10000*(1/(10000*(1-TBM(X(32),X(34)))+1))
5403 T(823)=1*10000*(1/(10000*(1-TBM(X(32),X(35)))+1))
5404 T(824)=2*10000*(1/(10000*(1-TBM(X(32),X(36)))+1))
5405 T(825)=5*10000*(1/(10000*(1-TBM(X(32),X(37)))+1))
5406 T(826)=10*10000*(1/(10000*(1-TBM(X(32),X(38)))+1))
5407 T(827)=1*10000*(1/(10000*(1-TBM(X(32),X(39)))+1))
5408 T(828)=6*10000*(1/(10000*(1-TBM(X(32),X(40)))+1))
5511 T(831)=2*10000*(1/(10000*(1-TBM(X(33),X(34)))+1))
5512 T(832)=2*10000*(1/(10000*(1-TBM(X(33),X(35)))+1))
5513 T(833)=10*10000*(1/(10000*(1-TBM(X(33),X(36)))+1))
5514 T(834)=2*10000*(1/(10000*(1-TBM(X(33),X(37)))+1))
5515 T(835)=10*10000*(1/(10000*(1-TBM(X(33),X(38)))+1))
5516 T(836)=2*10000*(1/(10000*(1-TBM(X(33),X(39)))+1))
5517 T(837)=4*10000*(1/(10000*(1-TBM(X(33),X(40)))+1))
5601 T(841)=6*10000*(1/(10000*(1-TBM(X(34),X(35)))+1))
5602 T(842)=0*10000*(1/(10000*(1-TBM(X(34),X(36)))+1))
5603 T(843)=0*10000*(1/(10000*(1-TBM(X(34),X(37)))+1))
5604 T(844)=1*10000*(1/(10000*(1-TBM(X(34),X(38)))+1))
5605 T(845)=4*10000*(1/(10000*(1-TBM(X(34),X(39)))+1))
5606 T(846)=4*10000*(1/(10000*(1-TBM(X(34),X(40)))+1))
5701 T(851)=2*10000*(1/(10000*(1-TBM(X(35),X(36)))+1))
5702 T(852)=5*10000*(1/(10000*(1-TBM(X(35),X(37)))+1))
5703 T(853)=5*10000*(1/(10000*(1-TBM(X(35),X(38)))+1))
5704 T(854)=2*10000*(1/(10000*(1-TBM(X(35),X(39)))+1))
5705 T(855)=1*10000*(1/(10000*(1-TBM(X(35),X(40)))+1))
5751 T(861)=0*10000*(1/(10000*(1-TBM(X(36),X(37)))+1))
5752 T(862)=0*10000*(1/(10000*(1-TBM(X(36),X(38)))+1))
5753 T(863)=10*10000*(1/(10000*(1-TBM(X(36),X(39)))+1))
5754 T(864)=2*10000*(1/(10000*(1-TBM(X(36),X(40)))+1))
5781 T(865)=6*10000*(1/(10000*(1-TBM(X(37),X(38)))+1))
5782 T(866)=2*10000*(1/(10000*(1-TBM(X(37),X(39)))+1))
5783 T(867)=2*10000*(1/(10000*(1-TBM(X(37),X(40)))+1))
5801 T(871)=6*10000*(1/(10000*(1-TBM(X(38),X(39)))+1))
5802 T(872)=0*10000*(1/(10000*(1-TBM(X(38),X(40)))+1))
5851 T(881)=5*10000*(1/(10000*(1-TBM(X(39),X(40)))+1))
5901 T(901)=2*10000*(1/(10000*(1-TBM(X(1),X(41)))+1))
5902 T(902)=10*10000*(1/(10000*(1-TBM(X(1),X(42)))+1))
5903 T(903)=5*10000*(1/(10000*(1-TBM(X(2),X(41)))+1))
5904 T(904)=4*10000*(1/(10000*(1-TBM(X(2),X(42)))+1))
5905 T(905)=1*10000*(1/(10000*(1-TBM(X(3),X(41)))+1))
5906 T(906)=5*10000*(1/(10000*(1-TBM(X(3),X(42)))+1))
5907 T(907)=6*10000*(1/(10000*(1-TBM(X(4),X(41)))+1))
5908 T(908)=5*10000*(1/(10000*(1-TBM(X(4),X(42)))+1))
5909 T(909)=0*10000*(1/(10000*(1-TBM(X(5),X(41)))+1))
5910 T(910)=2*10000*(1/(10000*(1-TBM(X(5),X(42)))+1))
5911 T(911)=0*10000*(1/(10000*(1-TBM(X(6),X(41)))+1))
5912 T(912)=0*10000*(1/(10000*(1-TBM(X(6),X(42)))+1))
5913 T(913)=1*10000*(1/(10000*(1-TBM(X(7),X(41)))+1))
5914 T(914)=2*10000*(1/(10000*(1-TBM(X(7),X(42)))+1))
5915 T(915)=5*10000*(1/(10000*(1-TBM(X(8),X(41)))+1))
5916 T(916)=0*10000*(1/(10000*(1-TBM(X(8),X(42)))+1))
5917 T(917)=5*10000*(1/(10000*(1-TBM(X(9),X(41)))+1))
5918 T(918)=0*10000*(1/(10000*(1-TBM(X(9),X(42)))+1))
5919 T(919)=2*10000*(1/(10000*(1-TBM(X(10),X(41)))+1))
5920 T(920)=5*10000*(1/(10000*(1-TBM(X(10),X(42)))+1))
5921 T(921)=5*10000*(1/(10000*(1-TBM(X(11),X(41)))+1))
5922 T(922)=0*10000*(1/(10000*(1-TBM(X(11),X(42)))+1))
5923 T(923)=10*10000*(1/(10000*(1-TBM(X(12),X(41)))+1))
5924 T(924)=0*10000*(1/(10000*(1-TBM(X(12),X(42)))+1))
5925 T(925)=1*10000*(1/(10000*(1-TBM(X(13),X(41)))+1))
5926 T(926)=4*10000*(1/(10000*(1-TBM(X(13),X(42)))+1))
5927 T(927)=2*10000*(1/(10000*(1-TBM(X(14),X(41)))+1))
5928 T(928)=0*10000*(1/(10000*(1-TBM(X(14),X(42)))+1))
5929 T(929)=0*10000*(1/(10000*(1-TBM(X(15),X(41)))+1))
5930 T(930)=0*10000*(1/(10000*(1-TBM(X(15),X(42)))+1))
5931 T(931)=2*10000*(1/(10000*(1-TBM(X(16),X(41)))+1))
5932 T(932)=0*10000*(1/(10000*(1-TBM(X(16),X(42)))+1))
5933 T(933)=5*10000*(1/(10000*(1-TBM(X(17),X(41)))+1))
5934 T(934)=5*10000*(1/(10000*(1-TBM(X(17),X(42)))+1))
5935 T(935)=0*10000*(1/(10000*(1-TBM(X(18),X(41)))+1))
5936 T(936)=10*10000*(1/(10000*(1-TBM(X(18),X(42)))+1))
5937 T(937)=10*10000*(1/(10000*(1-TBM(X(19),X(41)))+1))
5938 T(938)=0*10000*(1/(10000*(1-TBM(X(19),X(42)))+1))
5939 T(939)=0*10000*(1/(10000*(1-TBM(X(20),X(41)))+1))
5940 T(940)=0*10000*(1/(10000*(1-TBM(X(20),X(42)))+1))
5941 T(941)=2*10000*(1/(10000*(1-TBM(X(21),X(41)))+1))
5942 T(942)=1*10000*(1/(10000*(1-TBM(X(21),X(42)))+1))
5943 T(943)=0*10000*(1/(10000*(1-TBM(X(22),X(41)))+1))
5944 T(944)=5*10000*(1/(10000*(1-TBM(X(22),X(42)))+1))
5945 T(945)=1*10000*(1/(10000*(1-TBM(X(23),X(41)))+1))
5946 T(946)=0*10000*(1/(10000*(1-TBM(X(23),X(42)))+1))
5947 T(947)=0*10000*(1/(10000*(1-TBM(X(24),X(41)))+1))
5948 T(948)=0*10000*(1/(10000*(1-TBM(X(24),X(42)))+1))
5949 T(949)=2*10000*(1/(10000*(1-TBM(X(25),X(41)))+1))
5950 T(950)=1*10000*(1/(10000*(1-TBM(X(25),X(42)))+1))
5951 T(951)=0*10000*(1/(10000*(1-TBM(X(26),X(41)))+1))
5952 T(952)=5*10000*(1/(10000*(1-TBM(X(26),X(42)))+1))
5953 T(953)=10*10000*(1/(10000*(1-TBM(X(27),X(41)))+1))
5954 T(954)=0*10000*(1/(10000*(1-TBM(X(27),X(42)))+1))
5955 T(955)=1*10000*(1/(10000*(1-TBM(X(28),X(41)))+1))
5956 T(956)=0*10000*(1/(10000*(1-TBM(X(28),X(42)))+1))
5957 T(957)=10*10000*(1/(10000*(1-TBM(X(29),X(41)))+1))
5958 T(958)=10*10000*(1/(10000*(1-TBM(X(29),X(42)))+1))
5959 T(959)=0*10000*(1/(10000*(1-TBM(X(30),X(41)))+1))
5960 T(960)=2*10000*(1/(10000*(1-TBM(X(30),X(42)))+1))
5961 T(961)=1*10000*(1/(10000*(1-TBM(X(31),X(41)))+1))
5962 T(962)=1*10000*(1/(10000*(1-TBM(X(31),X(42)))+1))
5963 T(963)=0*10000*(1/(10000*(1-TBM(X(32),X(41)))+1))
5964 T(964)=0*10000*(1/(10000*(1-TBM(X(32),X(42)))+1))
5965 T(965)=3*10000*(1/(10000*(1-TBM(X(33),X(41)))+1))
5966 T(966)=5*10000*(1/(10000*(1-TBM(X(33),X(42)))+1))
5967 T(967)=10*10000*(1/(10000*(1-TBM(X(34),X(41)))+1))
5968 T(968)=0*10000*(1/(10000*(1-TBM(X(34),X(42)))+1))
5969 T(969)=5*10000*(1/(10000*(1-TBM(X(35),X(41)))+1))
5970 T(970)=0*10000*(1/(10000*(1-TBM(X(35),X(42)))+1))
5971 T(971)=2*10000*(1/(10000*(1-TBM(X(36),X(41)))+1))
5972 T(972)=1*10000*(1/(10000*(1-TBM(X(36),X(42)))+1))
5973 T(973)=1*10000*(1/(10000*(1-TBM(X(37),X(41)))+1))
5974 T(974)=3*10000*(1/(10000*(1-TBM(X(37),X(42)))+1))
5975 T(975)=10*10000*(1/(10000*(1-TBM(X(38),X(41)))+1))
5976 T(976)=5*10000*(1/(10000*(1-TBM(X(38),X(42)))+1))
5977 T(977)=0*10000*(1/(10000*(1-TBM(X(39),X(41)))+1))
5978 T(978)=10*10000*(1/(10000*(1-TBM(X(39),X(42)))+1))
5979 T(979)=0*10000*(1/(10000*(1-TBM(X(40),X(41)))+1))
5980 T(980)=1*10000*(1/(10000*(1-TBM(X(40),X(42)))+1))
5981 T(981)=4*10000*(1/(10000*(1-TBM(X(41),X(42)))+1))
6151 P1NEW=0
6152 FOR KAU7=1 TO 981
6153 P1NEW=P1NEW+T(KAU7)
6154 NEXT KAU7
6450 P=-P1NEW
6451 IF P<=M THEN 6670
6657 FOR KEW=1 TO 42
6658 A(KEW)=X(KEW)
6659 NEXT KEW
6661 M=P
6662 PP1=P1NEW
6669 GOTO 128
6670 NEXT I
6890 IF M>-2700000! THEN 6941 ELSE 6999
6941 PRINT A(1),A(2),A(3),A(4),A(5),A(6),A(7),A(8),A(9),A(10)
6942 PRINT A(11),A(12),A(13),A(14),A(15),A(16),A(17),A(18),A(19),A(20)
6943 PRINT A(21),A(22),A(23),A(24),A(25),A(26),A(27),A(28),A(29),A(30)
6944 PRINT A(31),A(32),A(33),A(34),A(35),A(36),A(37),A(38),A(39),A(40)
6945 PRINT A(41),A(42)
6951 PRINT JJJJ,M
6999 NEXT JJJJ

This BASIC computer program was run with the IBM basica/D interpreter, and its output through JJJJ=-21195 is presented below.

3 3 1 1 4
1 1 4 3 1
3 4 3 3 4
3 3 4 3 3
4 3 3 5 6
4 5 5 4 4
5 4 4 5 5
5 6 5 6 6
6 5
-30974 -2672059.794020599

1 3 1 3 3
4 4 1 3 4
3 3 1 1 4
3 3 4 3 3
4 3 4 5 4
5 6 5 6 5
5 5 5 5 6
4 4 4 6 6
5 5
-28948 -2692057.794220578

4 3 4 3 3
1 1 1 3 1
3 3 4 1 4
3 3 3 3 3
4 3 4 5 6
4 4 5 4 5
5 5 4 5 5
6 6 6 4 6
5 5
-24437 -2602066.793320668

1 3 1 3 3
4 4 1 3 4
3 3 1 1 4
3 3 3 3 3
4 3 4 6 4
5 6 5 6 5
5 5 5 5 5
4 4 6 5 6
4 4
-22969 -2652061.793820618

3 1 4 4 3
1 1 4 3 1
3 4 3 1 4
3 3 4 3 3
4 3 3 5 6
4 5 5 4 4
5 4 4 5 5
6 6 5 6 5
6 5
-21897 -2652061.793820619

4 3 1 3 3
1 1 4 3 1
3 3 4 1 4
3 3 3 3 3
4 3 4 5 6
4 6 5 6 4
5 4 4 5 5
5 6 5 4 6
4 5
-21195 -2582068.793120688

Copied manually from the computer monitor and interpreted in accordance with line 6941 through 6951, the output above was produced in 22 hours on a personal computer with an Intel 2.66 GHz. chip and the IBM interpreter, which is slower than the corresponding compiler.

References

[1] L. R. Arnold, R. E. Beckwith, C. M. Jones, Scheduling the 41st-ORSA-meeting sessions: the visiting-fireman problem, II, Operations Research 21 (1973) 1095-1103.

[2] R. C. Carlson, G. L. Nemhauser, Scheduling to minimize interaction cost, Operations Research 14 (1966) 52-58.

[3] F. E. Davis, M. D. Devine, R. P. Lutz, Scheduling activities among conflicting facilities to minimize conflict, Mathematical Programming 6 (1974) 224-228.

[4] R. P. Lutz, M. D. Devine, H. J. Kumin, W. C. Smith, An application of operations research to school desegregation, Management Science 19 (1972) P100-P109.

[5] C. E. Nugent, T. E. Vollmann, J. Ruml, An experimental comparison of techniques for the assignment of facilities to locations, Operations Research 16 (1968) 150-173.

[6] J. Skorin-Kapov, Tabu Search Applied to the Quadratic Assignment Problem, Journal on Computing 2 (1990) 33-45.