program grdeg
 
    integer, allocatable :: A(:,:), B(:,:), ecc(:)
    logical :: directed 
 
    read(*, *) directed, n
    print*, directed, n
    allocate(A(n,n), B(n,n), ecc(n))
 
    call rdmat(A, n, directed )
    	!call prmat (A, n, 2)
    call connect(A, n, B) 
 
    ecc = maxval(B, 1)
    irad =  minval(ecc, ecc > 0)
    if (directed) then 
        print*, 'directed graph'
    else
        print*, 'undirected graph'
    end if
 
1 format(A14, (T15, 4(I4:, ',')))    
    print 1, 'radius:        ',irad
    print 1, 'diameter:      ', maxval(ecc)
    print 1, 'center points: ', pack(  (/ (i, i = 1,n) /), ecc == irad)
 
   ! call prmat(B, n, 2)
 
contains
 
    subroutine prmat(A, n, d)
        character fmt*40
        integer, intent(in) :: n, d, A(n,n)
        write(fmt, '(A,I0,A,I0,A)') '(', n, '(I', d, ':,X))'
        write(*, fmt) ((A(j,i), j=1,n) ,i=1,n)
    end subroutine
 
    subroutine rdmat(A, n, directed)
        integer, intent(out) :: A(n,n)
        logical directed
        A = 0
        do
            read(*,*, end=1) i,j
            A(i,j) = A(i,j)+1
            if (.not. directed)  A(j,i) = A(j,i)+1
        end do
    1   continue  
    end subroutine
 
    subroutine connect(A, n,  B)
        integer, intent(in) :: A(n,n)
        integer, intent(out):: B(n,n)
        integer:: C(n,n)
 
        B = 0
        C = A
 
        do i=1,n
            where (B == 0 .and. C > 0) B = i
            C = MATMUL(A,C)
        end do
 
        do i=1,n
            B(i,i) = 0
        end do
 
    end subroutine
 
end program
 
 

CnByb2dyYW0gZ3JkZWcKCiAgICBpbnRlZ2VyLCBhbGxvY2F0YWJsZSA6OiBBKDosOiksIEIoOiw6KSwgZWNjKDopCiAgICBsb2dpY2FsIDo6IGRpcmVjdGVkIAogICAgCiAgICByZWFkKCosICopIGRpcmVjdGVkLCBuCiAgICBwcmludCosIGRpcmVjdGVkLCBuCiAgICBhbGxvY2F0ZShBKG4sbiksIEIobixuKSwgZWNjKG4pKQoKICAgIGNhbGwgcmRtYXQoQSwgbiwgZGlyZWN0ZWQgKQogICAgCSFjYWxsIHBybWF0IChBLCBuLCAyKQogICAgY2FsbCBjb25uZWN0KEEsIG4sIEIpIAogICAgCiAgICBlY2MgPSBtYXh2YWwoQiwgMSkKICAgIGlyYWQgPSAgbWludmFsKGVjYywgZWNjID4gMCkKICAgIGlmIChkaXJlY3RlZCkgdGhlbiAKICAgICAgICBwcmludCosICdkaXJlY3RlZCBncmFwaCcKICAgIGVsc2UKICAgICAgICBwcmludCosICd1bmRpcmVjdGVkIGdyYXBoJwogICAgZW5kIGlmCiAgICAKMSBmb3JtYXQoQTE0LCAoVDE1LCA0KEk0OiwgJywnKSkpICAgIAogICAgcHJpbnQgMSwgJ3JhZGl1czogICAgICAgICcsaXJhZAogICAgcHJpbnQgMSwgJ2RpYW1ldGVyOiAgICAgICcsIG1heHZhbChlY2MpCiAgICBwcmludCAxLCAnY2VudGVyIHBvaW50czogJywgcGFjayggICgvIChpLCBpID0gMSxuKSAvKSwgZWNjID09IGlyYWQpCiAgICAKICAgISBjYWxsIHBybWF0KEIsIG4sIDIpCgpjb250YWlucwoKICAgIHN1YnJvdXRpbmUgcHJtYXQoQSwgbiwgZCkKICAgICAgICBjaGFyYWN0ZXIgZm10KjQwCiAgICAgICAgaW50ZWdlciwgaW50ZW50KGluKSA6OiBuLCBkLCBBKG4sbikKICAgICAgICB3cml0ZShmbXQsICcoQSxJMCxBLEkwLEEpJykgJygnLCBuLCAnKEknLCBkLCAnOixYKSknCiAgICAgICAgd3JpdGUoKiwgZm10KSAoKEEoaixpKSwgaj0xLG4pICxpPTEsbikKICAgIGVuZCBzdWJyb3V0aW5lCiAgICAKICAgIHN1YnJvdXRpbmUgcmRtYXQoQSwgbiwgZGlyZWN0ZWQpCiAgICAgICAgaW50ZWdlciwgaW50ZW50KG91dCkgOjogQShuLG4pCiAgICAgICAgbG9naWNhbCBkaXJlY3RlZAogICAgICAgIEEgPSAwCiAgICAgICAgZG8KICAgICAgICAgICAgcmVhZCgqLCosIGVuZD0xKSBpLGoKICAgICAgICAgICAgQShpLGopID0gQShpLGopKzEKICAgICAgICAgICAgaWYgKC5ub3QuIGRpcmVjdGVkKSAgQShqLGkpID0gQShqLGkpKzEKICAgICAgICBlbmQgZG8KICAgIDEgICBjb250aW51ZSAgCiAgICBlbmQgc3Vicm91dGluZQogICAgCiAgICBzdWJyb3V0aW5lIGNvbm5lY3QoQSwgbiwgIEIpCiAgICAgICAgaW50ZWdlciwgaW50ZW50KGluKSA6OiBBKG4sbikKICAgICAgICBpbnRlZ2VyLCBpbnRlbnQob3V0KTo6IEIobixuKQogICAgICAgIGludGVnZXI6OiBDKG4sbikKCiAgICAgICAgQiA9IDAKICAgICAgICBDID0gQQogICAgCiAgICAgICAgZG8gaT0xLG4KICAgICAgICAgICAgd2hlcmUgKEIgPT0gMCAuYW5kLiBDID4gMCkgQiA9IGkKICAgICAgICAgICAgQyA9IE1BVE1VTChBLEMpCiAgICAgICAgZW5kIGRvCiAgICAgICAgCiAgICAgICAgZG8gaT0xLG4KICAgICAgICAgICAgQihpLGkpID0gMAogICAgICAgIGVuZCBkbwogICAgICAgIAogICAgZW5kIHN1YnJvdXRpbmUKCmVuZCBwcm9ncmFtCgo=

