A071403 -id:A071403 - cp's OEIS frontend

A005117 Squarefree numbers: numbers that are not divisible by a square greater than 1.

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 102, 103, 105, 106, 107, 109, 110, 111, 113

Offset: 1

N. J. A. Sloane

1 together with the numbers that are products of distinct primes.
Also smallest sequence with the property that a(m)*a(k) is never a square for k != m. - Ulrich Schimke (ulrschimke(AT)aol.com), Dec 12 2001
Numbers k such that there is only one Abelian group with k elements, the cyclic group of order k (the numbers such that A000688(k) = 1). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 25 2001
Numbers k such that A007913(k) > phi(k). - Benoit Cloitre, Apr 10 2002
a(n) is the smallest m with exactly n squarefree numbers <= m. - Amarnath Murthy, May 21 2002
k is squarefree <=> k divides prime(k)# where prime(k)# = product of first k prime numbers. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 30 2004
Numbers k such that omega(k) = Omega(k) = A072047(k). - Lekraj Beedassy, Jul 11 2006
The LCM of any finite subset is in this sequence. - Lekraj Beedassy, Jul 11 2006
This sequence and the Beatty Pi^2/6 sequence (A059535) are "incestuous": the first 20000 terms are bounded within (-9, 14). - Ed Pegg Jr, Jul 22 2008
Let us introduce a function D(n) = sigma_0(n)/2^(alpha(1) + ... + alpha(r)), sigma_0(n) number of divisors of n (A000005), prime factorization of n = p(1)^alpha(1) * ... * p(r)^alpha(r), alpha(1) + ... + alpha(r) is sequence (A001222). Function D(n) splits the set of positive integers into subsets, according to the value of D(n). Squarefree numbers (A005117) has D(n)=1, other numbers are "deviated" from the squarefree ideal and have 0 < D(n) < 1. For D(n)=1/2 we have A048109, for D(n)=3/4 we have A060687. - Ctibor O. Zizka, Sep 21 2008
Numbers k such that gcd(k,k')=1 where k' is the arithmetic derivative (A003415) of k. - Giorgio Balzarotti, Apr 23 2011
Numbers k such that A007913(k) = core(k) = k. - Franz Vrabec, Aug 27 2011
Numbers k such that sqrt(k) cannot be simplified. - Sean Loughran, Sep 04 2011
Indices m where A057918(m)=0, i.e., positive integers m for which there are no integers k in {1,2,...,m-1} such that k*m is a square. - John W. Layman, Sep 08 2011
It appears that these are numbers j such that Product_{k=1..j} (prime(k) mod j) = 0 (see Maple code). - Gary Detlefs, Dec 07 2011. - This is the same claim as Mohammed Bouayoun's Mar 30 2004 comment above. To see why it holds: Primorial numbers, A002110, a subsequence of this sequence, are never divisible by any nonsquarefree number, A013929, and on the other hand, the index of the greatest prime dividing any n is less than n. Cf. A243291. - Antti Karttunen, Jun 03 2014
Conjecture: For each n=2,3,... there are infinitely many integers b > a(n) such that Sum_{k=1..n} a(k)*b^(k-1) is prime, and the smallest such an integer b does not exceed (n+3)*(n+4). - Zhi-Wei Sun, Mar 26 2013
The probability that a random natural number belongs to the sequence is 6/Pi^2, A059956 (see Cesàro reference). - Giorgio Balzarotti, Nov 21 2013
Booker, Hiary, & Keating give a subexponential algorithm for testing membership in this sequence without factoring. - Charles R Greathouse IV, Jan 29 2014
Because in the factorizations into prime numbers these a(n) (n >= 2) have exponents which are either 0 or 1 one could call the a(n) 'numbers with a fermionic prime number decomposition'. The levels are the prime numbers prime(j), j >= 1, and the occupation numbers (exponents) e(j) are 0 or 1 (like in Pauli's exclusion principle). A 'fermionic state' is then denoted by a sequence with entries 0 or 1, where, except for the zero sequence, trailing zeros are omitted. The zero sequence stands for a(1) = 1. For example a(5) = 6 = 2^1*3^1 is denoted by the 'fermionic state' [1, 1], a(7) = 10 by [1, 0, 1]. Compare with 'fermionic partitions' counted in A000009. - Wolfdieter Lang, May 14 2014
From Vladimir Shevelev, Nov 20 2014: (Start)
The following is an Eratosthenes-type sieve for squarefree numbers. For integers > 1:
1) Remove even numbers, except for 2; the minimal non-removed number is 3.
2) Replace multiples of 3 removed in step 1, and remove multiples of 3 except for 3 itself; the minimal non-removed number is 5.
3) Replace multiples of 5 removed as a result of steps 1 and 2, and remove multiples of 5 except for 5 itself; the minimal non-removed number is 6.
4) Replace multiples of 6 removed as a result of steps 1, 2 and 3 and remove multiples of 6 except for 6 itself; the minimal non-removed number is 7.
5) Repeat using the last minimal non-removed number to sieve from the recovered multiples of previous steps.
Proof. We use induction. Suppose that as a result of the algorithm, we have found all squarefree numbers less than n and no other numbers. If n is squarefree, then the number of its proper divisors d > 1 is even (it is 2^k - 2, where k is the number of its prime divisors), and, by the algorithm, it remains in the sequence. Otherwise, n is removed, since the number of its squarefree divisors > 1 is odd (it is 2^k-1).
(End)
The lexicographically least sequence of integers > 1 such that each entry has an even number of proper divisors occurring in the sequence (that's the sieve restated). - Glen Whitney, Aug 30 2015
0 is nonsquarefree because it is divisible by any square. - Jon Perry, Nov 22 2014, edited by M. F. Hasler, Aug 13 2015
The Heinz numbers of partitions with distinct parts. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product_{j=1..r} prime(j) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] the Heinz number is 2*2*3*7*29 = 2436. The number 30 (= 2*3*5) is in the sequence because it is the Heinz number of the partition [1,2,3]. - Emeric Deutsch, May 21 2015
It is possible for 2 consecutive terms to be even; for example a(258)=422 and a(259)=426. - Thomas Ordowski, Jul 21 2015. [These form a subsequence of A077395 since their product is divisible by 4. - M. F. Hasler, Aug 13 2015]
There are never more than 3 consecutive terms. Runs of 3 terms start at 1, 5, 13, 21, 29, 33, ... (A007675). - Ivan Neretin, Nov 07 2015
a(n) = product of row n in A265668. - Reinhard Zumkeller, Dec 13 2015
Numbers without excess, i.e., numbers k such that A001221(k) = A001222(k). - Juri-Stepan Gerasimov, Sep 05 2016
Numbers k such that b^(phi(k)+1) == b (mod k) for every integer b. - Thomas Ordowski, Oct 09 2016
Boreico shows that the set of square roots of the terms of this sequence is linearly independent over the rationals. - Jason Kimberley, Nov 25 2016 (reference found by Michael Coons).
Numbers k such that A008836(k) = A008683(k). - Enrique Pérez Herrero, Apr 04 2018
The prime zeta function P(s) "has singular points along the real axis for s=1/k where k runs through all positive integers without a square factor". See Wolfram link. - Maleval Francis, Jun 23 2018
Numbers k such that A007947(k) = k. - Kyle Wyonch, Jan 15 2021
The Schnirelmann density of the squarefree numbers is 53/88 (Rogers, 1964). - Amiram Eldar, Mar 12 2021
Comment from Isaac Saffold, Dec 21 2021: (Start)
Numbers k such that all groups of order k have a trivial Frattini subgroup [Dummit and Foote].
Let the group G have order n. If n is squarefree and n > 1, then G is solvable, and thus by Hall's Theorem contains a subgroup H_p of index p for all p | n. Each H_p is maximal in G by order considerations, and the intersection of all the H_p's is trivial. Thus G's Frattini subgroup Phi(G), being the intersection of G's maximal subgroups, must be trivial. If n is not squarefree, the cyclic group of order n has a nontrivial Frattini subgroup. (End)
Numbers for which the squarefree divisors (A206778) and the unitary divisors (A077610) are the same; moreover they are also the set of divisors (A027750). - Bernard Schott, Nov 04 2022
0 = A008683(a(n)) - A008836(a(n)) = A001615(a(n)) - A000203(a(n)). - Torlach Rush, Feb 08 2023
From Robert D. Rosales, May 20 2024: (Start)
Numbers n such that mu(n) != 0, where mu(n) is the Möbius function (A008683).
Solutions to the equation Sum_{d|n} mu(d)*sigma(d) = mu(n)*n, where sigma(n) is the sum of divisors function (A000203). (End)
a(n) is the smallest root of x = 1 + Sum_{k=1..n-1} floor(sqrt(x/a(k))) greater than a(n-1). - Yifan Xie, Jul 10 2024
Number k such that A001414(k) = A008472(k). - Torlach Rush, Apr 14 2025
To elaborate on the formula from Greathouse (2018), the maximum of a(n) - floor(n*Pi^2/6 + sqrt(n)/17) is 10 at indices n = 48715, 48716, 48721, and 48760. The maximum is 11, at the same indices, if floor is taken individually for the two addends and the square root. If the value is rounded instead, the maximum is 9 at 10 indices between 48714 and 48765. - M. F. Hasler, Aug 08 2025

