help with Pari
Don Reble
djr at nk.ca
Sun Jun 22 00:55:30 CEST 2003
Seqfans:
I hate it, when programs disagree. Even worse, they sometimes
disagree with themselves.
Is the number "25 1097226443 0988869706 6444184023 1990201201
8513329120 2368754457 4398706423" prime, or not?
(It's the 81st cyclotomic polynomial, evaluated at 21.)
Sure enough, here's a "p-1" certificate. (My format is tighter than
Pari's.)
251097226443098886970664441840231990201201851332912023687544574398706423:
(2,3) (3^27,2) (7^27,2) (11,5) (19,2) (37,2) (199,2)
(421,2) (613,2) (5077,2) (17497,2) (7101932659132249,2)
7101932659132249: (2^3,17) (3^3,7) (7,2) (305479,2) (15376001,2)
And if I ask Pari-2.2.5 to do that,
isprime(..., 1)
it generates the same certificate.
But if I ask Pari-2.2.5 to use APRCL,
isprime(..., 2)
it says, "composite". What does your Pari say?
--
Don Reble djr at nk.ca
More information about the SeqFan
mailing list