Rev 4330 | Blame | Compare with Previous | Last modification | View Log | RSS feed
!if $wims_read_parm!=slib_header
!goto proc
!endif
slib_parms=4\
,n\
,a\
,b\
,m
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
slib_example=10, 3,7,11
slib_require=pari
!exit
:proc
slib_out=
!distribute item $wims_read_parm into slib_n,slib_a,slib_b,slib_m
slib_A=!randint $slib_m
slib_out=!exec pari {slib_L=vector($slib_n);slib_A=$slib_A;\
for(slib_k=1,$slib_n,\
slib_A=lift(Mod($slib_a*slib_A+$slib_b,$slib_m));\
slib_L[slib_k]=slib_A);print(slib_L);}
slib_out =!trim $slib_out