A375340 -id:A375340 - cp's OEIS frontend

A375341 The maximum exponent in the prime factorization of the numbers that have exactly one non-unitary prime factor.

2, 3, 2, 2, 4, 2, 2, 3, 2, 3, 2, 5, 3, 2, 2, 4, 2, 2, 2, 3, 3, 2, 2, 6, 2, 2, 2, 4, 4, 2, 3, 2, 2, 5, 2, 2, 3, 4, 2, 2, 3, 2, 2, 3, 2, 7, 2, 3, 3, 2, 2, 2, 2, 3, 2, 2, 5, 4, 2, 3, 2, 2, 2, 2, 4, 3, 2, 3, 6, 2, 2, 2, 4, 2, 2, 5, 2, 3, 2, 2, 4, 2, 5, 2, 2, 3, 3, 8, 2, 2, 3, 2, 3, 4, 2, 2, 2, 3, 2, 2, 2, 2, 3, 3

Offset: 1

Amiram Eldar, Aug 12 2024

The positive terms in A375339.

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

Cf. A005361, A048109, A051903, A060687, A082293, A154945, A190641, A375339, A375340.

s[n_] := Module[{e = Select[FactorInteger[n][[;; , 2]], # > 1 &]}, If[Length[e] == 1, e[[1]], Nothing]]; Array[s, 300]

lista(kmax) = {my(e); for(k = 1, kmax, e = select(x -> x > 1, factor(k)[,2]); if(#e == 1, print1(e[1], ", ")));}

a(n) = A051903(A190641(n)).
a(n) = A005361(A190641(n)).
a(n) = A375339(A190641(n)).
a(n) = A132349(A057521(A190641(n))).
a(n) = 2 if and only if A190641(n) is in A060687.
a(n) = 3 if and only if A190641(n) is in A048109.
a(n) <= 3 if and only if A190641(n) is in A082293.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} (2*p-1)/((p-1)*(p^2-1)) / Sum_{p prime} 1/(p^2-1) = A375340 / A154945 = 2.74622231282166656595... .
Asymptotic second raw moment:  = Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)^2 = Sum_{p prime} (4*p^2-3*p+1)/((p-1)^3*(p+1)) / Sum_{p prime} 1/(p^2-1) = 9.064902009520365378603... .
Asymptotic second central moment, or variance, is  - ^2 = 1.52316501808078192104... and the asymptotic standard deviation is sqrt( - ^2) = 1.23416571743051667098... .

A375339 If n has exactly one non-unitary prime factor then a(n) is the exponent of the highest power of this prime that divides n, otherwise a(n) = 0.

Original entry on oeis.org

0, 0, 0, 2, 0, 0, 0, 3, 2, 0, 0, 2, 0, 0, 0, 4, 0, 2, 0, 2, 0, 0, 0, 3, 2, 0, 3, 2, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 2, 2, 0, 0, 4, 2, 2, 0, 2, 0, 3, 0, 3, 0, 0, 0, 2, 0, 0, 2, 6, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 2, 2, 0, 0, 0, 4, 4, 0, 0, 2, 0, 0, 0, 3, 0, 2, 0, 2, 0, 0, 0, 5, 0, 2, 2, 0, 0, 0, 0, 3, 0

Offset: 1

Amiram Eldar, Aug 12 2024

First differs from A212172, A275812 and A372603 at n = 36.
If n = m * p^e, such that m is squarefree, p is a prime that does not divide m and e >= 2, then a(n) = e, otherwise a(n) = 0.
By definition all the positive terms are larger than 1.
The asymptotic density of 0's in this sequence is 1 - Sum_{p prime} (1/(p^2-1)) / zeta(2) = 1 - A059956 * A154945 = 0.66461069244308962639... .
The asymptotic density of the occurrences of k >= 2 in this sequence is Sum_{p prime} (1/(p^(k-1)*(p+1))) / zeta(2). E.g., 0.200755... (A271971) for k = 2, 0.0741777... for k = 3, and 0.0320652... for k = 4.

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

Cf. A005361, A048109, A051903, A060687, A082293, A190641, A359466, A375341 (the positive terms).
Cf. A013661, A059956, A154945, A271971, A375340.
Cf. A212172, A275812, A372603.

a[n_] := Module[{e = Select[FactorInteger[n][[;; , 2]], # > 1 &]}, If[Length[e] == 1, e[[1]], 0]]; Array[a, 100]

a(n) = {my(e = select(x -> x > 1, factor(n)[,2])); if(#e == 1, e[1], 0);}

a(n) = A051903(n) * A359466(n).
a(n) = A005361(n) * A359466(n).
a(A190641(n)) >= 2.
a(n) = 2 if and only if n is in A060687.
a(n) = 3 if and only if n is in A048109.
a(n) <= 3 if and only if n is in A082293.
Asymptotic mean:  = Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} (2*p-1)/((p-1)^2*(p+1)) / zeta(2) = A375340 / A013661 = 0.92105359989459565838... .
Asymptotic second raw moment:  = Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)^2 = Sum_{p prime} (4*p^2-3*p+1)/((p-1)^3*(p+1)) / zeta(2) = 3.04027120804428071157... .
The asymptotic second central moment, or variance, is  - ^2 = 2.19193147416548680815... and the asymptotic standard deviation is sqrt( - ^2) = 1.48051729951577627898... .

cp's OEIS Frontend

A375341 The maximum exponent in the prime factorization of the numbers that have exactly one non-unitary prime factor.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula