OFFSET
1,3
COMMENTS
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..8179 (terms up to 13!)
EXAMPLE
The first terms, alongside the corresponding factorial base representations, are:
n a(n) fac(a(n))
-- ---- ---------
1 0 (0)
2 1 (1)
3 2 (1,0)
4 4 (2,0)
5 5 (2,1)
6 10 (1,2,0)
7 13 (2,0,1)
8 14 (2,1,0)
9 19 (3,0,1)
10 20 (3,1,0)
11 22 (3,2,0)
12 23 (3,2,1)
13 46 (1,3,2,0)
14 67 (2,3,0,1)
MAPLE
b:= proc(n, i) local r; `if`(n<i, [n],
[b(iquo(n, i, 'r'), i+1)[], r])
end:
t:= n-> (l-> is(nops(l)=nops({l[]})))(b(n, 2)):
select(t, [$0..1000])[]; # Alois P. Heinz, Nov 16 2018
MATHEMATICA
q[n_] := Module[{k = n, m = 2, r, s = {}}, While[{k, r} = QuotientRemainder[k, m]; k != 0|| r != 0, AppendTo[s, r]; m++]; UnsameQ @@ s]; Select[Range[0, 720], q] (* Amiram Eldar, Feb 21 2024 *)
PROG
(PARI) is(n) = my (s=0); for (k=2, oo, if (n==0, return (1)); my (d=n%k); if (bittest(s, d), return (0), s+=2^d; n\=k))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Rémy Sigrist, Nov 16 2018
STATUS
approved