A368713 -id:A368713 - cp's OEIS frontend

A368711 The maximal exponent in the prime factorization of the exponentially odd numbers (A268335).

0, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1

Offset: 1

Amiram Eldar, Jan 04 2024

Differs from A368472 at n = 1, 154, 610, 707, 762, ... .

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

Cf. A033150, A051903, A268335, A368472.

Similar sequences: A368710, A368712, A368713.

f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, OddQ], Max @@ e, Nothing]]; f[1] = 0; Array[f, 150]

lista(kmax) = {my(e); print1(0, ", "); for(k = 2, kmax, e = factor(k)[,2]; if(vecprod(e)%2, print1(vecmax(e), ", ")));}

a(n) = A051903(A268335(n)).
a(n) is odd for n >= 2.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + 2 * Sum_{k>=1} (1 - Product_{p prime} (1 - 1/(p^(2*k-1)*(p^2+p-1)))) = 1.34877064483679975726... .

A368712 The maximal exponent in the prime factorization of the cubefree numbers.

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Jan 04 2024

The asymptotic density of occurrences of 1 is zeta(3)/zeta(2) = 0.730762... (A253905), and the asymptotic density of occurrences of 2 is 1 - zeta(3)/zeta(2) = 0.269237... .

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

Cf. A004709, A005117, A008966, A013929, A033150, A051903, A067259.
Cf. A002117, A013661, A253905.
Similar sequences: A368710, A368711, A368713.

