A019554 -id:A019554 - cp's OEIS frontend

A007947 Largest squarefree number dividing n: the squarefree kernel of n, rad(n), radical of n.

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, 2, 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, 2, 65, 66, 67, 34, 69, 70, 71, 6, 73, 74, 15, 38, 77, 78

Offset: 1

R. Muller, Mar 15 1996

Multiplicative with a(p^e) = p.
Product of the distinct prime factors of n.
a(k)=k for k=squarefree numbers A005117. - Lekraj Beedassy, Sep 05 2006
A note on square roots of numbers: we can write sqrt(n) = b*sqrt(c) where c is squarefree. Then b = A000188(n) is the "inner square root" of n, c = A007913(n), b*c = A019554(n) = "outer square root" of n, and a(n) = lcm(a(b),c). Unless n is biquadrateful (A046101), a(n) = lcm(b,c). [Edited by Jeppe Stig Nielsen, Oct 10 2021, and Andrey Zabolotskiy, Feb 12 2025]
a(n) = A128651(A129132(n-1) + 2) for n > 1. - Reinhard Zumkeller, Mar 30 2007
Also the least common multiple of the prime factors of n. - Peter Luschny, Mar 22 2011
The Mobius transform of the sequence generates the sequence of absolute values of A097945. - R. J. Mathar, Apr 04 2011
Appears to be the period length of k^n mod n. For example, n^12 mod 12 has period 6, repeating 1,4,9,4,1,0, so a(12)= 6. - Gary Detlefs, Apr 14 2013
a(n) differs from A014963(n) when n is a term of A024619. - Eric Desbiaux, Mar 24 2014
a(n) is also the smallest base (also termed radix) for which the representation of 1/n is of finite length.  For example a(12) = 6 and 1/12 in base 6 is 0.03, which is of finite length. - Lee A. Newberg, Jul 27 2016
a(n) is also the divisor k of n such that d(k) = 2^omega(n). a(n) is also the smallest divisor u of n such that n divides u^n. - Juri-Stepan Gerasimov, Apr 06 2017

			G.f. = x + 2*x^2 + 3*x^3 + 2*x^4 + 5*x^5 + 6*x^6 + 7*x^7 + 2*x^8 + 3*x^9 + ... - _Michael Somos_, Jul 15 2018

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from T. D. Noe)
Masum Billal, Divisible Sequence and its Characteristic Sequence, arXiv:1501.00609 [math.NT], 2015, theorem 11 page 5.
Henry Bottomley, Some Smarandache-type multiplicative sequences
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
Jarosław Grytczuk, Thue type problems for graphs, points and numbers, Discrete Math., 308 (2008), 4419-4429.
Neville Holmes, Integer Sequences [Broken link]
Serge Lang, Old and New Conjectured Diophantine Inequalities, Bull. Amer. Math. Soc., 23 (1990), 37-75. see p. 39.
Wolfdieter Lang, Cantor's List of Real Algebraic Numbers of Heights 1 to 7, arXiv:2307.10645 [math.NT], 2023.
D. H. Lehmer, Euler constants for arithmetical progressions, Collection of articles in memory of Juriĭ Vladimirovič Linnik. Acta Arith. 27 (1975), 125--142. MR0369233 (51 #5468). See N_k on page 131.
Ivar Peterson, The Amazing ABC Conjecture
Paul Tarau, Emulating Primality with Multiset Representations of Natural Numbers, in Theoretical Aspects of Computing, ICTAC 2011, Lecture Notes in Computer Science, 2011, Volume 6916/2011, 218-238
Paul Tarau, Towards a generic view of primality through multiset decompositions of natural numbers, Theoretical Computer Science, Volume 537, 5 June 2014, Pages 105-124.
Wikipedia, Radical of an integer.

See A007913, A062953, A000188, A019554, A003557, A066503, A087207 for other properties related to square and squarefree divisors of n.
More general factorization-related properties, specific to n: A020639, A028234, A020500, A010051, A284318, A000005, A001221, A005361, A034444, A014963, A128651, A267116.
Range of values is A005117.
Cf. A048803, A019565, A024619, A053462, A060735, A073355 (partial sums), A079277, A097945, A129132, A225546.
Bisections: A099984, A099985.
Sequences about numbers that have the same squarefree kernel: A065642, array A284311 (A284457).
A003961, A059896 are used to express relationship between terms of this sequence.
Cf. A001615, A008683, A076479, A078615.

a007947 = product . a027748_row  -- Reinhard Zumkeller, Feb 27 2012

[ &*PrimeDivisors(n): n in [1..100] ]; // Klaus Brockhaus, Dec 04 2008

with(numtheory); A007947 := proc(n) local i,t1,t2; t1 := ifactors(n)[2]; t2 := mul(t1[i][1],i=1..nops(t1)); end;
A007947 := n -> ilcm(op(numtheory[factorset](n))):
seq(A007947(i),i=1..69); # Peter Luschny, Mar 22 2011
A:= n -> convert(numtheory:-factorset(n),`*`):
seq(A(n),n=1..100); # Robert Israel, Aug 10 2014
seq(NumberTheory:-Radical(n), n = 1..78); # Peter Luschny, Jul 20 2021

rad[n_] := Times @@ (First@# & /@ FactorInteger@ n); Array[rad, 78] (* Robert G. Wilson v, Aug 29 2012 *)
Table[Last[Select[Divisors[n],SquareFreeQ]],{n,100}] (* Harvey P. Dale, Jul 14 2014 *)
a[ n_] := If[ n < 1, 0, Sum[ EulerPhi[d] Abs @ MoebiusMu[d], {d, Divisors[ n]}]]; (* Michael Somos, Jul 15 2018 *)
Table[Product[p, {p, Select[Divisors[n], PrimeQ]}], {n, 1, 100}] (* Vaclav Kotesovec, May 20 2020 *)

a(n) = factorback(factorint(n)[,1]); \\ Andrew Lelechenko, May 09 2014

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

from sympy import primefactors, prod
def a(n): return 1 if n < 2 else prod(primefactors(n))
[a(n) for n in range(1, 51)]  # Indranil Ghosh, Apr 16 2017

def A007947(n): return mul(p for p in prime_divisors(n))
[A007947(n) for n in (1..60)] # Peter Luschny, Mar 07 2017

(define (A007947 n) (if (= 1 n) n (* (A020639 n) (A007947 (A028234 n))))) ;; ;; Needs also code from A020639 and A028234. - Antti Karttunen, Jun 18 2017

If n = Product_j (p_j^k_j) where p_j are distinct primes, then a(n) = Product_j (p_j).
a(n) = Product_{k=1..A001221(n)} A027748(n,k). - Reinhard Zumkeller, Aug 27 2011
Dirichlet g.f.: zeta(s)*Product_{primes p} (1+p^(1-s)-p^(-s)). - R. J. Mathar, Jan 21 2012
a(n) = Sum_{d|n} phi(d) * mu(d)^2 = Sum_{d|n} |A097945(d)|. - Enrique Pérez Herrero, Apr 23 2012
a(n) = Product_{d|n} d^moebius(n/d) (see Billal link). - Michel Marcus, Jan 06 2015
a(n) = n/( Sum_{k=1..n} (floor(k^n/n)-floor((k^n - 1)/n)) ) = e^(Sum_{k=2..n} (floor(n/k) - floor((n-1)/k))*A010051(k)*M(k)) where M(n) is the Mangoldt function. - Anthony Browne, Jun 17 2016
a(n) = n/A003557(n). - Juri-Stepan Gerasimov, Apr 07 2017
G.f.: Sum_{k>=1} phi(k)*mu(k)^2*x^k/(1 - x^k). - Ilya Gutkovskiy, Apr 11 2017
From Antti Karttunen, Jun 18 2017: (Start)
a(1) = 1; for n > 1, a(n) = A020639(n) * a(A028234(n)).
a(n) = A019565(A087207(n)). (End)
Dirichlet g.f.: zeta(s-1) * zeta(s) * Product_{primes p} (1 + p^(1-2*s) - p^(2-2*s) - p^(-s)). - Vaclav Kotesovec, Dec 18 2019
From Peter Munn, Jan 01 2020: (Start)
a(A059896(n,k)) = A059896(a(n), a(k)).
a(A003961(n)) = A003961(a(n)).
a(n^2) = a(n).
a(A225546(n)) = A019565(A267116(n)). (End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = A065463/2. - Vaclav Kotesovec, Jun 24 2020
From Richard L. Ollerton, May 07 2021: (Start)
a(n) = Sum_{k=1..n} mu(n/gcd(n,k))^2.
a(n) = Sum_{k=1..n} mu(gcd(n,k))^2*phi(gcd(n,k))/phi(n/gcd(n,k)).
For n>1, Sum_{k=1..n} a(gcd(n,k))*mu(a(gcd(n,k)))*phi(gcd(n,k))/gcd(n,k) = 0.
For n>1, Sum_{k=1..n} a(n/gcd(n,k))*mu(a(n/gcd(n,k)))*phi(gcd(n,k))*gcd(n,k) = 0. (End)
a(n) = (-1)^omega(n) * Sum_{d|n} mu(d)*psi(d), where omega = A001221 and psi = A001615. - Ridouane Oudra, Aug 01 2025

More terms from several people including David W. Wilson

Definition expanded by Jonathan Sondow, Apr 26 2013

A007913 Squarefree part of n: a(n) is the smallest positive number m such that n/m is a square.

Original entry on oeis.org

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

Offset: 1

R. Muller, Mar 15 1996

Also called core(n). [Not to be confused with the squarefree kernel of n, A007947.]
Sequence read mod 4 gives A065882. - Philippe Deléham, Mar 28 2004
This is an arithmetic function and is undefined if n <= 0.
A note on square roots of numbers: we can write sqrt(n) = b*sqrt(c) where c is squarefree. Then b = A000188(n) is the "inner square root" of n, c = A007913(n), lcm(A007947(b),c) = A007947(n) = "squarefree kernel" of n and bc = A019554(n) = "outer square root" of n. [Corrected by M. F. Hasler, Mar 01 2018]
If n > 1, the quantity f(n) = log(n/core(n))/log(n) satisfies 0 <= f(n) <= 1; f(n) = 0 when n is squarefree and f(n) = 1 when n is a perfect square. One can define n as being "epsilon-almost squarefree" if f(n) < epsilon. - Kurt Foster (drsardonicus(AT)earthlink.net), Jun 28 2008
a(n) is the smallest natural number m such that product of geometric mean of the divisors of n and geometric mean of the divisors of m are integers. Geometric mean of the divisors of number n is real number b(n) = Sqrt(n). a(n) = 1 for infinitely many n. a(n) = 1 for numbers from A000290: a(A000290(n)) = 1. For n = 8; b(8) = sqrt(8), a(n) = 2 because b(2) = sqrt(2); sqrt(8) * sqrt(2) = 4 (integer). - Jaroslav Krizek, Apr 26 2010
Dirichlet convolution of A010052 with the sequence of absolute values of A055615. - R. J. Mathar, Feb 11 2011
Booker, Hiary, & Keating outline a method for bounding (on the GRH) a(n) for large n using L-functions. - Charles R Greathouse IV, Feb 01 2013
According to the formula a(n) = n/A000188(n)^2, the scatterplot exhibits the straight lines y=x, y=x/4, y=x/9, ..., i.e., y=x/k^2 for all k=1,2,3,... - M. F. Hasler, May 08 2014
The Dirichlet inverse of this sequence is A008836(n) * A063659(n). - Álvar Ibeas, Mar 19 2015
a(n) = 1 if n is a square, a(n) = n if n is a product of distinct primes. - Zak Seidov, Jan 30 2016
All solutions of the Diophantine equation n*x=y^2 or, equivalently, G(n,x)=y, with G being the geometric mean, are of the form x=k^2*a(n), y=k*sqrt(n*a(n)), where k is a positive integer. - Stanislav Sykora, Feb 03 2016
If f is a multiplicative function then Sum_{d divides n} f(a(d)) is also multiplicative. For example, A010052(n) = Sum_{d divides n} mu(a(d)) and A046951(n) = Sum_{d divides n} mu(a(d)^2). - Peter Bala, Jan 24 2024

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 1000 terms from T. D. Noe)
Krassimir T. Atanassov, On Some of Smarandache's Problems.
Krassimir T. Atanassov, On the 22nd, 23rd, and the 24th Smarandache Problems, Notes on Number Theory and Discrete Mathematics, Sophia, Bulgaria, Vol. 5 (1999), No. 2, 80-82.
Andrew Booker, Ghaith Hiary, and Jon Keating, Detecting squarefree numbers, CNTA XII (2012).
Henry Bottomley, Some Smarandache-type multiplicative sequences.
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
Vlad Copil and Laurenţiu Panaitopol, Properties of a sequence generated by positive integers, Bulletin mathématique de la Société des Sciences Mathématiques de Roumanie, Nouvelle Série, Vol. 50 (98), No. 2 (2007), pp. 131-137; alternative link.
Florentin Smarandache, Only Problems, Not Solutions!, Xiquan Publ., Phoenix-Chicago, 1993.
Eric Weisstein's World of Mathematics, Squarefree Part.

See A000188, A007947, A008833, A019554, A117811 for related information, specific to n.
See A027746, A027748, A124010 for factorization data for n.
Analogous sequences: A050985, A053165, A055231.
Cf. A002734, A005117 (range of values), A059897, A069891 (partial sums), A090699, A350389.
Related to A006519 via A225546.

a007913 n = product $
            zipWith (^) (a027748_row n) (map (`mod` 2) $ a124010_row n)
-- Reinhard Zumkeller, Jul 06 2012

[ Squarefree(n) : n in [1..256] ]; // N. J. A. Sloane, Dec 23 2006

A007913 := proc(n) local f,a,d; f := ifactors(n)[2] ; a := 1 ; for d in f do if type(op(2,d),'odd') then a := a*op(1,d) ; end if; end do: a; end proc: # R. J. Mathar, Mar 18 2011
# second Maple program:
a:= n-> mul(i[1]^irem(i[2], 2), i=ifactors(n)[2]):
seq(a(n), n=1..100);  # Alois P. Heinz, Jul 20 2015
seq(n / expand(numtheory:-nthpow(n, 2)), n=1..77);  # Peter Luschny, Jul 12 2022

data = Table[Sqrt[n], {n, 1, 100}]; sp = data /. Sqrt[] -> 1; sfp = data/sp /. Sqrt[x] -> x (* Artur Jasinski, Nov 03 2008 *)
Table[Times@@Power@@@({#[[1]],Mod[ #[[2]],2]}&/@FactorInteger[n]),{n,100}] (* Zak Seidov, Apr 08 2009 *)
Table[{p, e} = Transpose[FactorInteger[n]]; Times @@ (p^Mod[e, 2]), {n, 100}] (* T. D. Noe, May 20 2013 *)
Sqrt[#] /. (c_:1)*a_^(b_:0) -> (c*a^b)^2& /@ Range@100 (* Bill Gosper, Jul 18 2015 *)

PARI
```
a(n)=core(n)
```

from sympy import factorint, prod
def A007913(n):
    return prod(p for p, e in factorint(n).items() if e % 2)
# Chai Wah Wu, Feb 03 2015

[squarefree_part(n) for n in (1..77)] # Peter Luschny, Feb 04 2015

Multiplicative with a(p^k) = p^(k mod 2). - David W. Wilson, Aug 01 2001
a(n) modulo 2 = A035263(n); a(A036554(n)) is even; a(A003159(n)) is odd. - Philippe Deléham, Mar 28 2004
Dirichlet g.f.: zeta(2s)*zeta(s-1)/zeta(2s-2). - R. J. Mathar, Feb 11 2011
a(n) = n/( Sum_{k=1..n} floor(k^2/n)-floor((k^2 -1)/n) )^2. - Anthony Browne, Jun 06 2016
a(n) = rad(n)/a(n/rad(n)), where rad = A007947. This recurrence relation together with a(1) = 1 generate the sequence. - Velin Yanev, Sep 19 2017
From Peter Munn, Nov 18 2019: (Start)
a(k*m) = A059897(a(k), a(m)).
a(n) = n / A008833(n).
(End)
a(A225546(n)) = A225546(A006519(n)). - Peter Munn, Jan 04 2020
From Amiram Eldar, Mar 14 2021: (Start)
Theorems proven by Copil and Panaitopol (2007):
Lim sup_{n->oo} a(n+1)-a(n) = oo.
Lim inf_{n->oo} a(n+1)-a(n) = -oo.
Sum_{k=1..n} 1/a(k) ~ c*sqrt(n) + O(log(n)), where c = zeta(3/2)/zeta(3) (A090699). (End)
a(n) = A019554(n)^2/n. - Jianing Song, May 08 2022
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/30 = 0.328986... . - Amiram Eldar, Oct 25 2022
a(n) = A007947(A350389(n)). - Amiram Eldar, Jan 20 2024

More terms from Michael Somos, Nov 24 2001

Definition reformulated by Daniel Forgues, Mar 24 2009

A000188 (1) Number of solutions to x^2 == 0 (mod n). (2) Also square root of largest square dividing n. (3) Also max_{ d divides n } gcd(d, n/d).

Original entry on oeis.org

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

Offset: 1

N. J. A. Sloane

Shadow transform of the squares A000290. - Vladeta Jovovic, Aug 02 2002
Labos Elemer and Henry Bottomley independently proved that (2) and (3) define the same sequence. Bottomley also showed that (1) and (2) define the same sequence.
Proof that (2) = (3): Let max{gcd(d, n/d)} = K, then d = Kx, n/d = Ky so n = KKxy where xy is the squarefree part of n, otherwise K is not maximal. Observe also that g = gcd(K, xy) is not necessarily 1. Thus K is also the "maximal square-root factor" of n. - Labos Elemer, Jul 2000
We can write sqrt(n) = b*sqrt(c) where c is squarefree. Then b = A000188(n) is the "inner square root" of n, c = A007913(n) and b*c = A019554(n) = "outer square root" of n.

			a(8) = 2 because the largest square dividing 8 is 4, the square root of which is 2.
a(9) = 3 because 9 is a perfect square and its square root is 3.
a(10) = 1 because 10 is squarefree.

T. D. Noe, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Kevin A. Broughan, Restricted divisor sums, preprint.
Kevin A. Broughan, Restricted divisor sums, Acta Arithmetica, 101(2) (2002), 105-114.
Kevin A. Broughan, Relationship between the integer conductor and k-th root functions, Int. J. Pure Appl. Math. 5(3) (2003), 253-275.
Kevin A. Broughan, Relaxations of the ABC Conjecture using integer k'th roots, New Zealand J. Math. 35(2) (2006), 121-136.
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
Steven R. Finch and Pascal Sebah, Squares and Cubes Modulo n, arXiv:math/0604465 [math.NT], 2006-2016.
Vaclav Kotesovec, Graph - the asymptotic ratio (100000 terms)
Gerry Myerson, Trifectas in Geometric Progression, Australian Mathematical Society Gazette 35(3) (2008), 189-194
Andrew Reiter, On (mod n) spirals (2014), and posting to Number Theory Mailing List, Mar 23 2014.
N. J. A. Sloane, Transforms.
Florentin Smarandache, Collected Papers, Vol. II, Tempus Publ. Hse, Bucharest, 1996.
László Tóth, Counting solutions of quadratic congruences in several variables revisited, arXiv preprint arXiv:1404.4214 [math.NT], 2014.
László Tóth, Counting Solutions of Quadratic Congruences in Several Variables Revisited, J. Int. Seq. 17 (2014), # 14.11.6.

Cf. A019554 (outer square root), A053150 (inner 3rd root), A019555 (outer 3rd root), A053164 (inner 4th root), A053166 (outer 4th root), A015052 (outer 5th root), A015053 (outer 6th root).
Cf. A000189, A000190, A007947, A008833, A027748, A046951, A055210, A007913, A117811, A124010. For partial sums see A120486.
Cf. A240976 (Dgf at s=2).

a000188 n = product $ zipWith (^)
                      (a027748_row n) $ map (`div` 2) (a124010_row n)
-- Reinhard Zumkeller, Apr 22 2012

with(numtheory):A000188 := proc(n) local i: RETURN(op(mul(i,i=map(x->x[1]^floor(x[2]/2),ifactors(n)[2])))); end;

Array[Function[n, Count[Array[PowerMod[#, 2, n ] &, n, 0 ], 0 ] ], 100]
(* Second program: *)
nMax = 90; sList = Range[Floor[Sqrt[nMax]]]^2; Sqrt[#] &/@ Table[ Last[ Select[ sList, Divisible[n, #] &]], {n, nMax}] (* Harvey P. Dale, May 11 2011 *)
a[n_] := With[{d = Divisors[n]}, Max[GCD[d, Reverse[d]]]] (* Mamuka Jibladze, Feb 15 2015 *)
f[p_, e_] := p^Floor[e/2]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 18 2020 *)

PARI
```
a(n)=if(n<1,0,sum(i=1,n,i*i%n==0))
```

a(n)=sqrtint(n/core(n)) \\ Zak Seidov, Apr 07 2009

a(n)=core(n, 1)[2] \\ Michel Marcus, Feb 27 2013

from sympy.ntheory.factor_ import core
from sympy import integer_nthroot
def A000188(n): return integer_nthroot(n//core(n),2)[0] # Chai Wah Wu, Jun 14 2021

a(n) = n/A019554(n) = sqrt(A008833(n)).
a(n) = Sum_{d^2|n} phi(d), where phi is the Euler totient function A000010.
Multiplicative with a(p^e) = p^floor(e/2). - David W. Wilson, Aug 01 2001
Dirichlet series: Sum_{n >= 1} a(n)/n^s = zeta(2*s - 1)*zeta(s)/zeta(2*s), (Re(s) > 1).
Dirichlet convolution of A037213 and A008966. - R. J. Mathar, Feb 27 2011
Finch & Sebah show that the average order of a(n) is 3 log n/Pi^2. - Charles R Greathouse IV, Jan 03 2013
a(n) = sqrt(n/A007913(n)). - M. F. Hasler, May 08 2014
Sum_{n>=1} lambda(n)*a(n)*x^n/(1-x^n) = Sum_{n>=1} n*x^(n^2), where lambda() is the Liouville function A008836 (cf. A205801). - Mamuka Jibladze, Feb 15 2015
a(2*n) = a(n)*(A096268(n-1) + 1). - observed by Velin Yanev, Jul 14 2017, The formula says that a(2n) = 2*a(n) only when 2-adic valuation of n (A007814(n)) is odd, otherwise a(2n) = a(n). This follows easily from the definition (2). - Antti Karttunen, Nov 28 2017
Sum_{k=1..n} a(k) ~ 3*n*((log(n) + 3*gamma - 1)/Pi^2 - 12*zeta'(2)/Pi^4), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Dec 01 2020
Conjecture: a(n) = Sum_{k=1..n} A010052(n*k). - Velin Yanev, Jul 04 2021
G.f.: Sum_{k>=1} phi(k) * x^(k^2) / (1 - x^(k^2)). - Ilya Gutkovskiy, Aug 20 2021

Edited by M. F. Hasler, May 08 2014

A008833 Largest square dividing n.

Original entry on oeis.org

1, 1, 1, 4, 1, 1, 1, 4, 9, 1, 1, 4, 1, 1, 1, 16, 1, 9, 1, 4, 1, 1, 1, 4, 25, 1, 9, 4, 1, 1, 1, 16, 1, 1, 1, 36, 1, 1, 1, 4, 1, 1, 1, 4, 9, 1, 1, 16, 49, 25, 1, 4, 1, 9, 1, 4, 1, 1, 1, 4, 1, 1, 9, 64, 1, 1, 1, 4, 1, 1, 1, 36, 1, 1, 25, 4, 1, 1, 1, 16, 81, 1, 1, 4, 1, 1, 1, 4, 1, 9, 1, 4, 1, 1, 1, 16, 1

Offset: 1

N. J. A. Sloane

The Dirichlet generating function of the arithmetic function of the largest t-th power dividing n is zeta(s)*zeta(t*s-t)/zeta(s*t), here with t=2 and in A008834 and A008835 with t=3 and t=4, respectively. - R. J. Mathar, Feb 19 2011

Daniel Forgues, Table of n, a(n) for n = 1..100000
Henry Bottomley, Some Smarandache-type multiplicative sequences
R. J. Mathar, Survey of Dirichlet series of multiplicative arithmetic functions arXiv:1106.4038 [math.NT], 2011-2012, Remark 16.
Andrew Reiter, On (mod n) spirals, 2014, see also posting to Number Theory Mailing List, Mar 23 2014.
Eric Weisstein's World of Mathematics, Square part

Cf. A000188, A005563, A007913, A019554, A059897, A068310.

Cf. A008836, A358272, A001615.

a008833 n = head $ filter ((== 0) . (mod n)) $
   reverse $ takeWhile (<= n) $ tail a000290_list
-- Reinhard Zumkeller, Nov 13 2011

A008833 := proc(n)
    expand(numtheory:-nthpow(n,2)) ;
end proc:
seq(A008833(n), n=1..100) ;

a[n_] := First[ Select[ Reverse[ Divisors[n]], IntegerQ[Sqrt[#]]&, 1]]; Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Dec 12 2011 *)
f[p_, e_] := p^(2*Floor[e/2]); a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Jul 07 2020 *)

A008833(n)=n/core(n) \\ Michael B. Porter, Oct 17 2009

from sympy.ntheory.factor_ import core
def A008833(n): return n//core(n) # Chai Wah Wu, Dec 30 2021

a(n) = A000188(n)^2 = n/A007913(n). Cf. A019554.
Multiplicative with a(p^e) = p^(2[e/2]). - David W. Wilson, Aug 01 2001
Dirichlet g.f.: zeta(s)*zeta(2s-2)/zeta(2s). - R. J. Mathar, Oct 31 2011
a(n) = A005563(n-1) / A068310(n) for n > 1. - Reinhard Zumkeller, Nov 26 2011
Sum_{k=1..n} a(k) ~ Zeta(3/2) * n^(3/2) / (3*Zeta(3)). - Vaclav Kotesovec, Feb 01 2019
a(A059897(n,k)) = A059897(a(n), a(k)). - Peter Munn, Nov 30 2019
From Ridouane Oudra, May 11 2025: (Start)
a(n) = Sum_{d|n} lambda(d)*d*psi(n/d), where lambda = A008836 and psi = A001615.
a(n) = lambda(n) * Sum_{d|n} lambda(d)*d*phi(n/d).
a(n) = A008836(n) * A358272(n). (End)

A322483 The number of semi-unitary divisors of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 3, 2, 4, 2, 4, 2, 4, 4, 3, 2, 4, 2, 4, 4, 4, 2, 6, 2, 4, 3, 4, 2, 8, 2, 4, 4, 4, 4, 4, 2, 4, 4, 6, 2, 8, 2, 4, 4, 4, 2, 6, 2, 4, 4, 4, 2, 6, 4, 6, 4, 4, 2, 8, 2, 4, 4, 4, 4, 8, 2, 4, 4, 8, 2, 6, 2, 4, 4, 4, 4, 8, 2, 6, 3, 4, 2, 8, 4, 4, 4

Offset: 1

Amiram Eldar, Dec 11 2018

The notion of semi-unitary divisor was introduced by Chidambaraswamy in 1967.
A semi-unitary divisor of n is defined as the largest divisor d of n such that the largest divisor of d that is a unitary divisor of n/d is 1. In terms of the relation defined in A322482, d is the largest divisor of n such that T(d, n/d) = 1 (the largest divisor d that is semiprime to n/d).
The number of divisors of n that are exponentially odd numbers (A268335). - Amiram Eldar, Sep 08 2023

			The semi-unitary divisors of 8 are 1, 2, 8 (4 is not semi-unitary divisor since the largest divisor of 4 that is a unitary divisor of 8/4 = 2 is 2 > 1), and their number is 3, thus a(8) = 3.

J. Chidambaraswamy, Sum functions of unitary and semi-unitary divisors, J. Indian Math. Soc., Vol. 31 (1967), pp. 117-126.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Krishnaswami Alladi, On arithmetic functions and divisors of higher order, Journal of the Australian Mathematical Society, Vol. 23, No. 1 (1977), pp. 9-27.
Pentti Haukkanen, Basic properties of the bi-unitary convolution and the semi-unitary convolution, Indian J. Math, Vol. 40 (1998), pp. 305-315.
Vaclav Kotesovec, Graph - the asymptotic ratio (100000 terms).
D. Suryanarayana, The number of bi-unitary divisors of an integer, The theory of arithmetic functions. Springer, Berlin, Heidelberg, 1972, pp. 273-282.
D. Suryanarayana and V. Siva Rama Prasad, Sum functions of k-ary and semi-k-ary divisors, Journal of the Australian Mathematical Society, Vol. 15, No. 2 (1973), pp. 148-162.
D. Suryanarayana and R. Sita Rama Chandra Rao, The number of unitarily k-free divisors of an integer, Journal of the Australian Mathematical Society, Vol. 21, No. 1 (1976), pp. 19-35.
Laszlo Tóth, Sum functions of certain generalized divisors, Ann. Univ. Sci. Budap. Rolando Eötvös, Sect. Math., Vol. 41 (1998), pp. 165-180.

Cf. A000005, A005117, A019554, A034444, A268335, A286324, A322482.

Cf. A001620, A056671, A343443, A365488, A365491.

f[p_, e_] := Floor[(e+3)/2]; sud[n_] := If[n==1, 1, Times @@ (f @@@ FactorInteger[n])]; Array[sud, 100]

a(n) = {my(f = factor(n)); for (k=1, #f~, f[k,1] = (f[k,2]+3)\2; f[k,2] = 1;); factorback(f);} \\ Michel Marcus, Dec 14 2018

for(n=1, 100, print1(direuler(p=2, n, 1/(1-X) * 1/(1-X^2) * (1 + X - X^2))[n], ", ")) \\ Vaclav Kotesovec, Sep 06 2023

Multiplicative with a(p^e) = floor((e+3)/2).
a(n) <= A000005(n) with equality if and only if n is squarefree (A005117).
a(n) = Sum_{d|n} mu(d/gcd(d, n/d))^2. - Ilya Gutkovskiy, Feb 21 2020
a(n) = A000005(A019554(n)) (the number of divisors of the smallest number whose square is divisible by n). - Amiram Eldar, Sep 02 2023
From Vaclav Kotesovec, Sep 06 2023: (Start)
Dirichlet g.f.: zeta(s) * zeta(2*s) * Product_{p prime} (1 + 1/p^s - 1/p^(2*s)).
Dirichlet g.f.: zeta(s)^2 * zeta(2*s) * Product_{p prime} (1 - 2/p^(2*s) + 1/p^(3*s)).
Let f(s) = Product_{p prime} (1 - 2/p^(2*s) + 1/p^(3*s)).
Sum_{k=1..n} a(k) ~ Pi^2 * f(1) * n / 6 * (log(n) + 2*gamma - 1 + 12*zeta'(2)/Pi^2 + f'(1)/f(1)), where
f(1) = Product_{p prime} (1 - 2/p^2 + 1/p^3) = A065464 = 0.42824950567709444...,
f'(1) = f(1) * Sum_{p prime} (4*p-3) * log(p) / (p^3 - 2*p + 1) = 0.808661108949590913395... and gamma is the Euler-Mascheroni constant A001620. (End)

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

Henry Bottomley, Feb 28 2000

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

Antti Karttunen, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Kevin A. Broughan, Restricted divisor sums, Acta Arithmetica, 101(2) (2002), 105-114.
Kevin A. Broughan, Relationship between the integer conductor and k-th root functions, Int. J. Pure Appl. Math. 5(3) (2003), 253-275.
Kevin A. Broughan, Relaxations of the ABC Conjecture using integer k'th roots, New Zealand J. Math. 35(2) (2006), 121-136.
Vaclav Kotesovec, Graph - the asymptotic ratio.

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).

Cf. A000189, A000190, A008834, A008835, A015051, A061704.

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 *)

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

a(n) = my(f = factor(n)); for (k=1, #f~, f[k,2] = f[k,2]\3); factorback(f); \\ Michel Marcus, Jul 28 2017

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

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

More terms from Antti Karttunen, Jul 28 2017

A019555 Smallest number whose cube is divisible by n.

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 2, 3, 10, 11, 6, 13, 14, 15, 4, 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, 12, 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, 77, 78

Offset: 1

R. Muller

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

Peter Kagey, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Kevin A. Broughan, Restricted divisor sums, Acta Arithmetica, 101(2) (2002), 105-114.
Kevin A. Broughan, Relationship between the integer conductor and k-th root functions, Int. J. Pure Appl. Math. 5(3) (2003), 253-275.
Kevin A. Broughan, Relaxations of the ABC Conjecture using integer k'th roots, New Zealand J. Math. 35(2) (2006), 121-136.
Vaclav Kotesovec, Graph - the asymptotic ratio (1000000 terms).
Ana Rechtman, Mai 2021, 4e défi, Images des Mathématiques, CNRS, 2021 (in French).
Florentin Smarandache, Collected Papers, Vol. II, Tempus Publ. Hse, Bucharest, 1996.
Eric Weisstein's World of Mathematics, Smarandache Ceil Function.

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

Cf. A000189, A015050.

f:= n -> mul(t[1]^ceil(t[2]/3), t = ifactors(n)[2]):
map(f, [$1..100]); # Robert Israel, Sep 22 2015

cubes=Range[85]^3; Table[Position[Divisible[cubes,i],True,1,1][[1,1]],{i,85}] (* Harvey P. Dale, Jan 12 2011 *)
f[p_, e_] := p^Ceiling[e/3]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]  (* Amiram Eldar, Jan 06 2024 *)

a(n)=my(r=1);while(r^3%n!=0,r++);r \\ Anders Hellström, Sep 22 2015

for(n=1, 100, print1(direuler(p=2, n, (1 + p*X + p*X^2)/(1 - p*X^3))[n], ", ")) \\ Vaclav Kotesovec, Aug 30 2021

a(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i,1]^ceil(f[i,2]/3));} \\ Amiram Eldar, Jan 06 2024

from math import prod
from sympy import factorint
def A019555(n): return prod(p**((q%3 != 0)+(q//3)) for p, q in factorint(n).items()) # Chai Wah Wu, Aug 18 2021

[prod([t[0]^(ceil(t[1]/3)) for t in factor(n)]) for n in range(1,79)] # Danny Rorabaugh, Sep 22 2015

Replace any cubic factors in n by their cube roots.
a(n) = n/A000189(n).
Multiplicative with a(p^e) = p^ceiling(e/3). - R. J. Mathar, May 29 2011
From Vaclav Kotesovec, Aug 30 2021: (Start)
Dirichlet g.f.: zeta(3*s-1) * Product_{p prime} (1 + p^(1 - s) + p^(1 - 2*s)).
Dirichlet g.f.: zeta(3*s-1) * zeta(s-1) * Product_{p prime} (1 - p^(2 - 3*s) + p^(1 - 2*s) - p^(2 - 2*s)).
Sum_{k=1..n} a(k) ~ c * zeta(5) * n^2 / 2, where c = Product_{p prime} (1 - 1/p^2 + 1/p^3 - 1/p^4) = 0.684286924186862318141968725791218083472312736723163777284618226290055... (End)

Corrected and extended by David W. Wilson

cp's OEIS Frontend

A346597 Partial sums of A019554.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A346612 Moebius transform of A019554.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A346613 Inverse Moebius transform of A019554.

Views

Author

Keywords

Links

Crossrefs

A007947 Largest squarefree number dividing n: the squarefree kernel of n, rad(n), radical of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Sage

Scheme

Formula

Extensions

A007913 Squarefree part of n: a(n) is the smallest positive number m such that n/m is a square.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Sage

Formula

Extensions

A000188 (1) Number of solutions to x^2 == 0 (mod n). (2) Also square root of largest square dividing n. (3) Also max_{ d divides n } gcd(d, n/d).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Python

Formula