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)

cmd 2010-10-25 00:40

I
 
[QUOTE=cmd;234301][URL="http://4.bp.blogspot.com/_rvR3ouziO8g/TMTL4iDv-uI/AAAAAAAAAuA/NPWgKtjtflE/s1600/rep.PNG"]ep[/URL]

the right mathematics can not cheat[/QUOTE]


[url]http://www.xkcd.com/747/[/url]

[url]http://www.xkcd.com/759/[/url]

[URL="http://4.bp.blogspot.com/_rvR3ouziO8g/TLxMpwU62TI/AAAAAAAAAtI/Vsj4-gSD0lM/s1600/cmdugo.PNG"]third element[/URL] of the first group,
advance of 30 and then tell us where to find last

CRGreathouse 2010-10-25 02:46

[QUOTE=3.14159;234302]Hopefully it would be accurate..[/QUOTE]

It should be possible to calculate a reasonable range in which the answer lies. I'm not sure how precise it would be possible to make it.

3.14159 2010-10-28 00:29

I think I can sieve by hand for k * b^n ± 1 numbers; I just need to do some modular arithmetic to get rid of each prime p which divides a certain k;

Ex: k * 6^90 + 1;

Step 1; 6^90 mod 7 = 1.

Step 2; Multiply by 6; 6 * 6^90 = 6 mod 7;

6^91 + 1 mod 7 = 0.

7 divides k = {6, 13, 20, 27, 34, 41, 48, 55, 62, 69, 76, 83, 90, 97, 104, 111, ...}

Repeat for 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, and.. as far as you'll be willing to strain your hand..

Although even a computer is approx. 10^4 to 10^8 times faster than you.

3.14159 2010-10-28 00:42

I'll use that whenever working with no computer assistance; Unless there's an easier method to go by.

science_man_88 2010-11-01 13:19

[CODE](10:15) gp > for(i=1,#mersenne,print(isprime(floor(prime(i)^(2*Pi)))))
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
(10:16) gp > for(i=1,#mersenne,print(isprime(floor(mersenne[i]^(2*Pi)))))
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
(10:16) gp >[/CODE]

can we use this fact that all the Mersenne exponents up to #39 work but not all primes up to prime 39 could this help narrow down exponents if used ?

though it may be doubtful these continue and it took me a long time to find a 1 in the prime[i]^(5*Pi) one and they seem to work for Mersenne exponents so far as well.

CRGreathouse 2010-11-01 14:02

[QUOTE=science_man_88;235192]can we use this fact that all the Mersenne exponents up to #39 work but not all primes up to prime 39 could this help narrow down exponents if used ?[/QUOTE]

No.

First, we don't know that the pattern will continue -- in fact I doubt it does. Second, there are way too many primes for which this will work -- heuristically, [TEX]\frac{n}{2\pi\log^2n}[/TEX] up to n.

But if you think you're onto something, try testing
49982629, 49985387, 49985521, 49986509, 49986737, 49987409, 49987759, 49988693, 49989491, ...
for which the relation holds and are toward the small end of the unassigned range. (I don't know which of these may already be taken, possibly all.)

science_man_88 2010-11-01 14:07

[QUOTE=CRGreathouse;235200]No.

First, we don't know that the pattern will continue -- in fact I doubt it does. Second, there are way too many primes for which this will work -- heuristically, [TEX]\frac{n}{2\pi\log^2n}[/TEX] up to n.

But if you think you're onto something, try testing
49982629, 49985387, 49985521, 49986509, 49986737, 49987409, 49987759, 49988693, 49989491, ...
for which the relation holds and are toward the small end of the unassigned range. (I don't know which of these may already be taken, possibly all.)[/QUOTE]

okay but how many of those fit for the 5*Pi one ? or the sequence of x*Pi that are :

2 5 8 10 11 12 14 15 18 all Mersenne exponents so far tested by me work to be 0 ( false) for all of these but not every prime will work in these relations.

science_man_88 2010-11-01 14:15

[QUOTE=CRGreathouse;235200]No.

First, we don't know that the pattern will continue -- in fact I doubt it does. Second, there are way too many primes for which this will work -- heuristically, [TEX]\frac{n}{2\pi\log^2n}[/TEX] up to n.

But if you think you're onto something, try testing
49982629, 49985387, 49985521, 49986509, 49986737, 49987409, 49987759, 49988693, 49989491, ...
for which the relation holds and are toward the small end of the unassigned range. (I don't know which of these may already be taken, possibly all.)[/QUOTE]

okay but how many of those fit for the 5*Pi one ? or the sequence of x*Pi that are :

2 5 8 10 11 12 14 15 18 all Mersenne exponents so far tested by me work to be 0 ( false) for all of these but not every prime will work in these relations.

f each combo of isprime !isprime etc. knocks out a prime when the range is up to 1,000,000 then this can knock out over 60% of prime exponents in the range (though it might slump lower) without a single LL test being done. if confirmed it would help a lot I think. not quite sure anymore got to figure it out there are 623530 combinations of taking the numbers i've found to work so far off isprime I believe if so then if each one eliminated just one exponent we could eliminate 623530 primes at least from ever needing LL test.

CRGreathouse 2010-11-01 15:55

[QUOTE=science_man_88;235209]okay but how many of those fit for the 5*Pi one ?[/QUOTE]

You have Pari/GP, feel free to test. Heuristically, if I understand what you're asking, about [TEX]\frac{x}{10\pi^2\log^4x}[/TEX] numbers up to x will have this property and the other.

But I emphasize that this is likely to remove as many (as a fraction) of future Mersenne exponents as Mersenne non-exponents.

Edit: I can't decide on what exponent to use in the denominator, but not less than 2 or more than 4.

science_man_88 2010-11-01 16:03

[QUOTE=CRGreathouse;235222]You have Pari/GP, feel free to test. Heuristically, if I understand what you're asking, about [TEX]\frac{x}{10\pi^2\log^4x}[/TEX] numbers up to x will have this property and the other.

But I emphasize that this is likely to remove as many (as a fraction) of future Mersenne exponents as Mersenne non-exponents.[/QUOTE]

yeah scrap 18 lol for use for Mersenne exponents yeah well it's based on properties of the rest of them to predict the next lol.

3.14159 2010-11-11 01:47

Here's a nice challenge for you all;

Express the number 988741313296003[sup]47[/sup] as a sum of prime powers, which meet these conditions:

1. No more than 25 primes may be used.

Ex: 2[sup]30[/sup] + 3[sup]18[/sup] + 5[sup]12[/sup] + 19[sup]6[/sup] + 293[sup]3[/sup] + 577[sup]3[/sup] + 191[sup]3[/sup] + 89[sup]3[/sup] + 223[sup]2[/sup] + 37[sup]2[/sup] + 11[sup]2[/sup] + 73[sup]1[/sup] + 2[sup]1[/sup] = 7[sup]11[/sup]

Or, for a harder one; No primes smaller than 250000 may be used, no primes may be used more than once, all exponents must be prime and used once and only once, and no more than 30 primes may be used.


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

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