OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
f2nzdQ[n_]:=Module[{nf=n!}, Mod[nf/10^IntegerExponent[nf], 100]==16]; Select[ Range[1000], f2nzdQ] (* Harvey P. Dale, May 28 2015 *)
PROG
(PARI) last(n, k)=if(n*8\9-1<k, n=n!; return(n/10^valuation(n, 10)%10^k)); my(m=Mod(1, 5^k)); for(i=2, n, m*=i/10^valuation(i, 5)); lift(chinese(m, Mod(0, 2^k))) \\ Gives the last k decimal digits of n!.
is(n)=last(n, 2)==16 \\ Charles R Greathouse IV, Oct 22 2014
(PARI) is(n)=my(m=Mod(1, 25)); for(i=2, n, m*=i/10^valuation(i, 5)); m==16 \\ Charles R Greathouse IV, Oct 22 2014
(PARI) v=List(); m=Mod(1, 25); for(n=2, 1e6, m*=n/10^valuation(n, 5); if(m==16, listput(v, n); if(#v==100, return(Vec(v))))) \\ Charles R Greathouse IV, Oct 22 2014
(Python)
from math import factorial
def ok(n): return str(factorial(n)).rstrip("0").endswith("16")
print([k for k in range(1103) if ok(k)]) # Michael S. Branicky, Dec 29 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
STATUS
approved