A304098 Restricted growth sequence transform of A046523(A282291(n)).
1, 2, 3, 2, 4, 5, 4, 2, 3, 6, 3, 7, 8, 2, 3, 6, 3, 7, 4, 7, 3, 6, 9, 6, 10, 11, 12, 13, 12, 13, 12, 13, 10, 2, 3, 6, 3, 7, 4, 7, 3, 6, 9, 6, 8, 14, 3, 6, 9, 6, 8, 7, 8, 6, 9, 15, 9, 16, 12, 13, 10, 13, 10, 13, 17, 5, 4, 7, 4, 18, 19, 18, 4, 7, 8, 7, 20, 21, 4, 7, 8, 7, 20, 18, 20, 7, 8, 14, 8, 22, 23, 24, 17, 24, 17, 2, 3, 6, 3, 7, 4, 7, 3, 6, 9
Offset: 1
Links
- Antti Karttunen, Table of n, a(n) for n = 1..60824
Programs
-
PARI
\\ Needs also code from A282291: 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; }; A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ From A046523 write_to_bfile(1,rgs_transform(vector(60824,n,A046523(A282291(n)))),"b304098.txt");
Comments