VCAzNgoxMCAyCjI4IDIKMiAxMAoyIDQKMiAyOQoyIDE1CjIzIDI0CjIzIDI5CjE1IDI5CjE1IDE0CjE1IDM0CjcgNAo3IDI0CjE0IDIKMTQgNwoxNCAyOQoxNCAxMQoxNCA5CjE0IDE1CjM0IDE1CjM0IDE0CjM0IDI5CjM0IDI0CjM0IDExCjM0IDMzCjM0IDIwCjI5IDIzCjI5IDcKMjkgMgoyOSAxOAoyOSAyNwoyOSA0CjI5IDEzCjI5IDI0CjI5IDExCjI5IDIwCjI5IDkKMjkgMzQKMjkgMTQKMjkgMTUKMTggMjcKMTggMTMKMTggMTEKMTggMjkKMjcgMTgKMjcgNAoyNyAyNAo0IDIKNCAyNwo0IDEzCjQgMzUKNCAyNAo0IDIwCjQgMjkKMTMgMTgKMTMgMTYKMTMgMzAKMTMgMjAKMTMgMjkKMTMgNAoxMyAyCjI0IDQKMjQgMzAKMjQgNQoyNCAxOQoyNCAyMQoyNCAyMAoyNCAxMQoyNCAyOQoyNCA3CjExIDE4CjExIDI0CjExIDMwCjExIDMzCjExIDIwCjExIDM0CjExIDE0CjIwIDI5CjIwIDExCjIwIDQKMjAgMjQKMjAgMTMKMjAgMzMKMjAgMjEKMjAgMjYKMjAgMjIKMjAgMzQKMjIgMzQKMjIgMTEKMjIgMjAKOSAyOQo5IDIwCjIxIDkKMjEgMjAKMjEgMTkKMjEgNgozMyAyNAozMyAzNQozMyAyMAozMyAzNAozMyAxNAozMyAxMQozNSAzMwozNSA0CjM1IDMwCjM1IDE2CjM1IDE5CjM1IDEyCjM1IDI2CjMwIDEzCjMwIDE5CjMwIDM1CjMwIDExCjMwIDI0CjE2IDM2CjE2IDE5CjE2IDM1CjE2IDEzCjM2IDE2CjMxIDE2CjMxIDE5CjUgMTkKMTkgMzAKMTkgMTYKMTkgNQoxOSAzNQoxOSAzMwoxOSAyNAoxMiAzMwoxMiAzNQoxMiAzCjEyIDI2CjI2IDIxCjI2IDM1CjYgMjEKNiAxOQoxIDYKOCAzCjggNgozIDgKMyA2CjMgMTIKMyAzNQozMyAyOQoyOSAzMwoxNCAzMwoyOSAyMQo=

T 36
10 2
28 2
2 10
2 4
2 29
2 15
23 24
23 29
15 29
15 14
15 34
7 4
7 24
14 2
14 7
14 29
14 11
14 9
14 15
34 15
34 14
34 29
34 24
34 11
34 33
34 20
29 23
29 7
29 2
29 18
29 27
29 4
29 13
29 24
29 11
29 20
29 9
29 34
29 14
29 15
18 27
18 13
18 11
18 29
27 18
27 4
27 24
4 2
4 27
4 13
4 35
4 24
4 20
4 29
13 18
13 16
13 30
13 20
13 29
13 4
13 2
24 4
24 30
24 5
24 19
24 21
24 20
24 11
24 29
24 7
11 18
11 24
11 30
11 33
11 20
11 34
11 14
20 29
20 11
20 4
20 24
20 13
20 33
20 21
20 26
20 22
20 34
22 34
22 11
22 20
9 29
9 20
21 9
21 20
21 19
21 6
33 24
33 35
33 20
33 34
33 14
33 11
35 33
35 4
35 30
35 16
35 19
35 12
35 26
30 13
30 19
30 35
30 11
30 24
16 36
16 19
16 35
16 13
36 16
31 16
31 19
5 19
19 30
19 16
19 5
19 35
19 33
19 24
12 33
12 35
12 3
12 26
26 21
26 35
6 21
6 19
1 6
8 3
8 6
3 8
3 6
3 12
3 35
33 29
29 33
14 33
29 21