Subversion Repositories wimsdev

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

slib_author=Bernadette, Perrin-Riou

slib_example=6,3\
6,8\
6
slib_require=pari
!exit

:proc
!reset slib_m
!distribute item $wims_read_parm into slib_n, slib_m

!default slib_m =!randint $slib_n, $[min(floor(5*$slib_n/2), $slib_n*($slib_n-1)/2)]

slib_liste = !values 1 for x=1 to $slib_m
slib_liste0=!values 0 for x=1 to $slib_n*($slib_n-1)/2 -$slib_m
slib_liste =!shuffle $slib_liste,$slib_liste0
slib_G= !exec pari slib_n =$slib_n;slib_ liste =[$slib_liste];\
slib_M=matrix(slib_n,slib_n,slib_i,slib_j, if(slib_i < slib_j, slib_liste[slib_n*(slib_n-1)/2-(slib_n-slib_i+1)*(slib_n-slib_i)/2 +slib_j-slib_i]));\
matrix(slib_n,slib_n,slib_i,slib_j, if(slib_i<slib_j,slib_M[slib_i,slib_j],slib_M[slib_j,slib_i]))

!!On rend le graphe connexe
slib_debut=1

:debut

!read slib/graph/distance [$slib_G]
slib_D=$slib_out

!if -1 notin $slib_D
  !goto end
!endif
slib_D=!translate internal ; to $\
$  in $slib_D
!for slib_c=$slib_debut to $slib_n
  slib_ligne=!line $slib_c of $slib_D
  !if -1 isin $slib_ligne
    slib_J=!positionof item -1 in $slib_ligne
    !for slib_JJ in $slib_J
      slib_G=!exec pari slib_N=Mat([$slib_G]);slib_N[$slib_c,$slib_JJ]=1;slib_N[$slib_JJ,$slib_c]=1;slib_N
    !next slib_JJ
    slib_debut=$[$slib_c+1]
    !goto debut
  !endif
!next slib_c

:end

slib_out=$slib_G