Jean-Marie De Koninck, Ces nombres qui nous fascinent, Entry 165, p. 53, Ellipses, Paris, 2008.
David S. Dummit and Richard M. Foote, Abstract algebra. Vol. 1999. Englewood Cliffs, NJ: David S.Prentice Hall, 1991.
Ivan M. Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers. 2nd ed., Wiley, NY, 1966, p. 251.
Michael Pohst and Hans J. Zassenhaus, Algorithmic Algebraic Number Theory, Cambridge Univ. Press, page 432.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Daniel Forgues, Table of n, a(n) for n = 1..60794 (first 10000 terms from T. D. Noe)
Zenon B. Batang, Squarefree integers and the abc conjecture, arXiv:2109.10226 [math.GM], 2021.
Andrew R. Booker, Ghaith A. Hiary and Jon P. Keating, Detecting squarefree numbers, Duke Mathematical Journal, Vol. 164, No. 2 (2015), pp. 235-275; arXiv preprint, arXiv:1304.6937 [math.NT], 2013-2015.
Iurie Boreico, Linear independence of radicals, The Harvard College Mathematics Review 2(1), 87-92, Spring 2008.
Ernesto Cesàro, La serie di Lambert in aritmetica assintotica, Rendiconto della Reale Accademia delle Scienze di Napoli, Serie 2, Vol. 7 (1893), pp. 197-204.
Henri Cohen, Francois Dress, and Mohamed El Marraki, Explicit estimates for summatory functions linked to the Möbius μ-function, Functiones et Approximatio Commentarii Mathematici 37 (2007), part 1, pp. 51-63.
H. Gent, Letter to N. J. A. Sloane, Nov 27 1975.
Andrew Granville, ABC means we can count squarefrees, International Mathematical Research Notices 19 (1998), 991-1009.
Pentti Haukkanen, Mika Mattila, Jorma K. Merikoski and Timo Tossavainen, Can the Arithmetic Derivative be Defined on a Non-Unique Factorization Domain?, Journal of Integer Sequences, Vol. 16 (2013), Article 13.1.2.
Aaron Krowne, squarefree number, PlanetMath.org.
Louis Marmet, First occurrences of squarefree gaps and an algorithm for their computation.
Louis Marmet, First occurrences of square-free gaps and an algorithm for their computation, arXiv preprint arXiv:1210.3829 [math.NT], 2012.
Srinivasa Ramanujan, Irregular numbers, J. Indian Math. Soc., Vol. 5 (1913), pp. 105-106.
Kenneth Rogers, The Schnirelmann density of the squarefree integers, Proceedings of the American Mathematical Society, Vol. 15, No. 4 (1964), pp. 515-516.
J. A. Scott, Square-freedom revisited, The Mathematical Gazette, Vol. 90, No. 517 (2006), pp. 112-113.
Vladimir Shevelev, Set of all densities of exponentially S-numbers, arXiv preprint arXiv:1511.03860 [math.NT], 2015.
O. Trifonov, On the Squarefree Problem II, Math. Balkanica, Vol. 3 (1989), Fasc. 3-4.
Eric Weisstein's World of Mathematics, Squarefree.
Eric Weisstein's World of Mathematics, Prime Zeta Function.
Wikipedia, Squarefree integer.
Index entries for "core" sequences.

