OFFSET
2,1
COMMENTS
Based on function solution of: Solve[(Prime[n]+a*PrimePi[n])/(1+a)-Sqrt[Prime[n]*PrimePi[n]]=0,a] which gives a log type function a[n].
FORMULA
b(n) = Floor[Exp[(Prime[n] - Sqrt[(Prime[n]* PrimePi[n]]/(-PrimePi[n] + Sqrt[Prime[n]* PrimePi[n]]]]
MATHEMATICA
digits=200 b = Table[Floor[Exp[(Prime[n] - Sqrt[(Prime[n]* PrimePi[n]]/(-PrimePi[n] + Sqrt[Prime[n]* PrimePi[n]]]], {n, 2, digits}]
CROSSREFS
KEYWORD
nonn,uned
AUTHOR
Roger L. Bagula, Jan 05 2004
STATUS
approved