cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 23 results. Next

A073184 Number of cubefree divisors of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 3, 3, 4, 2, 6, 2, 4, 4, 3, 2, 6, 2, 6, 4, 4, 2, 6, 3, 4, 3, 6, 2, 8, 2, 3, 4, 4, 4, 9, 2, 4, 4, 6, 2, 8, 2, 6, 6, 4, 2, 6, 3, 6, 4, 6, 2, 6, 4, 6, 4, 4, 2, 12, 2, 4, 6, 3, 4, 8, 2, 6, 4, 8, 2, 9, 2, 4, 6, 6, 4, 8, 2, 6, 3, 4, 2, 12, 4, 4, 4, 6, 2, 12, 4, 6, 4, 4, 4, 6, 2, 6, 6, 9, 2, 8, 2
Offset: 1

Views

Author

Reinhard Zumkeller, Jul 19 2002

Keywords

Comments

a(n) = number of divisors of the cubefree kernel of n: a(n) = A000005(A007948(n)); [corrected by Amiram Eldar, Oct 08 2022]
Multiplicative because it is the Inverse Möbius transform of the characteristic function of cubefree numbers. a(n) is a prime signature sequence. a(p) = 2, a(p^e) = 3, e>1. - Christian G. Bower, May 18 2005

Examples

			The divisors of 56 are {1, 2, 4, 7, 8, 14, 28, 56}, 8=2^3 and 56=7*2^3 are not cubefree, therefore a(56) = 6.
		

Crossrefs