Complement of A013929. Subsequence of A072774 and A209061.
Characteristic function: A008966 (mu(n)^2, where mu = A008683).
Subsequences: A000040, A002110, A235488.
Cf. A001414, A008472, A076259 (first differences), A173143 (partial sums), A000688, A003277, A013928, A020753, A020754, A020755, A030059, A030229, A033197, A034444, A039956, A048672, A053797, A057918, A059956, A071403, A072284, A120992, A133466, A136742, A136743, A160764, A243289, A243347, A243348, A243351, A215366, A046660, A265668, A265675.
Cf. A027750, A077610, A206778.
Subsequences: numbers j such that j*a(k) is squarefree where k > 1: A056911 (k = 2), A261034 (k = 3), A274546 (k = 5), A276378 (k = 6).

a005117 n = a005117_list !! (n-1)
a005117_list = filter ((== 1) . a008966) [1..]
-- Reinhard Zumkeller, Aug 15 2011, May 10 2011

[ n : n in [1..1000] | IsSquarefree(n) ];

with(numtheory); a := [ ]; for n from 1 to 200 do if issqrfree(n) then a := [ op(a), n ]; fi; od:
t:= n-> product(ithprime(k),k=1..n): for n from 1 to 113 do if(t(n) mod n = 0) then print(n) fi od; # Gary Detlefs, Dec 07 2011
A005117 := proc(n) option remember; if n = 1 then 1; else for a from procname(n-1)+1 do if numtheory[issqrfree](a) then return a; end if; end do: end if; end proc:  # R. J. Mathar, Jan 09 2013

