OFFSET
1,4
COMMENTS
a(n)=1 iff n is squarefree.
LINKS
Robert Israel, Table of n, a(n) for n = 1..3315
FORMULA
a(n) = binomial(n, A007947(n)).
MAPLE
f:= proc(n) local k;
k:= convert(numtheory:-factorset(n), `*`);
binomial(n, k)
end proc:
map(f, [$1..60]); # Robert Israel, May 07 2021
MATHEMATICA
a[n_] := Binomial[n, Times @@ FactorInteger[n][[All, 1]]];
Table[a[n], {n, 1, 60}] (* Jean-François Alcover, May 11 2023 *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Reinhard Zumkeller, Jul 29 2002
STATUS
approved