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-10 of 48 results. Next

A055556 Differences of arithmetic progressions of which the terms give chains of n consecutive nonsquarefree numbers if started with terms of A045882(n) or A051681(n).

Original entry on oeis.org

4, 72, 58800, 1757711340, 54169994838960, 950645196424756293060600, 94630125612279498512066506747436400, 14119054639791549212725337736060964803166626000
Offset: 1

Views

Author

Labos Elemer, Jul 10 2000

Keywords

Examples

			For n=5, the difference of the relevant progression is 54169994838960.
		

Crossrefs

Formula

a(n) = LCM(x, x+1, ..., x+n-1), where x = A045882(n).

A013929 Numbers that are not squarefree. Numbers that are divisible by a square greater than 1. The complement of A005117.

Original entry on oeis.org

4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 49, 50, 52, 54, 56, 60, 63, 64, 68, 72, 75, 76, 80, 81, 84, 88, 90, 92, 96, 98, 99, 100, 104, 108, 112, 116, 117, 120, 121, 124, 125, 126, 128, 132, 135, 136, 140, 144, 147, 148, 150, 152, 153, 156, 160
Offset: 1

Views

Author

Keywords

Comments

Sometimes misnamed squareful numbers, but officially those are given by A001694.
This is different from the sequence of numbers k such that A007913(k) < phi(k). The two sequences differ at the values: 420, 660, 780, 840, 1320, 1560, 4620, 5460, 7140, ..., which is essentially A070237. - Ant King, Dec 16 2005
Numbers k such that Sum_{d|k} (d/phi(d))*mu(k/d) = 0. - Benoit Cloitre, Apr 28 2002
Also, k with at least one x < k such that A007913(x) = A007913(k). - Benoit Cloitre, Apr 28 2002
Numbers k for which there exists a partition into two parts p and q such that p + q = k and p*q is a multiple of k. - Amarnath Murthy, May 30 2003
Numbers k such that there is a solution 0 < x < k to x^2 == 0 (mod k). - Franz Vrabec, Aug 13 2005
Numbers k such that moebius(k) = 0.
a(n) = k such that phi(k)/k = phi(m)/m for some m < k. - Artur Jasinski, Nov 05 2008
Appears to be numbers such that when a column with index equal to a(n) in A051731 is deleted, there is no impact on the result in the first column of A054525. - Mats Granvik, Feb 06 2009
Numbers k such that the number of prime divisors of (k+1) is less than the number of nonprime divisors of (k+1). - Juri-Stepan Gerasimov, Nov 10 2009
Orders for which at least one non-cyclic finite abelian group exists: A000688(a(n)) > 1. This follows from the fact that not all exponents in the prime factorization of a(n) are 1 (moebius(a(n)) = 0). The number of such groups of order a(n) is A192005(n) = A000688(a(n)) - 1. - Wolfdieter Lang, Jul 29 2011
Subsequence of A193166; A192280(a(n)) = 0. - Reinhard Zumkeller, Aug 26 2011
It appears that terms are the numbers m such that Product_{k=1..m} (prime(k) mod m) <> 0. See Maple code. - Gary Detlefs, Dec 07 2011
A008477(a(n)) > 1. - Reinhard Zumkeller, Feb 17 2012
A057918(a(n)) > 0. - Reinhard Zumkeller, Mar 27 2012
A056170(a(n)) > 0. - Reinhard Zumkeller, Dec 29 2012
Numbers k such that A001221(k) != A001222(k). - Felix Fröhlich, Aug 13 2014
Numbers k such that A001222(k) > A001221(k), since in this case at least one prime factor of k occurs more than once, which implies that k is divisible by at least one perfect square > 1. - Carlos Eduardo Olivieri, Aug 02 2015
Lexicographically least sequence such that each term has a positive even number of proper divisors not occurring in the sequence, cf. the sieve characterization of A005117. - Glen Whitney, Aug 30 2015
There are arbitrarily long runs of consecutive terms. Record runs start at 4, 8, 48, 242, ... (A045882). - Ivan Neretin, Nov 07 2015
A number k is a term if 0 < min(A000010(k) + A023900(k), A000010(k) - A023900(k)). - Torlach Rush, Feb 22 2018
Every squareful number > 1 is nonsquarefree, but the converse is false and the nonsquarefree numbers that are not squareful (see first comment) are in A332785. - Bernard Schott, Apr 11 2021
Integers m where at least one k < m exists such that m divides k^m. - Richard R. Forberg, Jul 31 2021
Consider the Diophantine equation S(x,y) = (x+y) + (x-y) + (x*y) + (x/y) = z, when x and y are both positive integers with y | x. Then, there is a solution (x,y) iff z is a term of this sequence; in this case, if x = K*y, then z = S(K*y,y) = K*(y+1)^2 (see A351381, link and references Perelman); example: S(12,4) = 75 = a(28). The number of solutions for S(x,y) = a(n) is A353282(n). - Bernard Schott, Mar 29 2022
For each positive integer m, the number of unitary divisors of m = the number of squarefree divisors of m (see A034444); but only for the terms of this sequence does the set of unitary divisors differ from the set of squarefree divisors. Example: the set of unitary divisors of 20 is {1, 4, 5, 20}, while the set of squarefree divisors of 20 is {1, 2, 5, 10}. - Bernard Schott, Oct 15 2022