s[n_] := If[n == 1, 0, Max @@ Last /@ FactorInteger[n]]; s /@ Select[Range[100], Max[FactorInteger[#][[;; , 2]]] < 3 &]
(* or *)
f[n_] := Module[{e = Max @@ FactorInteger[n][[;; , 2]]}, If[e < 3, e, Nothing]]; f[1] = 0; Array[f, 100]

lista(kmax) = {my(e); print1(0, ", "); for(k = 2, kmax, e = vecmax(factor(k)[,2]); if(e < 3, print1(e, ", ")));}

from sympy import mobius, integer_nthroot, factorint
def A368712(n):
    def f(x): return n+x-sum(mobius(k)*(x//k**3) for k in range(1, integer_nthroot(x,3)[0]+1))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return max(factorint(m).values(),default=0) # Chai Wah Wu, Aug 12 2024

a(n) = A051903(A004709(n)).
a(n) = 2 - A008966(A004709(n)) for n >= 2.
Except for n = 1, a(n) = 1 or 2.
a(n) = 1 if and only if A004709(n) is squarefree (A005117).
a(n) = 2 if and only if A004709(n) > 1 and is nonsquarefree (A013929), i.e., A004709(n) is in A067259.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2 - zeta(3)/zeta(2) = 2 - A253905 = 1.269237030598... .

A368710 The maximal exponent in the prime factorization of the powerful numbers.

Original entry on oeis.org

0, 2, 3, 2, 4, 2, 3, 5, 2, 2, 6, 3, 4, 2, 3, 2, 3, 7, 4, 2, 2, 3, 3, 2, 5, 8, 5, 2, 4, 3, 2, 3, 4, 4, 2, 2, 3, 9, 2, 6, 4, 4, 3, 2, 6, 4, 5, 2, 5, 2, 2, 3, 5, 3, 10, 2, 3, 7, 2, 2, 4, 3, 3, 3, 2, 3, 2, 2, 5, 6, 2, 6, 2, 3, 2, 4, 5, 4, 4, 11, 2, 7, 3, 2, 8, 3, 4

Offset: 1

Amiram Eldar, Jan 04 2024

Amiram Eldar, Table of n, a(n) for n = 1..10000
Rafael Jakimczuk, Generalizations of the Niven constant and the Feller-Tornier constant, International Mathematical Forum, Vol. 13, No. 9 (2018), pp. 415-425.

Cf. A001694, A033150, A048250, A051903, A059956, A090699.

Similar sequences: A368711, A368712, A368713.

s[n_] := If[n == 1, 0, Max @@ Last /@ FactorInteger[n]]; s /@ Select[Range[3000], # == 1 || Min[FactorInteger[#][[;;, 2]]] > 1 &]
(* or *)
f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[n == 1, 0, If[Min[e] > 1, Max[e], Nothing]]]; Array[f, 3000]

lista(kmax) = {my(e); for(k = 1, kmax, e = factor(k)[,2]; if(k == 1, print1(0, ", "), if(vecmin(e) > 1, print1(vecmax(e), ", "))));}

a(n) = A051903(A001694(n)).
a(n) >= 2 for n >= 2.
Sum_{a(n)<=x} = D_{2,1} * sqrt(x) + O(sqrt(x)), where D_{2,1} = (6/Pi^2) * (2 + Sum_{k>=1} (A051903(k)+2)/(sqrt(k) * A048250(k))) (Jakimczuk, 2018; Theorem 2.1 and Remark 2.3).
Asymptotic mean (consequence of the formula above): Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = D_{2,1} * zeta(3)/zeta(3/2) = D_{2,1} / A090699.
The sum in the formula for D_{2,1} converges slowly: for k up to 10^8, 10^9 and 10^10 the sums are 14.845..., 14.908... and 14.938..., respectively. Thus, a lower bound for the value of this mean, calculated by summing over k=1..10^10, is 4.738... .

A372466 The maximal exponent in the prime factorization of the numbers whose number of divisors is a power of 2 (A036537).

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, May 01 2024

All the terms are of the form 2^k-1 (A000225).

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

Cf. A000225, A036537, A051903, A327839, A372467.

Similar sequences: A368710, A368711, A368712, A368713, A369933, A369935.

pow2Q[n_] := n == 2^IntegerExponent[n, 2]; f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, pow2Q[# + 1] &], Max @@ e, Nothing]]; f[1] = 0; Array[f, 150]

ispow2(n) = n >> valuation(n, 2) == 1;
lista(kmax) = {my(e); print1(0, ", "); for(k = 2, kmax, e = factor(k)[, 2]; if(ispow2(vecprod(apply(x -> x + 1, e))), print1(vecmax(e), ", "))); }

a(n) = A051903(A036537(n)).
a(n) = 2^A372467(n) - 1.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = (d(1) + Sum_{k>=2} (2^k-1) * (d(k) - d(k-1))) / A327839 = 1.25306367526166810834..., where d(k) = Product_{p prime} (1 - 1/p + Sum_{i=1..k} (1/p^(2^i-1)-1/p^(2^i))).

A373058 The sum of the aliquot coreful divisors of the nonsquarefree numbers.

Original entry on oeis.org

2, 6, 3, 6, 14, 6, 10, 18, 5, 12, 14, 30, 36, 30, 22, 15, 42, 7, 10, 26, 24, 42, 30, 21, 62, 34, 96, 15, 38, 70, 39, 42, 66, 30, 46, 90, 14, 33, 80, 78, 126, 98, 58, 39, 90, 11, 62, 30, 42, 126, 66, 60, 102, 70, 216, 21, 74, 30, 114, 51, 78, 150, 78, 82, 126, 13

Offset: 1

Amiram Eldar, May 21 2024

A coreful divisor d of n is a divisor that is divisible by every prime that divides n (see also A307958).
The positive terms of A336563: if k is a squarefree number (A005117) then the only coreful divisor of k is k itself, so k has no aliquot coreful divisors.
The number of the aliquot coreful divisors of the n-th nonsquarefree number is A368039(n).

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

Cf. A005117, A013661, A013929, A065487, A229099, A307958, A336563.

Cf. A084936, A174961, A275699, A368038, A368039, A368040, A368541, A368713.

f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; s[1] = 0; s[n_] := Times @@ f @@@ FactorInteger[n] - n; Select[Array[s, 300], # > 0 &]

lista(kmax) = {my(f); for(k = 1, kmax, f = factor(k); if(!issquarefree(f), print1(prod(i = 1, #f~, (f[i, 1]^(f[i, 2]+1) - 1)/(f[i, 1] - 1) - 1) - k, ", "))); }

from math import prod, isqrt
from sympy import mobius, factorint
def A373058(n):
    def f(x): return n+sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return prod((p**(e+1)-1)//(p-1)-1 for p, e in factorint(m).items())-m # Chai Wah Wu, Jul 22 2024

a(n) = A336563(A013929(n)).

Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = (A065487 - 1)/(1-1/zeta(2))^2 = 1.50461493205911656114... .

cp's OEIS Frontend

A368711 The maximal exponent in the prime factorization of the exponentially odd numbers (A268335).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A368712 The maximal exponent in the prime factorization of the cubefree numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A368710 The maximal exponent in the prime factorization of the powerful numbers.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A372466 The maximal exponent in the prime factorization of the numbers whose number of divisors is a power of 2 (A036537).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A373058 The sum of the aliquot coreful divisors of the nonsquarefree numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula