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