Select[ Range[ 113], SquareFreeQ] (* Robert G. Wilson v, Jan 31 2005 *)
Select[Range[150], Max[Last /@ FactorInteger[ # ]] < 2 &] (* Joseph Biberstine (jrbibers(AT)indiana.edu), Dec 26 2006 *)
NextSquareFree[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sf = n + sgn; While[c < Abs[k], While[ ! SquareFreeQ@ sf, If[sgn < 0, sf--, sf++]]; If[ sgn < 0, sf--, sf++]; c++]; sf + If[ sgn < 0, 1, -1]]; NestList[ NextSquareFree, 1, 70] (* Robert G. Wilson v, Apr 18 2014 *)
Select[Range[250], MoebiusMu[#] != 0 &] (* Robert D. Rosales, May 20 2024 *)

bnd = 1000; L = vector(bnd); j = 1; for (i=1,bnd, if(issquarefree(i),L[j]=i; j=j+1)); L

{a(n)= local(m,c); if(n<=1,n==1, c=1; m=1; while( cMichael Somos, Apr 29 2005 */

list(n)=my(v=vectorsmall(n,i,1),u,j); forprime(p=2,sqrtint(n), forstep(i=p^2, n, p^2, v[i]=0)); u=vector(sum(i=1,n,v[i])); for(i=1,n,if(v[i],u[j++]=i)); u \\ Charles R Greathouse IV, Jun 08 2012

for(n=1, 113, if(core(n)==n, print1(n, ", "))); \\ Arkadiusz Wesolowski, Aug 02 2016

S(n) = my(s); forsquarefree(k=1,sqrtint(n),s+=n\k[1]^2*moebius(k)); s;
a(n) = my(min=1, max=231, k=0, sc=0); if(n >= 144, min=floor(zeta(2)*n - 5*sqrt(n)); max=ceil(zeta(2)*n + 5*sqrt(n))); while(min <= max, k=(min+max)\2; sc=S(k); if(abs(sc-n) <= sqrtint(n), break); if(sc > n, max=k-1, if(sc < n, min=k+1, break))); while(!issquarefree(k), k-=1); while(sc != n, my(j=1); if(sc > n, j = -1); k += j; sc += j; while(!issquarefree(k), k += j)); k; \\ Daniel Suteu, Jul 07 2022

first(n)=my(v=vector(n),i); forsquarefree(k=1,if(n<268293,(33*n+30)\20,(n*Pi^2/6+0.058377*sqrt(n))\1), if(i++>n, return(v)); v[i]=k[1]); v \\ Charles R Greathouse IV, Jan 10 2023

A5117=[1..3]; A005117(n)={if(n>#A5117, my(N=#A5117); A5117=Vec(A5117, max(n+999, N*5\4)); iferr(forsquarefree(k=A5117[N]+1, #A5117*Pi^2\6+sqrtint(#A5117)\17+11, A5117[N++]=k[1]),E,)); A5117[n]} \\ M. F. Hasler, Aug 08 2025

from sympy.ntheory.factor_ import core
def ok(n): return core(n, 2) == n
print(list(filter(ok, range(1, 114)))) # Michael S. Branicky, Jul 31 2021

from itertools import count, islice
from sympy import factorint
def A005117_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:all(x == 1 for x in factorint(n).values()),count(max(startvalue,1)))
A005117_list = list(islice(A005117_gen(),20)) # Chai Wah Wu, May 09 2022

from math import isqrt
from sympy import mobius
def A005117(n):
    def f(x): return n+x-sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m # Chai Wah Wu, Jul 22 2024

Limit_{n->oo} a(n)/n = Pi^2/6 (see A013661). - Benoit Cloitre, May 23 2002
Equals A039956 UNION A056911. - R. J. Mathar, May 16 2008
A122840(a(n)) <= 1; A010888(a(n)) < 9. - Reinhard Zumkeller, Mar 30 2010
a(n) = A055229(A062838(n)) and a(n) > A055229(m) for m < A062838(n). - Reinhard Zumkeller, Apr 09 2010
A008477(a(n)) = 1. - Reinhard Zumkeller, Feb 17 2012
A055653(a(n)) = a(n); A055654(a(n)) = 0. - Reinhard Zumkeller, Mar 11 2012
A008966(a(n)) = 1. - Reinhard Zumkeller, May 26 2012
Sum_{n>=1} 1/a(n)^s = zeta(s)/zeta(2*s). - Enrique Pérez Herrero, Jul 07 2012
A056170(a(n)) = 0. - Reinhard Zumkeller, Dec 29 2012
A013928(a(n)+1) = n. - Antti Karttunen, Jun 03 2014
A046660(a(n)) = 0. - Reinhard Zumkeller, Nov 29 2015
Equals {1} UNION A000040 UNION A006881 UNION A007304 UNION A046386 UNION A046387 UNION A067885 UNION A123321 UNION A123322 UNION A115343 ... - R. J. Mathar, Nov 05 2016
|a(n) - n*Pi^2/6| < 0.058377*sqrt(n) for n >= 268293; this result can be derived from Cohen, Dress, & El Marraki, see links. - Charles R Greathouse IV, Jan 18 2018
From Amiram Eldar, Jul 07 2021: (Start)
Sum_{n>=1} (-1)^(a(n)+1)/a(n)^2 = 9/Pi^2.
Sum_{k=1..n} 1/a(k) ~ (6/Pi^2) * log(n).
Sum_{k=1..n} (-1)^(a(k)+1)/a(k) ~ (2/Pi^2) * log(n).
(all from Scott, 2006) (End)

A061399 Number of nonsquarefree integers between primes prime(n) and prime(n+1).

Original entry on oeis.org

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

Offset: 1

Labos Elemer, Jun 07 2001

nonn

			Between 113 and 127 the 7 numbers which are not squarefree are {116,117,120,121,124,125,126}, so a(30)=7.
From _Gus Wiseman_, Dec 07 2024: (Start)
The a(n) nonsquarefree numbers for n = 1..15:
   1   2   3   4   5   6   7   8   9  10  11  12  13  14  15
  ----------------------------------------------------------
   .   4   .   8  12  16  18  20  24   .  32  40   .  44  48
               9                  25      36          45  49
                                  27                      50
                                  28                      52
(End)

Harry J. Smith, Table of n, a(n) for n=1..1000

Zeros are A068361.
First differences of A378086, restriction of A057627 to the primes.
Other classes (instead of nonsquarefree):
- For composite we have A046933, first differences of A065890.
- For squarefree see A061398, A068360, A071403, A373197, A373198, A377431.
- For prime power we have A080101.
- For non prime power we have A368748, see A378616.
- For perfect power we have A377432, zeros A377436.
- For non perfect power we have A377433, A029707.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers, differences A076259.
A013929 lists the nonsquarefree numbers, differences A078147.
A120327 gives the least nonsquarefree number >= n.
Cf. A007674, A053806, A072284, A073247, A179211, A224363.
Cf. A013928, A067535, A070321, A112925, A112926, A240473, A240474.
Cf. A045882, A377049, A377783, A378032, A378618.

Count[Range[#[[1]]+1,#[[2]]-1],?(!SquareFreeQ[#]&)]&/@Partition[Prime[Range[120]],2,1] (* _Harvey P. Dale, Mar 31 2024 *)

{ n=0; q=2; forprime (p=3, prime(1001), a=0; for (i=q+1, p-1, a+=!issquarefree(i)); write("b061399.txt", n++, " ", a); q=p ) } \\ Harry J. Smith, Jul 22 2009

a(n) = my(p=prime(n)); sum(k=p, nextprime(p+1), ! issquarefree(k)); \\ Michel Marcus, Dec 09 2024

from sympy import mobius, prime
def A061399(n): return sum(not mobius(m) for m in range(prime(n)+1,prime(n+1))) # Chai Wah Wu, Jul 20 2024

A068361 Numbers n such that the number of squarefree numbers between prime(n) and prime(n+1) = prime(n+1)-prime(n)-1.

Original entry on oeis.org

1, 3, 10, 13, 26, 33, 60, 89, 104, 113, 116, 142, 148, 201, 209, 212, 234, 265, 268, 288, 313, 320, 332, 343, 353, 384, 398, 408, 477, 484, 498, 542, 545, 551, 577, 581, 601, 625, 636, 671, 719, 723, 726, 745, 794, 805, 815, 862, 864, 884, 944, 964, 995, 1054

Offset: 1

Benoit Cloitre, Feb 28 2002

nonn

Also numbers k such that all numbers from prime(k) to prime(k+1) are squarefree. All such primes are twins, so this is a subset of A029707. The other twin primes are A061368. - Gus Wiseman, Dec 11 2024

Amiram Eldar, Table of n, a(n) for n = 1..10000

A subset of A029707 (lesser index of twin primes).
Prime index of each (prime) term of A061351.
Positions of zeros in A061399.
For perfect power instead of squarefree we have A377436, zeros of A377432.
Positions of zeros in A377784.
The rest of the twin primes are at A378620, indices of A061368.
A000040 lists the primes, differences A001223, (run-lengths A333254, A373821).
A005117 lists the squarefree numbers, differences A076259.
A006562 finds balanced primes.
A013929 lists the nonsquarefree numbers, differences A078147.
A014574 is the intersection of A006093 and A008864.
A038664 locates the first prime gap of size 2n.
A046933 counts composite numbers between primes.
A061398 counts squarefree numbers between primes, zeros A068360.
A120327 gives the least nonsquarefree number >= n.
Cf. A000720, A007674, A013928, A057627, A070321, A071403, A072284, A073247, A112925, A122535, A155752, A224363, A240473, A251092.

Select[Range[100],And@@SquareFreeQ/@Range[Prime[#],Prime[#+1]]&] (* Gus Wiseman, Dec 11 2024 *)

isok(n) = for (k=prime(n)+1, prime(n+1)-1, if (!issquarefree(k), return (0))); 1; \\ Michel Marcus, Apr 29 2016

n such that A061398(n) = prime(n+1)-prime(n)-1.

prime(a(n)) = A061351(n). - Gus Wiseman, Dec 11 2024

A378086 Number of nonsquarefree numbers < prime(n).

Original entry on oeis.org

0, 0, 1, 1, 3, 4, 5, 6, 7, 11, 11, 13, 14, 14, 16, 20, 22, 23, 25, 26, 27, 29, 31, 33, 36, 39, 39, 40, 41, 42, 49, 50, 53, 53, 57, 58, 61, 63, 64, 68, 70, 71, 74, 75, 76, 77, 81, 84, 86, 87, 88, 90, 91, 97, 99, 101, 103, 104, 107, 109, 109, 113, 119, 120, 121

Offset: 1

Gus Wiseman, Dec 04 2024

nonn

			The nonsquarefree numbers counted under each term begin:
  n=1: n=2: n=3: n=4: n=5: n=6: n=7: n=8: n=9: n=10: n=11: n=12:
  --------------------------------------------------------------
   .    .    4    4    9    12   16   18   20   28    28    36
                       8    9    12   16   18   27    27    32
                       4    8    9    12   16   25    25    28
                            4    8    9    12   24    24    27
                                 4    8    9    20    20    25
                                      4    8    18    18    24
                                           4    16    16    20
                                                12    12    18
                                                9     9     16
                                                8     8     12
                                                4     4     9
                                                            8
                                                            4

For nonprime numbers we have A014689.
Restriction of A057627 to the primes.
First-differences are A061399 (zeros A068361), squarefree A061398 (zeros A068360).
For composite instead of squarefree we have A065890.
For squarefree we have A071403, differences A373198.
Greatest is A378032 (differences A378034), restriction of A378033 (differences A378036).
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147, seconds A376593.
A070321 gives the greatest squarefree number up to n.
A112925 gives the greatest squarefree number between primes, differences A378038.
A112926 gives the least squarefree number between primes, differences A378037.
A120327 gives the least nonsquarefree number >= n, first-differences A378039.
A377783 gives the least nonsquarefree > prime(n), differences A377784.
Cf. A013928, A046933, A053797, A053806, A072284, A076259, A224363, A337030, A377049, A378040, A378082, A378084.

Table[Length[Select[Range[Prime[n]],!SquareFreeQ[#]&]],{n,100}]

from math import isqrt
from sympy import prime, mobius
def A378086(n): return (p:=prime(n))-sum(mobius(k)*(p//k**2) for k in range(1,isqrt(p)+1)) # Chai Wah Wu, Dec 05 2024

a(n) = A057627(prime(n)).

A378040 Union of A377783(n) = least nonsquarefree number > prime(n).

Original entry on oeis.org

4, 8, 12, 16, 18, 20, 24, 32, 40, 44, 48, 54, 60, 63, 68, 72, 75, 80, 84, 90, 98, 104, 108, 112, 116, 128, 132, 140, 150, 152, 160, 164, 168, 175, 180, 184, 192, 196, 198, 200, 212, 224, 228, 232, 234, 240, 242, 252, 260, 264, 270, 272, 279, 284, 294, 308, 312

Offset: 1

Gus Wiseman, Nov 20 2024

nonn

Numbers k such that, if p is the greatest prime < k, all numbers from p to k (exclusive) are squarefree.

Harvey P. Dale, Table of n, a(n) for n = 1..1000

For squarefree we have A112926 (diffs A378037), opposite A112925 (diffs A378038).
For prime-power instead of nonsquarefree we have A345531, differences A377703.
Union of A377783 (diffs A377784), restriction of A120327 (diffs A378039).
Nonsquarefree numbers not appearing are A378084, see also A378082, A378083.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147, seconds A376593.
A061398 counts squarefree numbers between primes, zeros A068360.
A061399 counts nonsquarefree numbers between primes, zeros A068361.
A070321 gives the greatest squarefree number up to n.
A071403(n) = A013928(prime(n)) counts squarefree numbers up to prime(n).
A378086(n) = A057627(prime(n)) counts nonsquarefree numbers up to prime(n).
Cf. A378034 (differences of A378032), restriction of A378036 (differences A378033).
Cf. A000015, A053797, A053806, A072284, A224363, A337030, A377047, A377049, A377430, A377431.

Union[Table[NestWhile[#+1&,Prime[n],SquareFreeQ],{n,100}]]
lns[p_]:=Module[{k=p+1},While[SquareFreeQ[k],k++];k]; Table[lns[p],{p,Prime[Range[70]]}]//Union (* Harvey P. Dale, Jun 12 2025 *)

A378082 Terms appearing only once in A377783 = least nonsquarefree number > prime(n).

Original entry on oeis.org

12, 16, 18, 20, 24, 40, 48, 54, 60, 63, 68, 72, 75, 80, 84, 90, 98, 108, 112, 116, 128, 132, 150, 152, 160, 164, 168, 175, 180, 184, 192, 196, 198, 200, 212, 224, 228, 232, 234, 240, 242, 252, 260, 264, 270, 272, 279, 294, 308, 312, 315, 320, 332, 338, 348

Offset: 1

Gus Wiseman, Nov 20 2024

nonn

Nonsquarefree numbers k such that if p < q are the two greatest primes < k, there is at least one nonsquarefree number between p and q but all numbers between q and k are squarefree. - Robert Israel, Nov 20 2024

			The terms together with their prime indices begin:
   12: {1,1,2}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   24: {1,1,1,2}
   40: {1,1,1,3}
   48: {1,1,1,1,2}
   54: {1,2,2,2}
   60: {1,1,2,3}
   63: {2,2,4}
   68: {1,1,7}
   72: {1,1,1,2,2}
   75: {2,3,3}
   80: {1,1,1,1,3}
   84: {1,1,2,4}
   90: {1,2,2,3}
   98: {1,4,4}
  108: {1,1,2,2,2}
  112: {1,1,1,1,4}
  116: {1,1,10}
  128: {1,1,1,1,1,1,1}
  132: {1,1,2,5}

This is a transformation of A377783 (union A378040, differences A377784).
Note also A377783 restricts A120327 (differences A378039) to the primes.
Terms appearing twice are A378083.
Terms not appearing at all are A378084.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147, seconds A376593.
A061398 counts squarefree numbers between primes, zeros A068360.
A061399 counts nonsquarefree numbers between primes, zeros A068361.
A070321 gives the greatest squarefree number up to n.
A071403(n) = A013928(prime(n)) counts squarefree numbers < prime(n).
A378086(n) = A057627(prime(n)) counts nonsquarefree numbers < prime(n).
Cf. A112926 (diffs A378037), opposite A112925 (diffs A378038).
Cf. A378032 (diffs A378034), restriction of A378033 (diffs A378036).
Cf. A053797, A053806, A072284, A112929, A120992, A224363, A337030, A377430, A377431, A377703.

q:= 3: R:= NULL: flag:= false: count:= 0:
while count < 100 do
  p:= q; q:= nextprime(q);
  for k from p+1 to q-1 do
    found:= false;
    if not numtheory:-issqrfree(k) then
      if flag then
          count:= count+1; R:= R,k
      fi;
      found:= true; break
    fi;
   od;
   flag:= found;
od:
R; # Robert Israel, Nov 20 2024

y=Table[NestWhile[#+1&,Prime[n],SquareFreeQ],{n,100}];
Select[Most[Union[y]],Count[y,#]==1&]

A261144 Irregular triangle of numbers that are squarefree and smooth (row n contains squarefree p-smooth numbers, where p is the n-th prime).

Original entry on oeis.org

1, 2, 1, 2, 3, 6, 1, 2, 3, 5, 6, 10, 15, 30, 1, 2, 3, 5, 6, 7, 10, 14, 15, 21, 30, 35, 42, 70, 105, 210, 1, 2, 3, 5, 6, 7, 10, 11, 14, 15, 21, 22, 30, 33, 35, 42, 55, 66, 70, 77, 105, 110, 154, 165, 210, 231, 330, 385, 462, 770, 1155, 2310, 1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 21, 22, 26, 30, 33, 35, 39, 42

Offset: 1

Jean-François Alcover, Nov 26 2015

If we define a triangle whose n-th row consists of all squarefree numbers whose prime factors are all less than prime(k), we get this same triangle except starting with a row {1}, with offset 1. - Gus Wiseman, Aug 24 2021

			Triangle begins:
1, 2;                        squarefree and 2-smooth
1, 2, 3, 6;                  squarefree and 3-smooth
1, 2, 3, 5, 6, 10, 15, 30;
1, 2, 3, 5, 6,  7, 10, 14, 15, 21, 30, 35, 42, 70, 105, 210;
...

Jean-François Alcover, Table of n, a(n) for n = 1..2046 (first 10 rows)
A. Hildebrand and G. Tenenbaum, Integers without large prime factors, Journal de théorie des nombres de Bordeaux (1993) Volume:5, Issue:2, p. 411-484.
Eric Weisstein's MathWorld, Smooth number.
Wikipedia, Smooth number

Cf. A000079 (2-smooth), A003586 (3-smooth), A051037 (5-smooth), A002473 (7-smooth), A018336 (7-smooth & squarefree), A051038 (11-smooth), A087005 (11-smooth & squarefree), A080197 (13-smooth), A087006 (13-smooth & squarefree), A087007 (17-smooth & squarefree), A087008 (19-smooth & squarefree).
Row lengths are A000079.
Rightmost terms (or column k = 2^n) are A002110.
Rows are partial unions of rows of A019565.
Row n is A027750(A002110(n)), i.e., divisors of primorials.
Row sums are A054640.
Column k = 2^n-1 is A070826.
Multiplying row n by prime(n+1) gives A339195, row sums A339360.
A005117 lists squarefree numbers.
A056239 adds up prime indices, row sums of A112798.
A072047 counts prime factors of squarefree numbers.
A246867 groups squarefree numbers by Heinz weight, row sums A147655.
A329631 lists prime indices of squarefree numbers, sums A319246.
A339116 groups squarefree semiprimes by greater factor, sums A339194.
Cf. A000040, A001221, A006881, A071403, A209862, A319247.

b:= proc(n) option remember; `if`(n=0, [1],
      sort(map(x-> [x, x*ithprime(n)][], b(n-1))))
    end:
T:= n-> b(n)[]:
seq(T(n), n=1..7);  # Alois P. Heinz, Nov 28 2015

primorial[n_] := Times @@ Prime[Range[n]]; row[n_] := Select[ Divisors[ primorial[n]], SquareFreeQ]; Table[row[n], {n, 1, 10}] // Flatten

T(n-1,k) = A339195(n,k)/prime(n). - Gus Wiseman, Aug 24 2021

A372541 Least k such that the k-th squarefree number has exactly n ones in its binary expansion.

Original entry on oeis.org

1, 3, 6, 11, 20, 60, 78, 157, 314, 624, 1245, 3736, 4982, 9962, 19920, 39844, 79688, 239046, 318725, 956194, 1912371, 2549834, 5099650, 15298984, 20398664, 40797327, 81594626, 163189197, 326378284, 979135127, 1305513583, 2611027094, 5222054081, 10444108051

Offset: 0

Gus Wiseman, May 09 2024

			The squarefree numbers A005117(a(n)) together with their binary expansions and binary indices begin:
       1:                   1 ~ {1}
       3:                  11 ~ {1,2}
       7:                 111 ~ {1,2,3}
      15:                1111 ~ {1,2,3,4}
      31:               11111 ~ {1,2,3,4,5}
      95:             1011111 ~ {1,2,3,4,5,7}
     127:             1111111 ~ {1,2,3,4,5,6,7}
     255:            11111111 ~ {1,2,3,4,5,6,7,8}
     511:           111111111 ~ {1,2,3,4,5,6,7,8,9}
    1023:          1111111111 ~ {1,2,3,4,5,6,7,8,9,10}
    2047:         11111111111 ~ {1,2,3,4,5,6,7,8,9,10,11}
    6143:       1011111111111 ~ {1,2,3,4,5,6,7,8,9,10,11,13}
    8191:       1111111111111 ~ {1,2,3,4,5,6,7,8,9,10,11,12,13}
   16383:      11111111111111 ~ {1,2,3,4,5,6,7,8,9,10,11,12,13,14}
   32767:     111111111111111 ~ {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15}
   65535:    1111111111111111 ~ {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16}
  131071:   11111111111111111 ~ {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17}

Chai Wah Wu, Table of n, a(n) for n = 0..56
Wikipedia, Hamming weight.

Positions of firsts appearances in A372433.
Counting zeros instead of ones gives A372473, firsts in A372472.
For prime instead of squarefree we have A372517, firsts of A014499.
Counting bits (length) gives A372540, firsts of A372475, runs A077643.
A000120 counts ones in binary expansion (binary weight), zeros A080791.
A005117 lists squarefree numbers.
A030190 gives binary expansion, reversed A030308.
A048793 lists positions of ones in reversed binary expansion, sum A029931.
A145037, A097110 count ones minus zeros, for primes A372516, A177796.
A371571 lists positions of zeros in binary expansion, sum A359359.
A371572 lists positions of ones in binary expansion, sum A230877.
A372515 lists positions of zeros in reversed binary expansion, sum A359400.
Cf. A023416, A049093, A049094, A069010, A070939, A071403, A211997, A280296, A372474.

nn=10000;
spnm[y_]:=Max@@NestWhile[Most,y,Union[#]!=Range[0,Max@@#]&];
dcs=DigitCount[Select[Range[nn],SquareFreeQ],2,1];
Table[Position[dcs,i][[1,1]],{i,spnm[dcs-1]}]

from math import isqrt
from itertools import count
from sympy import factorint, mobius
from sympy.utilities.iterables import multiset_permutations
def A372541(n):
    if n==0: return 1
    for l in count(n):
        m = 1<Chai Wah Wu, May 10 2024

a(23)-a(33) from Chai Wah Wu, May 10 2024

A284311 Array T(n,k) read by antidiagonals (downward): T(1,k) = A005117(k+1) (squarefree numbers > 1); for n > 1, columns are nonsquarefree numbers (in ascending order) with exactly the same prime factors as T(1,k).

Original entry on oeis.org

2, 3, 4, 5, 9, 8, 6, 25, 27, 16, 7, 12, 125, 81, 32, 10, 49, 18, 625, 243, 64, 11, 20, 343, 24, 3125, 729, 128, 13, 121, 40, 2401, 36, 15625, 2187, 256, 14, 169, 1331, 50, 16807, 48, 78125, 6561, 512, 15, 28, 2197, 14641, 80, 117649, 54, 390625, 19683, 1024

Offset: 1

Bob Selcoe, Mar 24 2017

A permutation of the natural numbers > 1.

T(1,k)= A005117(m) with m > 1; terms in column k = T(1,k) * A162306(T(1,k)) only not bounded by T(1,k). Let T(1,k) = b. Terms in column k are multiples of b and numbers c such that c | b^e with e >= 0. Alternatively, terms in column k are multiples bc with c those numbers whose prime divisors p also divide b. - Michael De Vlieger, Mar 25 2017

			Array starts:
    2    3     5  6      7  10       11        13  14  15
    4    9    25 12     49  20      121       169  28  45
    8   27   125 18    343  40     1331      2197  56  75
   16   81   625 24   2401  50    14641    371293  98 135
   32  243  3125 36  16807  80   161051   4826809 112 225
   64  729 15625 48 117649 100  1771561  62748517 196 375
  128 2187 78125 54 823543 160 19487171 815730721 224 405
Column 6 is: T(1,6) = 2*5; T(2,6) = 2^2*5; T(3,6) = 2^3*5; T(4,6) = 2*5^2; T(5,6) = 2^4*5, etc.

Alois P. Heinz, Antidiagonals n = 1..141, flattened

Cf. A005117 (squarefree numbers), A033845 (column 4), columns 1,2,3,5 are powers of primes, A033846 (column 6), A033847 (column 9), A033849 (column 10).
The columns that are powers of primes have indices A071403(n) - 1. - Michel Marcus, Mar 24 2017
See also A007947; the k-th column of the array corresponds to the numbers with radical A005117(k+1). - Rémy Sigrist, Mar 24 2017
Cf. A284457 (this sequence read by antidiagonals upwards), A285321 (a similar array, but columns come in different order).
Cf. A065642.
Cf. A008479 (index of the row where n is located), A285329 (of the column).

f[n_, k_: 1] := Block[{c = 0, sgn = Sign[k], sf}, sf = n + sgn; While[c < Abs[k], While[! SquareFreeQ@ sf, If[sgn < 0, sf--, sf++]]; If[sgn < 0, sf--, sf++]; c++]; sf + If[sgn < 0, 1, -1]] (* after Robert G. Wilson v at A005117 *); T[n_, k_] := T[n, k] = Which[And[n == 1, k == 1], 2, k == 1, f@ T[n - 1, k], PrimeQ@ T[n, 1], T[n, 1]^k, True, Module[{j = T[n, k - 1]/T[n, 1] + 1}, While[PowerMod[T[n, 1], j, j] != 0, j++]; j T[n, 1]]]; Table[T[n - k + 1, k], {n, 10}, {k, n}] // Flatten (* Michael De Vlieger, Mar 25 2017 *)

(define (A284311 n) (A284311bi  (A002260 n) (A004736 n)))
(define (A284311bi row col) (if (= 1 row) (A005117 (+ 1 col)) (A065642 (A284311bi (- row 1) col))))
;; Antti Karttunen, Apr 17 2017

From Antti Karttunen, Apr 17 2017: (Start)
A(1,k) = A005117(1+k), A(n,k) = A065642(A(n-1,k)).
A(A008479(n), A285329(n)) = n for all n >= 2.
(End)

A378083 Nonsquarefree numbers appearing exactly twice in A377783 (least nonsquarefree number > prime(n)).

Original entry on oeis.org

4, 8, 32, 44, 104, 140, 284, 464, 572, 620, 644, 824, 860, 1232, 1292, 1304, 1484, 1700, 1724, 1880, 2084, 2132, 2240, 2312, 2384, 2660, 2732, 2804, 3392, 3464, 3560, 3920, 3932, 4004, 4220, 4244, 4424, 4640, 4724, 5012, 5444, 5480, 5504, 5660, 6092, 6200

Offset: 1

Gus Wiseman, Nov 23 2024

nonn

Warning: do not confuse with A377783.

			The terms together with their prime indices begin:
     4: {1,1}
     8: {1,1,1}
    32: {1,1,1,1,1}
    44: {1,1,5}
   104: {1,1,1,6}
   140: {1,1,3,4}
   284: {1,1,20}
   464: {1,1,1,1,10}
   572: {1,1,5,6}
   620: {1,1,3,11}
   644: {1,1,4,9}
   824: {1,1,1,27}
   860: {1,1,3,14}
  1232: {1,1,1,1,4,5}

Subset of A377783 (union A378040, diffs A377784), restriction of A120327 (diffs A378039).
Terms appearing once are A378082.
Terms not appearing at all are A378084.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A061398 counts squarefree numbers between primes, zeros A068360.
A061399 counts nonsquarefree numbers between primes, zeros A068361.
A071403(n) = A013928(prime(n)) counts squarefree numbers < prime(n).
A378086(n) = A057627(prime(n)) counts nonsquarefree numbers < prime(n).
Cf. A112926 (diffs A378037), opposite A112925 (diffs A378038).
Cf. A378032 (diffs A378034), restriction of A378033 (diffs A378036).
Cf. A053797, A053806, A070321, A072284, A112929, A120992, A224363, A337030, A377430, A377431, A377703.

y=Table[NestWhile[#+1&,Prime[n],SquareFreeQ[#]&],{n,1000}];
Select[Union[y],Count[y,#]==2&]

cp's OEIS Frontend

A005117 Squarefree numbers: numbers that are not divisible by a square greater than 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

PARI

PARI

Python

Python

Python

Formula

A061399 Number of nonsquarefree integers between primes prime(n) and prime(n+1).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

A068361 Numbers n such that the number of squarefree numbers between prime(n) and prime(n+1) = prime(n+1)-prime(n)-1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A378086 Number of nonsquarefree numbers < prime(n).

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Python

Formula

A378040 Union of A377783(n) = least nonsquarefree number > prime(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

A378082 Terms appearing only once in A377783 = least nonsquarefree number > prime(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

A261144 Irregular triangle of numbers that are squarefree and smooth (row n contains squarefree p-smooth numbers, where p is the n-th prime).

Views

Author