The OEIS is supported by the many generous donors to the OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A001694 Powerful numbers, definition (1): if a prime p divides n then p^2 must also divide n (also called squareful, square full, square-full or 2-powerful numbers). (Formerly M3325 N1335) 250
 1, 4, 8, 9, 16, 25, 27, 32, 36, 49, 64, 72, 81, 100, 108, 121, 125, 128, 144, 169, 196, 200, 216, 225, 243, 256, 288, 289, 324, 343, 361, 392, 400, 432, 441, 484, 500, 512, 529, 576, 625, 648, 675, 676, 729, 784, 800, 841, 864, 900, 961, 968, 972, 1000 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 COMMENTS In other words if the prime factorization of n is Product_k p_k^e_k then all e_k are greater than 1. Numbers n such that Sum_{d|n} phi(d)*phi(n/d)*mu(d) > 0; places of nonzero A300717. - Benoit Cloitre, Nov 30 2002 This sequence is closed under multiplication. The primitive elements are A168363. - Franklin T. Adams-Watters, May 30 2011 Complement of A052485; A112526(a(n)) = 1. - Reinhard Zumkeller, Sep 16 2011 a(n) = A224866(n) - 1. - Reinhard Zumkeller, Jul 23 2013 The number of terms less than or equal to 10^k beginning with k = 0: 1, 4, 14, 54, 185, 619, 2027, 6553, 21044, ...: A118896. - Robert G. Wilson v, Aug 11 2014 a(10^n): 1, 49, 3136, 253472, 23002083, 2200079025, 215523459072, 21348015504200, 2125390162618116, ... . - Robert G. Wilson v, Aug 15 2014 a(m) mod prime(n) > 0 for m < A258599(n); a(A258599(n)) = A001248(n) = prime(n)^2. - Reinhard Zumkeller, Jun 06 2015 From Des MacHale, Mar 07 2021: (Start) A number m is powerful if and only if |R/Z(R)| = m, for some finite non-commutative ring R. A number m is powerful if and only if |G/Z(G)| = m, for some finite nilpotent class two group G (Reference Aine Nishe). (End) Numbers n such that Sum_{k=1..n} phi(gcd(n,k))*mu(gcd(n,k)) > 0. - Richard L. Ollerton, May 09 2021 REFERENCES G. E. Hardy and M. V. Subbarao, Highly powerful numbers, Congress. Numer. 37 (1983), 277-307. A. Ivic, The Riemann Zeta-Function, Wiley, NY, 1985, see p. 407. R. A. Mollin, Quadratics, CRC Press, 1996, Section 1.6. Aine NiShe, Commutativity and Generalisations in Finite Groups, Ph.D. Thesis, University College Cork, 2000. Paulo Ribenboim, Meine Zahlen, meine Freunde, 2009, Springer, 9.1 Potente Zahlen, pp. 241-247. N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence). N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence). G. Tenenbaum, Introduction to analytic and probabilistic number theory, Cambridge University Press, 1995, p. 54, exercise 10 (in the third edition 2015, p. 63, exercise 70). LINKS G. C. Greubel, Table of n, a(n) for n = 1..5000 (terms 1..1000 from T. D. Noe) Paul T. Bateman and Emil Grosswald, On a theorem of Erdős and Szekeres, Illinois J. Math. 2:1 (1958), pp. 88-98. Valentin Blomer, Binary quadratic forms with large discriminants and sums of two squareful numbers II, Journal of the London Mathematical Society 71:1 (2005), pp. 69-84. C. K. Caldwell, Powerful Numbers J.-M. de Koninck, N. Doyon, and F. Luca, Powerful Values of Quadratic Polynomials, J. Int. Seq. 14 (2011) # 11.3.3. P. Erdős and G. Szekeres, Über die Anzahl der Abelschen Gruppen gegebener Ordnung und über ein verwandtes zahlentheoretisches Problem, Acta Sci. Math. (Szeged), 7 (1935), 95-102. [Zahlen i-ter Art, p. 101] S. W. Golomb, Powerful numbers, Amer. Math. Monthly, Vol. 77 (1970), 848-852. K. Schneider, PlanetMath.org, Squarefull Number V. Shevelev, S-exponential numbers, Acta Arithmetica, Vol. 175(2016), 385-395. D. Suryanarayana and R. Sita Rama Chandra Rao, The distribution of square-full integers, Ark. Mat., Volume 11, Number 1-2 (1973), 195-201. Eric Weisstein's World of Mathematics, Powerful Number Eric Weisstein's World of Mathematics, Squareful Wikipedia, Powerful number FORMULA Numbers of the form a^2*b^3, a >= 1, b >= 1. Bateman & Grosswald prove that there are zeta(3/2)/zeta(3) x^{1/2} + zeta(2/3)/zeta(2) x^{1/3} + O(x^{1/6}) terms up to x; see section 5 for a more precise error term. - Charles R Greathouse IV, Nov 19 2012 Sum_{n>=1} 1/a(n) = zeta(2)*zeta(3)/zeta(6). - Ivan Neretin, Aug 30 2015 EXAMPLE 1 is a term because for every prime p that divides 1, p^2 also divides 1. 2 is not a term since 2 divides 2 but 2^2 does not. 4 is a term because 2 is the only prime that divides 4 and 2^2 does divide 4. - N. J. A. Sloane, Jan 16 2022 MAPLE isA001694 := proc(n) for p in ifactors(n) do if op(2, p) = 1 then return false; end if; end do; return true; end proc: A001694 := proc(n) option remember; if n = 1 then 1; else for a from procname(n-1)+1 do if isA001694(a) then return a; end if; end do; end if; end proc: seq(A001694(n), n=1..20) ; # R. J. Mathar, Jun 07 2011 MATHEMATICA Join[{1}, Select[ Range@ 1250, Min@ FactorInteger[#][[All, 2]] > 1 &]] (* Harvey P. Dale, Sep 18 2011; modified by Robert G. Wilson v, Aug 11 2014 *) max = 10^3; Union@ Flatten@ Table[a^2*b^3, {b, max^(1/3)}, {a, Sqrt[max/b^3]}] (* Robert G. Wilson v, Aug 11 2014 *) nextPowerfulNumber[n_] := Block[{r = Range[ Floor[1 + n^(1/3)]]^3}, Min@ Select[ Sort[ r*Floor[1 + Sqrt[n/r]]^2], # > n &]]; NestList[ nextPowerfulNumber, 1, 55] (* Robert G. Wilson v, Aug 16 2014 *) PROG (PARI) isA001694(n)=n=factor(n)[, 2]; for(i=1, #n, if(n[i]==1, return(0))); 1 \\ Charles R Greathouse IV, Feb 11 2011 (PARI) list(lim, mn=2)=my(v=List(), t); for(m=1, sqrtnint(lim\1, 3), t=m^3; for(n=1, sqrtint(lim\t), listput(v, t*n^2))); Set(v) \\ Charles R Greathouse IV, Jul 31 2011; edited Sep 22 2015 (PARI) is=ispowerful \\ Charles R Greathouse IV, Nov 13 2012 (Haskell) a001694 n = a001694_list !! (n-1) a001694_list = filter ((== 1) . a112526) [1..] -- Reinhard Zumkeller, Nov 30 2012 (Python) from sympy import factorint A001694 = +[n for n in range(2, 10**6) if min(factorint(n).values()) > 1] # Chai Wah Wu, Aug 14 2014 (Sage) sloane.A001694.list(54) # Peter Luschny, Feb 08 2015 CROSSREFS Cf. A007532 (Powerful numbers, definition (2)), A005934, A005188, A003321, A014576, A023052 (Powerful numbers, definition (3)), A046074, A013929, A076871, A258599, A001248, A261883. Cf. A076446 (first differences). Sequence in context: A339497 A348121 A080366 * A317102 A157985 A001597 Adjacent sequences:  A001691 A001692 A001693 * A001695 A001696 A001697 KEYWORD nonn,nice,easy AUTHOR EXTENSIONS More terms from Henry Bottomley, Mar 16 2000 Definition expanded by Jonathan Sondow, Jan 03 2016 STATUS approved

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

Last modified September 27 09:48 EDT 2022. Contains 357054 sequences. (Running on oeis4.)