mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   PARI/GP (https://www.mersenneforum.org/forumdisplay.php?f=155)
-   -   PARI's commands (https://www.mersenneforum.org/showthread.php?t=13636)

science_man_88 2010-09-02 16:17

yeah on the kbb(a,x,m) idea I apparently forgot I was doing something in Pari and I didn't realize it even though by the time it gives it should of been at 3 in the morning when I was asleep this started.

a=1;for(m=1,10,for(x=1,m^m-1,kbb(a,x,m)))

*** print1: user interrupt after 10h, 40mn, 54,500 ms.
(13:11) gp > a=1;for(m=1,10,for(x=1,m^m-1,kbb(a,x,m)))


oh and the highest I see in the last few lines before I stopped it:

65390961287 which is over 65 billion

stand corrected:

98418329759 over 98 billion found.

3.14159 2010-09-03 00:43

I'm currently looking for a record k-b-b. (b = 28657), and that will also be the largest prime I ever found,. :smile: (Expected digits: 127737)

3.14159 2010-09-03 13:34

I have now commenced testing for k * 2[sup]328750[/sup] + 1 and k * 28657[sup]28657[/sup] + 1.

Here are the odds for each:

For the base 2 search: 1 in 4283.
For the base 28657 search: 1 in 5812.

Time taken for each test: 3 minutes for base 2.
Expected time for base 2 completion: Approx. 9 days.
Expected time for base 28657 completion: .. 81 days.

science_man_88 2010-09-03 19:02

[QUOTE=3.14159;228275]I have now commenced testing for k * 2[sup]328750[/sup] + 1 and k * 28657[sup]28657[/sup] + 1.

Here are the odds for each:

For the base 2 search: 1 in 4283.
For the base 28657 search: 1 in 5812.

Time taken for each test: 3 minutes for base 2.
Expected time for base 2 completion: Approx. 9 days.
Expected time for base 28657 completion: .. 81 days.[/QUOTE]

well you're k range for your expected length is 10^114477-10^114478 how many candidates in that range ?

3.14159 2010-09-04 00:26

[QUOTE=science_man_88]well you're k range for your expected length is 10^114477-10^114478 how many candidates in that range ?
[/QUOTE]

For b = 2, there are approximately 4880 candidates. For b = 28657, there are approx. 9289 candidates.

Therefore, I will expect only one prime for each, in 9 and 81 days, respectively.

3.14159 2010-09-04 02:11

I'll report back in November for b = 28657.

science_man_88 2010-09-06 16:49

[CODE](13:45) gp > ifprime(x,p)=if(isprime(x),p)
%11 = (x,p)->if(isprime(x),p)
(13:47) gp > ifprime(100,print(100))
100[/CODE]

I was trying something like this but I can't get it to work. why ? is my best question.

axn 2010-09-06 17:16

[QUOTE=science_man_88;228683]
I was trying something like this but I can't get it to work. why ? is my best question.[/QUOTE]

[CODE]ifprime(p, s)=if(isprime(p), eval(s));
ifprime(100, "print(p)");
ifprime(101, "print(p)");
101
[/CODE]

science_man_88 2010-09-06 17:25

[QUOTE=axn;228687][CODE]ifprime(p, s)=if(isprime(p), eval(s));
ifprime(100, "print(p)");
ifprime(101, "print(p)");
101
[/CODE][/QUOTE]

I should of known this by now lol thank you though in a way this is useless but it is useful in the fact that it could create shorter code when typing it in. as by the look of it it saves 10 characters I was also thinking of one where you can do a && part in the if I think I'll try that next.

science_man_88 2010-09-06 17:32

[CODE](14:28) gp > ifprime(p,u,s)=if(isprime(p) && eval(u), eval(s))
%42 = (p,u,s)->if(isprime(p)&&eval(u),eval(s))
(14:29) gp > ifprime(100,"p%9==1","print(100)")
(14:29) gp > ifprime(101,"p%9==1","print(100)")
(14:29) gp > ifprime(101,"p%9==2","print(100)")
100
(14:29) gp > ifprime(101,"p%9==2","print(p)")
101[/CODE]

thank you axn

science_man_88 2010-09-06 17:42

[QUOTE=science_man_88;228690]I should of known this by now lol thank you though in a way this is useless but it is useful in the fact that it could create shorter code when typing it in. as by the look of it it saves 10 characters I was also thinking of one where you can do a && part in the if I think I'll try that next.[/QUOTE]

okay maybe not but it's at least as small as what it equals.


All times are UTC. The time now is 23:15.

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.