A304732 Restricted growth sequence transform of A046523(A304531(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, 12, 13, 22, 13, 21, 10, 24, 4, 7, 19, 7, 18, 9, 25, 9, 24, 15, 16, 15, 10, 21, 10, 26, 27, 28, 15, 10, 21, 10, 29, 30, 31, 32, 31

Offset: 1

Author

Antti Karttunen, May 19 2018

Keywords

• nonn

Comments

For all i, j: a(i) = a(j) => A304536(i-1) = A304536(j-1).

Links

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

Crossrefs

Cf. A046523, A304531, A304536.

Differs from related A304729 for the first time at n=66, where a(66) = 12, while A304729(66) = 15.

Compare also to scatter plots of A304535, A304098, A304729.

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
v304732 = rgs_transform(vector(74431,n,A046523(A304531(n)))); \\ Needs also code from A304531
A304732(n) = v304732[n];