A331109 -id:A331109 - cp's OEIS frontend

A115063 Natural numbers of the form p^F(n_p)q^F(n_q)r^F(n_r)...z^F(n_z), where p,q,r,... are distinct primes and F(n) is a Fibonacci number.

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67

Offset: 1

Giovanni Teofilatto, Mar 01 2006

The complementary sequence is 16, 48, 64, 80, 81, 112, 128, 144, 162, 176, 192, 208, 240, 272, 304, 320, 324, 336, 368, 384, 400, ... - R. J. Mathar, Apr 22 2010
Or exponentially Fibonacci numbers. - Vladimir Shevelev, Nov 15 2015
Sequences A004709, A005117, A046100 are subsequences. - Vladimir Shevelev, Nov 16 2015
Let h_k be the density of the subsequence of A115063 of numbers whose prime power factorization has the form Product_i p_i^e_i where the e_i all squares <= k^2. Then for every k>1 there exists eps_k>0 such that for any x from the interval (h_k-eps_k, h_k) there is no sequence S of positive integers such that x is the density of numbers whose prime power factorization has the form Product_i p_i^e_i where the e_i are all in S. For a proof, see [Shevelev], the second link. - Vladimir Shevelev, Nov 17 2015
Numbers whose sets of unitary divisors (A077610) and Zeckendorf-infinitary divisors (see A318465) coincide. Also, numbers whose sets of unitary divisors and dual-Zeckendorf-infinitary divisors (see A331109) coincide. - Amiram Eldar, Aug 09 2024

			12 is a term, since 12=2^2*3^1 and the exponents 2 and 1 are terms of Fibonacci sequence (A000045). - _Vladimir Shevelev_, Nov 15 2015

Amiram Eldar, Table of n, a(n) for n = 1..10000
Vladimir Shevelev, Exponentially S-numbers, arXiv:1510.05914 [math.NT], 2015.
Vladimir Shevelev, Set of all densities of exponentially S-numbers, arXiv:1511.03860 [math.NT], 2015.
Vladimir Shevelev, S-exponential numbers, Acta Arithmetica, Vol. 175 (2016), 385-395.

Cf. A004709, A005117, A046100, A077610, A115975, A197680, A209061, A318465, A331109.

fibQ[n_] := IntegerQ @ Sqrt[5 n^2 - 4] || IntegerQ @ Sqrt[5 n^2 + 4]; aQ[n_] := AllTrue[FactorInteger[n][[;;, 2]], fibQ]; Select[Range[100], aQ] (* Amiram Eldar, Oct 06 2019 *)

Sum_{i<=x, i is in A115063} 1 = h*x + O(sqrt(x)*log x*e^(c*sqrt(log x)/(log(log x))), where c = 4*sqrt(2.4/log 2) = 7.44308... and h = Product_{prime p}(1 + Sum_{i>=2} (u(i)-u(i-1))/p^i) = 0.944335905... where u(n) is the characteristic function of sequence A000045. The calculations of h over the formula were done independently by Juan Arias-de-Reyna and Peter J. C. Moses.

For a proof of the formula, see [Shevelev], the first link. - Vladimir Shevelev, Nov 17 2015

a(35) inserted by Amiram Eldar, Oct 06 2019

A375430 The maximum exponent in the unique factorization of n in terms of distinct terms of A115975 using the dual Zeckendorf representation of the exponents in the prime factorization of n; a(1) = 0.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 3, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 3, 3, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 3, 1, 2, 2, 2, 1, 1, 1, 2, 1

Offset: 1

Amiram Eldar, Aug 15 2024

First differs from A299090 at n = 128. Differs from A046951 and A159631 at n = 1, 36, 64, 72, ... .

When the exponents in the prime factorization of n are expanded as sums of distinct Fibonacci numbers using the dual Zeckendorf representation (A104326), we get a unique factorization of n in terms of distinct terms of A115975, i.e., n is represented as a product of prime powers (A246655) whose exponents are Fibonacci numbers. a(n) is the maximum exponent of these prime powers. Thus all the terms are Fibonacci numbers.

			For n = 8 = 2^3, the dual Zeckendorf representation of 3 is 11, i.e., 3 = Fibonacci(2) + Fibonacci(3) = 1 + 2. Therefore 8 = 2^(1+2) = 2^1 * 2^2, and a(8) = 2.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A000045, A095791, A104326, A115975, A246655, A331109, A368781, A375428, A375431.

Cf. A046951, A159631, A299090.

A130312[n_] := Module[{k = 0}, While[Fibonacci[k] <= n, k++]; Fibonacci[k-2]]; a[n_] := A130312[1 + Max[FactorInteger[n][[;;, 2]]]]; a[1] = 0; Array[a, 100]

A130312(n) = {my(k = 0); while(fibonacci(k) <= n, k++); fibonacci(k-2);}
a(n) = if(n == 1, 0, A130312(1 + vecmax(factor(n)[,2])));

a(n) = A130312(1 + A051903(n)).
a(n) = A000045(A375431(n)).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=4} Fibonacci(k) * (1 - 1/zeta(Fibonacci(k)-1)) = 1.48543763231328442311... .

