A304744 Restricted growth sequence transform of A046523(A052330(n)).

1, 2, 2, 3, 4, 5, 6, 7, 2, 3, 3, 8, 6, 7, 9, 10, 2, 3, 3, 8, 6, 7, 9, 10, 3, 8, 8, 11, 9, 10, 12, 13, 4, 6, 5, 7, 14, 15, 15, 16, 6, 9, 7, 10, 17, 18, 18, 19, 6, 9, 7, 10, 17, 18, 18, 19, 9, 12, 10, 13, 20, 21, 21, 22, 2, 3, 3, 8, 6, 7, 9, 10, 3, 8, 8, 11, 9, 10, 12, 13, 3, 8, 8, 11, 9, 10, 12, 13, 8, 11, 11, 23, 12, 13, 24, 25, 6, 9, 7, 10, 17, 18, 18, 19

Offset

0,2

Comments

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

Links

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

Index entries for sequences related to binary expansion of n

Prog

(PARI)
up_to_e = 17; \\ Good for computing up to n = (2^up_to_e)-1
v050376 = vector(up_to_e);
ispow2(n) = (n && !bitand(n, n-1));
i = 0; for(n=1, oo, if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to_e, break));
A050376(n) = v050376[n];
A052330(n) = { my(p=1, i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); };
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
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; };
v304744 = rgs_transform(vector(65538, n, A046523(A052330(n-1))));
A304744(n) = v304744[1+n];

Crossrefs

Cf. A046523, A052330.

Cf. also A286622 (A278222), A302791, A304535, A304745.

Sequence in context: A038668 A251629 A279033 * A327225 A071754 A266113

Adjacent sequences: A304741 A304742 A304743 * A304745 A304746 A304747

Keyword

nonn

Author

Antti Karttunen, May 27 2018

Status

approved