mersenneforum.org  

Go Back   mersenneforum.org > Extra Stuff > Blogorrhea > enzocreti

Reply
 
Thread Tools
Old 2018-06-26, 13:42   #45
Dr Sardonicus
 
Dr Sardonicus's Avatar
 
Feb 2017
Nowhere

10010001000112 Posts
Default

Quote:
Originally Posted by enzocreti View Post
I am astonished that nobody knows something about this sequence.
Au contraire, there is a lot known about both the whole sequence, and the primes embedded therein. For example: none of the terms are divisible either by 2 or 3. The terms for which n is congruent to 1 (mod 4) are divisible by 5.

For any prime p > 5, the proportion of terms divisible by p can be found (though which terms are divisible by p is not easily predicted).

The rate of growth of the sequence is fairly clear.

There is a formula to calculate the terms of the sequence.

True, the sequence has no obvious divisibility properties like a^n - 1.

And, alas, the sequence has no obvious theoretical importance, so there is little impetus to throw a lot of resources at determining which terms are prp's.

I have no sense of how many prp's one might "expect" to find for n <= X for large X.

Quote:
Only a person told me that this sequence is NOT random, but he didn't explain why it is NOT random...I am very angry for that
I suggest you reread this post to this thread. Long story short-- the given sequence isn't random. The sequence of primes isn't random. So their intersection isn't random.

BTW, since the constant log(2)/log(10) = log102 plays a role in defining the terms of this sequence, I will mention that this number is not only irrational (easily proved by invoking unique factorization), but transcendental (easily proved by invoking the Gelfond-Schneider Theorem).

Last fiddled with by Dr Sardonicus on 2018-06-26 at 13:44 Reason: fixing typos
Dr Sardonicus is offline   Reply With Quote
Old 2018-06-26, 13:56   #46
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

597910 Posts
Default

Quote:
Originally Posted by enzocreti View Post
Only a person told me that this sequence is NOT random, but he didn't explain why it is NOT random...I am very angry for that
What do you think that it means to be random?
CRGreathouse is offline   Reply With Quote
Old 2018-06-26, 14:10   #47
enzocreti
 
Mar 2018

2×5×53 Posts
Default

Quote:
Originally Posted by Dr Sardonicus View Post
Au contraire, there is a lot known about both the whole sequence, and the primes embedded therein. For example: none of the terms are divisible either by 2 or 3. The terms for which n is congruent to 1 (mod 4) are divisible by 5.

For any prime p > 5, the proportion of terms divisible by p can be found (though which terms are divisible by p is not easily predicted).

The rate of growth of the sequence is fairly clear.

There is a formula to calculate the terms of the sequence.

True, the sequence has no obvious divisibility properties like a^n - 1.

And, alas, the sequence has no obvious theoretical importance, so there is little impetus to throw a lot of resources at determining which terms are prp's.

I have no sense of how many prp's one might "expect" to find for n <= X for large X.


I suggest you reread this post to this thread. Long story short-- the given sequence isn't random. The sequence of primes isn't random. So their intersection isn't random.

BTW, since the constant log(2)/log(10) = log102 plays a role in defining the terms of this sequence, I will mention that this number is not only irrational (easily proved by invoking unique factorization), but transcendental (easily proved by invoking the Gelfond-Schneider Theorem).
What is the formula for calculating the terms of the sequence?
enzocreti is offline   Reply With Quote
Old 2018-06-26, 14:19   #48
enzocreti
 
Mar 2018

10228 Posts
Default

log(2)/log(10) = log102 how does it play a role in the terms of the sequence?
enzocreti is offline   Reply With Quote
Old 2018-06-26, 14:24   #49
enzocreti
 
Mar 2018

2×5×53 Posts
Default

Why this sequence there wasn't in Oeis if you say that it is yet known?
enzocreti is offline   Reply With Quote
Old 2018-06-26, 14:36   #50
enzocreti
 
Mar 2018

2·5·53 Posts
Thumbs down

Quote:
Originally Posted by CRGreathouse View Post
What do you think that it means to be random?
These numbers seem to be not random at all...why residue 6 didn't yet occur? Why five in a row probable primes with residue 5? Are that all coincidences???
enzocreti is offline   Reply With Quote
Old 2018-06-26, 14:38   #51
enzocreti
 
Mar 2018

2·5·53 Posts
Default

Why we found 3 probable primes with residue 1 when the residue 1 occurs with half the frequency of the residue 6?
enzocreti is offline   Reply With Quote
Old 2018-06-26, 14:47   #52
Dr Sardonicus
 
