X-Git-Url: https://git.distorted.org.uk/~mdw/catacomb/blobdiff_plain/0f00dc4c8eb47e67bc0f148c2dd109f73a451e0a..1519ef6617fd7351bd40026b9eeb47c7f7f42234:/math/t/pgen

diff --git a/math/t/pgen b/math/t/pgen
index 8d52d197..be26c183 100644
--- a/math/t/pgen
+++ b/math/t/pgen
@@ -14,6 +14,29 @@ pgen {
   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;
@@ -24,6 +47,30 @@ primep {
   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 {