OFFSET
1,4
COMMENTS
a(n) = 1 iff n belongs to A059590.
a(a(n)) <= n for any n > 0.
LINKS
FORMULA
a(n) = A286820(n) / n.
EXAMPLE
The first terms, alongside n*a(n) in decimal and factorial base, are:
n a(n) n*a(n) n*a(n) in factorial base
-- ---- ------ ------------------------
1 1 1 1
2 1 2 1,0
3 1 3 1,1
4 2 8 1,1,0
5 5 25 1,0,0,1
6 1 6 1,0,0
7 1 7 1,0,1
8 1 8 1,1,0
9 1 9 1,1,1
10 3 30 1,1,0,0
11 3 33 1,1,1,1
12 2 24 1,0,0,0
13 2 26 1,0,1,0
14 9 126 1,0,1,0,0
15 2 30 1,1,0,0
16 2 32 1,1,1,0
17 9 153 1,1,1,1,1
18 7 126 1,0,1,0,0
19 8 152 1,1,1,1,0
20 6 120 1,0,0,0,0
PROG
(PARI) isA059590(n) = my (r=2); while (n, if (n%r > 1, return (0), n\=r; r++)); return (1)
a(n) = for (k=1, oo, if (isA059590(k*n), return (k)))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Rémy Sigrist, Jun 25 2017
STATUS
approved