A330596 -id:A330596 - cp's OEIS frontend

A057723 Sum of positive divisors of n that are divisible by every prime that divides n.

1, 2, 3, 6, 5, 6, 7, 14, 12, 10, 11, 18, 13, 14, 15, 30, 17, 24, 19, 30, 21, 22, 23, 42, 30, 26, 39, 42, 29, 30, 31, 62, 33, 34, 35, 72, 37, 38, 39, 70, 41, 42, 43, 66, 60, 46, 47, 90, 56, 60, 51, 78, 53, 78, 55, 98, 57, 58, 59, 90, 61, 62, 84, 126, 65, 66, 67, 102, 69, 70

Offset: 1

Leroy Quet, Oct 27 2000

			The divisors of 12 that are divisible by both 2 and 3 are 6 and 12. So a(12) = 6 + 12 = 18.

Ivan Neretin, Table of n, a(n) for n = 1..10000
József Sándor and Krassimir Atanassov, Some new arithmetic functions, Notes on Number Theory and Discrete Mathematics, Volume 30, 2024, Number 4, Pages 851-856. See sigma- function.

Row sums of triangle A284318.
Cf. A000203 (sigma), A007947 (rad), A005361 (number of these divisors).
Cf. A049060 and A060640 (other sigma-like functions).
Cf. A065487, A330596.

[&*PrimeDivisors(n)*SumOfDivisors(n div &*PrimeDivisors(n)): n in [1..70]]; // Vincenzo Librandi, May 14 2015

seq(mul(f[1]*(f[1]^f[2]-1)/(f[1]-1), f = ifactors(n)[2]), n = 1 .. 100); # Robert Israel, May 13 2015

Table[(b = Times @@ FactorInteger[n][[All, 1]])*DivisorSigma[1, n/b], {n, 70}] (* Ivan Neretin, May 13 2015 *)
f[p_, e_] := (p^(e+1)-1)/(p-1) - 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 15 2023 *)

