PARI's commands - Page 81

CRGreathouse · 2010-08-25, 19:25

Quote:

Originally Posted by 3.14159

Code:

41441994149199491949199941949494914919499419411441141114499941994199

68 digits, and only has square numbers for digits!

Update:

Code:

41449199414911994941149149491999144449499419114949194149949499414499419449919999191449199941949494914919499419411441141114499941994199

134 digits, and only composed of square numbers! HA!

Hmm, I wonder how rare such numbers are. Heuristics suggest about
$\frac43\cdot\frac{3^n}{n\log10}$
n-digit examples; how well does this hold in practice?

There are 0, 3, 11, 12, 43, 94, 239, 566, 1710 such primes with 1, 2, ..., 9 digits. The formula predicts 2, 3, 5, 12, 28, 70, 181, 475, 1266 primes; not too bad, perhaps, though biased on the low side. Maybe with the next correction term, the -1 in the denominator? That gives 3, 3, 6, 13, 31, 76, 193, 502, 1331.

This is an asymptotic density of something like n^0.477, so maybe somewhat rarer than palindromatic primes.

CRGreathouse · 2010-08-25, 19:26

Quote:

Originally Posted by science_man_88

okay well since these other seem to fit the same relation and are all of form 2x+1,2(2x+1)+1,etc they should all deviate the same way in my eyes so what would the odds for an exception to give you a sequence that works in the same way as the original being inserted into 2^p-1 with p prime will actually be prime ?. can you give those odds ?

I don't understand.

CRGreathouse · 2010-08-25, 19:32

Quote:

Originally Posted by science_man_88

if only there was a pattern you'd find a 266 digit one.

Here's a 266-digit one:

Code:

11111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111494999

science_man_88 · 2010-08-25, 19:32

Code:

(14:37) gp > findrec([11,23,47,95,191,383,767])
Recurrence relation is a(n) = 3a(n-1) - 2a(n-2).
3 d.f.
%63 = [3, -2]~
(14:39) gp > findrec([29,59,119,239,479])
Recurrence relation is a(n) = 3a(n-1) - 2a(n-2).
1 d.f.
%66 = [3, -2]~
(14:39) gp > findrec([1,3,7,15,31,63,127])
Recurrence relation is a(n) = 3a(n-1) - 2a(n-2).

what are the odds that 2^ (the primes in the first 2)-1 are prime (which would put them in the primes of the last sequence). answer this and not only can we predict Mersenne primes but double Mersenne primes.

science_man_88 · 2010-08-25, 19:34

Quote:

Originally Posted by CRGreathouse

Here's a 266-digit one:

Code:

11111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111494999

is there a 530 digit one ? if so it would possibly be part of this idea i have lol.

CRGreathouse · 2010-08-25, 19:49

Quote:

Originally Posted by science_man_88

is there a 530 digit one ? if so it would possibly be part of this idea i have lol.

If you look at post #881, we expect lots of 530-digit ones -- about 8 × 10²⁴⁹. The first is

Code:

11111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111491141

CRGreathouse · 2010-08-25, 19:53

Quote:

Originally Posted by science_man_88

Code:

(14:37) gp > findrec([11,23,47,95,191,383,767])
Recurrence relation is a(n) = 3a(n-1) - 2a(n-2).
3 d.f.
%63 = [3, -2]~
(14:39) gp > findrec([29,59,119,239,479])
Recurrence relation is a(n) = 3a(n-1) - 2a(n-2).
1 d.f.
%66 = [3, -2]~
(14:39) gp > findrec([1,3,7,15,31,63,127])
Recurrence relation is a(n) = 3a(n-1) - 2a(n-2).

what are the odds that 2^ (the primes in the first 2)-1 are prime (which would put them in the primes of the last sequence). answer this and not only can we predict Mersenne primes but double Mersenne primes.

I don't understand your question at all. If I wanted to test if a small number was the exponent of a Mersenne prime, I would check if it's in the listed members of A000043 (by a binary search). If I wanted to test if a large number was a Mersenne exponent, I'd check if it was a prime, trial-divide, look for some small factors with ECM, and then do a Lucas-Lehmer test on it. I don't have a special test for numbers of a certain form unless I can show an algebraic factor, but that would be caught by my tests anyway. It wouldn't speed up the search at all.

science_man_88 · 2010-08-25, 19:56

