Displaying 1-5 of 5 results found.
page
1
a(n) is the total number of terms (1-digits) in Zeckendorf representation of all divisors of n.
+10
11
1, 2, 2, 4, 2, 5, 3, 5, 4, 5, 3, 10, 2, 6, 5, 7, 4, 9, 4, 10, 5, 6, 3, 13, 5, 5, 7, 11, 3, 13, 4, 10, 8, 6, 6, 16, 3, 8, 5, 14, 4, 12, 4, 11, 10, 8, 3, 18, 6, 11, 9, 10, 5, 16, 5, 14, 7, 6, 4, 23, 4, 8, 9, 13, 6, 16, 5, 10, 7, 14, 4, 23, 4, 8, 12, 12, 8, 13, 4, 20, 10, 9, 5, 23, 9, 9, 8, 17, 2, 22, 6, 12, 8, 6, 8, 24, 3, 12, 13, 19, 5, 15, 4, 14, 13
EXAMPLE
For n=12, its divisors are 1, 2, 3, 4, 6 and 12. Zeckendorf-representations ( A014417) of these numbers are 1, 10, 100, 101, 1001 and 10101. Total number of 1's present is 10 (ten), thus a(12) = 10.
PROG
(PARI)
A072649(n) = { my(m); if(n<1, 0, m=0; until(fibonacci(m)>n, m++); m-2); }; \\ From A072649
A007895(n) = { my(s=0); while(n>0, s++; n -= fibonacci(1+ A072649(n))); (s); };
1, 2, 2, 6, 2, 30, 2, 60, 10, 42, 2, 4200, 2, 126, 70, 660, 2, 9240, 2, 13860, 210, 330, 2, 5082000, 14, 78, 220, 32760, 2, 3783780, 2, 42900, 550, 780, 294, 924924000, 2, 1092, 130, 41621580, 2, 3898440, 2, 112200, 60060, 306, 2, 28078050000, 42, 235620, 1300, 92820, 2, 200119920, 770, 128648520, 1820, 1122, 2, 424964656116000, 2, 3366
PROG
(PARI)
A072649(n) = { my(m); if(n<1, 0, m=0; until(fibonacci(m)>n, m++); m-2); }; \\ From A072649
A003714(n) = { my(s=0, w); while(n>2, w = A072649(n); s += 2^(w-1); n -= fibonacci(w+1)); (s+n); }
A019565(n) = {my(j, v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
Restricted growth sequence transform of A304102, a filter sequence related to the proper divisors of n expressed in Fibonacci number system.
+10
5
1, 2, 2, 3, 2, 4, 2, 5, 3, 4, 2, 6, 2, 5, 4, 7, 2, 8, 2, 6, 5, 5, 2, 9, 3, 4, 10, 11, 2, 12, 2, 6, 5, 13, 5, 14, 2, 13, 4, 9, 2, 15, 2, 11, 8, 10, 2, 16, 17, 18, 13, 6, 2, 19, 5, 20, 13, 5, 2, 21, 2, 13, 6, 22, 4, 23, 2, 24, 10, 25, 2, 26, 2, 10, 18, 27, 28, 12, 2, 29, 30, 13, 2, 31, 13, 32, 5, 33, 2, 34, 5, 35, 13, 5, 13, 21, 2, 36, 37, 38, 2, 39, 2, 9, 15
PROG
(PARI)
up_to = 65537;
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(start_offset, vec, bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A304102(n) = { my(m=1); fordiv(n, d, if(d<n, m *= prime( A304101(d)-1))); (m); };
write_to_bfile(1, rgs_transform(vector(up_to, n, A304102(n))), "b304103.txt");
CROSSREFS
Cf. A305793 (analogous filter for base 2).
1, 2, 2, 3, 2, 4, 2, 4, 5, 6, 2, 7, 2, 8, 9, 8, 2, 10, 2, 11, 12, 13, 2, 14, 15, 16, 12, 14, 2, 17, 2, 18, 19, 20, 21, 22, 2, 23, 24, 25, 2, 26, 2, 27, 28, 29, 2, 30, 9, 31, 32, 33, 2, 34, 24, 35, 36, 37, 2, 38, 2, 39, 40, 39, 41, 42, 2, 43, 44, 45, 2, 46, 2, 47, 48, 49, 50, 51, 2, 52, 53, 54, 2, 55, 56, 57, 58, 59, 2, 60, 61, 62, 63, 64, 65, 66, 2
PROG
(PARI)
up_to = 65537;
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; };
A019565(n) = {my(j, v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A318834(n) = { my(m=1); fordiv(n, d, if(d < n, m *= A019565(eulerphi(d)))); m; };
v318835 = rgs_transform(vector(up_to, n, A318834(n)));
Filter sequence combining primorial base representations of the proper divisors of n; Restricted growth sequence transform of A319708.
+10
2
1, 2, 2, 3, 2, 4, 2, 5, 6, 7, 2, 8, 2, 9, 10, 11, 2, 12, 2, 13, 14, 15, 2, 16, 4, 17, 18, 19, 2, 20, 2, 21, 22, 23, 18, 24, 2, 25, 26, 27, 2, 28, 2, 29, 30, 31, 2, 32, 33, 34, 12, 35, 2, 36, 37, 38, 39, 40, 2, 41, 2, 42, 43, 44, 45, 46, 2, 47, 48, 49, 2, 50, 2, 51, 52, 53, 45, 54, 2, 55, 56, 57, 2, 58, 59, 60, 61, 62, 2, 63, 64, 65, 66, 67, 68, 69, 2, 70, 71
PROG
(PARI)
up_to = 65537;
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; };
A276086(n) = { my(i=0, m=1, pr=1, nextpr); while((n>0), i=i+1; nextpr = prime(i)*pr; if((n%nextpr), m*=(prime(i)^((n%nextpr)/pr)); n-=(n%nextpr)); pr=nextpr); m; };
A319708(n) = { my(m=1); fordiv(n, d, if(d<n, m *= A276086(d))); (m); };
v319709 = rgs_transform(vector(up_to, n, A319708(n)));
Search completed in 0.005 seconds
|