A331110 The sum of dual-Zeckendorf-infinitary divisors of n = Product_{i} p(i)^r(i): divisors d = Product_{i} p(i)^s(i), such that the dual Zeckendorf expansion (A104326) of each s(i) contains only terms that are in the dual Zeckendorf expansion of r(i).

Original entry on oeis.org

1, 3, 4, 5, 6, 12, 8, 15, 10, 18, 12, 20, 14, 24, 24, 27, 18, 30, 20, 30, 32, 36, 24, 60, 26, 42, 40, 40, 30, 72, 32, 45, 48, 54, 48, 50, 38, 60, 56, 90, 42, 96, 44, 60, 60, 72, 48, 108, 50, 78, 72, 70, 54, 120, 72, 120, 80, 90, 60, 120, 62, 96, 80, 135, 84, 144

Offset: 1

Amiram Eldar, Jan 09 2020

First differs from A188999 at n = 32.

			a(32) = 45 since 32 = 2^5 and the dual Zeckendorf expansion of 5 is 110, i.e., its dual Zeckendorf representation is a set with 2 terms: {2, 3}. There are 4 possible exponents of 2: 0, 2, 3 and 5, corresponding to the subsets {}, {2}, {3} and {2, 3}. Thus 32 has 4 dual-Zeckendorf-infinitary divisors: 2^0 = 1, 2^2 = 4, 2^3 = 8, and 2^5 = 32, and their sum is 1 + 4 + 8 + 32 = 45.

Amiram Eldar, Table of n, a(n) for n = 1..10000

The number of dual-Zeckendorf-infinitary divisors of n is in A331109.

Cf. A049417, A049418, A074847, A097863, A104326, A188999, A331107.

