OFFSET
1,2
COMMENTS
Restricted growth sequence transform of A291750, which means that this is the lexicographically least sequence a, such that for all i, j: a(i) = a(j) <=> A291750(i) = A291750(j) <=> A003557(i) = A003557(j) and A048250(i) = A048250(j). That this is equal to the definition given in the title follows because any such lexicographically least sequence satisfying relation <=> is also the least sequence satisfying relation => with the same parameters.
Sigma (A000203) and psi (A001615) are functions of this sequence. See comments in A291750 for the reason. For example, to find the value of A001615(n) when we know just a(n), but without knowing n, let m be the least i for which a(i) = a(n); then A001615(n) = A003991(A291750(m)) = A003557(m) * A048250(m).
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
PROG
(PARI)
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om, invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om, invec[i], i); outvec[i] = u; u++ )); outvec; };
A048250(n) = if(n<1, 0, sumdiv(n, d, if(core(d)==d, d)));
v291751 = rgs_transform(vector(65537, n, A291750(n)));
A291751(n) = v291751[n];
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 06 2017
EXTENSIONS
Name changed and comments added by Antti Karttunen, Nov 24 2018
STATUS
approved