A304729 Restricted growth sequence transform of A046523(A303751(n)).

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

Offset: 1

Author

Antti Karttunen, May 19 2018

Keywords

• nonn

Links

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

Crossrefs

Cf. A046523, A303751.

Cf. also A304728, A304098, A304732.

Programs

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; };
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
v304729 = rgs_transform(vector(76503,n,A046523(A303751(n)))); \\ Needs also code from A303751
A304729(n) = v304729[n];