fibTerms[n_] := Module[{k = Ceiling[Log[GoldenRatio, n*Sqrt[5]]], t = n, fr = {}}, While[k > 1, If[t >= Fibonacci[k], AppendTo[fr, 1]; t = t - Fibonacci[k], AppendTo[fr, 0]]; k--]; fr];
dualZeck[n_] := Module[{v = fibTerms[n]}, nv = Length[v]; i = 1; While[i <= nv - 2, If[v[[i]] == 1 && v[[i + 1]] == 0 && v[[i + 2]] == 0, v[[i]] = 0; v[[i + 1]] = 1; v[[i + 2]] = 1; If[i > 2, i -= 3]]; i++]; i = Position[v, _?(# > 0 &)]; If[i == {}, {}, v[[i[[1, 1]] ;; -1]]]];
f[p_, e_] := p^Fibonacci[1 + Position[Reverse@dualZeck[e], _?(# == 1 &)]];
a[1] = 1; a[n_] := Times @@ (Flatten@(f @@@ FactorInteger[n]) + 1); Array[a, 100]

Multiplicative with a(p^e) = Product_{i} (p^s(i) + 1), where s(i) are the terms in the dual Zeckendorf representation of e (A104326).

A375272 The number of factors of n of the form p^Fibonacci(k), where p is a prime and k >= 2, when the factorization is uniquely done using the dual Zeckendorf representation of the exponents in the prime factorization of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 2, 2, 1, 2, 2, 3, 1, 3, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 2, 3, 2, 2, 1, 3, 1, 2, 2, 3, 2, 3, 1, 2, 2, 3, 1, 3, 1, 2, 2, 2, 2, 3, 1, 3, 2, 2, 1, 3, 2, 2, 2

Offset: 1

Amiram Eldar, Aug 09 2024

First differs from A086435 at n = 36. Differs from A266226 at n = 1, 36, ... .
The number of dual-Zeckendorf-infinitary divisors of n (defined in A331109) that are prime powers (A246655).
a(n) depends only on the prime signature of n.
Analogous to A064547 (binary representation) and A318464 (Zeckendorf representation).

			For n = 8 = 2^3, the dual Zeckendorf representation of 3 is 11, i.e., 3 = Fibonacci(2) + Fibonacci(3) = 1 + 2. Therefore 8 = 2^(1+2) = 2^1 * 2^2, and a(8) = 2.
For n = 256 = 2^8, the dual Zeckendorf representation of 8 is 1011, i.e., 8 = Fibonacci(2) + Fibonacci(3) + Fibonacci(5) = 1 + 2 + 5. Therefore 256 = 2^(1+2+5) = 2^1 * 2^2 * 2^5, and a(256) = 3.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A000045, A064547, A077761, A086435, A112309, A112310, A246655, A266226, A318464, A331109.

toDualZeck[n_] := Module[{s = 0, v = 0, i = 0, f}, While[s < n, s += Fibonacci[i + 2]; v += 2^i; i++]; i--; While[i >= 0, f = Fibonacci[i + 2]; If[s - f >= n, s -= f; v -= 2^i]; i--]; v]; (* A003754, after Rémy Sigrist's PARI code in A112309 *)
f[p_, e_] := DigitCount[toDualZeck[e], 2, 1]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]

todualzeck(n) = {my (s=0, v=0); for (i=0, oo, if (s>=n, forstep (j=i-1, 0, -1, if (s-fibonacci(2+j)>=n, s-=fibonacci(2+j); v-=2^j;);); return (v);); s+=fibonacci(2+i); v+=2^i;);} \\ A003754, Rémy Sigrist's code in A112309
a(n) = vecsum(apply(x -> hammingweight(todualzeck(x)), factor(n)[, 2]));

Additive with a(p^e) = A112310(e).
a(n) = log_2(A331109(n)).
Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761), C = Sum_{k>=2} (A112310(k)-A112310(k-1)) * P(k) = 0.18790467121403662496..., and P(s) is the prime zeta function.

A376887 The number of divisors of n that are products of factors of the form p^(k!) with multiplicities not larger than their multiplicity in n, where p is a prime and k >= 1, when the factorization of n is uniquely done using the factorial-base representation of the exponents in the prime factorization of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 4, 3, 2, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 2, 6, 4, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 2, 6, 2, 4, 4, 4, 2, 8, 4, 8, 4, 4, 2, 8, 2, 4, 4, 2, 4, 8, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 4, 8, 2, 6, 3, 4, 2, 8, 4, 4, 4

Offset: 1

Amiram Eldar, Oct 08 2024

See A376885 for details about this factorization.
If n = Product p_i^e_i is the canonical prime factorization of n, then the divisors that are counted by this function are d = Product p_i^s_i, where all the digits of s_i in factorial base are not larger than the corresponding digits of e_i.
The sum of these divisors is given by A376888(n).

			For n = 12 = 2^2 * 3^1, the representation of 2 in factorial base is 10, i.e., 2 = 2!, so 12 = (2^(2!))^1 * (3^(1!))^1 and a(12) = (1+1) * (1+1) = 4, corresponding to the 4 divisors 1, 3, 4 and 12.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A007623, A034968, A376885, A376886, A376888.

Similar sequences: A037445, A038148, A074848, A318465, A331109.

fdigprod[n_] := Module[{k = n, m = 2, r, s = 1}, While[{k, r} = QuotientRemainder[k, m]; k != 0 || r != 0, s *= (r+1); m++]; s]; f[p_, e_] := fdigprod[e]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

fdigprod(n) = {my(k = n, m = 2, r, s = 1); while([k, r] = divrem(k, m); k != 0 || r != 0, s *= (r+1); m++); s;}
a(n) = {my(e = factor(n)[, 2]); prod(i = 1, #e, fdigprod(e[i]));}

Multiplicative with a(p^e) = A227154(e).

A372381 The number of divisors of the largest divisor of n whose number of divisors is a power of 2.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 4, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 2, 8, 2, 4, 4, 4, 2, 8, 4, 8, 4, 4, 2, 8, 2, 4, 4, 4, 4, 8, 2, 4, 4, 8, 2, 8, 2, 4, 4, 4, 4, 8, 2, 8, 4, 4, 2, 8, 4, 4, 4

Offset: 1

Amiram Eldar, Apr 29 2024

First differs from A286324 at n = 32, and from A331109 at n = 64.

Also, the number of infinitary divisors of the largest divisor of n whose number of divisors is a power of 2.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A000005, A004709, A036537, A037445, A372379, A372380.

Cf. A286324, A331109.

f[p_, e_] := 2^Floor[Log2[e + 1]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

a(n) = vecprod(apply(x -> 2^exponent(x+1), factor(n)[, 2]));

from math import prod
from sympy import factorint
def A372381(n): return prod(1<<(e+1).bit_length()-1 for e in factorint(n).values()) # Chai Wah Wu, Apr 30 2024

Multiplicative with a(p^e) = 2^floor(log_2(e+1)).
a(n) = A000005(A372379(n)).
a(n) = A037445(A372379(n)).
a(n) = A000005(n) if and only if n is in A036537.
a(n) <= A372380(n), with equality if and only if n is cubefree (A004709).

A375269 Partial products of A115975.

Original entry on oeis.org

1, 2, 6, 24, 120, 840, 6720, 60480, 665280, 8648640, 147026880, 2793510720, 64250746560, 1606268664000, 43369253928000, 1257708363912000, 38988959281272000, 1247646697000704000, 46162927789026048000, 1892680039350067968000, 81385241692052922624000, 3825106359526487363328000

Offset: 1

Amiram Eldar, Aug 09 2024

nonn

First differs from A334395 at n = 42.
Numbers with a record number of dual-Zeckendorf-infinitary divisors (A331109). Also, indices of records in A375272.
a(n) is the least number k such that A375272(k) = n-1 and A331109(k) = 2^(n-1).

			A115975 begins with 1, 2, 3, 4, 5, 7, ..., so, a(1) = 1, a(2) = 1 * 2 = 2, a(3) = 1 * 2 * 3 = 6, ..., a(6) = 1 * 2 * 3 * 4 * 5 * 7 = 840.

Amiram Eldar, Table of n, a(n) for n = 1..358

Cf. A037992 (analogous with "Fermi-Dirac primes", A050376), A115975, A331109, A334395, A375271, A375272.

Subsequence of A025487.

fib[lim_] := Module[{s = {}, f = 1, k = 2}, While[f <= lim, AppendTo[s, f]; k++; f = Fibonacci[k]]; s];
seq[max_] := Module[{s = {}, p = 2, e = 1, f = {}}, While[e > 0, e = Floor[Log[p, max]]; If[f == {}, f = fib[e], f = Select[f, # <= e &]]; s = Join[s, p^f]; p = NextPrime[p]]; FoldList[Times, 1, Sort[s]]]; seq[250]

fib(lim) = {my(s = List(), f = 1, k = 2); while(f <= lim, listput(s, f); k++; f = fibonacci(k)); Vec(s);}
lista(pmax) = {my(s = [1], p = 2, e = 1, f = [], r = 1); while(e > 0, e = logint(pmax, p); if(#f == 0, f = fib(e), f = select(x -> x <= e, f)); s = concat(s, apply(x -> p^x, f)); p = nextprime(p+1)); s = vecsort(s); for(i = 1, #s, r *= s[i]; print1(r, ", "));}

a(n) = Product_{k=1..n} A115975(k).

cp's OEIS Frontend

A115063 Natural numbers of the form p^F(n_p)*q^F(n_q)*r^F(n_r)*...*z^F(n_z), where p,q,r,... are distinct primes and F(n) is a Fibonacci number.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A375430 The maximum exponent in the unique factorization of n in terms of distinct terms of A115975 using the dual Zeckendorf representation of the exponents in the prime factorization of n; a(1) = 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A331110 The sum of dual-Zeckendorf-infinitary divisors of n = Product_{i} p(i)^r(i): divisors d = Product_{i} p(i)^s(i), such that the dual Zeckendorf expansion (A104326) of each s(i) contains only terms that are in the dual Zeckendorf expansion of r(i).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A375272 The number of factors of n of the form p^Fibonacci(k), where p is a prime and k >= 2, when the factorization is uniquely done using the dual Zeckendorf representation of the exponents in the prime factorization of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A372381 The number of divisors of the largest divisor of n whose number of divisors is a power of 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A375269 Partial products of A115975.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

A115063 Natural numbers of the form p^F(n_p)q^F(n_q)r^F(n_r)...z^F(n_z), where p,q,r,... are distinct primes and F(n) is a Fibonacci number.