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

Author

Antti Karttunen, May 17 2018

Keywords

• nonn

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

Crossrefs

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

Cf. also A304113.

Programs

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");