Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header
  !goto proc
!endif

slib_author=Bernadette, Perrin-Riou; Dominique, Bernardi
slib_example=[0,0,1,1,1,1,0,0;0,0,1,1,1,1,1,0;0,0,0,0,0,0,0,1;0,0,0,0,0,0,1,0;0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0;0,0,0,0,0,0,0,0],3,8\
[0,0,1,1,0,1;0,0,1,1,0,0;0,0,0,0,1,0;0,0,0,0,1,1;0,0,0,0,0,0;0,0,0,0,0,0],1,5
!exit

:proc
!reset slib_out

!distribute item $wims_read_parm into slib_A,slib_s,slib_t
slib_A=!declosing $slib_A
slib_A=!lines2rows $slib_A

Directed Acyclic Graph

slib_res=!exec pari (allpaths(M,a,b)={my(n=#M,res=[],res1);if(a==b,res=[[a]],for(i=1,n,if(M[a,i],res1=allpaths(M,i,b);for(j=1,#res1,res=concat(res,[concat(a,res1[j])]))))); res});A=[$slib_A]; allpaths(A,$slib_s,$slib_t)

slib_rescnt=!itemcnt $slib_res
!for slib_ch=1 to $slib_rescnt
  slib_chh=!declosing $(slib_res[$slib_ch])
  slib_out=!append line $slib_chh to $slib_out
!next

!exit

slib_A=!replace internal ; by $\
$ in $slib_A

slib_A= !declosing $slib_A
slib_A=!replace internal ; by $\
$ in $slib_A
slib_n = !linecnt $slib_A
slib_chemin=$slib_t

!for slib_i=1 to $slib_n
  !for slib_k=1 to $slib_n
    !for slib_j=1 to $slib_n
      slib_r=!line $slib_k of $slib_A
      slib_r=!item $slib_j of $slib_r
      !if $slib_r != 0
        slib_cnt=!linecnt $slib_chemin
        !for slib_u = 1 to $slib_cnt
          slib_che=!line $slib_u of $slib_chemin
          slib_che1=!item 1 of $slib_che
          !if $slib_che1 = $slib_j and $slib_k,$slib_che notlineof $slib_chemin
            slib_chemin=$slib_chemin\
$slib_k,$slib_che
          !endif
        !next slib_u
      !endif
    !next slib_j
  !next slib_k
!next slib_i

slib_chemin=!select $slib_chemin where column 1=$slib_s

slib_out=$slib_chemin