login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A178969 Last nonzero decimal digit of (10^10^n)!. 1

%I #39 Dec 08 2023 09:36:31

%S 8,2,6,4,2,2,6,2,6,4,2,2

%N Last nonzero decimal digit of (10^10^n)!.

%C It is possible to find a(10) using the program from the second Bomfim link, or a similar GMP program. Reserve 312500001 words instead of 31250001. Use a computer with at least 6 GB of RAM. - _Washington Bomfim_, Jan 06 2011

%H W. Bomfim, <a href="http://oeis.org/w/images/6/61/Proof3.txt">A property of the last non-zero digit of factorials</a>.

%H W. Bomfim, <a href="http://oeis.org/w/images/c/cb/Program.txt">C program</a>.

%F From _Washington Bomfim_, Jan 06 2011: (Start)

%F a(n) = A008904(10^(10^n)).

%F a(n) = A008904(2^(10^n)), if n > 0.

%F For n >= 1, with N = 10^n, 2^N represented in base 5 as (a_h, ..., a_0)_5, t = Sum_{i = h, h-1, ..., 0} (a_i even), x = Sum_{i = h, h-1, ..., 1}(Sum_{k = h, h-1, ..., i} (a_i)), z = (x + t/2) mod 4, and y = 2^z, a(n) = 6*(y mod 2) + y*(1 - (y mod 2)).

%F (End)

%t f[n_] := Mod[6Times @@ (Rest[FoldList[{ 1 + #1[[1]], #2!2^(#1[[1]]#2)} &, {0, 0}, Reverse[IntegerDigits[n, 5]]]]), 10][[2]]; (* _Jacob A. Siehler_ *) Table[ f[10^10^n], {n, 0, 4}]

%o (PARI)\\ L is the list of the N digits of 2^(10^n) in base 5.

%o \\ With 2^(10^n) in base 5 as (a_h, ... , a_0)5,

%o \\ L[1] = a_0, ... ,L[N] = a_h.

%o convert(n)={n=2^(10^n); x=n; N=floor(log(n)/log(5)) + 1;

%o L = listcreate(N);

%o while(x, n=floor(n/5); r=x-5*n; listput(L, r); x=n; );

%o L; N

%o };

%o print("0 8"); for(n=1, 5, print1(n, " "); convert(n); q=0; t=0; x=0; forstep(i=N, 2, -1, a_i=L[i]; q+=a_i; x+=q; t+=a_i*(1-a_i%2); ); a_i=L[1]; t+=a_i*(1-a_i%2); z=(x+t/2)%4; y=2^z; an=6*(y%2)+y*(1-(y%2)); print(an)); \\ _Washington Bomfim_, Jan 06 2011

%o (Python)

%o from functools import reduce

%o from sympy.ntheory.factor_ import digits

%o def A178969(n): return reduce(lambda x,y:x*y%10,(((6,2,4,8,6,2,4,8,2,4,8,6,6,2,4,8,4,8,6,2)[(a<<2)|(i*a&3)] if i*a else (1,1,2,6,4)[a]) for i, a in enumerate(sympydigits(1<<10**n,5)[-1:0:-1])),6) if n else 8 # _Chai Wah Wu_, Dec 07 2023

%Y Cf. A008904, A173670

%K nonn,base,more

%O 0,1

%A _Robert G. Wilson v_, Jan 01 2011

%E a(9) from _Washington Bomfim_, Jan 06 2011

%E a(10) from _Chai Wah Wu_, Dec 07 2023

%E a(11) from _Chai Wah Wu_, Dec 08 2023

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 00:26 EDT 2024. Contains 371264 sequences. (Running on oeis4.)