if we prove all exceptions that aren't knocked out by others have a sequence of the same form and that series of that form starting with an exception can't have a 2^p-1 prime then we can knock out all primes that have form of that have sequence that start at an exception.

CRGreathouse · 2010-08-25, 20:11

Quote:

Originally Posted by science_man_88

if we prove all exceptions that aren't knocked out by others

What's an exception? What does "knocked out" mean? What are "others"?

I assume that "knocked out" means that, in some sense, the number is guaranteed to not be the exponent of a Mersenne prime. But what determines if it's knocked out or not?

Quote:

Originally Posted by science_man_88

have a sequence of the same form and that series of that form starting with an exception can't have a 2^p-1 prime

The same form as what?

Quote:

Originally Posted by science_man_88

then we can knock out all primes that have form of that have sequence that start at an exception.

Without definitions for the terms, this doesn't tell me much of anything.

But let me try something. Let's say that a number is "knocked out" if it is a member of one of the sequences S1, S2, S3, ..., Sk. Further, let's suppose that the sequences are, like A055010, exponential integer sequences, where each term is roughly twice the one before it. And let's look at large numbers, say exponents at least a million.

Suppose we want to find all Mersenne exponents between 1,000,000 and 1,100,000 by this method. (We know from Slowinski & Gage that there are none, but we're testing this!) Now, by our assumptions, a given sequence Sj will remove about lg (1100000/1000000) ≈ 0.13 primes from the list. Since there are about 100,000 members, we'll need a minimum of > 720,000 sequences to clear the list.

Not one sequence or three sequences; not even a hundred. Three-quarters of a million sequences! And that's just to clear a small part of a small range; a larger one would take many more. Further, if there's overlap, you'll need more.

So if you don't want to use millions of sequences, make sure your method doesn't follow my assumptions.

science_man_88 · 2010-08-25, 21:00

I realized that my idea takes a lot but if we can come up with something it's like a sieve.

3.14159 · 2010-08-25, 22:13

Another try at predicting the next Mersenne? I think there are figures that place it at around 19M to 20M digits. I believe this was stated somewhere in the Prime Pages.

2010-08-25, 19:32	#884
science_man_88 "Forget I exist" Jul 2009 Dumbassville 2⁶×131 Posts	Code: (14:37) gp > findrec([11,23,47,95,191,383,767]) Recurrence relation is a(n) = 3a(n-1) - 2a(n-2). 3 d.f. %63 = [3, -2]~ (14:39) gp > findrec([29,59,119,239,479]) Recurrence relation is a(n) = 3a(n-1) - 2a(n-2). 1 d.f. %66 = [3, -2]~ (14:39) gp > findrec([1,3,7,15,31,63,127]) Recurrence relation is a(n) = 3a(n-1) - 2a(n-2). what are the odds that 2^ (the primes in the first 2)-1 are prime (which would put them in the primes of the last sequence). answer this and not only can we predict Mersenne primes but double Mersenne primes. Last fiddled with by science_man_88 on 2010-08-25 at 19:32

2010-08-25, 22:13	#891
3.14159 May 2010 Prime hunting commission. 2⁴·3·5·7 Posts	Another try at predicting the next Mersenne? I think there are figures that place it at around 19M to 20M digits. I believe this was stated somewhere in the Prime Pages. Last fiddled with by 3.14159 on 2010-08-25 at 22:13

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Why do I sometimes see all the <> formatting commands when I quote or edit?	cheesehead	Forum Feedback	3	2013-05-25 12:56
Passing commands to PARI on Windows	James Heinrich	Software	2	2012-05-13 19:19
Ubiquity commands	Mini-Geek	Aliquot Sequences	1	2009-09-22 19:33
64-bit Pari?	CRGreathouse	Software	2	2009-03-13 04:22
Are these commands correct?	jasong	Linux	2	2007-10-18 23:40

2010-08-25, 19:56	#888
science_man_88 "Forget I exist" Jul 2009 Dumbassville 2⁶·131 Posts	if we prove all exceptions that aren't knocked out by others have a sequence of the same form and that series of that form starting with an exception can't have a 2^p-1 prime then we can knock out all primes that have form of that have sequence that start at an exception.

2010-08-25, 21:00	#890
science_man_88 "Forget I exist" Jul 2009 Dumbassville 10000011000000₂ Posts	I realized that my idea takes a lot but if we can come up with something it's like a sieve.