OFFSET
1,4
COMMENTS
For n>1: a(n)=1 iff n is prime.
LINKS
Robin Visser, Table of n, a(n) for n = 1..3000
FORMULA
a(n) = binomial(n, A006530(n)). - Michel Marcus, May 06 2020
MATHEMATICA
Table[Binomial[n, FactorInteger[n][[-1, 1]]], {n, 50}] (* Harvey P. Dale, Feb 25 2015 *)
PROG
(PARI) a(n) = if (n==1, 1, binomial(n, vecmax(factor(n)[, 1]))); \\ Michel Marcus, May 06 2020
(Sage)
def a(n):
if n==1: return 1
return binomial(n, factor(n)[-1][0]) # Robin Visser, Nov 25 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Feb 06 2003
STATUS
approved