A304736 Restricted growth sequence transform of A048675(A064413(n)).

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

Offset: 1

Author

Antti Karttunen, May 18 2018

Keywords

• nonn

Links

• Antti Karttunen, Table of n, a(n) for n = 1..65539
• Index entries for sequences related to EKG sequence

Crossrefs

Cf. A048675, A064413, A304735.

Cf. also A064301, A304113.

Programs

PARI

A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
A304735(n) = A048675(A064413(n)); \\ Needs also code for A064413.
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; };
v304736 = rgs_transform(vector(65539,n,A304735(n)));
A304736(n) = v304736[n];