login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A007913 Squarefree part of n: a(n) = smallest positive number m such that n/m is a square. 170

%I

%S 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,

%T 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,

%U 6,55,14,57,58,59,15,61,62,7,1,65,66,67,17,69,70,71,2,73,74,3,19,77

%N Squarefree part of n: a(n) = smallest positive number m such that n/m is a square.

%C Also called core(n). [Not to be confused with the squarefree kernel of n, A007947.]

%C Sequence read mod 4 gives A065882. - _Philippe Deléham_, Mar 28 2004

%C This is an arithmetic function and is undefined if n <= 0.

%C 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]

%C 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

%C a(n) = the smallest natural numbers 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

%C Dirichlet convolution of A010052 with the sequence of absolute values of A055615. - _R. J. Mathar_, Feb 11 2011

%C 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

%C 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

%C The Dirichlet inverse of this sequence is A008836(n) * A063659(n). - _Álvar Ibeas_, Mar 19 2015

%C a(n) = 1 if n is a square, a(n) = n if n is a product of distinct primes. - _Zak Seidov_, Jan 30 2016

%C 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

%D K. 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.

%H T. D. Noe and Daniel Forgues, <a href="/A007913/b007913.txt">Table of n, a(n) for n = 1..100000</a> (first 1000 terms from T. D. Noe)

%H K. Atanassov, <a href="http://www.gallup.unm.edu/~smarandache/Atanassov-SomeProblems.pdf">On Some of Smarandache's Problems</a>

%H Andrew Booker, Ghaith Hiary, and Jon Keating, <a href="http://www.cs.uleth.ca/~cnta2012/slides-cnta12/Andrew-Booker.pdf">Detecting squarefree numbers</a>, CNTA XII (2012).

%H H. Bottomley, <a href="http://fs.gallup.unm.edu/Bottomley-Sm-Mult-Functions.htm">Some Smarandache-type multiplicative sequences</a>

%H John M. Campbell, <a href="http://arxiv.org/abs/1105.3399">An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences</a>, arXiv preprint arXiv:1105.3399 [math.GM], 2011.

%H F. Smarandache, <a href="http://www.gallup.unm.edu/~smarandache/OPNS.pdf">Only Problems, Not Solutions!</a>, Xiquan Publ., Phoenix-Chicago, 1993.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/SquarefreePart.html">Squarefree Part</a>

%F Multiplicative with a(p^k) = p^(k mod 2). - _David W. Wilson_, Aug 01 2001

%F a(n) modulo 2 = A035263(n); a(A036554(n)) is even; a(A003159(n)) is odd. - _Philippe Deléham_, Mar 28 2004

%F Dirichlet g.f.: zeta(2s)*zeta(s-1)/zeta(2s-2). - _R. J. Mathar_, Feb 11 2011

%F a(n) = n/( Sum_{k=1..n} floor(k^2/n)-floor((k^2 -1)/n) )^2. - _Anthony Browne_, Jun 06 2016

%F 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

%p 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

%p # second Maple program:

%p a:= n-> mul(i[1]^irem(i[2], 2), i=ifactors(n)[2]):

%p seq(a(n), n=1..100); # _Alois P. Heinz_, Jul 20 2015

%t data = Table[Sqrt[n], {n, 1, 100}]; sp = data /. Sqrt[_] -> 1; sfp = data/sp /. Sqrt[x_] -> x (* _Artur Jasinski_, Nov 03 2008 *)

%t Table[Times@@Power@@@({#[[1]],Mod[ #[[2]],2]}&/@FactorInteger[n]),{n,100}] (* _Zak Seidov_, Apr 08 2009 *)

%t Table[{p, e} = Transpose[FactorInteger[n]]; Times @@ (p^Mod[e, 2]), {n, 100}] (* _T. D. Noe_, May 20 2013 *)

%t √#/.(c_:1)*a_^(b_:0)->(c*a^b)^2&/@Range@100 (* _Bill Gosper_, Jul 18 2015 *)

%o (MAGMA) [ Squarefree(n) : n in [1..256] ]; // _N. J. A. Sloane_, Dec 23 2006

%o (PARI) a(n)=core(n)

%o (Haskell)

%o a007913 n = product $

%o zipWith (^) (a027748_row n) (map (`mod` 2) $ a124010_row n)

%o -- _Reinhard Zumkeller_, Jul 06 2012

%o (Python)

%o from operator import mul

%o from functools import reduce

%o from sympy import factorint

%o def A007913(n):

%o ....return reduce(mul,[1]+[p for p,e in factorint(n).items() if e % 2])

%o # _Chai Wah Wu_, Feb 03 2015

%o (Sage)

%o [squarefree_part(n) for n in (1..77)] # _Peter Luschny_, Feb 04 2015

%Y Cf. A000188, A002734, A117811, A007947, A019554, A027748, A124010.

%K nonn,easy,mult,nice

%O 1,2

%A R. Muller, Mar 15 1996

%E More terms from _Michael Somos_, Nov 24 2001

%E Definition reformulated by _Daniel Forgues_, Mar 24 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 26 04:32 EDT 2019. Contains 321481 sequences. (Running on oeis4.)