A318839 Restricted growth sequence transform of A318838.

1, 1, 2, 2, 3, 2, 4, 5, 6, 3, 7, 8, 9, 4, 10, 11, 12, 6, 13, 10, 14, 7, 15, 16, 17, 9, 18, 14, 19, 10, 20, 21, 22, 12, 23, 24, 25, 13, 26, 27, 28, 14, 29, 22, 30, 15, 31, 32, 33, 17, 34, 26, 35, 18, 36, 37, 38, 19, 39, 40, 41, 20, 42, 43, 44, 22, 45, 34, 46, 23, 47, 48, 49, 25, 50, 38, 51, 26, 52, 53, 54, 28, 55, 56, 57, 29, 58, 59, 60, 30, 61, 46, 62, 31, 63

Offset

1,3

Comments

For all i, j: a(i) = a(j) => A000010(i) = A000010(j).

Links

Antti Karttunen, Table of n, a(n) for n = 1..65537

Prog

(PARI)
up_to = 65537;
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; };
A007431(n) = sumdiv(n, d, moebius(n/d)*eulerphi(d));
A318838(n) = { my(m=1); fordiv(n, d, if((A007431(d)!=0), m *= prime(A007431(d)))); (m); };
v318839 = rgs_transform(vector(up_to, n, A318838(n)));
A318839(n) = v318839[n];

Crossrefs

Cf. A000010, A007431, A318838.

Cf. also A318837.

Sequence in context: A324343 A328474 A369065 * A347664 A304105 A305978

Adjacent sequences: A318836 A318837 A318838 * A318840 A318841 A318842

Keyword

nonn

Author

Antti Karttunen, Sep 05 2018

Status

approved