Rev 4330 | Details | Compare with Previous | Last modification | View Log | RSS feed
| Rev | Author | Line No. | Line |
|---|---|---|---|
| 20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
| 13586 | bpr | 2 | !goto proc |
| 20 | reyssat | 3 | !endif |
| 4161 | bpr | 4 | |
| 20 | reyssat | 5 | slib_parms=4\ |
| 6 | ,n\ |
||
| 7 | ,a\ |
||
| 8 | ,b\ |
||
| 9 | ,m |
||
| 4330 | bpr | 10 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
| 4161 | bpr | 11 | |
| 20 | reyssat | 12 | slib_example=10, 3,7,11 |
| 3265 | bpr | 13 | slib_require=pari |
| 20 | reyssat | 14 | !exit |
| 15 | |||
| 16 | :proc |
||
| 17 | slib_out= |
||
| 18 | |||
| 19 | !distribute item $wims_read_parm into slib_n,slib_a,slib_b,slib_m |
||
| 20 | slib_A=!randint $slib_m |
||
| 21 | |||
| 22 | slib_out=!exec pari {slib_L=vector($slib_n);slib_A=$slib_A;\ |
||
| 23 | for(slib_k=1,$slib_n,\ |
||
| 24 | slib_A=lift(Mod($slib_a*slib_A+$slib_b,$slib_m));\ |
||
| 25 | slib_L[slib_k]=slib_A);print(slib_L);} |
||
| 26 | |||
| 27 | slib_out =!trim $slib_out |