Programs

  • Haskell
    a073184 = sum . map a212793 . a027750_row
    -- Reinhard Zumkeller, May 27 2012
    
  • Mathematica
    a[1] = 1; a[p_?PrimeQ] = 2; a[n_] := Times @@ (If[#[[2]] == 1, 2, 3] & /@ FactorInteger[n]); Table[a[n], {n, 1, 103}] (* Jean-François Alcover, May 24 2012, after Christian G. Bower *)
  • PARI
    a(n) = {my(e = factor(n)[,2]); prod(i = 1, #e, if(e[i] == 1, 2, 3))}; \\ Amiram Eldar, Oct 08 2022

Formula

a(n) <= A073182(n).
Dirichlet g.f.: zeta(s)^2/zeta(3*s). Dirichlet convolution of the characteristic function of cubefree numbers by A000012. - R. J. Mathar, Apr 12 2011
a(n) = Sum_{k = 1..A000005(n)} A212793(A027750(n,k)). - Reinhard Zumkeller, May 27 2012
Sum_{k=1..n} a(k) ~ n / zeta(3) * (log(n) - 1 + 2*gamma - 3*zeta'(3)/zeta(3)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Jan 31 2019

A360539 a(n) is the cubefree part of n: the largest unitary divisor of n that is a cubefree number (A004709).

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Feb 11 2023

Keywords

Comments

Equivalently, a(n) is the least divisor d of n such that n/d is a cubefull number (A036966).

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e < 3, p^e, 1]; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n)); prod(i=1, #f~, if(f[i, 2] < 3, f[i, 1]^f[i, 2], 1));}
    
  • Python
    from math import prod
    from sympy import factorint
    def A360539(n): return prod(p**e for p,e in factorint(n).items() if e<=2) # Chai Wah Wu, Aug 06 2024

Formula

a(n) = 1 if and only if n is a cubefull number (A036966).
a(n) = n if and only if n is a cubefree number (A004709).
a(n) >= A055231(n) with equality if and only if n is in A337050.
a(n) = n/A360540(n).
Multiplicative with a(p^e) = p^e if e <= 2, and 1 otherwise.
Dirichlet g.f.: zeta(s) * Product_{p prime} (1 + p^(1-s) - p^(-s) + p^(2-2*s) - p^(1-2*s) - p^(2-3*s) + p^(-3*s)).
Sum_{k=1..n} a(k) ~ c * n^2, where c = (Pi^2/12) * Product_{p prime} (1 - 1/p^2 - 1/p^3 + 1/p^5 + 1/p^6 - 1/p^7) = 0.4213813264... .

A073185 Sum of cubefree divisors of n.

Original entry on oeis.org

1, 3, 4, 7, 6, 12, 8, 7, 13, 18, 12, 28, 14, 24, 24, 7, 18, 39, 20, 42, 32, 36, 24, 28, 31, 42, 13, 56, 30, 72, 32, 7, 48, 54, 48, 91, 38, 60, 56, 42, 42, 96, 44, 84, 78, 72, 48, 28, 57, 93, 72, 98, 54, 39, 72, 56, 80, 90, 60, 168, 62, 96, 104, 7, 84, 144, 68, 126, 96, 144, 72
Offset: 1

Views

Author

Reinhard Zumkeller, Jul 19 2002

Keywords

Comments

Sum of divisors of the cubefree kernel of n (see first formula).

Examples

			The divisors of 56 are {1, 2, 4, 7, 8, 14, 28, 56}, 8=2^3 and 56=7*2^3 are not cubefree, therefore a(56) = 1 + 2 + 4 + 7 + 14 + 28 = 56.
		

Crossrefs

Programs

  • Haskell
    a073185 = sum . filter ((== 1) . a212793) . a027750_row
    -- Reinhard Zumkeller, May 27 2012
    
  • Maple
    charFfree := proc(n,t) local f; for f in ifactors(n)[2] do if op(2,f) >= t then return 0 ; end if; end do: return 1 ; end proc:
    A073185 := proc(n) add( d*charFfree(d,3),d =numtheory[divisors](n) ); end proc: # R. J. Mathar, Apr 12 2011
  • Mathematica
    nn = 71;f[list_, i_] := list[[i]]; a =Table[If[Max[FactorInteger[n][[All, 2]]] <= 2, n, 0], {n, 1, nn}]; b = Table[1, {nn}]; Select[Table[DirichletConvolve[f[a, n], f[b, n], n, m], {m, 1, nn}], # > 0 &] (* Geoffrey Critzer, Mar 22 2015 *)
    f[p_, e_] := 1 + p + If[e > 1, p^2, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 14 2020 *)
  • PARI
    a(n) = {my(f=factor(n)); for (i=1, #f~, p = f[i,1]; if ((e=f[i,2]) == 1, f[i,1] = 1+p, f[i,1] = 1+p+p^2); f[i,2] = 1;); factorback(f);} \\ Michel Marcus, Feb 06 2015

Formula

a(n) = A000203(A007948(n)).
a(n) <= A073183(n).
Multiplicative with a(p) = 1+p, a(p^e) = 1 + p + p^2, e>1. - Christian G. Bower, May 18 2005
a(n) = sum(A212793(A027750(n,k)) * A027750(n,k): k=1..A000005(n)). - Reinhard Zumkeller, May 27 2012
Dirichlet g.f.: zeta(s)*zeta(s-1)/zeta(3s-3). - R. J. Mathar, Apr 12 2011
Sum_{k=1..n} a(k) ~ Pi^2 * n^2 / (12*Zeta(3)). - Vaclav Kotesovec, Feb 01 2019

Extensions

Incorrect comment removed by Álvar Ibeas, Feb 06 2015

A197863 Smallest powerful number that is a multiple of n.

Original entry on oeis.org

1, 4, 9, 4, 25, 36, 49, 8, 9, 100, 121, 36, 169, 196, 225, 16, 289, 36, 361, 100, 441, 484, 529, 72, 25, 676, 27, 196, 841, 900, 961, 32, 1089, 1156, 1225, 36, 1369, 1444, 1521, 200, 1681, 1764, 1849, 484, 225, 2116, 2209, 144, 49, 100, 2601, 676, 2809, 108, 3025
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Haskell
    a197863 n = product $
       zipWith (^) (a027748_row n) (map (max 2) $ a124010_row n)
    -- Reinhard Zumkeller, Jan 06 2012
  • Mathematica
    With[{pwrnos=Join[{1},Select[Range[5000],Min[Transpose[ FactorInteger[#]] [[2]]]>1&]]},Flatten[Table[Select[pwrnos,Divisible[#,n]&,1],{n,60}]]] (* Harvey P. Dale, Aug 14 2012 *)
    f[p_, e_] := p^Max[e, 2]; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jul 09 2022 *)
  • PARI
    a(n)=local(fm=factor(n));prod(k=1,matsize(fm)[1],fm[k,1]^max(fm[k,2],2))
    

Formula

Multiplicative with a(p^e) = p^max(e,2).
Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + (2*p-1)/(p^2*(p-1))) = 2.71098009471568319328... . - Amiram Eldar, Jul 29 2022
Sum_{k=1..n} a(k) ~ c * n^3, where c = (Pi^2/18) * Product_{p prime} (1 - 2/p^2 + 2/p^4 - 1/p^5) = 0.2165355664... . - Amiram Eldar, Nov 19 2022
a(n) = n * A055231(n). - Amiram Eldar, Sep 01 2023

A007424 a(n) = 1 if n is squarefree, otherwise 2.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 2
Offset: 1

Views

Author

Keywords

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    Table[If[SquareFreeQ[n],1,2],{n,100}] (* Harvey P. Dale, Jul 09 2014 *)
  • MuPAD
    func(2-abs(numlib::moebius(n)), n):
    
  • PARI
    A007424(n) = (2-issquarefree(n)); \\ Antti Karttunen, Nov 18 2017

Formula

a(n) = 2 - A008966(n). - Antti Karttunen, Nov 18 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2 - 6/Pi^2 = 1 + A229099. - Amiram Eldar, Sep 26 2022
a(n) = A051903(A007948(n)) = A130130(A051903(n)) for n >= 2. - Amiram Eldar, May 07 2024

A062378 n divided by largest cubefree factor of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 16, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 9, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 2
Offset: 1

Views

Author

Henry Bottomley, Jun 18 2001

Keywords

Comments

Numerator of n/rad(n)^2, where rad is the squarefree kernel of n (A007947), denominator: A055231. - Reinhard Zumkeller, Dec 10 2002

Crossrefs

Cf. A000189, A000578, A007948, A008834, A019555, A048798, A050985, A053149, A053150, A056551, A056552. See A003557 for squares and A062379 for 4th powers.
Differs from A073753 for the first time at n=90, where a(90) = 1, while A073753(90) = 3.

Programs

Formula

a(n) = n / A007948(n).
a(n) = A003557(A003557(n)). - Antti Karttunen, Nov 28 2017
Multiplicative with a(p^e) = p^max(e-2, 0). - Amiram Eldar, Sep 07 2020
Dirichlet g.f.: zeta(s-1) * Product_{p prime} (1 - 1/p^(s-1) + 1/p^s - 1/p^(2*s-1) + 1/p^(2*s)). - Amiram Eldar, Dec 07 2023

A071773 a(n) = gcd(rad(n), n/rad(n)), where rad(n) = A007947(n) is the squarefree kernel of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 1, 1, 1, 2, 5, 1, 3, 2, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 2, 7, 5, 1, 2, 1, 3, 1, 2, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 5, 2, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 1, 1, 1, 2, 1, 7, 3, 10, 1, 1, 1, 2
Offset: 1

Views

Author

Reinhard Zumkeller, Jun 24 2002

Keywords

Comments

n is squarefree iff a(n)=1.
Product of primes dividing n more than once. - Charles R Greathouse IV, Aug 08 2013
Squarefree kernel of the square part of n. - Peter Munn, Jun 12 2020

Crossrefs

Cf. A003415, A003557, A005117, A007947, A007948, A008833, A057521, A166486 (parity of terms), A359433 (Dirichlet inverse).
Cf. A065464.

Programs

Formula

a(n) = gcd(A007947(n), A003557(n)).
Multiplicative with p^e -> p^ceiling((e-1)/e), p prime.
a(n) = rad(n/rad(n)) = A007947(A003557(n)). - Velin Yanev, Antti Karttunen, Aug 20 2017, Nov 28 2017
a(n) = A007947(A057521(n)). - Antti Karttunen, Nov 28 2017
a(n) = A007947(A008833(n)). - Peter Munn, Jun 12 2020
a(n) = gcd(A003415(n), A007947(n)). - Antti Karttunen, Jan 02 2023
Dirichlet g.f.: zeta(s) * Product_{p prime} (1 + 1/p^(2*s-1) - 1/p^(2*s)). - Amiram Eldar, Nov 09 2023
From Vaclav Kotesovec, May 06 2025: (Start)
Let f(s) = Product_{p prime} (1 - 1/p^(2*s) + 1/p^(4*s-1) - 1/p^(4*s-2)).
Dirichlet g.f.: zeta(s) * zeta(2*s-1) * f(s).
Sum_{k=1..n} a(k) ~ f(1) * n * (log(n) + 3*gamma - 1 + f'(1)/f(1)) / 2, where
f(1) = A065464 = Product_{p prime} (1 - 2/p^2 + 1/p^3) = 0.428249505677094440218765707581823546121298513355936144031901379532123...
f'(1) = f(1) * Sum_{p prime} 2*(3*p-2)*log(p) / (p^3-2*p+1) = f(1) * 2.939073481649229666406787986900328729326669597518287791424059647447664...
and gamma is the Euler-Mascheroni constant A001620. (End)

A078615 a(n) = rad(n)^2, where rad is the squarefree kernel of n (A007947).

Original entry on oeis.org

1, 4, 9, 4, 25, 36, 49, 4, 9, 100, 121, 36, 169, 196, 225, 4, 289, 36, 361, 100, 441, 484, 529, 36, 25, 676, 9, 196, 841, 900, 961, 4, 1089, 1156, 1225, 36, 1369, 1444, 1521, 100, 1681, 1764, 1849, 484, 225, 2116, 2209, 36, 49, 100, 2601, 676, 2809, 36, 3025, 196
Offset: 1

Views

Author

Reinhard Zumkeller, Dec 10 2002

Keywords

Comments

It is conjectured that only 1 and 1782 satisfy a(k) = sigma(k). See Broughan link. - Michel Marcus, Feb 28 2019

Crossrefs

Programs

  • Maple
    a := n -> mul(f,f=map(x->x^2,select(isprime,divisors(n))));
    seq(a(n), n=1..56);  # Peter Luschny, Mar 30 2014
  • Mathematica
    a[n_] := Times @@ FactorInteger[n][[All, 1]]^2; Array[a, 60] (* Jean-François Alcover, Jun 04 2019 *)
  • PARI
    a(n)=my(f=factor(n)[,1]);prod(i=1,#f,f[i])^2 \\ Charles R Greathouse IV, Aug 06 2013

Formula

Multiplicative with a(p^e) = p^2. - Mitch Harris, May 17 2005
G.f.: Sum_{k>=1} mu(k)^2*J_2(k)*x^k/(1 - x^k), where J_2() is the Jordan function. - Ilya Gutkovskiy, Nov 06 2018
Sum_{k=1..n} a(k) ~ c * n^3, where c = (zeta(3)/3) * Product_{p prime} (1 - 1/p^2 - 1/p^3 + 1/p^4) = A002117 * A330523 / 3 = 0.214725... . - Amiram Eldar, Oct 30 2022
a(n) = Sum_{1 <= i, j <= n} ( mobius(n/gcd(i, j, n)) )^2. - Peter Bala, Jan 28 2024
a(n) = Sum_{d|n} mu(d)^2*J_2(d), where J_2 = A007434. - Ridouane Oudra, Jul 24 2025
a(n) = (-1)^omega(n) * Sum_{d|n} mu(d)*Psi_2(d), where omega = A001221 and Psi_2 = A065958. - Ridouane Oudra, Aug 01 2025

A058035 Largest 4th-power-free number dividing n.

Original entry on oeis.org

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

Views

Author

Henry Bottomley, Nov 16 2000

Keywords

Examples

			a(96) = 24 since the factors of 96 are {1,2,3,4,6,8,12,16,24,32,48,96} but 32, 48 and 96 all contain a 4th power factor (16).
		

Crossrefs

Programs

  • Haskell
    a058035 n = product $
       zipWith (^) (a027748_row n) (map (min 3) $ a124010_row n)
    -- Reinhard Zumkeller, Jan 06 2012
    
  • Mathematica
    f[p_, e_] := p^Min[e, 3]; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jul 09 2022 *)
  • PARI
    a(n) = my(f=factor(n)); for(k=1,#f~,f[k,2]=min(3,f[k,2])); factorback(f); \\ Michel Marcus, Sep 13 2017

Formula

Multiplicative with a(p^e) = p ^ min(e,3), p prime, e > 0. - Reinhard Zumkeller, Jan 06 2012
Sum_{k=1..n} a(k) ~ (1/2) * c * n^2, where c = Product_{p prime} (1 - 1/(p^3*(p+1))) = 0.947733... (A065466). - Amiram Eldar, Oct 13 2022

A073182 Number of divisors of n which are not greater than the cubefree kernel of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 3, 3, 4, 2, 6, 2, 4, 4, 3, 2, 6, 2, 6, 4, 4, 2, 7, 3, 4, 3, 6, 2, 8, 2, 3, 4, 4, 4, 9, 2, 4, 4, 7, 2, 8, 2, 6, 6, 4, 2, 7, 3, 6, 4, 6, 2, 6, 4, 7, 4, 4, 2, 12, 2, 4, 6, 3, 4, 8, 2, 6, 4, 8, 2, 11, 2, 4, 6, 6, 4, 8, 2, 8, 3, 4, 2, 12, 4, 4, 4, 7, 2, 12, 4, 6, 4, 4, 4, 7, 2, 6, 6, 9, 2, 8
Offset: 1

Views

Author

Reinhard Zumkeller, Jul 19 2002

Keywords

Comments

a(n) >= A073184(n).

Examples

			The cubefree kernel of 56 = 7*2^3 is 28 = 7*2^2 and the divisors <= 28 of 56 are {1, 2, 4, 7, 8, 14, 28}, therefore a(56) = 7.
		

Crossrefs

Programs

  • Mathematica
    Table[Function[k, DivisorSum[n, 1 &, # <= k &]]@ Apply[Times, FactorInteger[n] /. {p_, e_} /; p > 0 :> p^Min[e, 2]], {n, 102}] (* Michael De Vlieger, Jul 18 2017 *)
  • PARI
    a007948(n) = my(f=factor(n)); for (i=1, #f~, f[i, 2] = min(f[i, 2], 2)); factorback(f);
    a(n) = my(cfk = a007948(n)); sumdiv(n, d, d<=cfk); \\ Michel Marcus, May 14 2015
Showing 1-10 of 23 results. Next