PARI's commands

[CRGreathouse]

Originally Posted by 3.14159

I guess I am unable to do much uploading, as the network is failing yet again.

Not a big deal -- you're the one doing the calculation at the moment, and as long as that work is able to continue there's no problem in waiting a day, or even a week, to update.

[3.14159]

Originally Posted by CRGreathouse

Not a big deal -- you're the one doing the calculation at the moment, and as long as that work is able to continue there's no problem in waiting a day, or even a week, to update.

Any programs that you would recommend me using for the larger primes?

My current idea is NewPGen + PFGW + Proving it prime with Proth.exe

And:

912646 * 798336²⁰¹⁶⁰ + 1 is prime! (Approx. 119000 digits)

[CRGreathouse]

Originally Posted by 3.14159

Any programs that you would recommend me using for the larger primes?

Not particularly. As long as your basic setup is sieve + probable prime test + primality proving, you're pretty well set.

[3.14159]

Now, that find lifted a load of CPU usage off the computer. It reduced it by 40% !

Originally Posted by CRGreathouse

Not particularly. As long as your basic setup is sieve + probable prime test + primality proving, you're pretty well set.

Sieve: NewPGen or Vk.

PRP Test: PFGW or LLR.

Primality proof: Proth.exe (Items 1, 2, 3, 4, 8, 9, 10, 11, 12, 13, 14, 15)

[3.14159]

Okay: Past 250, programs will be put into use; The primes reach 600 digits by b = 250. Or begin using "isPRP", for PARI.

[CRGreathouse]

Originally Posted by 3.14159

Or begin using "isPRP", for PARI.

Probably better to use ispseudoprime(), which does a small amount of trial division then a BPSW test. It's coded directly in C which means you get about the same speed as isPRP but a much stronger test.

[3.14159]

Okay: Completed the file for b = 60 to 250.

Now: For 251 to 500.

[CRGreathouse]

Originally Posted by 3.14159

Other members only:

20. n-1 analogues of Proths.
21. n-1 analogues of k-b-b, (k * b^b - 1)
22. Twins.

I stay within 1-19. Ban me for a week if I defy that rule.

Since Pi's not allowed to participate and no one else has stepped forward, I'll submit 1 * 2^2 - 1 = 3 for #20 and #21 and 3 = 5 - 2 for #22. Primality certificates available on request.

[science_man_88]

obviously even though you told me to post again you don't care of my idea.

[CRGreathouse]

Originally Posted by science_man_88

obviously even though you told me to post again you don't care of my idea.

I can't make heads or tails of it.

[science_man_88]

turn this into vectors:

for(s=1,50,for(x=s,20000000*s,if(isprime(x),print1(x","));x=x*2);print("<--"s))

then compare with the exponents of Mersenne primes.

some don't fit them so we could say they may never hit them. if we can figure this out we can use it like a sieve.

3,7,31,127,8191,131071,524287,<--1
2,5,11,23,47,191,383,6143,786431,<--2
3,7,31,127,8191,131071,524287,<--3
19,79,1279,5119,20479,81919,1310719,<--4
5,11,23,47,191,383,6143,786431,<--5
13,223,3583,917503,14680063,<--6
7,31,127,8191,131071,524287,<--7
17,71,1151,73727,294911,18874367,<--8
19,79,1279,5119,20479,81919,1310719,<--9
43,<--10
11,23,47,191,383,6143,786431,<--11
103,1663,109051903,<--12
13,223,3583,917503,14680063,<--13
29,59,239,479,15359,245759,1966079,<--14
31,127,8191,131071,524287,<--15
67,271,1087,1114111,17825791,<--16
17,71,1151,73727,294911,18874367,<--17
37,151,607,39845887,<--18
19,79,1279,5119,20479,81919,1310719,<--19
41,83,167,2687,21503,172031,5505023,<--20
43,<--21
367,1471,94207,<--22
23,47,191,383,6143,786431,<--23
199,12799,51199,3276799,209715199,<--24
103,1663,109051903,<--25
53,107,431,863,6911,27647,442367,<--26
223,3583,917503,14680063,<--27