OFFSET
1,3
COMMENTS
a(p_i) < a(p_j) for i<j.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1000
MATHEMATICA
f[n_] := Length[ Select[ Divisors[n! ], GCD[ #, n] == 1 &]]; Table[ f[n], {n, 33}]
PROG
(PARI) a(n)={if(n<=2, 1, numdiv(Mat(select(f->n%f[1], Col(factor(n!))))))} \\ Andrew Howroyd, Nov 05 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Jul 19 2004, based on a suggestion from Leroy Quet, Jun 18 2004
EXTENSIONS
Terms a(34) and beyond from Andrew Howroyd, Nov 05 2019
STATUS
approved