OFFSET
1,1
FORMULA
Let a(1)=2, S(1)=log(2), n=1; (*) n=n+1 a(n)=S(n-1)+(Pi)*i^(Pi/6) ( Pi=3.14159....), S(i)=S(i-1)+log(a(i)), a(i)=int(a(i)), goto (*).
EXAMPLE
a(i)=p(i) for i= 1,81,101,109,142,169,286,334,446,479,624,846,959,1019,1034....
pi(10) = 4 whereas sequence gives 3
pi(100) = 25 whereas sequence gives 24
pi(1000) = 168 whereas sequence gives 167
pi(10000) = 1229 whereas sequence gives 1226
pi(100000) = 9592 whereas sequence gives 9586
pi(1000000) = 78498 whereas sequence gives 78524
pi(10000000) = 664579 whereas sequence gives 664668
pi(100000000) = 5761455 whereas sequence gives 5761527
pi(1000000000) = 50847534 whereas sequence gives 50847381
pi(10000000000) = 455052511 whereas sequence gives 455050443
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Jun 02 2004
STATUS
approved