OFFSET
1,1
COMMENTS
First term prime(2) = 3 is placed on 0th row.
If prime(n-1) = +1 mod 4 is on k-th row then we put prime(n) on (k-1)-st row.
If prime(n-1) = -1 mod 4 is on k-th row then we put prime(n) on (k+1)-st row.
This process makes an array of prime numbers:
3, 7, 19, 43, ....0th row
5, 11, 17, 23, 31, 41, 47, 59, 67, 103, 127, ....first row
13, 29, 37, 53, 61, 71, 79, 101, 107, 113 ....2nd row
73, 83, 97, 109, ....3rd row
89, ....4th row
LINKS
Robert Israel, Table of n, a(n) for n = 1..811
Andrew Granville and Greg Martin, Prime number races, Amer. Math. Monthly, 113 (No. 1, 2006), 1-33.
MAPLE
c1:= 0; c3:= 0: p:= 2: count:= 0: Res:= NULL:
while count < 100 do
p:= nextprime(p);
if c1 = c3 - 2 then
count:= count+1;
Res:= Res, p;
fi;
if p mod 4 = 1 then c1:=c1+1
else c3:= c3+1
fi
od:
Res; # Robert Israel, Nov 07 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Yasutoshi Kohmoto, Aug 12 2004
EXTENSIONS
More terms from Joshua Zucker, May 03 2006
STATUS
approved