A304115 Restricted growth sequence transform of A046523(A098550(n)), a filter sequence based on the prime signature of the Yellowstone permutation.

1, 2, 2, 3, 3, 4, 5, 5, 2, 5, 3, 6, 5, 7, 2, 5, 5, 6, 4, 5, 5, 2, 2, 5, 5, 6, 6, 5, 8, 2, 6, 5, 9, 5, 5, 5, 10, 5, 11, 3, 5, 6, 2, 11, 5, 6, 5, 9, 5, 6, 2, 12, 5, 6, 6, 5, 7, 9, 5, 5, 2, 2, 5, 5, 13, 6, 6, 5, 9, 5, 14, 5, 11, 5, 15, 6, 5, 11, 2, 11, 5, 6, 6, 12, 5, 5, 2, 2, 5, 5, 9, 9, 6, 4, 11, 5, 14, 5, 5, 11, 2, 14, 5, 16, 5

Offset

1,2

Comments

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

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

Links

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

Prog

(PARI)
\\ Needs also code for A098550:
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
Aux304115(n) = A046523(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, Aux304115(n))), "b304115.txt");

Crossrefs

Cf. A046523, A098550, A251138, A251140, A304114.

Cf. also A304113.

Sequence in context: A349122 A076918 A162620 * A338365 A194249 A285759

Adjacent sequences: A304112 A304113 A304114 * A304116 A304117 A304118

Keyword

nonn

Author

Antti Karttunen, May 17 2018

Status

approved