Dr Sardonicus's Avatar
 
Feb 2017
Nowhere

4,643 Posts
Default

Quote:
Originally Posted by enzocreti View Post
log(2)/log(10) = log102 how does it play a role in the terms of the sequence?
The sequence in which you're looking for primes, as given here, is

f(k) = 10^m * (2^k - 1) + 2^(k-1) - 1,

where m is the number of decimal digits in 2^(k-1).

The value of m is

m = 1 + floor((k-1)*log10(2)), i.e.

m = 1 + floor((k-1)*log(2)/log(10)).
Dr Sardonicus is offline   Reply With Quote
Old 2018-06-26, 18:52   #53
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3·1,993 Posts
Default

Quote:
Originally Posted by enzocreti View Post
These numbers seem to be not random at all...why residue 6 didn't yet occur? Why five in a row probable primes with residue 5? Are that all coincidences???
Probably.
http://mathworld.wolfram.com/StrongL...llNumbers.html
https://www.ime.usp.br/~rbrito/docs/2322249.pdf
CRGreathouse is offline   Reply With Quote
Old 2018-06-27, 13:06   #54
Dr Sardonicus
 
Dr Sardonicus's Avatar
 
Feb 2017
Nowhere

4,643 Posts
Default (non) periodicity mod p, p > 5

Quote:
Originally Posted by science_man_88 View Post
Mod 7 it repeats a pattern every 60 entries. [3,3,6,28,7,6,7] is the distribution mod 7
The sequence

f(k) = 10^m*(2^k - 1) + 2^(k-1) - 1, where

m = 1 + floor((k-1)*log(2)/log(10)), the number of decimal digits in 2^(k-1),

is 1 (mod 2) for n > 1, and 1 (mod 3) for all n.

Since 10 is divisible by 5, the sequence is also periodic (mod 5) with period 5-1 = 4, the sequence of remainders (mod 5) being 0, 1, 3, 2, 0, 1, 3, 2,...

For p > 5, the sequence is not periodic (mod p), but there is one specific remainder, namely (p - 1)/2, that does recur when n is divisible by the multiplicative order of 2 (mod p), which multiplies the factor 10^m (mod p) by 0 (mod p), taking it out of consideration. (I am too lazy to work out whether this remainder can also occur at other values of n.)

The reason the sequence is not periodic (mod p) for any p > 5 is, first of all, the terms 2^k - 1 and 2^(k-1) - 1 are periodic (mod p), with period equal to the multiplicative order h of 2 (mod p). The non-periodicity of f(k) comes from the non-periodicity of m (mod h). (This problem disappears when p divides 2^k - 1, as noted above.)

To show that m is not periodic modulo any integer N > 1, we show that, for any given positive positive irrational constant c, and any given positive integer L,

floor((n + L)*c) - floor(n*c) is not constant. We use frac(x) to denote x - floor(x).

Clearly, (n + L)*c = floor(n*c) + frac(n*c) + floor(L*c) + frac(L*c), so that

floor((n+L)*c) = floor(n*c) + floor(L*c) if frac(n*c) + frac(L*c) < 1, and floor(n*c) + floor(L*c) + 1, if frac(n*c) + frac(L*c) > 1.

Now it is a well known result that frac(n*c) is uniformly distributed in (0, 1), so both possibilities occur (and, given the value of frac(L*c), we can say how often). Thus,

floor((n + L)*c) - floor(n*c) is not constant, so m is not periodic modulo any integer N greater than 1. Hence f(k) is not periodic (mod p) for any prime p > 5.
Dr Sardonicus is offline   Reply With Quote
Old 2018-06-27, 14:09   #55
enzocreti
 
Mar 2018

2×5×53 Posts
Default

And so residue 6 mod 7 cannot occur?
enzocreti is offline   Reply With Quote
Reply



Similar Threads
Thread Thread Starter Forum Replies Last Post
so who got the $25,000 cut of the $100k EFF prize? ixfd64 Lounge 10 2012-09-22 15:55
Name the trolls and win a prize. Xyzzy Lounge 42 2010-03-08 07:06
EFF prize and error rate S485122 PrimeNet 15 2009-01-16 11:27
EFF Prize? Unregistered Information & Answers 73 2007-08-11 11:38
GIMPS' Prize. T.Rex Math 8 2007-03-13 10:59

All times are UTC. The time now is 04:28.


Sat Jul 17 04:28:15 UTC 2021 up 50 days, 2:15, 1 user, load averages: 2.05, 2.25, 2.31

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

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.