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-5 of 5 results.

A053150 Cube root of largest cube dividing n.

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, 3, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1
Offset: 1

Views

Author

Henry Bottomley, Feb 28 2000

Keywords

Comments

This can be thought as a "lower 3rd root" of a positive integer. Upper k-th roots were studied by Broughan (2002, 2003, 2006). The sequence of "upper 3rd root" of positive integers is given by A019555. - Petros Hadjicostas, Sep 15 2019

Crossrefs

Cf. A000188 (inner square root), A019554 (outer square root), A019555 (outer third root), A053164 (inner 4th root), A053166 (outer 4th root), A015052 (outer 5th root), A015053 (outer 6th root).

Programs

  • Mathematica
    f[list_] := list[[1]]^Quotient[list[[2]], 3]; Table[Apply[Times, Map[f,FactorInteger[n]]], {n, 1, 81}] (* Geoffrey Critzer, Jan 21 2015 *)
    Table[SelectFirst[Reverse@ Divisors@ n, IntegerQ[#^(1/3)] &]^(1/3), {n, 105}] (* Michael De Vlieger, Jul 28 2017 *)
    f[p_, e_] := p^Floor[e/3]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 18 2020 *)
  • PARI
    A053150(n) = { my(f = factor(n), m = 1); for (k=1, #f~, m *= (f[k, 1]^(f[k, 2]\3)); ); m; } \\ Antti Karttunen, Jul 28 2017
    
  • PARI
    a(n) = my(f = factor(n)); for (k=1, #f~, f[k,2] = f[k,2]\3); factorback(f); \\ Michel Marcus, Jul 28 2017
    
  • Python
    from math import prod
    from sympy import factorint
    def A053150(n): return prod(p**(q//3) for p, q in factorint(n).items()) # Chai Wah Wu, Aug 18 2021

Formula

Multiplicative with a(p^e) = p^[e/3]. - Mitch Harris, Apr 19 2005
a(n) = A008834(n)^(1/3) = sqrt(A000189(n)/A000188(A050985(n))).
Dirichlet g.f.: zeta(3s-1)*zeta(s)/zeta(3s). - R. J. Mathar, Apr 09 2011
Sum_{k=1..n} a(k) ~ Pi^2 * n / (6*zeta(3)) + 3*zeta(2/3) * n^(2/3) / Pi^2. - Vaclav Kotesovec, Jan 31 2019
a(n) = Sum_{d^3|n} phi(d). - Ridouane Oudra, Dec 30 2020
G.f.: Sum_{k>=1} phi(k) * x^(k^3) / (1 - x^(k^3)). - Ilya Gutkovskiy, Aug 20 2021

Extensions

More terms from Antti Karttunen, Jul 28 2017

A053166 Smallest positive integer for which n divides a(n)^4.

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 2, 3, 10, 11, 6, 13, 14, 15, 2, 17, 6, 19, 10, 21, 22, 23, 6, 5, 26, 3, 14, 29, 30, 31, 4, 33, 34, 35, 6, 37, 38, 39, 10, 41, 42, 43, 22, 15, 46, 47, 6, 7, 10, 51, 26, 53, 6, 55, 14, 57, 58, 59, 30, 61, 62, 21, 4, 65, 66, 67, 34, 69, 70, 71, 6, 73, 74, 15, 38
Offset: 1

Views

Author

Henry Bottomley, Feb 29 2000

Keywords

Comments

According to Broughan (2002, 2003, 2006), a(n) is the "upper 4th root of n". The "lower 4th root of n" is sequence A053164. - Petros Hadjicostas, Sep 15 2019

Crossrefs

Cf. A000188 (inner square root), A019554 (outer square root), A053150 (inner 3rd root), A019555 (outer 3rd root), A053164 (inner 4th root), A015052 (outer 5th root), A015053 (outer 6th root).

Programs

  • Mathematica
    f[p_, e_] := p^Ceiling[e/4]; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Sep 08 2020 *)
  • PARI
    a(n) = my(f=factor(n)); for (i=1, #f~, f[i,2] = ceil(f[i,2]/4)); factorback(f); \\ Michel Marcus, Jun 09 2014

Formula

a(n) = n/A000190(n) = A019554(n)/(A008835(A019554(n)^2))^(1/4).
If n is 5th-power-free (i.e., not 32, 64, 128, 243, ...) then a(n) = A007947(n).
Multiplicative with a(p^e) = p^(ceiling(e/4)). - Christian G. Bower, May 16 2005
Sum_{k=1..n} a(k) ~ c * n^2, where c = (zeta(7)/2) * Product_{p prime} (1 - 1/p^2 + 1/p^3 - 1/p^4 + 1/p^5 - 1/p^6) = 0.3528057925... . - Amiram Eldar, Oct 27 2022

A015049 Let m = A013929(n); then a(n) = smallest k such that m divides k^2.

Original entry on oeis.org

2, 4, 3, 6, 4, 6, 10, 12, 5, 9, 14, 8, 6, 20, 22, 15, 12, 7, 10, 26, 18, 28, 30, 21, 8, 34, 12, 15, 38, 20, 9, 42, 44, 30, 46, 24, 14, 33, 10, 52, 18, 28, 58, 39, 60, 11, 62, 25, 42, 16, 66, 45, 68, 70, 12, 21, 74, 30, 76, 51, 78, 40, 18, 82, 84, 13, 57, 86
Offset: 1

Views

Author

R. Muller

Keywords

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := p^Ceiling[e/2]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; s /@ Select[Range[200], !SquareFreeQ[#] &] (* Amiram Eldar, Feb 10 2021 *)
  • PARI
    lista(kmax) = {my(f); for(k = 2, kmax, f = factor(k); if(!issquarefree(f), print1(k/core(f, 1)[2], ", ")));} \\ Amiram Eldar, Jan 06 2024

Formula

a(n) = A019554(A013929(n)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(2)*(zeta(3)-1)/(zeta(2)-1)^2 = 0.799082... . - Amiram Eldar, Jan 06 2024

Extensions

Description corrected by Diego Torres (torresvillarroel(AT)hotmail.com), Jun 23 2002
Offset corrected by Amiram Eldar, Feb 10 2021

A015050 Let m = A013929(n); then a(n) = smallest k such that m divides k^3.

Original entry on oeis.org

2, 2, 3, 6, 4, 6, 10, 6, 5, 3, 14, 4, 6, 10, 22, 15, 12, 7, 10, 26, 6, 14, 30, 21, 4, 34, 6, 15, 38, 20, 9, 42, 22, 30, 46, 12, 14, 33, 10, 26, 6, 28, 58, 39, 30, 11, 62, 5, 42, 8, 66, 15, 34, 70, 12, 21, 74, 30, 38, 51, 78, 20, 18, 82, 42, 13, 57, 86
Offset: 1

Views

Author

R. Muller

Keywords

Crossrefs

Programs

  • Maple
    isA013929 := proc(n)
        not numtheory[issqrfree](n) ;
    end proc:
    A013929 := proc(n)
        option remember;
        local a;
        if n = 1 then
            4;
        else
            for a from procname(n-1)+1 do
                if isA013929(a) then
                    return a;
                end if;
            end do:
        end if;
    end proc:
    A015050 := proc(n)
        local m ;
        m := A013929(n) ;
        for k from 1 do
            if modp(k^3,m) = 0 then
                return k;
            end if;
        end do:
    end proc:
  • Mathematica
    f[p_, e_] := p^Ceiling[e/3]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; s /@ Select[Range[200], !SquareFreeQ[#] &] (* Amiram Eldar, Feb 09 2021 *)
  • PARI
    lista(kmax) = {my(f); for(k = 2, kmax, f = factor(k); if(!issquarefree(f), print1(prod(i = 1, #f~, f[i,1]^ceil(f[i,2]/3)), ", ")));} \\ Amiram Eldar, Jan 06 2024

Formula

a(n) = A019555(A013929(n)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(2) * (zeta(2) * zeta(5) * Product_{p prime} (1-1/p^2+1/p^3-1/p^4) - 1)/(zeta(2)-1)^2 = 0.6611256641303... . - Amiram Eldar, Jan 06 2024

Extensions

Description corrected by Diego Torres (torresvillarroel(AT)hotmail.com), Jun 23 2002

A015047 Numbers k such that S(k) = S(k-1) + S(k-2) where S is the Kempner function A002034.

Original entry on oeis.org

3, 11, 121, 4902, 26245, 32112, 64010, 368140, 415664, 2091206, 2519648, 4573053, 7783364, 12610170, 16739482, 20526651, 30744490, 31137773, 44893057, 79269727, 99681870, 112803824, 136193976, 213879072, 321022289, 445810543, 559199345, 617868162, 624292242
Offset: 1

Views

Author

R. Muller (Research37(AT)aol.com)

Keywords

References

  • Charles Ashbacher and Mike Mudge, Personal Computer World, London, #10, 1995, p. 302.

Crossrefs

Extensions

Missing terms inserted by Amiram Eldar, Jul 07 2021
Showing 1-5 of 5 results.