mersenneforum.org Primes in order in pi
 Register FAQ Search Today's Posts Mark Forums Read

2022-02-04, 15:54   #1
robert44444uk

Jun 2003
Suva, Fiji

204010 Posts
Primes in order in pi

I am a little surprised I can't find any mention on the web, of the sequence of primes, in order, that might be found in the decimal expansion of pi (3.1415926535..)

The first instance of 2 is at the 6th decimal position, the first 3 thereafter at the 9th position, the first 5 thereafter at the 10th, etc.

For definition's sake, the first 71 is not at the 1757th position following 67 at the 1756th, as the 1757th position was also used to ascertain 67. Instead, 71 is found at the 1923rd.

A file of the decimal positions of all primes up to 8713 which are to be found in order in the first 10,000,000 decimal positions of pi is attached.

A graph is also shown, showing the positions of each prime, supporting the hypothesis that pi is normal. Some of the results are shown below

Code:
14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196442881097566593344612847564823378678316527...
The calculation is trivial. My average 7 year old desktop carried out the following perl program in less than 10 seconds.

Code:
#!/usr/bin/env perl
use warnings;
use strict;
use Math::Prime::Util qw/:all/;
use feature ':5.10';

$|=1; my$pie = substr(Pi(10_000_000),2);
my $maxprime = 100_000; my$prevplace = 0;
forprimes {
my $i = 1+index($pie,$_,$prevplace);
say $_, " ",$i;
$prevplace =$i+length($_)-1; if ($i==0){
last;
}
} 2,\$maxprime;
So...the challenge (as this is the puzzle subsection) is to expand the series, and perhaps find quite large gaps or close instances..and post to OEIS
Attached Thumbnails

Attached Files
 pieres.txt (14.3 KB, 61 views)

 2022-02-04, 16:29 #2 Uncwilly 6809 > 6502     """"""""""""""""""" Aug 2003 101×103 Posts 24×33×52 Posts That inflection point around 1000 looks interesting. Any insight into that?
2022-02-04, 17:04   #3
robert44444uk

Jun 2003
Suva, Fiji

111111110002 Posts

Quote:
 Originally Posted by Uncwilly That inflection point around 1000 looks interesting. Any insight into that?
That's the point at which the next prime has 4 digits. There would be similar inflections at 2 to 3 and 1 to 2 digits as well.

When I inferred that pi is normal in the OP, I meant pi behaves as a normal number.

Last fiddled with by robert44444uk on 2022-02-04 at 17:48

2022-02-05, 01:48   #4
newalex

Feb 2019

24 Posts

Here is my result for 100b digits of Pi.
All primes up to 1042333 were found in order in these 100b digits.
Attached Files
 primesinpi.zip (463.5 KB, 62 views)

 2022-02-05, 14:50 #5 Dr Sardonicus     Feb 2017 Nowhere 2×11×277 Posts Assume pi is a normal number to base ten. It occurred to me to wonder how many times the block of decimal digits representing a prime p can reasonably be expected to occur before it occurs in its proper place in the sequence of primes.
2022-02-05, 14:59   #6
kruoli

"Oliver"
Sep 2017
Porta Westfalica, DE

23068 Posts

Quote:
 Originally Posted by newalex Here is my result for 100b digits of Pi. All primes up to 1042333 were found in order in these 100b digits.
Would you mind sharing/PM'ing your program? I could run it on 1t digits I have. And maybe some other constants.

2022-02-05, 17:05   #7
newalex

Feb 2019

24 Posts

Quote:
 Originally Posted by kruoli Would you mind sharing/PM'ing your program? I could run it on 1t digits I have. And maybe some other constants.
Yes, of course. But I am not sure if it would be suitable for you. It is written in Java and stores all digits in memory. So you need 1 Tb of RAM to run it for 1t. It can be fixed easily but need some work.
Here is the source code
Code:
import java.io.*;

