A304728 Restricted growth sequence transform of A246277(A303751(n)).

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

Offset

1,2

Comments

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

Links

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

Index entries for sequences computed from indices in prime factorization

Prog

(PARI)
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};
A246277(n) = { if(1==n, 0, while((n%2), n = A064989(n)); (n/2)); };
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; };
v304728 = rgs_transform(vector(76503, n, A246277(A303751(n)))); \\ Needs also code from A303751
A304728(n) = v304728[n];

Crossrefs

Cf. A246277, A303751, A304729, A304731.

Sequence in context: A297161 A227842 A304731 * A284318 A322791 A361255

Adjacent sequences: A304725 A304726 A304727 * A304729 A304730 A304731

Keyword

nonn

Author

Antti Karttunen, May 19 2018

Status

approved