A305296 Restricted growth sequence transform of A305295(n), formed from 1-digits in ternary representation of A245612(n).
1, 2, 3, 3, 4, 5, 2, 6, 7, 8, 9, 1, 3, 2, 8, 2, 10, 11, 2, 12, 13, 14, 2, 1, 4, 15, 16, 7, 17, 18, 8, 1, 19, 20, 21, 22, 3, 23, 16, 2, 24, 25, 26, 27, 3, 28, 2, 29, 7, 30, 28, 30, 31, 32, 33, 2, 34, 31, 20, 18, 17, 16, 2, 8, 35, 36, 37, 38, 39, 8, 40, 41, 4, 42, 43, 44, 31, 45, 16, 46, 47, 48, 49, 50, 51, 4, 52, 53, 4, 1, 54, 15, 3
Offset: 0
Keywords
Links
- Antti Karttunen, Table of n, a(n) for n = 0..65537
Crossrefs
Programs
-
PARI
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961 A048673(n) = (A003961(n)+1)/2; A254049(n) = A048673((2*n)-1); A245612(n) = if(n<2,1+n,if(!(n%2),(3*A245612(n/2))-1,A254049(A245612((n-1)/2)))); A289813(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); } \\ From A289813 A305295(n) = A289813(A245612(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; }; v305296 = rgs_transform(vector(65538,n,A305295(n-1))); A305296(n) = v305296[1+n];
Comments