public class PrimesInPi {

public static void main(String[] args) throws IOException {
long piSize = 100000000000L;
final int maxArraySize = 1000000000;
int chunkSize = (piSize>maxArraySize) ? maxArraySize : (int) piSize;
char[][] pi;
pi = new char[(int) (piSize/chunkSize)][chunkSize];
for (int i=0; i<piSize/chunkSize; i++) {
}
fr.close();
FileWriter result = new FileWriter("output.txt");
long pos = 0;
int i=0, j=0;
while (pos<piSize) {
if (isCurrentPrime(pi, prime, i, j)) {
result.write(prime+" "+(pos+1)+"\n");
j+=prime.length();
pos+=prime.length();
if (j>=chunkSize) {
j-=chunkSize;
i++;
}
} else {
pos++;
j++;
if (j==chunkSize) {
j=0;
i++;
}
}
}
result.close();
}

private static boolean isCurrentPrime(char[][] pi, String prime, int i, int j) {
boolean result = true;
int i1 = i, j1 = j;
for (int k=0; k<prime.length(); k++) {
if (prime.charAt(k) == pi[i1][j1]) {
j1++;
if (j1 == pi[i1].length) {
j1 = 0;
i1++;
}
if (i1 == pi.length) {
result = false;
break;
}
} else {
result = false;
break;
}
}
return result;
}
}
The program expects to have pi.txt and primes.txt files in the same folder. It ignores first 2 symbols in pi.txt file to start from fraction part(assuming the file starts with 3. ). primes.txt is a list of primes(one line - one prime).
piSize variable should set to count of pi digits which we want to process.

 2022-02-05, 17:16 #8 kruoli     "Oliver" Sep 2017 Porta Westfalica, DE 23068 Posts Thank you! For some reason I was expecting that the straightforward approach would be inefficient but you are correct that this approach can be applied to a file that is not fully stored in memory. I'll code up something later. I like your idea to have the primes precalculated in a file. This makes the code way less complex.
2022-02-05, 17:59   #9
uau

Jan 2017

23·19 Posts

Quote:
 Originally Posted by Dr Sardonicus It occurred to me to wonder how many times the block of decimal digits representing a prime p can reasonably be expected to occur before it occurs in its proper place in the sequence of primes.
Each sequence of n digits should occur equally often. Thus if you sequentially search for 50 6-digit sequences, you should on average expect to have seen 50 of the last sequence - on average encountering it once before finding each of the earlier 49 ones.

So the expected number should be (number of earlier equally long target sequences) + (number of one earlier one digit shorter target sequences)/10 + (two shorter)/100 and so on.

 2022-02-07, 09:15 #10 robert44444uk     Jun 2003 Suva, Fiji 204010 Posts You have to love the projected result of this exercise - there is exists within the decimal representation of pi, (which is itself just one of an infinite number of irrationals), all of the infinity of primes, in order. I would expect that this projected result will be the same for all integer bases, i.e. the representations of the primes in any integer base can be found in order within the representation of pi in that base. I can almost visualise this must be true for base 2. Can we go further to project or hypothesise that the primes, as depicted in any integer base are in order in the representation of any irrational number in the same integer base? As integers represented in irrational bases appear themselves irrational, then we can safely say that the primes in order, represented in irrational base form, cannot be found in the representation of the irrational in that base, as the irrational in its own base will have a simple non-irrational form. Because the representation of the irrational form is infinite, I would hypothesise that the primes represented in irrational base form (which I think are irrational representations themselves) are not included in order in the representation of an irrational number presented in an different rational or irrational base. Have I covered all the permutations here? PS: My teacher at school did day "don't mess with infinity, it bites"
2022-02-07, 12:11   #11
uau

Jan 2017

100110002 Posts

Quote:
 Originally Posted by robert44444uk You have to love the projected result of this exercise - there is exists within the decimal representation of pi, (which is itself just one of an infinite number of irrationals), all of the infinity of primes, in order.
I don't find it too interesting myself. Any random sequence would eventually produce a match. So would a sequence like 1 2 3 4 5 ... (you would get more matches arbitrarily late as parts of larger numbers).

Quote:
 Can we go further to project or hypothesise that the primes, as depicted in any integer base are in order in the representation of any irrational number in the same integer base?
No, you can easily construct an irrational number whose decimal representation consists of only 0 and 1.

 Similar Threads Thread Thread Starter Forum Replies Last Post bur Miscellaneous Math 9 2021-05-26 13:53 bhelmes Math 5 2018-09-13 16:23 Brain Miscellaneous Math 1 2010-12-08 01:00 Raman Puzzles 6 2010-09-05 17:43 Xyzzy Forum Feedback 5 2007-01-28 10:35

All times are UTC. The time now is 00:11.

Thu Dec 1 00:11:06 UTC 2022 up 104 days, 21:39, 0 users, load averages: 0.98, 1.02, 0.99