Carmichael 3 strong pseudoprimes: Difference between revisions

Content deleted Content added
GP
m →‎{{header|PARI/GP}}: fix + reformat
Line 280: Line 280:
Mod[p2 p3, p1 - 1] == 1 :> Print[p1, "*", p2, "*", p3]]</lang>
Mod[p2 p3, p1 - 1] == 1 :> Print[p1, "*", p2, "*", p3]]</lang>
=={{header|PARI/GP}}==
=={{header|PARI/GP}}==
<lang parigp>f(p)={
<lang parigp>f(p)=my(v=List(),q,r);for(h=2,p-1,for(d=1,h+p-1,if((h+p)*(p-1)%d==0&&Mod(-p,h)^2==d,if(isprime(q=(p-1)*(h+p)/d+1)&&isprime(r=p*q\h+1)&&q*r%(p-1)!=1,listput(v,p*q*r)))));Set(v)
my(v=List(),q,r);
for(h=2,p-1,
for(d=1,h+p-1,
if((h+p)*(p-1)%d==0 && Mod(p,h)^2==-d && isprime(q=(p-1)*(h+p)/d+1) && isprime(r=p*q\h+1)&&q*r%(p-1)==1,
listput(v,p*q*r)
)
)
);
Set(v)
};
forprime(p=3,67,v=f(p); for(i=1,#v,print1(v[i]", ")))</lang>
forprime(p=3,67,v=f(p); for(i=1,#v,print1(v[i]", ")))</lang>
{{out}}
{{out}}
<pre>561, 1105, 2465, 10585, 1729, 2821, 6601, 8911, 15841, 52633, 29341, 46657, 115921, 314821, 530881, 162401, 7207201, 334153, 1024651, 1615681, 3581761, 5444489, 399001, 512461, 1193221, 1857241, 5049001, 5481451, 471905281, 294409, 488881, 1461241, 8134561, 36765901, 252601, 410041, 1909001, 5148001, 7519441, 434932961, 2489462641, 1152271, 3057601, 5968873, 6868261, 11972017, 14913991, 15829633, 43331401, 67902031, 368113411, 564651361, 104569501, 902645857, 958762729, 2508013, 4335241, 17316001, 178837201, 6189121, 9439201, 15247621, 35703361, 60957361, 99036001, 101649241, 329769721, 824389441, 1574601601, 10267951, 163954561, 7991602081,</pre>
<pre>8695, 540793, 552563, 21607, 5814901, 8209657, 500867, 2158003, 6301153, 1480111, 30832211, 1363178701, 3729863, 636757889, 1061464331, 2546740153, 11674241, 4396516867, 673819, 7113859, 22226111, 43370641, 74183807,</pre>


=={{header|Perl 6}}==
=={{header|Perl 6}}==