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

A072774 Powers of squarefree numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97
Offset: 1

Views

Author

Reinhard Zumkeller, Jul 10 2002

Keywords

Comments

Essentially the same as A062770. - R. J. Mathar, Sep 25 2008
Numbers m such that in canonical prime factorization all prime exponents are identical: A124010(m,k) = A124010(m,1) for k = 2..A000005(m). - Reinhard Zumkeller, Apr 06 2014
Heinz numbers of uniform partitions. An integer partition is uniform if all parts appear with the same multiplicity. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 16 2018

Links

Crossrefs

Complement of A059404.
Cf. A072775, A072776, A072777 (subsequence), A005117, A072778, A124010, A329332 (tabular arrangement), A384667 (characteristic function).
A subsequence of A242414.
Cf. A000005, A000009, A000837, A007916, A013661, A047966, A052409, A052410, A072774, A078374, A289023, A289509, A300486, A302491, A302796, A302979.

Programs

  • Haskell
    import Data.Map (empty, findMin, deleteMin, insert)
import qualified Data.Map.Lazy as Map (null)
a072774 n = a072774_list !! (n-1)
(a072774_list, a072775_list, a072776_list) = unzip3 $
   (1, 1, 1) : f (tail a005117_list) empty where
   f vs'@(v:vs) m
    | Map.null m || xx > v = (v, v, 1) :
                             f vs (insert (v^2) (v, 2) m)
    | otherwise = (xx, bx, ex) :
                  f vs' (insert (bx*xx) (bx, ex+1) $ deleteMin m)
    where (xx, (bx, ex)) = findMin m
-- Reinhard Zumkeller, Apr 06 2014
  • Maple
    isA := n -> n=1 or is(1 = nops({seq(p[2], p in ifactors(n)[2])})):
select(isA, [seq(1..97)]);  # Peter Luschny, Jun 10 2025
  • Mathematica
    Select[Range[100], Length[Union[FactorInteger[#][[All, 2]]]] == 1 &] (* Geoffrey Critzer, Mar 30 2015 *)
  • PARI
    is(n)=ispower(n,,&n); issquarefree(n) \\ Charles R Greathouse IV, Oct 16 2015
  • Python
    from math import isqrt
from sympy import mobius, integer_nthroot
def A072774(n):
    def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))-1
    def f(x): return n-2+x-sum(g(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length()))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 19 2024

Formula

a(n) = A072775(n)^A072776(n).
Sum_{n>=1} 1/a(n)^s = 1 + Sum_{k>=1} (zeta(k*s)/zeta(2*k*s)-1) for s > 1. - Amiram Eldar, Mar 20 2025
a(n)/n ~ Pi^2/6 (A013661). - Friedjof Tellkamp, Jun 09 2025

A052410 Write n = m^k with m, k integers, k >= 1, then a(n) is the smallest possible choice for m.

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein

Keywords

Comments

Value of m in m^p = n, where p is the largest possible power (see A052409).
For n > 1, n is a perfect power iff a(n) <> n. - Reinhard Zumkeller, Oct 13 2002
a(n)^A052409(n) = n. - Reinhard Zumkeller, Apr 06 2014
Every integer root of n is a power of a(n). All entries (except 1) belong to A007916. - Gus Wiseman, Sep 11 2017

Links

Crossrefs

Cf. A001597, A025478, A007916, A027748, A052409, A072775, A124010, A175781, A278028, A288636, A289023.

Programs

  • Haskell
    a052410 n = product $ zipWith (^)
                      (a027748_row n) (map (`div` (foldl1 gcd es)) es)
            where es = a124010_row n
-- Reinhard Zumkeller, Jul 15 2012
  • Maple
    a:= n-> (l-> (t-> mul(i[1]^(i[2]/t), i=l))(
         igcd(seq(i[2], i=l))))(ifactors(n)[2]):
seq(a(n), n=1..74);  # Alois P. Heinz, Jul 22 2024
  • Mathematica
    Table[If[n==1, 1, n^(1/(GCD@@(Last/@FactorInteger[n])))], {n, 100}]
  • PARI
    a(n) = if (ispower(n,,&r), r, n); \\ Michel Marcus, Jul 19 2017
  • Python
    def upto(n):
    list = [1] + [0] * (n - 1)
    for i in range(2, n + 1):
        if not list[i - 1]:
            j = i
            while j <= n:
                list[j - 1] = i
                j *= i
    return list
# M. Eren Kesim, Jun 03 2021
  • Python
    from math import gcd
from sympy import integer_nthroot, factorint
def A052410(n): return integer_nthroot(n,gcd(*factorint(n).values()))[0] if n>1 else 1 # Chai Wah Wu, Mar 02 2024

Formula

a(A001597(k)) = A025478(k).
a(n) = A007916(A278028(n,1)). - Gus Wiseman, Sep 11 2017

Extensions

Definition edited (in a complementary form to A052409) by Daniel Forgues, Mar 14 2009
Corrected by Charles R Greathouse IV, Sep 02 2009
Definition edited by N. J. A. Sloane, Sep 03 2010

A072776 Exponents of powers of squarefree numbers.

Original entry on oeis.org

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

Views

Author

Reinhard Zumkeller, Jul 10 2002

Keywords

Comments

A072774(n) = A072775(n)^a(n);
A072774(n) is squarefree iff a(n)=1.

Links

Crossrefs

Cf. A052409.

Programs

  • Haskell
    a072776 n = a072776_list !! (n-1)  -- a072776_list defined in A072774.
-- Reinhard Zumkeller, Apr 06 2014
  • Python
    from math import isqrt
from sympy import mobius, integer_nthroot, perfect_power
def A072776(n):
    def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))-1
    def f(x): return n-2+x-sum(g(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length()))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return 1 if not (p:=perfect_power(kmax)) else p[1] # Chai Wah Wu, Aug 19 2024
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.