a(n) = {my(f = factor(n)); for (i=1, #f~, f[i,2]=1); my(pp = factorback(f)); sumdiv(n, d, if (! (d % pp), d, 0));} \\ Michel Marcus, May 14 2015

If n = Product p_i^e_i then a(n) = Product (p_i + p_i^2 + ... + p_i^e_i).
a(n) = rad(n)*sigma(n/rad(n)) = A007947(n)*A000203(A003557(n)). - Ivan Neretin, May 13 2015
Dirichlet g.f.: zeta(s) * zeta(s-1) * Product(p prime, 1 - p^(-s) + p^(1-2*s)). - Robert Israel, May 13 2015
Sum_{k=1..n} a(k) ~ c * Pi^2 * n^2 / 12, where c = A330596 = Product_{primes p} (1 - 1/p^2 + 1/p^3) = 0.7485352596823635646442150486379106016416403430053244045... - Vaclav Kotesovec, Dec 18 2019
a(n) = Sum_{d|n, rad(d)=rad(n)} d. - R. J. Mathar, Jun 02 2020
Lim_{n->oo} (1/n) * Sum_{k=1..n} a(k)/k = Product_{p prime}(1 + 1/(p*(p^2-1))) = 1.231291... (A065487). - Amiram Eldar, Jun 10 2020
a(n) = Sum_{d|n, gcd(d, n/d) = 1} (-1)^omega(n/d) * sigma(d). - Ilya Gutkovskiy, Apr 15 2021

A337050 Numbers without an exponent 2 in their prime factorization.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 8, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 24, 26, 27, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 46, 47, 48, 51, 53, 54, 55, 56, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 80, 81, 82, 83, 85, 86, 87

Offset: 1

Amiram Eldar, Aug 12 2020

Numbers k such that the powerful part (A057521) of k is a cubefull number (A036966).
Numbers k such that A003557(k) = k/A007947(k) is a powerful number (A001694).
The asymptotic density of this sequence is Product_{primes p} (1 - 1/p^2 + 1/p^3) = 0.748535... (A330596).
A304364 is apparently a subsequence.
These numbers were named semi-2-free integers by Suryanarayana (1971). - Amiram Eldar, Dec 29 2020

			6 = 2^1 * 3^1 is a term since none of the exponents in its prime factorization is equal to 2.
9 = 3^2 is not a term since it has an exponent 2 in its prime factorization.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Ertan Elma and Greg Martin, Distribution of the number of prime factors with a given multiplicity, Canadian Mathematical Bulletin, Vol. 67, No. 4 (2024), pp. 1107-1122; arXiv preprint, arXiv:2406.04574 [math.NT], 2024.
D. Suryanarayana, Semi-k-free integers, Elemente der Mathematik, Vol. 26 (1971), pp. 39-40.
D. Suryanarayana and R. Sitaramachandra Rao, Distribution of semi-k-free integers, Proceedings of the American Mathematical Society, Vol. 37, No. 2 (1973), pp. 340-346.

Complement of A038109.
Cf. A003557, A007947, A057521, A304364, A330596.
A005117, A036537, A036966, A048109, A175496, A268335 and A336590 are subsequences.
Numbers without an exponent k in their prime factorization: A001694 (k=1), this sequence (k=2), A386799 (k=3), A386803 (k=4), A386807 (k=5).
Numbers that have exactly m exponents in their prime factorization that are equal to 2: this sequence (m=0), A386796 (m=1), A386797 (m=2), A386798 (m=3).

q:= n-> andmap(i-> i[2]<>2, ifactors(n)[2]):
select(q, [$1..100])[];  # Alois P. Heinz, Aug 12 2020

Select[Range[100], !MemberQ[FactorInteger[#][[;;, 2]], 2] &]

is(n) = {my(f = factor(n)); for(i = 1, #f~, if(f[i, 2] == 2, return(0))); 1; } \\ Amiram Eldar, Oct 21 2023

Sum_{n>=1} 1/a(n)^s = zeta(s) * Product_{p prime} (1 - 1/p^(2*s) + 1/p^(3*s)), for s > 1. - Amiram Eldar, Oct 21 2023

A384050 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is a powerful number.

Original entry on oeis.org

1, 1, 2, 4, 4, 2, 6, 8, 9, 4, 10, 8, 12, 6, 8, 16, 16, 9, 18, 16, 12, 10, 22, 16, 25, 12, 27, 24, 28, 8, 30, 32, 20, 16, 24, 36, 36, 18, 24, 32, 40, 12, 42, 40, 36, 22, 46, 32, 49, 25, 32, 48, 52, 27, 40, 48, 36, 28, 58, 32, 60, 30, 54, 64, 48, 20, 66, 64, 44

Offset: 1

Amiram Eldar, May 18 2025

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

Unitary analog of A384039.
Cf. A000010, A001694, A055231, A057521, A330596, A384046, A384047.
The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is: A047994 (1), A384048 (squarefree), A384049 (cubefree), this sequence (powerful), A384051 (cubefull), A384052 (square), A384053 (cube), A384054 (exponentially odd), A384055 (odd), A384056 (power of 2), A384057 (3-smooth), A384058 (5-rough).

f[p_, e_] := p^e - If[e < 2, 1, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a,100]

a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^f[i,2] - if(f[i,2] == 1, 1, 0));}

Multiplicative with a(p) = p-1, and p^e if e >= 2.
a(n) = n * A047994(n) / A384048(n).
a(n) = A047994(A055231(n)) * A057521(n) = A000010(A055231(n)) * A057521(n).
Dirichlet g.f.: zeta(s-1) * Product_{p prime} (1 - 1/p^s + 1/p^(2*s-1)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{p prime} (1 - 1/p^2 + 1/p^3) = 0.748535... (A330596).

A038109 Divisible exactly by the square of a prime.

Original entry on oeis.org

4, 9, 12, 18, 20, 25, 28, 36, 44, 45, 49, 50, 52, 60, 63, 68, 72, 75, 76, 84, 90, 92, 98, 99, 100, 108, 116, 117, 121, 124, 126, 132, 140, 144, 147, 148, 150, 153, 156, 164, 169, 171, 172, 175, 180, 188, 196, 198, 200, 204, 207, 212, 220, 225, 228, 234, 236, 242

Offset: 1

Felice Russo

Numbers for which at least one prime factor exponent is exactly 2.
Sometimes called squarefull numbers, although that term is usually reserved for A001694. - N. J. A. Sloane, Jul 22 2012
The asymptotic density of this sequence is 1 - A330596 = 0.2514647... - Amiram Eldar, Aug 12 2020

			20=5*2*2 is divisible by 2^2.

R. J. Mathar, Table of n, a(n) for n = 1..1247

Cf. A001694, A013929, A284017, A284018.

isA038109 := proc(n)
    local p;
    for p in ifactors(n)[2] do
        if op(2,p) = 2 then
            return true;
        end if;
    end do:
    false ;
end proc: # R. J. Mathar, Dec 08 2015
# second Maple program:
q:= n-> ormap(i-> i[2]=2, ifactors(n)[2]):
select(q, [$1..300])[];  # Alois P. Heinz, Aug 12 2020

Select[Range[250],MemberQ[Transpose[FactorInteger[#]][[2]],2]&] (* Harvey P. Dale, Sep 24 2012 *)

is(n)=#select(n->n==2, Set(factor(n)[,2])) \\ Charles R Greathouse IV, Sep 17 2015

Corrected and extended by Erich Friedman

A053149 Smallest cube divisible by n.

Original entry on oeis.org

1, 8, 27, 8, 125, 216, 343, 8, 27, 1000, 1331, 216, 2197, 2744, 3375, 64, 4913, 216, 6859, 1000, 9261, 10648, 12167, 216, 125, 17576, 27, 2744, 24389, 27000, 29791, 64, 35937, 39304, 42875, 216, 50653, 54872, 59319, 1000, 68921, 74088, 79507, 10648

Offset: 1

Henry Bottomley, Feb 28 2000

Amiram Eldar, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.

Cf. A000188, A000189, A007913, A008834, A019554, A048798, A050985.

Cf. A002117, A013667, A330596.

a[n_] := For[k = 1, True, k++, If[ Divisible[c = k^3, n], Return[c]]]; Table[a[n], {n, 1, 44}] (* Jean-François Alcover, Sep 03 2012 *)
f[p_, e_] := p^(e + Mod[3 - e, 3]); a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Aug 29 2019 *)
scdn[n_]:=Module[{c=Ceiling[Surd[n,3]]},While[!Divisible[c^3,n],c++];c^3]; Array[scdn,50] (* Harvey P. Dale, Jun 13 2020 *)

a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^(f[i,2] + (3-f[i,2])%3));} \\ Amiram Eldar, Oct 27 2022

a(n) = (n/A000189(n))^3 = A008834(n)*A019554(A050985(n))^3 = n*A050985(n)^2/A000188(A050985(n))^3.
a(n) = n * A048798(n). - Franklin T. Adams-Watters, Apr 08 2009
From Amiram Eldar, Jul 29 2022: (Start)
Multiplicative with a(p^e) = p^(e + ((3-e) mod 3)).
Sum_{n>=1} 1/a(n) = Product_{p prime} ((p^3+2)/(p^3-1)) = 1.655234386560802506... . (End)
Sum_{k=1..n} a(k) ~ c * n^4, where c = (zeta(9)/(4*zeta(3))) * Product_{p prime} (1 - 1/p^2 + 1/p^3) = A013667*A330596/(4*A002117) = 0.1559906... . - Amiram Eldar, Oct 27 2022

A330595 Decimal expansion of Product_{primes p} (1 + 1/p^2 + 1/p^3).

Original entry on oeis.org

1, 7, 4, 8, 9, 3, 2, 9, 9, 7, 8, 4, 3, 2, 4, 5, 3, 0, 3, 0, 3, 3, 9, 0, 6, 9, 9, 7, 6, 8, 5, 1, 1, 4, 8, 0, 2, 2, 5, 9, 8, 8, 3, 4, 9, 3, 5, 9, 5, 4, 8, 0, 8, 9, 7, 2, 7, 3, 6, 6, 2, 1, 4, 4, 0, 8, 4, 8, 4, 9, 7, 9, 1, 3, 0, 0, 1, 0, 1, 3, 1, 4, 0, 6, 8, 1, 7, 8, 1, 3, 0, 2, 6, 4, 5, 5, 1, 0, 8, 9, 7, 0, 5, 9, 1

Offset: 1

Vaclav Kotesovec, Dec 19 2019

			1.748932997843245303033906997685114802259883493595480897273662144084849...

Cf. A065464, A065473, A065476, A160889, A328017, A330594, A330596, A338325.

Do[Print[N[Exp[-Sum[q = Expand[(-p^2 - p^3)^j]; Sum[PrimeZetaP[Exponent[q[[k]], p]] * Coefficient[q[[k]], p^Exponent[q[[k]], p]], {k, 1, Length[q]}]/j, {j, 1, t}]], 110]], {t, 20, 200, 20}]

prodeulerrat(1 + 1/p^2 + 1/p^3) \\ Vaclav Kotesovec, Sep 19 2020

Equals Sum_{n>=1} 1/A338325(n). - Amiram Eldar, Oct 26 2020

A386796 Numbers that have exactly one exponent in their prime factorization that is equal to 2.

Original entry on oeis.org

4, 9, 12, 18, 20, 25, 28, 44, 45, 49, 50, 52, 60, 63, 68, 72, 75, 76, 84, 90, 92, 98, 99, 108, 116, 117, 121, 124, 126, 132, 140, 144, 147, 148, 150, 153, 156, 164, 169, 171, 172, 175, 188, 198, 200, 204, 207, 212, 220, 228, 234, 236, 242, 244, 245, 260, 261, 268

Offset: 1

Amiram Eldar, Aug 02 2025

First differs from its subsequence A060687 at n = 16: a(16) = 72 is not a term of A060687.
Differs from A286228 by having the terms 60, 72, 84, 90, ..., and not having the term 1.
Numbers k such that A369427(k) = 1.
The asymptotic density of this sequence is Product_{p primes} (1 - 1/p^2 + 1/p^3) * Sum_{p prime} (p-1)/(p^3 - p + 1) = 0.22661832022705616779... (the product is A330596) (Elma and Martin, 2024).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Ertan Elma and Greg Martin, Distribution of the number of prime factors with a given multiplicity, Canadian Mathematical Bulletin, Vol. 67, No. 4 (2024), pp. 1107-1122; arXiv preprint, arXiv:2406.04574 [math.NT], 2024.
Index entries for sequences computed from indices in prime factorization.

A060687 is a subsequence.
Cf. A286228, A330596, A369427.
Numbers that have exactly one exponent in their prime factorization that is equal to k: A119251 (k=1), this sequence (k=2), A386800 (k=3), A386804 (k=4), A386808 (k=5).
Numbers that have exactly m exponents in their prime factorization that are equal to 2: A337050 (m=0), this sequence (m=1), A386797 (m=2), A386798 (m=3).

f[p_, e_] := If[e == 2, 1, 0]; s[1] = 0; s[n_] := Plus @@ f @@@ FactorInteger[n]; Select[Range[300], s[#] == 1 &]

isok(k) = vecsum(apply(x -> if(x == 2, 1, 0), factor(k)[, 2])) == 1;

A386797 Numbers that have exactly two exponents in their prime factorization that are equal to 2.

Original entry on oeis.org

36, 100, 180, 196, 225, 252, 300, 396, 441, 450, 468, 484, 588, 612, 676, 684, 700, 828, 882, 980, 1044, 1089, 1100, 1116, 1156, 1225, 1260, 1300, 1332, 1444, 1452, 1476, 1521, 1548, 1575, 1692, 1700, 1800, 1900, 1908, 1980, 2028, 2100, 2116, 2124, 2156, 2178, 2196

Offset: 1

Amiram Eldar, Aug 02 2025

First differs from its subsequence A375144 at n = 38: a(38) = 1800 = 2^3 * 3^2 * 5^2 is not a term of A375144.
Numbers k such that A369427(k) = 2.
The asymptotic density of this sequence is Product_{p primes} (1 - 1/p^2 + 1/p^3) * ((Sum_{p prime} (p-1)/(p^3 - p + 1))^2 - Sum_{p prime} ((p-1)^2/(p^3 - p + 1)^2)) / 2 = 0.023701044250873975412... (the product is A330596) (Elma and Martin, 2024).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Ertan Elma and Greg Martin, Distribution of the number of prime factors with a given multiplicity, Canadian Mathematical Bulletin, Vol. 67, No. 4 (2024), pp. 1107-1122; arXiv preprint, arXiv:2406.04574 [math.NT], 2024.
Index entries for sequences computed from indices in prime factorization.

A375144 is a subsequence.	
Cf. A330596, A369427.
Numbers that have exactly two exponents in their prime factorization that are equal to k: this sequence (k=2), A386801 (k=3), A386805 (k=4), A386809 (k=5).
Numbers that have exactly m exponents in their prime factorization that are equal to 2: A337050 (m=0), A386796 (m=1), this sequence (m=2), A386798 (m=3).

f[p_, e_] := If[e == 2, 1, 0]; s[1] = 0; s[n_] := Plus @@ f @@@ FactorInteger[n]; Select[Range[2200], s[#] == 2 &]

isok(k) = vecsum(apply(x -> if(x == 2, 1, 0), factor(k)[, 2])) == 2;

A386798 Numbers that have exactly three exponents in their prime factorization that are equal to 2.

Original entry on oeis.org

900, 1764, 4356, 4900, 6084, 6300, 8820, 9900, 10404, 11025, 11700, 12100, 12996, 14700, 15300, 16900, 17100, 19044, 19404, 20700, 21780, 22050, 22932, 23716, 26100, 27225, 27900, 28900, 29988, 30276, 30420, 30492, 33124, 33300, 33516, 34596, 36100, 36300, 36900, 38025, 38700

Offset: 1

Amiram Eldar, Aug 02 2025

Numbers k such that A369427(k) = 2.

The asymptotic density of this sequence is Product_{p primes} (1 - 1/p^2 + 1/p^3) * (s(1)^3 + 3*s(1)*s(2) + 2*s(3)) / 6 = 0.0011175284878980531468... (the product is A330596), where s(m) = (-1)^(m-1) * Sum_{p prime} (1/(p^3/(p-1)-1))^m (Elma and Martin, 2024).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Ertan Elma and Greg Martin, Distribution of the number of prime factors with a given multiplicity, Canadian Mathematical Bulletin, Vol. 67, No. 4 (2024), pp. 1107-1122; arXiv preprint, arXiv:2406.04574 [math.NT], 2024.
Index entries for sequences computed from indices in prime factorization.

Cf. A330596, A369427.
Numbers that have exactly three exponents in their prime factorization that are equal to k: this sequence (k=2), A386802 (k=3), A386806 (k=4), A386810 (k=5).
Numbers that have exactly m exponents in their prime factorization that are equal to 2: A337050 (m=0), A386796 (m=1), A386797 (m=2), this sequence (m=3).

f[p_, e_] := If[e == 2, 1, 0]; s[1] = 0; s[n_] := Plus @@ f @@@ FactorInteger[n]; Select[Range[40000], s[#] == 3 &]

isok(k) = vecsum(apply(x -> if(x == 2, 1, 0), factor(k)[, 2])) == 3;

A304411 If n = Product (p_j^k_j) then a(n) = Product ((p_j + 1)*k_j).

Original entry on oeis.org

1, 3, 4, 6, 6, 12, 8, 9, 8, 18, 12, 24, 14, 24, 24, 12, 18, 24, 20, 36, 32, 36, 24, 36, 12, 42, 12, 48, 30, 72, 32, 15, 48, 54, 48, 48, 38, 60, 56, 54, 42, 96, 44, 72, 48, 72, 48, 48, 16, 36, 72, 84, 54, 36, 72, 72, 80, 90, 60, 144, 62, 96, 64, 18, 84, 144, 68, 108, 96, 144, 72, 72

Offset: 1

Ilya Gutkovskiy, May 12 2018

			a(24) = a(2^3*3) = (2 + 1)*3 * (3 + 1)*1 = 36.

Andrew Howroyd, Table of n, a(n) for n = 1..1000
Ilya Gutkovskiy, Scatter plot of a(n) up to n=50000 .
Index entries for sequences computed from exponents in factorization of n.
Index entries for sequences computed from indices in prime factorization.

Cf. A000005, A000026, A000040, A000203, A003959, A001615, A003961, A005117, A005361, A007947, A008864, A045965, A048250, A064478, A081294 (numbers n such that a(n) is odd), A181797, A304407, A304408, A304409, A304412.

Cf. A072691, A330596.

a[n_] := Times @@ ((#[[1]] + 1) #[[2]] & /@ FactorInteger[n]); a[1] = 1; Table[a[n], {n, 72}]
Table[Total[Select[Divisors[n], SquareFreeQ]] DivisorSigma[0, n/Last[Select[Divisors[n], SquareFreeQ]]], {n, 72}]

a(n)={my(f=factor(n)); prod(i=1, #f~, my(p=f[i,1], e=f[i,2]); (p+1)*e)} \\ Andrew Howroyd, Jul 24 2018

a(n) = A005361(n)*A048250(n) = A000005(n/A007947(n))*A000203(A007947(n)).
a(p^k) = (p + 1)*k where p is a prime and k > 0.
a(n) = Product_{p|n} (p + 1) if n is a squarefree (A005117).
Sum_{k=1..n} a(k) ~ c * n^2, where c = (Pi^2/12) * Product_{p prime} (1 - 1/p^2 + 1/p^3) = A072691 * A330596 = 0.6156455744... . - Amiram Eldar, Nov 30 2022

cp's OEIS Frontend

A057723 Sum of positive divisors of n that are divisible by every prime that divides n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A337050 Numbers without an exponent 2 in their prime factorization.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A384050 The number of integers k from 1 to n such that the greatest divisor of k that is a unitary divisor of n is a powerful number.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A038109 Divisible exactly by the square of a prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A053149 Smallest cube divisible by n.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A330595 Decimal expansion of Product_{primes p} (1 + 1/p^2 + 1/p^3).

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

PARI

Formula

A386796 Numbers that have exactly one exponent in their prime factorization that is equal to 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI