The fastest primality test for Fermat numbers.

[Arkadiusz]

The test states that for n > 2,
F(n) is prime iff 5^((F(n)-1)/4) == sqrt(F(n)-1) (mod F(n)).


<Thread posted on April 01, 2011.

[science_man_88]

Originally Posted by Arkadiusz

The test states that for n > 2,
F(n) is prime iff 5^((F(n)-1)/4) == sqrt(F(n)-1) (mod F(n)).


<Thread posted on April 01, 2011.

it says the date already:

so you say:

5^(2^(2^n-2)) mod (2^(2^n)+1) = 2^((2^n)/2) according to my research. I'll have to think on this more, I'm not that advanced.

[science_man_88]

Originally Posted by science_man_88

it says the date already:

so you say:

5^(2^(2^n-2)) mod (2^(2^n)+1) = 2^((2^n)/2) according to my research. I'll have to think on this more, I'm not that advanced.

5^(2^(n+1)-4) mod (2^(2^n)+1) = 2^(2^(n-1)) is what I've worked it down to.

[science_man_88]

Originally Posted by science_man_88

5^(2^(n+1)-4) mod (2^(2^n)+1) = 2^(2^(n-1)) is what I've worked it down to.

which I believe simplifies to 5^(2n-2) mod (2^(2^n)+1) = 2^(2n-2). though I'm not great at remembering math when I want it.

[science_man_88]

Originally Posted by science_man_88

which I believe simplifies to 5^(2n-2) mod (2^(2^n)+1) = 2^(2n-2). though I'm not great at remembering math when I want it.

okay I have to remember things better.

[science_man_88]

(13:56)>for(n=1,100,print1(isprime(F(n))","))
1,1,1,1,0,0,0,0,0,0,0,0,0,
  *** isprime: user interrupt after 15,468 ms.
(13:57)>for(n=1,100,print1(5^((F(n)-1)/4)%(F(n)) == sqrt(F(n)-1)","))
0,0,1,1,
  *** _^_: user interrupt after 12,782 ms.

F(n)= 2^(2^n)+1

it fails twice in the first 4.

[science_man_88]

Originally Posted by science_man_88

(13:56)>for(n=1,100,print1(isprime(F(n))","))
1,1,1,1,0,0,0,0,0,0,0,0,0,
  *** isprime: user interrupt after 15,468 ms.
(13:57)>for(n=1,100,print1(5^((F(n)-1)/4)%(F(n)) == sqrt(F(n)-1)","))
0,0,1,1,
  *** _^_: user interrupt after 12,782 ms.

F(n)= 2^(2^n)+1

it fails twice in the first 4.

forgot the n>2 part