OFFSET
2,1
MATHEMATICA
(*first entropy constant to n=10000*)
p0[n_]=PrimePi[n]/n
a=Table[ -p0[n]*Log[p0[n]], {n, 2, 10000}];
w0=N[Apply[Plus, a], 200]/10000
(*Second entropy constant to n=10000*)
z0=Sum[PrimePi[n]/(n*Prime[n]), {n, 1, 10000}];
(* semi-empirical prime-like function based on prime entropy sums*)
p[n_]=Sum[ -Log[PrimePi[i]/i], {i, 2, n}]
a0=Table[If[PrimeQ[Floor[p[n]*(2+E+z0/n+w0)]]==True, Floor[p[n]*(2+E+z0/n+w0)], 0], {n, 2, 1000}];
c0=Delete[Union[a0], 1]
CROSSREFS
KEYWORD
nonn,uned,obsc,less
AUTHOR
Roger L. Bagula, Dec 03 2003
STATUS
approved