книги / Надежность и живучесть систем связи
..pdfИсходные массивы
1. МРТ[\ : КРТ]. Отображает матрицу смежности транзитных узлов друг с другом. Для его подготовки составляется таблица из двух столбцов. В каждую строку записываются номера соседних узлов в порядке их возрастания: сначала для узла аи затем для а2 и т. д. При этом если записана пара (а,-, а,), i< /, то в даль нейшем пара (aj, а.,-) также записывается в таблицу. Число таких пар в массиве равно 2NBt.
По окончании записи матрицы смежности в левый столбец таблицы записывается номер оконечного узла той пары, относи тельно которой начинается анализ системы связи. В правый стол бец таблицы записываются номера транзитных узлов, смежных с данным оконечным. Если rmax=max{r(a°i)}f то rmax+ l последних
строк таблицы МРТ заполняются нулями. Массив МРТ представ ляет собой правый столбец подготовленной указанным способом таблицы.
2. MSB[ 1 : KSB] — адресный массив, предназначенный для ускорения поиска нужной информации в массиве МРТ. Номер ячейки массива соответствует номеру транзитного узла. Число соответствует номеру ячейки массива МРТ, с которой на чинается запись информации о смежности узла ай В последней ячейке массива MSB записан номер первой ячейки массива МРТ,
вкоторой записан нуль.
3.Массив МР[ 1 :КМР]. Содержит список соседних узлов сис темы связи. Для подготовки массива составляется таблица, в ко торую записываются пары соседних узлов по порядку возраста ния их номеров: сначала для транзитных, затем для оконечных.
При этом если записана пара {аи а3), *</. то пара (а3) сц) в дальнейшем не записывается.
Таким образом, в первой части таблицы номер узла в левом столбце таблицы всегда меньше номера, записанного в правом столбце. Во второй части таблицы (матрица смежности оконечных
узлов |
с транзитными) в левый столбец записываются номера по |
||||
люсов, |
в правый— номера смежных с ними транзитных |
узлов. |
|||
Массив МР представляет собой правый столбец таблицы. |
|
||||
4. Массив |
: KSB1] — адресный массив, |
предназначен |
|||
ный для ускорения |
поиска нужной |
информации в |
массиве МР. |
||
Номер ячейки массива MSB1 соответствует номеру |
транзитного |
||||
узла или полюса |
системы связи. |
Число MSBl[i] ~ это |
номер |
ячейки массива МР, с которой начинается запись информации о
смежности узла (полюса) а,-. Могут |
быть случаи, когда для како |
|
го-либо транзитного узла ai при |
af, |
1 в массиве нет информации |
(сведения об узлах, соседних с |
записаны ранее). Тогда в |
|
MSBl\[i] записывается число из |
(i-H )-ft ячейки и т. д. |
|
5. Массив VP[[ : MNAL]. Записываются вероятности исправ |
||
ности элементов системы связи |
(за |
исключением полюсов). |
181
6. Массив MPUZ{\ :LP, 1 \KPUZ]. Первые два столбца мас сива содержат номера тяготеющих друг к другу полюсов. В столб цах 3,..., KPVZ записаны номера транзитных узлов, смежных с полюсом, номер которого записан во втором столбце. Часть этих столбцов заполнена нулями из-за неравномерной связности сис
тем связи. |
|
|
|
7. |
Массив YPUZ[ 1 : LP] — интенсивности потоков |
вызовов |
|
между тяготеющими парами полюсов. |
|
|
|
8. |
Массив МК\[\:КМК, 1 :2 ] — массив координат полюсов |
и |
|
У К системы связи. В первом столбце записывается абсцисса, |
во |
||
втором — ордината. В первых MK.UZ строках записаны |
координа |
||
ты У К |
а ^ А т, в остальных — координаты полюсов. |
|
|
Для подготовки массива на карту наносится произвольная ко ординатная сетка и стороны координатных прямоугольников ну меруются слева направо и снизу вверх. Таким образом, первый координатный прямоугольник находится в левом нижнем углу листа. Номера квадрата, где находится вершина, записываются в массив МК.
11.2. Программа формирования множества путей передачи информации между полюсами системы связи
Программа предназначена для формирования множества воз можных путей между оконечными узлами системы связи* при че тырех заданных ограничениях. Объем программы составляет 277 перфокарт.
Исходные данные:
а) переменные KPT, MKUZ, KPUZ, ККТ, N, KPRAL, /Z, LO, KSB, KSB1, КМР, LP, KMPR, NLP-
б) массивы MPUZ, МР, MSB1, MPT, MSB.
Для каждой пары оконечных узлов на печать выводятся сле дующие результаты: -
1. Сообщение об окончании подготовки исходных данных. В случаях некорректности исходных данных, сбоях в работе внеш них устройств ЭВМ, недостаточности размера массива путей про грамма выдает необходимые диагностические сообщения.
2.Номера оконечных узлов и число возможных путей переда чи информации между ними с учетом ограничений по пропускной способности и рангу путей.
3.Номера оконечных узлов и число путей передачи информа ции между ними с учетом дополнительного ограничения по числу исходящих направлений из транзитных узлов.
Для нормальной работы программы необходимо иметь 100 кбайт оперативной 'памяти: 27 кбайт для объектного модуля про граммы; 13 кбайт для исходных данных; 60 кбайт для окончатель ных и промежуточных результатов.
162
|
CCniGM/KTrC/'rfT П 00 >1 2 4 ) I MSHP(100, 20) |
|
|
OtnENSlON h ° UZ( 5 0 0 , 5 1 , HPT( 1 6 0 ) | |
|
|
ZHSatSO) , МР( Л00, 1» ,Hf,ai (3001 , |
|
|
1ЯТ |400) IHA 1 1 fZ » 4 яр. I |
|
|
READ |
1 l»Kp T »M4IJZ , KPUZ I KKT >N , |
|
l KPRAL. IZ. LO. KSDiKsBl,KHP. LP, |
|
|
2KMPP.NUP |
|
|
2iKp TiHKUZ, KPUZ,KKTI N, |
|
|
IKPRAL•1Z , Lfl.KS0.KSB1iHMP. LP, |
|
|
2KHPR,HUP |
|
|
Nil I =N- 1 |
|
|
MNALBMKUZ*KMP |
|
|
READ |
l i t i ( H P U Z t l i J ) , J=1,KPUZ)i |
|
11 = 1. LP) |
|
|
2 2 ■ ( (HPUZ ( I , J 1 i J = 1 i KPUZI i |
|
|
l t = l. l-P) |
|
|
ПЕдО |
3 1 1 ( u p ( j , i ) , j = j , КНР) |
|
2 , ( м р ( j , D , i s i , KHP) |
|
|
REAO 1. ( H5 8 l ( I ) • I s 11KS 3 1) |
|
|
2 i (H Sq J ( I I i I ж1 i К SВ 1 > |
|
|
РЕдО 3 1 , ( HP T ( t ) . 1 = 1 , KPT) |
|
|
2 , (HP T ( I ) , 1 = 1 , KPT) |
|
|
READ |
1 , ( M s n u i . I = 1 I KSB ) |
|
2, (MSQlU 11 = 1, К S 0 ) |
|
101 |
L = l Я |
|
1=1*1 |
|
J=2
102J = J M
IF (HPUZ I L, J ) 1110, 110, 103
tj) |
I ?=HPUZ( L, J) *l |
|
I = h S 3 (12) |
|
I 1=KPT |
134 |
11 = 11-1 |
' |
ЦРТ ( 11♦ 1 1*MPT( I 1 > |
IF ( I 1 - J )90 1 , 9C 1 . 104 |
110 IF (J 1-J) П 6 , n o , j , » 119 0=0*1
1 F ( K T ( ] , J I - м к и г > 1 2 0 , 1 2 0 , 1 2 3 1 2 0 I F ( J - K P R A L 1 1 2 1 , 1 2 j i 122
121H S K TM ,J)
CO TO 113
122КТ ( 1 ,.))s0
0=0-1
CO TO 127
123MORPH I 1»=J
KCI l =HCRP( I , i ) - l 00 124. J 1 =2 , KCM
INO=KT(I,11)
124HAl l ( l , I ND) =i
IF(I*KKT> 1 125, 1 124 Ц t?4
1124 |
5 1, hPUZ (L , 1) i MPL)Z (L , 2) |
|||
1125 |
CO |
TO |
500 |
|
1=1*1 |
|
|
||
|
KCI l = J-2 |
|
||
|
IM 1 = 1-1 |
|
||
125 |
00 |
125 |
11 = 1 , КС I 1 |
|
KT( I , I 1)=KT( I M1 , I j ) |
||||
|
0= J - 2 |
|
|
|
|
CO |
TO |
1129 |
126IA=KTU, J)
1271NDT=MT(1A)-MSB(1A*1)*1 I F ( I NOT »120, 129, 129
120K=h T d A l « l CO TO, 115
129I F ( J - l l 130, 131, 13э
130KT( l i J ) =0
0 =0 - 1
1129 l F ( J r i H 2 6 , i i 3 0 , l 26 1130 1A=MKUZ*1
nsassHifcAUi
KCI 1SMS81( I NO*1) NCI1*I 1*1
10511=11*1
INOT=MP(I 1, 1)-. HPUZ(U, J) IF 11NflJ 1 106, 10H, l p6
106I F ( 1 1“ KCI1>105, 1071107
137 |
PRIHT |
50, HPUZ( Li O) *l «J i NCU» |
1KCI 1 |
259 |
|
U B |
CO TO |
|
MP T <I I =HKUZ * I I |
||
|
KCI1=HKUZ*2 |
|
|
00 109 |
1 1=12, KCM |
109 NSg(111=058( 111*1 |
||
110 |
l F ( J r * p u Z ) 1 3 2 , 1 1 0 , 1 1 0 |
|
00 111 |
I 1=1 , KKT |
|
, |
00 t i l J 1 = 1 , N |
|
И1 KT( I 1, J 1) =0 |
||
|
40 |
|
|
1 = 1 |
|
|
J=1 |
|
|
K T ( I , J l = MpUZ(L, 1) |
|
|
00 112 |
I 1 = 1 , HNAL |
112 |
HA 11 ( 1 , Г11=0 |
|
HT(I 11=0 |
||
|
ГАДКИЕ*! |
|
из K=nsa<m-i |
||
114 |
K=K*l |
|
И5 |
Jl = i |
|
116J1*J1*1
If'OT=MPT (Kl-KT ( I I 01) 1F ( INDr ) 1 1 0 , 1 1 7 * 1 1 0
131 KT(I, O)=0
500 1=1-1
PRINT 1? iHPUZ (1, 1) iHP’uZ (L ,2) , 1 IF ( 1 1259, 259, 132
132I Ac I
1=0
1331=1*1
J=0
134J=j*2
IF 10-2) 136, 135, 136 135 1N0 = KT(1 , J - U .
I N0i =KTU, J |
I 1xhSB 1 ( IND1-1 CO TO 130
136INOTSKTU i O- l > - KT( l , J ) IF(INDT)1 3 5 , 137»137
•137 IND=KT(1,J)
INDt=KT(I,J. -l>
l l * H S B l ( I N O l - l . 130 11=11*1
IF (HP ( 1111) - 1N01 1 130, 139, 130
139I2=HSRP( 1, 11*1
14012=12-1
кт(Ы2*п=кт1Ыг>
I F ( 1 2 - 0) 902, 982, 140
902HCRP ( I , 11=H3R®( 1 , j 1 * 1 KT(I,0)»HKUZ*1 1
I F [ J O- n CRP ( I , 11)134, П А , , A . 141- 1 F ( 1 - U ) 133, 142, 142
1420=0
14300 144 t=l,MNM.
П7 I F | K * H S D ( U n l * U U 4 i l 2 9 l l 2 9 |
HMl l i i l ) = o |
Н « |
МТ I | ) =в |
|
J = J*1 |
|
i f i*, * m i i J M ) H 6 i H 5 i H 4 ' |
1*5 |
lN [,T *«$ filJ* l) -M SB (J)-|Z |
146 |
IF(tNOT) 14 6, 14 6, M7 |
IF (J-4KUZ) И З . И 6 , 186 |
|
M7 |
K=0 |
|
J2*0 |
146K=K*1
01= 1
149O l » J l * l
l F ( K T ( K i J n - J ) 150 * 956 i 150
150IF (J 1* 2 -H C R P| 4, 1 ),1 4 9 , 1 5 1 ,1 5 1
151IF (H-IA) 140, 152 ,15 2
152I F ( H A I I I 2 , J | - I Z ) 15 3,1 53, 959
15300 1 54 Jl«t,HKUZ
154M l n l 2 , J l l « 0
00 155 J 1 *1 ,J2
155 KT(О 1 , Nh1> =0
30 TO 146
956J1*J1*1
IN0=KT(KtJl)
IF MAI 1 ( 2 ,1 4 0 ) ) 9 5 7 , 9 5 7 1958
957H*i 1 (2,J)*HA11 ( 2 , J ) * l
996M4l l(2,!ND)=MAll(2»IND)*l
J2*J2 +1
KT(J 2 I NM1)=K CO TO 151
959 12*0
1561=0
IC*99999
1571=1+1
IFt K TU ,N H l)) 156, 160,156 156 IND=KT(I,NK1)
I F ( MCRP ( I N0, 1 ) • 1 C) 159, 160,160
159IC*MCRP( INO, | ) K*KT(I,NH1)
160IF 1I - J2 ) 157, 161, 161
1611*1
1621*1*1
JF (K T(K ,I» -J1 162 ,1 63, 162
16312*12*1
MT11 2 ) *KTIК, I ♦ 1) 1*0
1641=1*1
INQ»KTII,NM1)
IF <INO»164 ,1 6 5 ,9 6 4
964 11*1
16511*11+1
lN O T * K T ( I N D ,l l) - n f (12)
IF( IN O T )16 7 , 9 4 5 ,1 6 7
965 IF(KT11 NO,I 1 - 1 ) - J ) 167, 166, 107
146K T (I|4 H 1) *0 CO TO 166
167IN O T all ^n C R P il N O ,] )*! I F ( I N 0 T ) 1 6 5 , l 4 6 , l 6 6
166 I F ( I - J 2 ) 1 6 4 , l 4 9 , U 9
169I F M 2 - I D 1 5 6 , 1 7 0Ц 70
1701*0
JA*0
171I * | * l
IF (K T( I ( N i l ) j 17 2 , 1 7 4 ,1 7 2
172K*K T ( I »NH 1)
JAB JA+1
KCI1*HCRP(K, 1)
DO 173 1 1 * 1 , K711
173К Г «К , 1 11*0 Ч5Я Р(К ,1) *0
174I F ( r - J 2 ) l 7 l , i 7 5 , 1 7 5
175I=e
1761=1*1
I F ( p C R p ( I , | ) ) 1 7 7 , 1 & 0 , j 7 7
177I F ( | M A) 1 7 6 , 1 7 6 ,1 7 6
K?( f I ? N M l " l | * * IA»fA-JA
CO TO 146
160(<= I - 1
161K=K*1
I F t n G R P J K , I I ) 1 63 , 1 6 2 , 1 6 3 162 IF(K -IA) 1 6 Ы 7 6 , 176
183KCI1=MCRP(K, 1)
00 164 1 1*1 , КСЦ
K T ( I , t l ) * K T ( K , I l )
184K T|K ,I 1 )* 0
riCRPI!, D=nCRP(K, i) HCRP«K,r))*0
186 |
CO TO |
i i r |
|
CONTINUE |
Ц |
||
|
1 2 , H P U Z < L ,1 ) , HPU2 (U ,2 I |
||
259 I F ( L - L P ) 1 6 7 , 9 9 9 ,9 9 9 |
|
||
167 |
J=2 |
|
|
166J»J*1
I F ( MPU2 ( l , J I ) 169, 192 , 1 6 9
189I 2 BMPUZ11* J >+ 1 l =H5B(12 I " l KCI1*KPT-1
00 |
190 ! 1 * 1 ,KCI 1 |
190НРТ(И)=ПРТ<11*1)
KCI1*HKUZ*2
00 191 I 1*12, KCI I
191NSB (I11* HSB (!1) -1
IF (J- K P U Z )1 0 0 , 1 9 2 , 1 9 2
1921N0T*NPUZ(U, l l - N P u Z ( L * 1, l ) IF ( I NOT ) 193,10.1 , 19 3
193NCIt =MS8( MKUZ*1> KCI 1*MSB(MKU2*2>
03 191 11 =NC11«КС 1 1
1 94 f l P r M 1 ) = 0
0 =MSB(NKUZ+1) IND=HPUZ(L* t , 1 ) J l*hSR 1 1INO)
195Ц Р Т ( Л * И Р ( Л , 1 )
) N O T * J l- flS B l( tN O * i) * i
IF (IN O T )1 9 6 , 1 9 7 , 1 9 7
196 J 1s J 1 + 1
0=j*l
CO TO 195
197 HSB (HKUZ*2 I= J * l CO TO 101
1FORMAT(2613)
2FORMAT(1 X,2 4 I 5)
11FORMAT(2014)
12FORMAT( IX, 3 0 | 4)
13FORHAT(1X,40I3)
21 FORMAT(1615)
22 F0RMAT(1X,2(1I4)
31 |
FORMAT(4012) |
|
|
|
40 |
FOR KATdX,'ОКОНЧЕН |
ЧИКЛ ПОДГО'* |
||
1 |
'т о в к и исходных |
Дан ных*) |
||
50 |
FORMAT(2 X, '9 |
HP |
HfT |
* , 2 Х , 90151 |
я |
FORMA T (2 X , * B |
кт |
м*ло с т р о к », |
1'ДЛЯ ПАРЫ*,216)
999CONTINUE STOP
ENp
184
11.3. Программа расчета надежности (живучести) двухполюсной сети связи корреляционным методом
Программа предназначена для расчета вероятности сохране ния связи каждой из заданных пар оконечных узлов.
Объем программы составляет 485 перфокарт. Исходные данные:
а) переменные KPT, MKUZ, KPUZ, ККТ, N, KVMGRP, L0, KSB, KSB1, КМР, LP, KVKT, KLRG, KPRAL, IZ, KMPR, NLP-V
б) массивы MPUZ, МР, MSB1, MPT, MSB, VP\
в) ключ KLRG.
Если KLRG = О, программа работает только в режиме POISK.
Если KLRG= 1, то дополнительно выполняется и расчет вероятно сти.
Для каждой пары оконечных узлов на печать выводятся сле дующие результаты:
1. Номера оконечных узлов и число путей передачи информа ции между ними с учетом ограничений по пропускной способно сти и рангу путей.
2.Номера оконечных узлов и число путей передачи информа ции между ними с учетом дополнительного ограничения по числу исходящих направлений из транзитных узлов.
3.Трассы путей (левая часть массива КТ) и ранги путей (ле
вый столбец массива MGRP). |
. |
4. При KLRG—1 печатаются вероятности Р{Е) |
исправности |
хотя бы одного пути множества и вероятность а.(е*), i=l,h, ис правности каждого пути.
О ходе работы программы могут быть выданы диагностиче ские сообщения в случае некорректности исходных данных или сбоя в работе машины, а также при достижении граничных разме ров массивов.
Для нормальной работы программы необходимо иметь I6S кбайт оперативной памяти: 63 кбайт для объектного модуля про граммы; 28 кбайт для исходных данных; 75 кбайт для промежу точных и окончательных результатов.
8-32
|
|
ПРОГРАММА |
POISK R |
|
|
COMMON - K T l l f l l i l Z l ) i H C R P f 1 0 0 | 2 0 > |
|||
|
DIMENSION J1PUZ(5'e0 j5 ) ,MPTI350> , H T U 5 5 0 ) , |
|||
|
1MSB ( 6J0 ) »MP (15001 1) I N S | 1 11000> »HCRD (1 3 0 ) • |
|||
|
2HA1 1 (2.1 1550) ,VP( 155») , N J ( 2 , 2 » ) iPHJ < 2 • » « Р Е Л 2 0 ) |
|||
|
READ |
ll i K P T iMKU? I KPUZ tK.KT »N | |
||
|
1 KVMCRPI 1 0 I KSBI K SB 1I KMPI |
|||
|
■2LP,KVKT,KURB,KPRAL » IZ ,K MPR»NLP |
|||
|
PRINT 2|KPT,MKUZ,KPUZ,KKT.N, |
|||
|
1MVMC RP,L 0I KSBI K SBI , KMP» |
|||
|
2LP|KVKT,KLRC,KPRAUIZ.KMPRI NLP |
|||
|
MNAL*MKUZ«-KMP |
|
||
|
KNJ«HVMCRP |
|
||
|
READ |
t i l l ( M P U Z ( I I J ) i J s 1 1KPUZ) . I = 1 »LP > |
||
|
2 | |
( (MPUZ ( 11 J) i J e t , KPUZ) . I s l . L P ) |
||
|
READ |
3 1 , |
IMP ( | , 1) I I M I KMPR) |
|
|
READ 3 1 , (ИР f I •1>» IaNLPiKMP) |
|||
|
SPRINT 2, |
(MP f 1,1) .I*1,KMP) |
||
|
READ |
1 , ( M S B 1 ( ( ) , I > 1 I HKUZ) |
||
|
!NDeMKUZ*l |
|
||
|
REAO |
1 1 , (MSaiU » • 1«1 NO I KSs 1) |
||
|
2, |
IMSB1U) , I * 1 , KS B 1) |
||
|
R E AD |
3 1 , |
I M P T ( I ) , I S I,KPT) |
|
|
2 , |
| M P T | I ) , 1 * 1 , KPT) |
||
|
READ |
1, IMSB(I) , I » i , K S B ) |
||
|
2 , (MSB|! ) i I a l , K S B ) |
|||
|
READ |
B, ( VP * I ) .I«I,MKUZ) |
||
|
1N01«MKUZ*KMPR |
|
||
|
READ |
4 , I V P ( | ) , Д с I MD, I NO 1) |
||
|
J N D » MKUZ*N L P |
|
||
|
READ |
9 , «VP( I ) ,I«lNO,MNAL) |
||
<999 |
9 , I V P ( I | , I » 1 , MNAL| |
|||
L»L0 |
|
|
|
|
111 |
L*L*1 |
|
|
|
|
J« 2 |
|
|
|
i | 2 |
J * J * 1 |
|
|
|
|
i F i M P u z a . a p i i b i i e . l e j |
|||
l f 3 I 2 » H P U Z ( L ,J ) * 1 |
|
|||
|
I - M SB I1 2) |
|
|
|
|
I1«KPT |
|
|
|
114 |
I l e i 1*1 |
|
|
|
|
MPT( I 1 ♦ 1 ) *MPT 11 11 |
|||
|
IF ( I 1 • I ) 9 1 1 , 9 0 1 , |
104 |
||
9 1 1 |
INDeMPUZ(L,2) |
|
||
|
1 1«MS'* 1 ( I NO ( . 1 |
|
||
|
K C li e M S B l( I N O + l) |
|
||
..'if5 |
NCI 1*1 i |
|
|
|
11»11♦ 1 |
|
|
||
106 |
i P ( M P ( I l | l | « N P U Z ( L i J H 106.1,00,106 |
|||
JF 11 1-KCl 1) 10 9 ,1 0 7 , 1 0 7 |
||||
1 1 7 P R I N T |
5 0 , M P U Z C L , J ) , L , J , N c U , K C I l |
|||
|
CO TO |
259 |
|
|
10R |
M P T ( I >*MKUZ + 1 1 |
|
||
|
KCl1eMKUZ *2 |
|
||
|
DO 109 |
I1 > I 2 , КС? 1 |
109MSB«Il)*MS6 II 1)M
J F ( J»KPUZ) J 0 2 1 11 0 , 1 IP
110 |
DO |
111 |
Il iJ,K K T |
111 |
DO |
111 |
J 1 * |, N |
K T I I l i J l l M |
|||
|
I»1 |
|
|
J * i
186
К Т < г I J ) s U P U Z ( L i t )
00 11» I1=1»MNAL
им к ы и ^
1 1 2 , Ц Т ( П » * Я
1 UrMKUZM
113K » H S B ( U ! - 1
114K«K*1
115Л = 1
ЙЙJ 1*01 ♦ 1
IF( M PT ( К ) - к г ( r . O i l ) 1 15 i 117 I 118
117 |
IF(K"MS D ( 1 A* l ) ♦ 1 ) 1 1 4 » 1 2 9 1 1 2 9 |
И в J F ( J l - J ) 1 16 , 1 1 9 . 1 1 9 |
|
119 |
J =J M |
|
кт(I ,xJ) = npr(K> |
MT(IA)=K
| F ( K T < ! | J ) » M KU Z > 1 ? 4 I ! 2 0 I 123
120 I F ( J - K P R A L ) 1 2 1 1 12 l 1 122 l i t IA = K T ( J , J )
СО TO П З
122 КT( 1• 0 | *j0
J*J- 1
|
CO |
TO |
127 |
125 |
MCRp( I , I ) s j |
||
|
KCI 1*M GRP( I , 1 ) - 1 |
||
|
DO |
124 |
11 = ? , К С П |
|
IND = KT ( 1 , 1 1 ) |
||
124 |
МАЦ (1 ,I M p ) = i |
||
|
I F ( I - K K T ) 1 1 2 5 , 1 12 4 i 11 24 |
||
1124 |
5 l , H P U Z < L , l ) , MP(JZ(L,2> |
||
|
CO |
TO |
500 |
1125 |
1=1*1 |
|
|
|
KCI |
l = J - 2 |
|
|
I M1= I • 1 |
||
125 |
DO |
125 |
I 1=1 , KCI 1 |
К Т ( I , I 1 ) с к т ( IM l , I 1) |
|||
|
J * J - 2 |
|
|
|
CO |
TO |
1129 |
126I A e K T ( I , J )
127| F ( M T ( l A ) - h S B ( l A * l ) * l ) 1 2 a , 1 2 9 i l 2 9
126 |
K=MT(IAJ*1 |
||
129 |
CO |
TO |
115 |
J F |
( J M ) |
1 3 0 , 13 1 , 13 0 |
130КТ( I 1-11*0 J =J r 1
1129 IF ( J * i ) 1 2 6 , Ц 3 0 , 126
1130 |
I Aс MKUZ * 1 |
||
|
I F ( Ц Т ( l A ) - r t S 0 ( M K U 2 * 2 I * 1 ) 1 2 8 , 1 3 1 , 1 3 1 |
||
131 |
К Т (I , J ) « 0 |
||
50 0 |
1= 1-1 |
|
|
I Ae I |
|
||
|
2 , M P U Z ( l , 1 J, MP U?(L » 2 | , l |
||
|
IF ( 1 ) 2 5 9 , 2 5 9 , 1 1 3 2 |
||
1132 |
DO |
132 |
I * I i I A |
|
DO |
132 |
J « 2 , КVMGRP |
132 |
MCRP(! , J ) = 0 |
||
|
1*0 |
|
|
133 |
1*1* 1 |
|
|
|
J * 0 |
|
|
1 3 4 |
J = J * 2 |
|
|
|
I F ( J - 2 1 |
136 , 13 5 , 13^ |
|
1 3 5 |
I N D = K T | I , J - 1 I |
||
|
I Mo 1 >*KT ( 11 o ) |
ItB n S e ltlN O !" t
1 3 6 |
I F i K T ( I i J - l > - K T M , J ) ) I 3 5 i l 3 7 » l 3 ? |
1 3 7 |
l MO=KT| r »J ) |
|
J N D l = K T ( I « 0 - 1 ) |
|
1 I s H S B 1 ( I ND) -1 |
13 8 |
11=11*1 |
|
I f (мР U 111 1 - IND1 ) 13 0 1 l 3 9 » 138 |
13 9 |
I2 =HCRP.(I . 1) + J |
340 |
12=12*1 |
|
KT(1 » l 2 * l ) = K T U i l 2 > |
|
IF { I 2 - J ) 9 0 2 , 9 0 2 ? 1 + 0 |
* 0 2 |
H G R P i 1 ) =MCRP11 , 1 ) + 1 |
|
КТ ( I , J ) = M K U 2* I1 |
|
IF ( J O - n C R P t l , i ) > i 3 4 t i 3 4 i M l |
1+1 |
I F ( i - l A ) 1 3 3, 1+2 ?1+2 |
1+2 |
J = 0 |
343 |
DO H 4 I = l?nNAL |
|
M i l l 1 2 , 1 ) = 0 |
1+4 |
MT( 11=0 |
|
J = J * 1 |
1+5 |
l F ( h A l l ( l ? J ) - l ) l 4 6 ? l 4 5 . 1 + 6 |
l F ( MS f l ( J * l ) - MSB [ J ) - l Z ) 1 4 6 ? 1+6*147 |
|
1+6 I F<J - HKUZ) 1 4 3 ? 1 0 6 , 1 0 6 |
|
1+7 |
K =0 |
|
J2=J* |
1+B |
K =K*1 |
1+9 S l s i l M
I F ( K T ( K , J l l * J ) 1 5 0 , 9 3 6 , 150
150I F ( J l * 2 - M C R P ( K ? D ) 1 + 9 , 1 5 1 , 1 5 1
151I F t K - I A ) 1 + 8 , 1 5 2 , 1 5 2
152 |
I F ( H A U ( 2 , J ) - t Z > 153 , 1 5 3 , 9 5 9 |
|
153 |
00 |
15+ J l s l . h K U i |
15+ |
MAI 1 1 1 , J 1 ) = 0 |
|
|
DO |
155 J 1 = 1 , J 2 |
155KT ( J 1 ? NHl >=J0 CO TO 1+6
9 S 6 |
J 1 = J 1 * 1 |
|
1ND*KT(K. J l ) |
|
I F 1ИА11 1 2 , IND>1 9 5 7 , 9 5 7 , 9 5 8 |
* 5 7 И A 11 ( ? , J ) A i 1 ( 2 . . DM |
|
958 |
flA 1 1 I 2 , I МО) =>,Л 11 I 2 * I KOI * 1 |
|
J 2 x J 2 M |
|
KT(J2, NM1)=K |
|
CO TO 15 Г |
959 |
12=0 |
156 |
1=0 |
|
10=99999 |
1 5 7 |
1=1*1 |
|
I F ( КТ( I , N M 1 ) j 1 5 8 , 1 6 3 , 1 5 8 |
1 58 |
J N0 =К T ( I ,N К 1J |
|
I F ( r CRP( I ND, 1 I - I C J 1 5 9 , 160, 161 |
159I Cs NCRPHNO, i)
КT ( I , N N 1 J
16 0 |
I f |
I I - J 2 H 5 7 , 161, |
161 |
161 |
1=1 |
|
|
1 6 2 |
1=1*1 |
* |
|
|
I F ( КТ I К, П - М ) 162 , 16 3 , 1 6 2 |
||
163 |
12 г 12 * 1 |
|
|
|
M T ( I 2 ) = K T ( K , 1*1) |
|
|
|
1 |
= 0 |
|
14+ |
1=1*1 |
|
|
|
J NDBKTI I | NH| ) |
|
'188
J Г ( I NO) 1 6 5 , 1 6 8 , 9 6 + 964 11=1
165I l s I l M
IF (KT t IND»I 1 ) —MT ( I2TJ J 167I ?45 J f $ 7
965 |
J F ( K T U N D , ! l - l ) - J j 1 67 , |
16tf*I(57 |
|
166 |
КТ ( I , NM1)= 0 |
* |
|
|
CO TO 160 |
|
|
167 |
I F ( I 1 * 2 « M C R P C 1 H 0 , U ) 1 6 ? r T 6 ( f t l 5 B |
||
165 |
I F ( ! " J 2 ) 1 6 4 , 1 6 9 , 1 6 9 |
' * |
iV |
169IF 11 2 - Ш 1 56 , 170 , i?fl
1701=0
171 1=1*1 |
v |
I F ( K T M i N H l l ) |
1 7 2 , 1 T 4 1 172 |
172K =KT H , N M1 ) J A=J A*1
KCI 1 = MCRP I K, D
173 |
DO |
173 |
1 1 = 1 , KCM |
|
|
||
КТ I K , 1 1) =0 |
|
|
|
|
|||
|
MCRPI K, 1 ) =0 |
|
|
|
|
||
174 |
1F ( I - J 2) 1 7 1 , |
175, |
175 |
|
|||
175 |
1=0 |
|
|
|
|
|
|
176 |
1=1*1 |
|
|
|
|
|
|
|
IF <MCRP( I , 1 ) ) 177, |
100, 177 - |
|
||||
177 |
I Г ( I - I A ) 1 7 6 , |
175, 170 |
|
||||
175 |
DO |
179 |
1 1 =1 , IA |
|
|
|
|
179 |
К Т ( I 1 , NM 1) = 0 |
|
|
|
|
||
|
I » = I A - J A |
|
|
|
|
||
|
GO |
TO |
146 |
|
|
|
|
150 |
K=I - 1 |
|
|
|
|
|
|
181 |
K=K* 1 |
|
|
|
|
|
|
15? |
IF(MCPP(K,1 )) 1 5 3 , 1 6 2 , 1 5 3 |
|
|||||
IF ( K- l |
A H 5 1 , |
178, 175 |
|
||||
153 |
KCI 1*WCRP IK, l) |
|
|
|
|||
|
00 |
10 4 |
I 1 = 1 , KCI 1 |
|
|
||
|
К Т ( I ,1 1)=KT (К, I 1) |
|
|
||||
154 K T ( K , l l ) S 0 |
|
|
|
|
|||
|
MGRP(I , 1)=MCRP IK, 1) |
|
|||||
|
HCRP(K,1)=0 |
|
|
|
|
||
|
GO |
TO |
177 |
|
|
|
|
106 |
I F ( 1 A ) 2 5 9 , 2 5 9 , 190 |
|
|||||
195 |
KPuT=IA |
|
|
|
, |
||
|
2,MPUZ IL , 1) , HPUZ IL , 2 ) |
, I A |
|||||
|
12, H K T U |
1 , 0 1) , - H= l * N) |
, I 1 =1| J A> |
||||
|
2, (HGRpl I |
, 1) , I = 1 , I A) |
|
||||
|
) F ( K L R G ) 2 5 9 , 2 5 9 , 2 6 0 |
|
|||||
260 |
DO |
199 |
1 = 1 , KN^ |
|
|
|
|
|
N J ( 1 , I ) =0 |
|
|
|
|
||
|
HJ(2 ,I)=0 |
|
|
|
|
||
|
PH'J |
I I ) = 1 . 0 |
|
|
|
|
|
199 |
PFJ I I ) * 1 . 0 |
|
|
|
|
||
|
.1=1 |
• |
|
|
|
|
|
|
NJ (l i 1)=2 |
|
|
|
|
||
|
KCI=MGRP ( 1 , 1 ) |
|
|
|
|
||
|
00 |
20 0 |
1=2, КС I |
|
|
|
|
л |
INQ i = KT( 1, I ) |
|
n , |
Mnll |
|
||
2 00 |
PEJ M ) =Pf J ( 1) *VP 11 NO 11 |
|
PRINT 3 , PEJ fJ )
201IC=0 1 0 = 0
I N Q = N J I l , J )
I F ( 7 - 1 ) 2 0 2 , 2 0 3 , 2 0 2
2 0 2 I 2 = MC R P ( 1 N D , J - 1 ) * 1
1€9
|
СО т о 2<j4 |
2 0 3 |
1 2 * 2 |
2 0 4 DO 2 0 5 1 = 1 i 2 |
|
2*5 |
tJO 205 ll=l.nMM. |
MA1111d 1> =0 |
|
|
KCI1SMCRP(IMDI )> |
|
00 206 I 1=I 2 i КС 11 |
|
I ND1»KT( 1ND, U) |
2 0 6 M A l H b l N D U a l |
|
|
1 =0 |
207 |
HAX=1 |
1=1*1 |
I F ( J » 1 ) 2 0 8 , 2 0 9 , 2 0 8 20 8 I l = M C R P ( I i J - l )
GO TO 210
20911=1
210I2*MGRP<1*J>
21111*11*1 1М01 =К Т ( Ы П
I F ( M A l l ( I » I ND1)“ 1 ) 2 1 2 , 3 1 1 , 2 1 2
З Ц U = 1
CO TO 213
21212*12*1
IF r i 2 - K V K M U 1213 , 1 2 1 3 1 1 2 1 2
1212 |
PRINT. 52 |
|
СО TO258 |
1213 |
KTII Н 2 1 = К Т ( Ь П ) |
|
liCRP (1 , J + 1) = 12 |
213 |
IFtl l-HGRP(I ,J D2l li 2l4 ,2 K |
2 1 4 IF( MAX- I 2*HGRP<I . J ) 1 2 1 5 , 2 1 0 1 2 1 6 |
|
215 |
MAX=I2-NGRP(1iJ> |
2lfl |
I F u i C R p t l , J M ) - HGRP( 1 , 0 l " t ) 2 2 2 » 2 l 9 i 2 2 i |
2191C=IC*1
IN£) 1 =KT СI » 12)
IF (MA 11 ( 2 , 1 N D D - 1 ) 2 2 0 , 2 2 1 ( 2 2 0
220M A I l ( 2 i I N D l ) s l CO TO 222
221I C = I C - 1
MCRPU , J * 1 ) = MCRPU , J + l J t l
222 |
I F d - l N O + D 2 0 7 ( 2 2 3 1 2 2 3 |
|
223 |
K=I NO- 1 |
|
|
1 F ( I J - 1 ) 2 1 6 , 2 1 7 , 2 1 6 |
|
21 6 |
J=J *1 |
248 |
|
CO TO |
|
217 |
12=0 |
|
|
I F ( I C 1 1 3 1 1 , 1 3 1 1 ( 2 2 4 |
|
224 |
1=0 |
|
2 2 5 |
1= 1 * 1 |
|
1251 |
IF{MCQP( I . , JM ) ) 2 3 n 2 3 l , 1 2 5 i |
|
I F ( MCRPM , J * U - M C R P ( l , J J M ) 2 2 6 , 2 3 1 , 2 2 6 |
||
226 |
1 1 aftCRP( I , J ) |
|
2 2 7 |
11*11*1 |
|
|
I N0 l = K T ( J , I 1 ) |
|
228 |
IF (MA 1112, I N D U - 1 ) 2 2 8 , 2 2 9 , 2 2 8 |
|
1F ( I 1 - MC R P ( I , J * 1 ) ) 2 2 7 , 2 3 1 , 2 3 1 |
||
229 |
12=12*1 |
|
|
КГ-f I 2 • KVKT) =1 |
|
|
NC11 =flCRP 1 I , Jl ♦ l |
|
|
КС I 1=HGRP( I , J * 1) |
|
|
DO 230 |
j i=b/Cj 1 ,KCl 1 |
230 |
KT(J,!1)=0 |
«C»P(I ,J * 1 )= 0