A305301 Restricted growth sequence transform of A278222(A304760(n)), constructed from runlengths of 1-digits in base-3 representation of A254103(n).
1, 2, 1, 2, 2, 3, 1, 1, 3, 2, 2, 2, 2, 4, 2, 5, 4, 3, 1, 2, 5, 1, 1, 3, 3, 3, 3, 2, 5, 2, 2, 5, 6, 2, 2, 5, 2, 6, 3, 1, 7, 7, 2, 5, 2, 7, 2, 2, 4, 7, 2, 2, 7, 3, 1, 2, 7, 3, 1, 3, 5, 2, 2, 2, 8, 3, 1, 5, 5, 2, 2, 5, 3, 4, 4, 2, 7, 4, 2, 7, 9, 5, 5, 2, 5, 2, 2, 7, 3, 5, 5, 5, 5, 1, 1, 5, 6, 5, 5, 5, 5, 1, 1, 5, 10, 2
Offset: 0
Links
- Antti Karttunen, Table of n, a(n) for n = 0..65537
Programs
-
PARI
A254103(n) = if(!n,n,if(!(n%2),(3*A254103(n/2))-1,(3*(1+A254103((n-1)/2)))\2)); A289813(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); } \\ From A289813 A304760(n) = A289813(A254103(n)); A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t }; A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ From A046523 A278222(n) = A046523(A005940(1+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; }; v305301 = rgs_transform(vector(65538,n,A278222(A304760(n-1)))); A305301(n) = v305301[1+n];
Comments