166359567317705838255275971708060308423814413741683015010175247351623188739655446196925981468626681882384215574706593049022467680136399439302347043107836749816290369600677730213469006507173065402294688841278559283358390567733443050775707749725690534182003442070447739085348456478911335969765393755383551520173 166359567317705838255275971708060308423814413741683015010175247351623188739655446196925981468626681882384215574706593049022467680136399439302347043107836749816290369600677730213469006507173065402294688841278559283358390567733443050775707749725690534182003442070447739085348456478911335969765393755383551520257;
}
+pgen-granfrob {
+ 5 0 0 -1;
+ 7 0 0 4;
+ 15 0 0 3;
+ 5777 1 -1 4; # pseudoprime
+ 40301809 0 0 4;
+ 86059163416987297647409667483582114939806237974424324409828198660056356336227 1 5 4;
+ 102508420970861015999300753620309481186457893679971500520427161277511389396803 1 5 4;
+ 72291866454056552194087337607224612505157525245486245416393486917859196707519 1 5 4;
+ 72291866454056552194087337607224612505157525245486265416393486917859196707519 1 5 3;
+
+ ## A large Frobenius pseudoprime: call the first number p_1; then p_2 = 31
+ ## (p_1 + 1) - 1 and p_3 = 43 (p_1 + 1) - 1. These three are all prime.
+ ## Their product is a strong Lucas, and Frobenius, pseudoprime.
+ ##
+ ## See `Prime and Prejudice' by Martin R. Albrecht, Jake Massimo, Kenneth
+ ## G. Paterson, and Juraj Somorovsky.
+ 3690125385954346893658786222051913500627130245213169388019826598097107079718295481926241398412699320815932808015860263240282855670239765686869973444864115322609857375876438922226372746215468824202413623127 0 0 4;
+ 114393886964584753703422372883609318519441037601608251028614624541010319471267159939713483350793678945293917048491668160448768525777432736292969176790787575000905578652169606589017555132679533550274822316967 0 0 4;
+ 158675391596036916427327807548232280526966600544166283684852543718175604427886705722828380131746070795085110744681991319332162793820309924535408858129156958872223867162686873655734028087265159440703785794503 0 0 4;
+ 66981291792500223036804182765508448534715465524671325885174850970812009004775815201151227900130153990294748113034471984909912807896550069799856170439734910206802409847773026240559371480115711600866989845251707737806461503879250232804362190067578216069266197879151809743235261582813331022213587929425243163096486125825510076936556242805690400001899138503900919499414951069309064408305196756524628693684938044145785145327821174180933033293089394794328963673467918652042794300291355500468079109432376296868174257674548727592142782202898031102246775544402811199608266683925072825828225074019194302318324623049819212337927 0 0 4;
+}
+
primep {
-5 0;
-1 0;
4 0;
40301809 1;
40301811 0;
+
+ ## A small Lucas pseudoprime: 5777 = 53*109.
+ 5777 0;
+
+ ## A large strong pseudoprime: this is the product of
+ ##
+ ## p_1 = 142445387161415482404826365418175962266689133006163
+ ## p_2 = 5840260873618034778597880982145214452934254453252643
+ ## p_3 = 14386984103302963722887462907235772188935602433622363
+ ##
+ ## See `Prime and Prejudice' by Martin R. Albrecht, Jake Massimo, Kenneth
+ ## G. Paterson, and Juraj Somorovsky.
+ 142445387161415482404826365418175962266689133006163 1;
+ 5840260873618034778597880982145214452934254453252643 1;
+ 14386984103302963722887462907235772188935602433622363 1;
+ 11968794224604718293549908104759518204343930652759288592987578098131927050572705181539873293848476235393230314654912729920657864630317971562727057595285667 0;
+
+ ## A large Lucas pseudoprime: call the first number p_1; then p_2 = 31 (p_1
+ ## + 1) - 1 and p_3 = 43 (p_1 + 1) - 1. These three are all prime. Their
+ ## product is a strong Lucas pseudoprime.
+ 3690125385954346893658786222051913500627130245213169388019826598097107079718295481926241398412699320815932808015860263240282855670239765686869973444864115322609857375876438922226372746215468824202413623127 1;
+ 114393886964584753703422372883609318519441037601608251028614624541010319471267159939713483350793678945293917048491668160448768525777432736292969176790787575000905578652169606589017555132679533550274822316967 1;
+ 158675391596036916427327807548232280526966600544166283684852543718175604427886705722828380131746070795085110744681991319332162793820309924535408858129156958872223867162686873655734028087265159440703785794503 1;
+ 66981291792500223036804182765508448534715465524671325885174850970812009004775815201151227900130153990294748113034471984909912807896550069799856170439734910206802409847773026240559371480115711600866989845251707737806461503879250232804362190067578216069266197879151809743235261582813331022213587929425243163096486125825510076936556242805690400001899138503900919499414951069309064408305196756524628693684938044145785145327821174180933033293089394794328963673467918652042794300291355500468079109432376296868174257674548727592142782202898031102246775544402811199608266683925072825828225074019194302318324623049819212337927 0;
}
primeiter {