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