A065487 -id:A065487 - cp's OEIS frontend

A003557 n divided by largest squarefree divisor of n; if n = Product p(k)^e(k) then a(n) = Product p(k)^(e(k)-1), with a(1) = 1.

1, 1, 1, 2, 1, 1, 1, 4, 3, 1, 1, 2, 1, 1, 1, 8, 1, 3, 1, 2, 1, 1, 1, 4, 5, 1, 9, 2, 1, 1, 1, 16, 1, 1, 1, 6, 1, 1, 1, 4, 1, 1, 1, 2, 3, 1, 1, 8, 7, 5, 1, 2, 1, 9, 1, 4, 1, 1, 1, 2, 1, 1, 3, 32, 1, 1, 1, 2, 1, 1, 1, 12, 1, 1, 5, 2, 1, 1, 1, 8, 27, 1, 1, 2, 1, 1, 1, 4, 1, 3, 1, 2, 1, 1, 1, 16, 1, 7

Offset: 1

Marc LeBrun

a(n) is the size of the Frattini subgroup of the cyclic group C_n - Ahmed Fares (ahmedfares(AT)my-deja.com), Jun 07 2001.
Also of the Frattini subgroup of the dihedral group with 2*n elements. - Sharon Sela (sharonsela(AT)hotmail.com), Jan 01 2002
Number of solutions to x^m==0 (mod n) provided that n < 2^(m+1), i.e. the sequence of sequences A000188, A000189, A000190, etc. converges to this sequence. - Henry Bottomley, Sep 18 2001
a(n) is the number of nilpotent elements in the ring Z/nZ. - Laszlo Toth, May 22 2009
The sequence of partial products of a(n) is A085056(n). - Peter Luschny, Jun 29 2009
The first occurrence of n in this sequence is at A064549(n). - Franklin T. Adams-Watters, Jul 25 2014
From Hal M. Switkay, Jul 03 2025: (Start)
For n > 1, a(n) is a proper divisor of n. Thus the sequence n, a(n), a(a(n)), ... eventually becomes 1. This yields a minimal factorization of n as a product of squarefree numbers (A005117), each factor dividing all larger factors, in a factorization that is conjugate to the minimal factorization of n as a product of prime powers (A000961), as follows.
Let f(n,0) = n, and let f(n,k) = a(f(n,k-1)) for k > 0. A051903(n) is the minimal value of k such that f(n,k) = 1. A051903(n) <= log(n)/log(2). Since n/a(n) = A007947(n) is always squarefree by definition, n is a product of squarefree factors in the form Product_{i=1..A051903(n)} [f(n,i-1)/f(n,i)].
The two factorizations correspond to conjugate partitions of bigomega(n) = A001222(n). (End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Steven R. Finch, Idempotents and Nilpotents Modulo n, arXiv:math/0605019 [math.NT], 2006-2017.

Cf. A007947, A062378, A062379, A064549, A300717 (Möbius transform), A326306 (inv. Möbius transf.), A328572.
Sequences that are multiples of this sequence (the other factor of a pointwise product is given in parentheses): A000010 (A173557), A000027 (A007947), A001615 (A048250), A003415 (A342001), A007434 (A345052), A057521 (A071773).
Cf. A082695 (Dgf at s=2), A065487 (Dgf at s=3).
Cf. A000961, A001222, A005117, A051903.

a003557 n = product $ zipWith (^)
                      (a027748_row n) (map (subtract 1) $ a124010_row n)
-- Reinhard Zumkeller, Dec 20 2013

using Nemo
function A003557(n)
    n < 4 && return 1
    q = prod([p for (p, e) ∈ Nemo.factor(fmpz(n))])
    return n == q ? 1 : div(n, q)
end
[A003557(n) for n in 1:90] |> println  # Peter Luschny, Feb 07 2021

[(&+[(Floor(k^n/n)-Floor((k^n-1)/n)): k in [1..n]]): n in [1..100]]; // G. C. Greubel, Nov 02 2018

A003557 := n -> n/ilcm(op(numtheory[factorset](n))):
seq(A003557(n), n=1..98); # Peter Luschny, Mar 23 2011
seq(n / NumberTheory:-Radical(n), n = 1..98); # Peter Luschny, Jul 20 2021

Prepend[ Array[ #/Times@@(First[ Transpose[ FactorInteger[ # ] ] ])&, 100, 2 ], 1 ] (* Olivier Gérard, Apr 10 1997 *)

a(n)=n/factorback(factor(n)[,1]) \\ Charles R Greathouse IV, Nov 17 2014

for(n=1, 100, print1(direuler(p=2, n, (1 - p*X + X)/(1 - p*X))[n], ", ")) \\ Vaclav Kotesovec, Jun 20 2020

from sympy.ntheory.factor_ import core
from sympy import divisors
def a(n): return n / max(i for i in divisors(n) if core(i) == i)
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Apr 16 2017

from math import prod
from sympy import primefactors
def A003557(n): return n//prod(primefactors(n)) # Chai Wah Wu, Nov 04 2022

def A003557(n) : return n*mul(1/p for p in prime_divisors(n))
[A003557(n) for n in (1..98)] # Peter Luschny, Jun 10 2012

Multiplicative with a(p^e) = p^(e-1). - Vladeta Jovovic, Jul 23 2001
a(n) = n/rad(n) = n/A007947(n) = sqrt(J_2(n)/J_2(rad(n))), where J_2(n) is A007434. - Enrique Pérez Herrero, Aug 31 2010
a(n) = (J_k(n)/J_k(rad(n)))^(1/k), where J_k is the k-th Jordan Totient Function: (J_2 is A007434 and J_3 A059376). - Enrique Pérez Herrero, Sep 03 2010
Dirichlet convolution of A000027 and A097945. - R. J. Mathar, Dec 20 2011
a(n) = A000010(n)/|A023900(n)|. - Eric Desbiaux, Nov 15 2013
a(n) = Product_{k = 1..A001221(n)} (A027748(n,k)^(A124010(n,k)-1)). - Reinhard Zumkeller, Dec 20 2013
a(n) = Sum_{k=1..n}(floor(k^n/n)-floor((k^n-1)/n)). - Anthony Browne, May 11 2016
a(n) = e^[Sum_{k=2..n} (floor(n/k)-floor((n-1)/k))*(1-A010051(k))*Mangoldt(k)] where Mangoldt is the Mangoldt function. - Anthony Browne, Jun 16 2016
a(n) = Sum_{d|n} mu(d) * phi(d) * (n/d), where mu(d) is the Moebius function and phi(d) is the Euler totient function (rephrases formula of Dec 2011). - Daniel Suteu, Jun 19 2018
G.f.: Sum_{k>=1} mu(k)*phi(k)*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Nov 02 2018
Dirichlet g.f.: Product_{primes p} (1 + 1/(p^s - p)). - Vaclav Kotesovec, Jun 24 2020
From Richard L. Ollerton, May 07 2021: (Start)
a(n) = Sum_{k=1..n} mu(n/gcd(n,k))*gcd(n,k).
a(n) = Sum_{k=1..n} mu(gcd(n,k))*(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). (End)
a(n) = A001615(n)/A048250(n) = A003415/A342001(n) = A057521(n)/A071773(n). - Antti Karttunen, Jun 08 2021

Secondary definition added to the name by Antti Karttunen, Jun 08 2021

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

Original entry on oeis.org

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

A335988 Cubefull exponentially odd numbers: numbers whose prime factorization contains only odd exponents that are larger than 1.

Original entry on oeis.org

1, 8, 27, 32, 125, 128, 216, 243, 343, 512, 864, 1000, 1331, 1944, 2048, 2187, 2197, 2744, 3125, 3375, 3456, 4000, 4913, 6859, 7776, 8192, 9261, 10648, 10976, 12167, 13824, 16000, 16807, 17496, 17576, 19683, 24389, 25000, 27000, 29791, 30375, 31104, 32768, 35937

Offset: 1

Amiram Eldar, Jul 03 2020

nonn

This sequence is a permutation of A355038.
This sequence is also a permutation of the exponentially odd numbers (A268335) multiplied by the square of their squarefree kernel (A007947).
a(n)/rad(a(n)) is a permutation of the squares.
a(n)/rad(a(n))^2 is a permutation of the exponentially odd numbers.

			8 = 2^3 is a term since the exponent of its prime factor 2 is 3 which is odd and larger than 1.

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

Intersection of A001694 and A268335.
Intersection of A036966 and A268335.
A355038 in ascending order.
A030078, A050997, A092759, A179665, A079395 and A138031 are subsequences.
Cf. A007947, A065487.

Join[{1}, Select[Range[10^5], AllTrue[Last /@ FactorInteger[#], #1 > 1 && OddQ[#1] &] &]]

from math import isqrt, prod
from sympy import factorint
def afind(N): # all terms up to limit N
    cands = (n**2*prod(factorint(n**2)) for n in range(1, isqrt(N//2)+2))
    return sorted(c for c in cands if c <= N)
print(afind(4*10**4)) # Michael S. Branicky, Jun 16 2022

Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + 1/(p*(p^2-1))) = 1.2312911... (A065487).

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

Original entry on oeis.org

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

Offset: 0

Vaclav Kotesovec, Dec 19 2019

The asymptotic density of A337050. - Amiram Eldar, Aug 13 2020

			0.748535259682363564644215048637910601641640343005324404515852793925925...

Cf. A057723, A065464, A065473, A065476, A065487, A328017, A330594, A330595, A337050, A372636.

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) \\ Amiram Eldar, Mar 17 2021

Equals (6/Pi^2) * A065487. - Amiram Eldar, Jun 10 2020

A127917 Product of three numbers: n-th prime, previous number, and following number.

Original entry on oeis.org

6, 24, 120, 336, 1320, 2184, 4896, 6840, 12144, 24360, 29760, 50616, 68880, 79464, 103776, 148824, 205320, 226920, 300696, 357840, 388944, 492960, 571704, 704880, 912576, 1030200, 1092624, 1224936, 1294920, 1442784, 2048256, 2247960, 2571216, 2685480, 3307800

Offset: 1

Artur Jasinski, Feb 06 2007

a(n) is the order of the matrix group SL(2,prime(n)). - Tom Edgar, Sep 28 2015

G. C. Greubel, Table of n, a(n) for n = 1..10000
J. B. Marshall, On the extension of Fermat's theorem to matrices of order n, Proceedings of the Edinburgh Mathematical Society 6 (1939) 85-91. See (11) page 90-91 when p=2.
Index to sequences related to prime powers.

Cf. A036689, A034953, A084920, A117762.

Cf. A065470, A065487.

[6] cat [NthPrime(n)*(NthPrime(n)^2-1): n in [2..40]]; // Vincenzo Librandi, Sep 29 2015

Table[(Prime[n] + 1) Prime[n](Prime[n] - 1), {n, 1, 100}]
Table[p(p^2-1),{p,Prime[Range[40]]}] (* Harvey P. Dale, Apr 26 2025 *)

forprime(p=2,1e3,print1(6*binomial(p+1,3)", ")) \\ Charles R Greathouse IV, Jun 16 2011

a(n) = prime(n)*(prime(n)^2-1);
vector(40, n, a(n)) \\ Altug Alkan, Sep 28 2015

a(n) = prime(n)*(prime(n)^2-1). - Tom Edgar, Sep 28 2015
a(n) = 2 * A117762(n), for n > 1. - Altug Alkan, Sep 28 2015
From Amiram Eldar, Nov 22 2022: (Start)
Product_{n>=1} (1 + 1/a(n)) = A065487.
Product_{n>=1} (1 - 1/a(n)) = A065470. (End)

A325837 The number of coreful divisors of n that are exponentially odd numbers (A268335).

Original entry on oeis.org

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

Offset: 1

Amiram Eldar, Sep 07 2019

First differs from A050361 at n = 64.
From Amiram Eldar, Sep 08 2023: (Start)
The number of exponentially odd divisors of n is A322483(n), and their sum is A033634(n).
A coreful divisor d of a number n is a divisor with the same set of distinct prime factors as n. (End)
Also, the number of divisors of n that are cubefull exponentially odd numbers (A335988). - Amiram Eldar, Feb 11 2024

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

Cf. A033634, A050361, A065487, A268335, A295316, A335988, A350390.

Cf. A003557, A005361 (number of coreful divisors), A046951, A268335.

fun[p_,e_] := Floor[(e+1)/2]; a[n_] := Times@@(fun@@@FactorInteger[n]); Array[a, 100]

a(n) = vecprod(apply(x -> (x+1)\2, factor(n)[, 2])); \\ Amiram Eldar, Sep 01 2023

Multiplicative with a(p^e) = floor((e+1)/2).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} (1 + 1/(p*(p^2-1))) = 1.231291... (A065487). - Amiram Eldar, Sep 10 2022
a(n) = A046951(A350390(n)) (the number of squares dividing the largest exponentially odd divisor of n). - Amiram Eldar, Sep 01 2023
From Amiram Eldar, Sep 08 2023: (Start)
a(n) = A046951(A003557(n)).
Dirichlet g.f.: zeta(s) * zeta(2*s) * Product_{p prime} (1 - 1/p^(2*s) + 1/p^(3*s)). (End)

Name corrected by Amiram Eldar, Sep 08 2023

A336563 Sum of proper divisors of n that are divisible by every prime that divides n.

Original entry on oeis.org

0, 0, 0, 2, 0, 0, 0, 6, 3, 0, 0, 6, 0, 0, 0, 14, 0, 6, 0, 10, 0, 0, 0, 18, 5, 0, 12, 14, 0, 0, 0, 30, 0, 0, 0, 36, 0, 0, 0, 30, 0, 0, 0, 22, 15, 0, 0, 42, 7, 10, 0, 26, 0, 24, 0, 42, 0, 0, 0, 30, 0, 0, 21, 62, 0, 0, 0, 34, 0, 0, 0, 96, 0, 0, 15, 38, 0, 0, 0, 70, 39, 0, 0, 42, 0, 0, 0, 66, 0, 30, 0, 46, 0, 0, 0, 90

Offset: 1

Antti Karttunen, Jul 27 2020

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10800
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537

Cf. A001065, A003557, A007947, A057723, A033879, A065487, A335341, A336564, A336565, A336566, A336567.

Cf. A005117 (positions of zeros).

f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; a[1] = 0; a[n_] := Times @@ f @@@ FactorInteger[n] - n; Array[a, 100] (* Amiram Eldar, May 06 2023 *)

A007947(n) = factorback(factorint(n)[, 1]);
A057723(n) = { my(r=A007947(n)); (r*sigma(n/r)); };
A336563(n) = (A057723(n)-n);
\\ Or just as:
A336563(n) = { my(x=A007947(n),y = n/x); (x*(sigma(y)-y)); };

a(n) = A057723(n) - n.
a(n) = A007947(n) * A336567(n) = A007947(n) * A001065(A003557(n)).
a(n) = A336564(n) - A033879(n).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = A065487 - 1 = 0.231291... . - Amiram Eldar, Dec 07 2023

A308135 Sum of non-coreful divisors of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 6, 1, 1, 1, 8, 1, 10, 1, 10, 9, 1, 1, 15, 1, 12, 11, 14, 1, 18, 1, 16, 1, 14, 1, 42, 1, 1, 15, 20, 13, 19, 1, 22, 17, 20, 1, 54, 1, 18, 18, 26, 1, 34, 1, 33, 21, 20, 1, 42, 17, 22, 23, 32, 1, 78, 1, 34, 20, 1, 19, 78, 1, 24, 27, 74, 1, 27, 1, 40

Offset: 1

Amiram Eldar and Paolo P. Lava, May 14 2019

nonn

Non-coreful divisor d of a number k is a divisor such that rad(d) != rad(k), where rad(k) is the largest squarefree divisor of k (A007947).

			a(15) = 9. Prime factors of 15 are 3, 5 and its divisors are 1, 3, 5, 15. The non-coreful divisors are 1, 3, 5 and their sum is 9.

Amiram Eldar, Table of n, a(n) for n = 1..10000
G. E. Hardy and M. V. Subbarao, Highly powerful numbers, Congress. Numer., Vol. 37 (1983), pp. 277-307. (Annotated scanned copy)

Cf. A000203, A005101, A007947, A057723, A307888, A307958, A308029, A308127.

Cf. A013661, A065487.

with(numtheory): P:=proc(k) local a,n; a:=mul(n,n=factorset(k));
sigma(k)-a*sigma(k/a); end: seq(P(i),i=1..74);

f[p_, e_] := (p^(e + 1) - 1)/(p - 1); fc[p_, e_] := f[p, e] - 1; a[1] = 0; a[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (fc @@@ FactorInteger[n]); Array[a, 100]

a(n) = A000203(n) - A057723(n).

Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = A013661 - A065487 = 0.413642... . - Amiram Eldar, Dec 08 2023

A336564 a(n) = n - A308135(n), where A308135(n) is the sum of non-coreful divisors of n.

Original entry on oeis.org

1, 1, 2, 3, 4, 0, 6, 7, 8, 2, 10, 2, 12, 4, 6, 15, 16, 3, 18, 8, 10, 8, 22, 6, 24, 10, 26, 14, 28, -12, 30, 31, 18, 14, 22, 17, 36, 16, 22, 20, 40, -12, 42, 26, 27, 20, 46, 14, 48, 17, 30, 32, 52, 12, 38, 34, 34, 26, 58, -18, 60, 28, 43, 63, 46, -12, 66, 44, 42, -4, 70, 45, 72, 34, 41, 50, 58, -12, 78, 44, 80, 38, 82, -14, 62

Offset: 1

Antti Karttunen, Jul 27 2020

sign

Antti Karttunen, Table of n, a(n) for n = 1..20000

Cf. A000203, A033879, A057723, A308135, A336563, A336565, A336566.

Cf. A013661, A065487.

f[p_, e_] := (p^(e + 1) - 1)/(p - 1); fc[p_, e_] := f[p, e] - 1; a[1] = 1; a[n_] := n - Times @@ f @@@ (fct = FactorInteger[n]) + Times @@ fc @@@ fct; Array[a, 100] (* Amiram Eldar, Dec 08 2023 *)

A007947(n) = factorback(factorint(n)[, 1]);
A057723(n) = { my(r=A007947(n)); (r*sigma(n/r)); };
A308135(n) = (sigma(n)-A057723(n));
A336564(n) = (n - A308135(n));

a(n) = n - A308135(n) = n - (sigma(n) - A057723(n)).
a(n) = A336563(n) + A033879(n). [Corrected by Georg Fischer, Dec 13 2022]
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = A065487 - A013661 + 1 = 0.586357... . - Amiram Eldar, Dec 08 2023

A371413 Dedekind psi function applied to the cubefull numbers (A036966).

Original entry on oeis.org

1, 12, 24, 36, 48, 96, 108, 150, 192, 432, 324, 384, 392, 864, 768, 750, 1296, 972, 1728, 1800, 1536, 2592, 1452, 3456, 3888, 3600, 3072, 2916, 2366, 2744, 5184, 4704, 3750, 5400, 6912, 7776, 7200, 6144, 5202, 9000, 10368, 9408, 11664, 8748, 7220, 13824, 15552

Offset: 1

Amiram Eldar, Mar 22 2024

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

Cf. A001615, A036966, A065464.

Similar sequences: A323332, A371412, A371415.

psi[n_] := n * Times @@ (1 + 1/FactorInteger[n][[;; , 1]]); psi[1] = 1; Join[{1}, psi /@ Select[Range[20000], AllTrue[Last /@ FactorInteger[#], #1 > 2 &] &]]
(* or *)
f[n_] := Module[{f = FactorInteger[n], p, e}, If[n == 1, 1, p = f[[;;, 1]]; e = f[[;;, 2]]; If[Min[e] > 2, Times @@ ((p+1) * p^(e-1)), Nothing]]]; Array[f, 20000]

dedpsi(f) = prod(i = 1, #f~, (f[i, 1] + 1) * f[i, 1]^(f[i, 2]-1));
lista(max) = {my(f); print1(1, ", "); for(k = 2, max, f = factor(k); if(vecmin(f[, 2]) > 2, print1(dedpsi(f), ", "))); }

a(n) = A001615(A036966(n)).

Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + 1/((p^2-1)*p)) = 1.231291... (A065487).

cp's OEIS Frontend

A003557 n divided by largest squarefree divisor of n; if n = Product p(k)^e(k) then a(n) = Product p(k)^(e(k)-1), with a(1) = 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Julia

Magma

Maple

Mathematica

PARI

PARI

Python

Python

Sage

Formula

Extensions

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

A335988 Cubefull exponentially odd numbers: numbers whose prime factorization contains only odd exponents that are larger than 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Formula

A127917 Product of three numbers: n-th prime, previous number, and following number.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A325837 The number of coreful divisors of n that are exponentially odd numbers (A268335).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI