OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
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)==28 \\ 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==28 \\ 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==28, listput(v, n); if(#v==100, return(Vec(v))))) \\ Charles R Greathouse IV, Oct 22 2014
CROSSREFS
KEYWORD
nonn,base
AUTHOR
EXTENSIONS
More terms from Robert G. Wilson v, Dec 14 2000
STATUS
approved