login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A202021
The leading digit of (10^n)!.
1
1, 3, 9, 4, 2, 2, 8, 1, 1, 9, 2, 3, 1, 2, 1, 1, 1, 1, 5, 2, 1, 5, 1, 1, 3, 5, 3, 9, 1, 1, 6, 7, 7, 6, 3, 4, 1, 9, 9, 3, 2, 1, 2, 6, 6, 1, 2, 3, 5, 1, 5, 2, 5, 1, 1, 5, 8, 2, 7, 3, 4, 1, 1, 5, 5, 2, 3, 1, 8, 1, 8, 9, 1, 6, 3, 1, 4, 6, 4, 1, 8, 1, 1, 9, 1, 4, 8, 8, 8, 9, 1, 3, 3, 2, 1, 5, 4, 2, 3, 3, 1, 1, 4, 6, 6
OFFSET
0,2
COMMENTS
I employed R. Wm. Gosper's approximation (A090583).
LINKS
Eric Weisstein's World of Mathematics, Stirling's Approximation.
EXAMPLE
(10^1)! = 3628800 begins with 3.
(10^6)! begins with 8 and (10^100)! begins with 1.
MATHEMATICA
f[n_] := IntegerPart[ 10^FractionalPart[ N[(n*Log[n] - n + (1/2) Log[2 Pi*n + 1/3])/Log[10], 150]]]; f[1] = 1; Table[ f[10^n], {n, 0, 104}]
PROG
(PARI) a(n)=my(g=lngamma(10^n+1)/log(10)); g-=g\1; 10^g\1 \\ Charles R Greathouse IV, Jan 09 2013
CROSSREFS
Sequence in context: A125143 A200012 A130701 * A342220 A197507 A264992
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Jan 09 2013
STATUS
approved