Examples

			For the terms up to 20, we compute the squares of primes up to floor(sqrt(20)) = 4. Those squares are 4 and 9. For every such square s, put the terms s*k^2 for k = 1 to floor(20 / s). This gives after sorting and removing duplicates the list 4, 8, 9, 12, 16, 18, 20. - _David A. Corneth_, Oct 25 2017
		

References

  • I. Perelman, L'Algèbre récréative, Deux nombres et quatre opérations, Editions en langues étrangères, Moscou, 1959, pp. 101-102.
  • Ya. I. Perelman, Algebra can be fun, Two numbers and four operations, Mir Publishers Moscow, 1979, pp. 131-132.

Crossrefs

Complement of A005117. Subsequences: A130897, A190641, A332785.
Partitions into: A114374, A256012.

Programs

  • Haskell
    a013929 n = a013929_list !! (n-1)
    a013929_list = filter ((== 0) . a008966) [1..]
    -- Reinhard Zumkeller, Apr 22 2012
    
  • Magma
    [ n : n in [1..1000] | not IsSquarefree(n) ];
    
  • Maple
    a := n -> `if`(numtheory[mobius](n)=0,n,NULL); seq(a(i),i=1..160); # Peter Luschny, May 04 2009
    t:= n-> product(ithprime(k),k=1..n): for n from 1 to 160 do (if t(n) mod n <>0) then print(n) fi od; # Gary Detlefs, Dec 07 2011
    with(NumberTheory): isQuadrateful := n -> irem(Radical(n), n) <> 0:
    select(isQuadrateful, [`$`(1..160)]);  # Peter Luschny, Jul 12 2022
  • Mathematica
    Union[ Flatten[ Table[ n i^2, {i, 2, 20}, {n, 1, 400/i^2} ] ] ]
    Select[ Range[2, 160], (Union[Last /@ FactorInteger[ # ]][[ -1]] > 1) == True &] (* Robert G. Wilson v, Oct 11 2005 *)
    Cases[Range[160], n_ /; !SquareFreeQ[n]] (* Jean-François Alcover, Mar 21 2011 *)
    Select[Range@160, ! SquareFreeQ[#] &] (* Robert G. Wilson v, Jul 21 2012 *)
    Select[Range@160, PrimeOmega[#] > PrimeNu[#] &] (* Carlos Eduardo Olivieri, Aug 02 2015 *)
    Select[Range[200], MoebiusMu[#] == 0 &] (* Alonso del Arte, Nov 07 2015 *)
  • PARI
    {a(n)= local(m,c); if(n<=1,4*(n==1), c=1; m=4; while( cMichael Somos, Apr 29 2005 */
    
  • PARI
    for(n=1, 1e3, if(omega(n)!=bigomega(n), print1(n, ", "))) \\ Felix Fröhlich, Aug 13 2014
    
  • PARI
    upto(n)=my(res = List()); forprime(p = 2, sqrtint(n), for(k = 1, n \ p^2, listput(res, k * p^2))); listsort(res, 1); res \\ David A. Corneth, Oct 25 2017
    
  • Python
    from sympy.ntheory.factor_ import core
    def ok(n): return core(n, 2) != n
    print(list(filter(ok, range(1, 161)))) # Michael S. Branicky, Apr 08 2021
    
  • Python
    from math import isqrt
    from sympy import mobius
    def A013929(n):
        def f(x): return n+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 20 2024

Formula

A008966(a(n)) = 0. - Reinhard Zumkeller, Apr 22 2012
Sum_{n>=1} 1/a(n)^s = (zeta(s)*(zeta(2*s)-1))/zeta(2*s). - Enrique Pérez Herrero, Jul 07 2012
a(n) ~ n/k, where k = 1 - 1/zeta(2) = 1 - 6/Pi^2 = A229099. - Charles R Greathouse IV, Sep 13 2013
A001222(a(n)) > A001221(a(n)). - Carlos Eduardo Olivieri, Aug 02 2015
phi(a(n)) > A003958(a(n)). - Juri-Stepan Gerasimov, Apr 09 2019

Extensions

More terms from Erich Friedman
More terms from Franz Vrabec, Aug 13 2005

A053797 Lengths of successive gaps between squarefree numbers.

Original entry on oeis.org

1, 2, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 3, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 2, 3, 1, 1, 2, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 2, 2, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 3, 1, 3, 1, 2, 2, 2, 1
Offset: 1

Views

Author

N. J. A. Sloane, Apr 07 2000

Keywords

Comments

From Gus Wiseman, Jun 11 2024: (Start)
Also the length of the n-th maximal run of nonsquarefree numbers. These runs begin:
4
8 9
12
16
18
20
24 25
27 28
32
36
40
44 45
48 49 50
(End)

Examples

			The first gap is at 4 and has length 1; the next starts at 8 and has length 2 (since neither 8 nor 9 are squarefree).
		

Crossrefs

Gaps between terms of A005117.
For squarefree runs we have A120992, antiruns A373127 (firsts A373128).
For composite runs we have A176246 (rest of A046933), antiruns A373403.
For prime runs we have A251092 (rest of A175632), antiruns A027833.
Position of first appearance of n is A373199(n).
For antiruns instead of runs we have A373409.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.

Programs

  • Maple
    SF:= select(numtheory:-issqrfree,[$1..1000]):
    map(`-`,select(`>`,SF[2..-1]-SF[1..-2],1),1); # Robert Israel, Sep 22 2015
  • Mathematica
    ReplaceAll[Differences[Select[Range@384, SquareFreeQ]] - 1, 0 -> Nothing] (* Michael De Vlieger, Sep 22 2015 *)

Extensions

Offset set to 1 by Peter Kagey, Sep 29 2015

A053806 Numbers where a gap begins in the sequence of squarefree numbers (A005117).

Original entry on oeis.org

4, 8, 12, 16, 18, 20, 24, 27, 32, 36, 40, 44, 48, 52, 54, 56, 60, 63, 68, 72, 75, 80, 84, 88, 90, 92, 96, 98, 104, 108, 112, 116, 120, 124, 128, 132, 135, 140, 144, 147, 150, 152, 156, 160, 162, 164, 168, 171, 175, 180, 184, 188, 192, 196, 198, 200, 204, 207, 212
Offset: 1

Views

Author

N. J. A. Sloane, Apr 07 2000

Keywords

Examples

			The first gap is at 4 and has length 1; the next starts at 8 and has length 2 (since neither 8 nor 9 are squarefree).
		

Crossrefs

Programs

  • PARI
    is(n)=!issquarefree(n) && issquarefree(n-1) \\ Charles R Greathouse IV, Nov 05 2017
    
  • PARI
    list(lim)=my(v=List(),t); forfactored(n=4,lim\1, if(vecmax(n[2][,2])>1, if(!t, listput(v,n[1])); t=1, t=0)); Vec(v) \\ Charles R Greathouse IV, Nov 05 2017

A112925 Largest squarefree integer < the n-th prime.

Original entry on oeis.org

1, 2, 3, 6, 10, 11, 15, 17, 22, 26, 30, 35, 39, 42, 46, 51, 58, 59, 66, 70, 71, 78, 82, 87, 95, 97, 102, 106, 107, 111, 123, 130, 134, 138, 146, 149, 155, 161, 166, 170, 178, 179, 190, 191, 195, 197, 210, 222, 226, 227, 231, 238, 239, 249, 255, 262, 267, 269, 274, 278
Offset: 1

Views

Author

Leroy Quet, Oct 06 2005

Keywords

Examples

			6 is the largest squarefree less than the 4th prime, 7. So a(4) = 6.
		

Crossrefs

For prime powers instead of squarefree numbers we have A065514, opposite A345531.
Restriction of A070321 (differences A378085) to the primes; see A378619.
The opposite is A112926, differences A378037.
Subtracting each term from prime(n) gives A240473, opposite A240474.
For nonsquarefree numbers we have A378033, differences A378036, see A378034, A378032.
For perfect powers we have A378035.
First differences are A378038.
A000040 lists the primes, differences A001223, seconds A036263.
A005117 lists the squarefree numbers, differences A076259.
A013928 counts squarefree numbers up to n - 1.
A013929 lists the nonsquarefree numbers, differences A078147.
A061398 counts squarefree numbers between primes, zeros A068360.
A061399 counts nonsquarefree numbers between primes, zeros A068361.
A112929 counts squarefree numbers up to prime(n).

Programs

  • Maple
    with(numtheory): a:=proc(n) local p,B,j: p:=ithprime(n): B:={}: for j from 1 to p-1 do if abs(mobius(j))>0 then B:=B union {j} else B:=B fi od: B[nops(B)] end: seq(a(m),m=1..75); # Emeric Deutsch, Oct 14 2005
  • Mathematica
    With[{k = 120}, Table[SelectFirst[Range[Prime@ n - 1, Prime@ n - Min[Prime@ n - 1, k], -1], SquareFreeQ], {n, 60}]] (* Michael De Vlieger, Aug 16 2017 *)
  • PARI
    a(n,p=prime(n))=while(!issquarefree(p--),); p \\ Charles R Greathouse IV, Aug 16 2017

Formula

a(n) = prime(n) - A240473(n). - Gus Wiseman, Jan 10 2025

Extensions

More terms from Emeric Deutsch, Oct 14 2005

A068781 Lesser of two consecutive numbers each divisible by a square.

Original entry on oeis.org

8, 24, 27, 44, 48, 49, 63, 75, 80, 98, 99, 116, 120, 124, 125, 135, 147, 152, 168, 171, 175, 188, 207, 224, 242, 243, 244, 260, 275, 279, 288, 296, 315, 324, 332, 342, 343, 350, 351, 360, 363, 368, 375, 387, 404, 423, 424, 440, 459, 475, 476, 495, 507, 512
Offset: 1

Views

Author

Robert G. Wilson v, Mar 04 2002

Keywords

Comments

Also numbers m such that mu(m)=mu(m+1)=0, where mu is the Moebius-function (A008683); A081221(a(n))>1. - Reinhard Zumkeller, Mar 10 2003
The sequence contains an infinite family of arithmetic progressions like {36a+8}={8,44,80,116,152,188,...} ={4(9a+2)}. {36a+9} provides 2nd nonsquarefree terms. Such AP's can be constructed to any term by solution of a system of linear Diophantine equation. - Labos Elemer, Nov 25 2002
1. 4k^2 + 4k is a member for all k; i.e., 8 times a triangular number is a member. 2. (4k+1) times an odd square - 1 is a member. 3. (4k+3) times odd square is a member. - Amarnath Murthy, Apr 24 2003
The asymptotic density of this sequence is 1 - 2/zeta(2) + Product_{p prime} (1 - 2/p^2) = 1 - 2 * A059956 + A065474 = 0.1067798952... (Matomäki et al., 2016). - Amiram Eldar, Feb 14 2021
Maximum of the n-th maximal anti-run of nonsquarefree numbers (A013929) differing by more than one. For runs instead of anti-runs we have A376164. For squarefree instead of nonsquarefree we have A007674. - Gus Wiseman, Sep 14 2024

Examples

			44 is in the sequence because 44 = 2^2 * 11 and 45 = 3^2 * 5.
From _Gus Wiseman_, Sep 14 2024: (Start)
Splitting nonsquarefree numbers into maximal anti-runs gives:
  (4,8)
  (9,12,16,18,20,24)
  (25,27)
  (28,32,36,40,44)
  (45,48)
  (49)
  (50,52,54,56,60,63)
  (64,68,72,75)
  (76,80)
  (81,84,88,90,92,96,98)
  (99)
The maxima are a(n). The corresponding pairs are (8,9), (24,25), (27,28), (44,45), etc.
(End)
		

Crossrefs

Subsequence of A261869.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.
A053797 gives lengths of runs of nonsquarefree numbers, firsts A373199.

Programs

  • Haskell
    a068781 n = a068781_list !! (n-1)
    a068781_list = filter ((== 0) . a261869) [1..]
    -- Reinhard Zumkeller, Sep 04 2015
    
  • Mathematica
    Select[ Range[2, 600], Max[ Transpose[ FactorInteger[ # ]] [[2]]] > 1 && Max[ Transpose[ FactorInteger[ # + 1]] [[2]]] > 1 &]
    f@n_:= Flatten@Position[Partition[SquareFreeQ/@Range@2000,n,1], Table[False,{n}]]; f@2 (* Hans Rudolf Widmer, Aug 30 2022 *)
    Max/@Split[Select[Range[100], !SquareFreeQ[#]&],#1+1!=#2&]//Most (* Gus Wiseman, Sep 14 2024 *)
  • PARI
    isok(m) = !moebius(m) && !moebius(m+1); \\ Michel Marcus, Feb 14 2021

Formula

A261869(a(n)) = 0. - Reinhard Zumkeller, Sep 04 2015

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

Views

Author

Labos Elemer, Jun 07 2001

Keywords

Examples

			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)
		

Crossrefs

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.

Programs

A378037 First differences of A112926 (smallest squarefree integer > prime(n)).

Original entry on oeis.org

2, 1, 4, 3, 1, 5, 2, 5, 4, 3, 5, 4, 4, 5, 4, 6, 1, 7, 4, 1, 8, 3, 6, 10, 1, 3, 4, 1, 4, 15, 4, 5, 3, 10, 3, 4, 7, 5, 4, 7, 1, 11, 1, 5, 2, 12, 13, 3, 1, 5, 6, 5, 7, 5, 7, 6, 2, 5, 4, 3, 10, 14, 4, 1, 4, 16, 5, 10, 4, 1, 8, 8, 4, 7, 4, 5, 8, 4, 8, 11, 1, 11, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 04 2024

Keywords

Crossrefs

First differences of A112926, restriction of A067535, differences A378087.
For prime powers we have A377703.
The nonsquarefree version is A377784 (differences of A377783), restriction of A378039.
The nonsquarefree opposite is A378034, first differences of A378032.
The opposite is A378038, differences of A112925.
The unrestricted opposite is A378085 except first term, differences of A070321.
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.

Programs

  • Mathematica
    Differences[Table[NestWhile[#+1&,Prime[n]+1,!SquareFreeQ[#]&],{n,100}]]

A378038 First differences of A112925 = greatest squarefree number < prime(n).

Original entry on oeis.org

1, 1, 3, 4, 1, 4, 2, 5, 4, 4, 5, 4, 3, 4, 5, 7, 1, 7, 4, 1, 7, 4, 5, 8, 2, 5, 4, 1, 4, 12, 7, 4, 4, 8, 3, 6, 6, 5, 4, 8, 1, 11, 1, 4, 2, 13, 12, 4, 1, 4, 7, 1, 10, 6, 7, 5, 2, 5, 4, 4, 9, 14, 5, 1, 3, 16, 5, 11, 1, 2, 9, 8, 5, 6, 5, 4, 9, 4, 8, 11, 1, 11, 1, 7
Offset: 1

Views

Author

Gus Wiseman, Dec 04 2024

Keywords

Crossrefs

First differences of A112925, restriction of A070321, differences A378085.
For prime powers we have A377781, opposite A377703.
The nonsquarefree opposite is A377784 (differences of A377783), restriction of A378039.
The nonsquarefree version is A378034, first differences of A378032.
The opposite is A378037, differences of A112926.
The unrestricted opposite is A378087, differences of A067535.
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.

Programs

  • Mathematica
    Differences[Table[NestWhile[#-1&,Prime[n]-1,!SquareFreeQ[#]&],{n,100}]]

A373199 Least k such that the k-th maximal run of nonsquarefree numbers has length n. Position of first appearance of n in A053797.

Original entry on oeis.org

1, 2, 13, 68, 241, 6278, 61921, 311759, 2530539
Offset: 1

Views

Author

Gus Wiseman, Jun 08 2024

Keywords

Comments

A run of a sequence (in this case A013929) is an interval of positions at which consecutive terms differ by one. The a(n)-th run of nonsquarefree numbers begins with A045882 = A051681, subset of A053806.

Examples

			The maximal runs of nonsquarefree numbers begin:
   4
   8   9
  12
  16
  18
  20
  24  25
  27  28
  32
  36
  40
  44  45
  48  49  50
  52
  54
  56
  60
  63  64
The a(n)-th rows are:
     4
     8     9
    48    49    50
   242   243   244   245
   844   845   846   847   848
For example, (48, 49, 50) is the first maximal run of 3 nonsquarefree numbers, so a(3) = 13.
		

Crossrefs

For composite instead of nonsquarefree we have A073051.
The version for squarefree runs is A373128.
For prime instead of nonsquarefree we have A373400.
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.

Programs

  • Mathematica
    seq=Length/@Split[Select[Range[10000],!SquareFreeQ[#]&],#1+1==#2&];
    spna[y_]:=Max@@Select[Range[Length[y]],SubsetQ[y,Range[#]]&];
    Table[Position[seq,i][[1,1]],{i,spna[seq]}]
Showing 1-10 of 48 results. Next