Collatz Sequences for
p
n
,
n!
, and
♯(n)

Paul Abbott
School of Physics, Mathematics and Computing
The University of Western Australia
paul.c.abbott@uwa.edu.au

Collatz sequence


Primes


Factorial


Primorial

The primorial
♯(n)
is defined as the product of first
n
prime numbers:
In[]:=
♯(n_):=
n
∏
k=1
p
k
Compute
♯(n)
for
1≤n≤10
:
In[]:=
♯/@Range[10]
Out[]=
{2,6,30,210,2310,30030,510510,9699690,223092870,6469693230}
By construction, each primorial is divisible by only one power
2
.
Compute the sequence length for
♯(n)
versus
log(♯(n))
for
1≤n≤1000
:
In[]:=
data=(n{Log[♯[n]],Length@cs[♯[n]]})/@Range[1000];
Display the data:
In[]:=
plp=ListPlot[data,ImageSizeLarge]
Out[]=
2000
4000
6000
10000
20000
30000
40000
50000
A linear fit works well:
In[]:=
fit=Fit[data,x,x]
Out[]=
6.95456x
In[]:=
Show[plp,Plot[fit,{x,0,8000},PlotStyleRed]]
Out[]=
2000
4000
6000
10000
20000
30000
40000
50000