OFFSET
1,2
COMMENTS
The number of zeros of n! base 3 is approaching n/2 as n grows. Similarly, the number of zeros of n! base 4 is approaching n/2 as n grows. Consequently, this sequence is expected to have high density.
From Robert Israel, Jan 19 2017: (Start)
Since typically A000120(n) ~ log_2(n) while typically A053735(n) ~ log_3(n), the density of this sequence should go to 0, contrary to the previous comment. (End)
Comment from N. J. A. Sloane, Dec 06 2019: (Start)
Appears to be the same as the list of positive numbers n such that the last nonzero digit of n! in base 12 belongs to the set [1, 2, 5, 7, 10, 11].
The first footnote in Deshouillers et al. (2016) says: "if the last nonzero digit of n! in base 12 belongs to {1, 2, 5, 7, 10, 11} then |(digit-sum of n in base 3) - (digit-sum of n in base 2)| is <= 1; this seems to occur infinitely many times." Compare A096288. (End)
REFERENCES
Jean-Marc Deshouillers, Laurent Habsieger, Shanta Laishram, Bernard Landreau, Sums of the digits in bases 2 and 3, arXiv:1611.08180, 2016
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
EXAMPLE
6! is 222200 in base 3 and 23100 in base 4, both of them have 2 zeros at the end, so 6 is in the sequence.
MAPLE
s2:= n -> convert(convert(n, base, 2), `+`):
s3:= n -> convert(convert(n, base, 3), `+`):
select(n -> s2(n) + (n+s2(n) mod 2) = s3(n), [$1..1000]); # Robert Israel, Jan 19 2017
MATHEMATICA
sntzQ[n_]:=Module[{f=n!}, Last[Split[IntegerDigits[f, 3]]]==Last[ Split[ IntegerDigits[ f, 4]]]]; Select[Range[200], sntzQ] (* Harvey P. Dale, Jul 11 2020 *)
PROG
(PARI) is(n)=my(L=log(n+1)); sum(k=1, L\log(3), n\3^k)==sum(k=1, L\log(2), n>>k)\2 \\ Charles R Greathouse IV, Oct 04 2012
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Tanya Khovanova, Oct 03 2012
EXTENSIONS
More terms from Alois P. Heinz, Oct 03 2012
STATUS
approved