A304114 Restricted growth sequence transform of A246277(A098550(n)), related to indices in the prime factorization of Yellowstone permutation.

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

Offset

1,2

Links

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

Formula

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

Prog

(PARI)
\\ Needs also code for A098550:
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)); };
Aux304114(n) = A246277(A098550(n));
write_to_bfile(start_offset, vec, bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
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; };
write_to_bfile(1, rgs_transform(vector(65539, n, Aux304114(n))), "b304114.txt");

Crossrefs

Cf. A098550, A246277, A304115.

Cf. also A304112.

Sequence in context: A359320 A029144 A173244 * A031234 A030585 A030565

Adjacent sequences: A304111 A304112 A304113 * A304115 A304116 A304117

Keyword

nonn

Author

Antti Karttunen, May 17 2018

Status

approved