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
a(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
