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 64 results. Next

A078147 First differences of sequence of nonsquarefree numbers, A013929.

Original entry on oeis.org

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

Views

Author

Labos Elemer, Nov 26 2002

Keywords

Comments

Run lengths in A132345, apart from initial run of zeros. - Reinhard Zumkeller, Apr 22 2012
The asymptotic density of the occurrences of 1 in this sequence is density(A068781)/density(A013929) = (1 - 2 * A059956 + A065474)/A229099 = 0.272347... - Amiram Eldar, Mar 09 2021

Examples

			a(1) = 4 = 8 - 4.
		

Crossrefs

Programs

  • Haskell
    a078147 n = a078147_list !! (n-1)
    a078147_list = zipWith (-) (tail a013929_list) a013929_list
    -- Reinhard Zumkeller, Apr 22 2012
    
  • Mathematica
    t=Flatten[Position[Table[MoebiusMu[w], {w, 1, 1000}], 0]]; t1=Delete[RotateLeft[t]-t, -1]
    Differences[Select[Range[300],!SquareFreeQ[#]&]] (* Harvey P. Dale, May 07 2012 *)
  • PARI
    lista(nn) = {my(prec=0); for (n=1, nn, if (!issquarefree(n), if (prec, print1(n-prec, ", ")); prec = n;););} \\ Michel Marcus, Mar 26 2020
    
  • Python
    from math import isqrt
    from sympy import mobius, factorint
    def A078147(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 next(i for i in range(1,5) if any(d>1 for d in factorint(m+i).values())) # Chai Wah Wu, Sep 10 2024

Formula

a(n) = A013929(n+1) - A013929(n).
a(n) = 1, 2, 3 or 4 since n = 4*k is always nonsquarefree.
Asymptotic mean: lim_{n->oo} (1/n) Sum_{k=1..n} a(k) = Pi^2/(Pi^2-6) = 2.550546... - Amiram Eldar, Oct 21 2020

Extensions

Offset fixed by Reinhard Zumkeller, Apr 22 2012

A045882 Smallest term of first run of (at least) n consecutive integers which are not squarefree.

Original entry on oeis.org

4, 8, 48, 242, 844, 22020, 217070, 1092747, 8870024, 221167422, 221167422, 47255689915, 82462576220, 1043460553364, 79180770078548, 3215226335143218, 23742453640900972, 125781000834058568
Offset: 1

Views

Author

Keywords

Comments

Solution for n=10 is same as for n=11.
This sequence is infinite and each term initiates a suitable arithmetic progression with large differences like squares of primorials or other suitable products of primes from prime factors being on power 2 in terms and in chains after. Proof includes solution of linear Diophantine equations and math. induction. See also A068781, A070258, A070284, A078144, A049535, A077640, A077647, A078143 of which first terms are recollected here. - Labos Elemer, Nov 25 2002

Examples

			a(3) = 48 as 48, 49 and 50 are divisible by squares.
n=5 -> {844=2^2*211; 845=5*13^2; 846=2*3^2*47; 847=7*11^2; 848=2^4*53}.
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 242, p. 67, Ellipses, Paris 2008.

Crossrefs

Cf. A013929, A053806, A049535, A077647, A078143. Also A069021 and A051681 are different versions.

Programs

  • Mathematica
    cnt = 0; k = 0; Table[While[cnt < n, k++; If[! SquareFreeQ[k], cnt++, cnt = 0]]; k - n + 1, {n, 7}]
  • PARI
    a(n)=my(s);for(k=1,9^99,if(issquarefree(k),s=0,if(s++==n,return(k-n+1)))) \\ Charles R Greathouse IV, May 29 2013

Formula

a(n) = 1 + A020754(n+1). - R. J. Mathar, Jun 25 2010
Correction from Jeppe Stig Nielsen, Mar 05 2022: (Start)
a(n) = 1 + A020754(n+1) for 1 <= n < 11.
a(n) = 1 + A020754(n) for 11 <= n < N where N is unknown.
Possibly a(n) = 1 + A020754(n-d) for some higher n, depending on how many repeated terms the sequence has. (End)
a(n) <= A061742(n) = A002110(n)^2 is the trivial bound obtained from the CRT. - Charles R Greathouse IV, Sep 06 2022

Extensions

a(9)-a(11) from Patrick De Geest, Nov 15 1998, Jan 15 1999
a(12)-a(15) from Louis Marmet (louis(AT)marmet.org) and David Bernier (ezcos(AT)yahoo.com), Nov 15 1999
a(16) was obtained as a result of a team effort by Z. McGregor-Dorsey et al. [Louis Marmet (louis(AT)marmet.org), Jul 27 2000]
a(17) was obtained as a result of a team effort by E. Wong et al. [Louis Marmet (louis(AT)marmet.org), Jul 13 2001]
a(18) was obtained as a result of a team effort by L. Marmet et al.

A120327 Smallest nonsquarefree number >= n.

Original entry on oeis.org

4, 4, 4, 4, 8, 8, 8, 8, 9, 12, 12, 12, 16, 16, 16, 16, 18, 18, 20, 20, 24, 24, 24, 24, 25, 27, 27, 28, 32, 32, 32, 32, 36, 36, 36, 36, 40, 40, 40, 40, 44, 44, 44, 44, 45, 48, 48, 48, 49, 50, 52, 52, 54, 54, 56, 56, 60, 60, 60, 60, 63, 63, 63, 64, 68, 68, 68, 68, 72, 72, 72, 72
Offset: 1

Views

Author

Zak Seidov, Aug 16 2006

Keywords

Crossrefs

For squarefree instead of nonsquarefree we have A067535, differences A378087.
The opposite for squarefree is A070321, differences A378085.
The run-lengths are A078147 if we prepend 4, differences A376593.
The restriction to primes is A377783 (union A378040), differences A377784.
The opposite is A378033 (differences A378036), for prime powers A031218.
First differences are A378039 if we assume that a(1) = 1.
A005117 lists the squarefree numbers.
A013929 lists the nonsquarefree numbers.
A061398 counts squarefree numbers between primes, zeros A068360.
A061399 counts nonsquarefree numbers between primes, zeros A068361.

Programs

  • Mathematica
    Table[NestWhile[ #+1&,n,SquareFreeQ],{n,100}] (* simplified by Harvey P. Dale, Apr 08 2014 *)

A373673 First element of each maximal run of powers of primes (including 1).

Original entry on oeis.org

1, 7, 11, 13, 16, 19, 23, 25, 27, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 125, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239
Offset: 1

Views

Author

Gus Wiseman, Jun 15 2024

Keywords

Comments

A run of a sequence (in this case A000961) is an interval of positions at which consecutive terms differ by one.
The last element of the same run is A373674.
Consists of all powers of primes k such that k-1 is not a power of primes.

Examples

			The maximal runs of powers of primes begin:
   1   2   3   4   5
   7   8   9
  11
  13
  16  17
  19
  23
  25
  27
  29
  31  32
  37
  41
  43
  47
  49
		

Crossrefs

For composite antiruns we have A005381, max A068780, length A373403.
For prime antiruns we have A006512, max A001359, length A027833.
For composite runs we have A008864, max A006093, length A176246.
For prime runs we have A025584, max A067774, length A251092 or A175632.
For runs of prime-powers:
- length A174965
- min A373673 (this sequence)
- max A373674
- sum A373675
For runs of non-prime-powers:
- length A110969 (firsts A373669, sorted A373670)
- min A373676
- max A373677
- sum A373678
For antiruns of prime-powers:
- length A373671
- min A120430
- max A006549
- sum A373576
For antiruns of non-prime-powers:
- length A373672
- min A373575
- max A255346
- sum A373679
A000961 lists all powers of primes (A246655 if not including 1).
A025528 counts prime-powers up to n.
A057820 gives first differences of consecutive prime-powers, gaps A093555.
A361102 lists all non-prime-powers (A024619 if not including 1).

Programs

  • Mathematica
    pripow[n_]:=n==1||PrimePowerQ[n];
    Min/@Split[Select[Range[100],pripow],#1+1==#2&]//Most

A373676 First element of each maximal run of non-prime-powers.

Original entry on oeis.org

1, 6, 10, 12, 14, 18, 20, 24, 26, 28, 30, 33, 38, 42, 44, 48, 50, 54, 60, 62, 65, 68, 72, 74, 80, 82, 84, 90, 98, 102, 104, 108, 110, 114, 122, 126, 129, 132, 138, 140, 150, 152, 158, 164, 168, 170, 174, 180, 182, 192, 194, 198, 200, 212, 224, 228, 230, 234
Offset: 1

Views

Author

Gus Wiseman, Jun 16 2024

Keywords

Comments

We consider 1 to be a power of a prime and a non-prime-power, but not a prime-power.
A run of a sequence (in this case A000961) is an interval of positions at which consecutive terms differ by one.
The last element of the same run is A373677.
Consists of 1 and all non-prime-powers k such that k-1 is a power of a prime.

Examples

			The maximal runs of non-prime-powers begin:
   1
   6
  10
  12
  14  15
  18
  20  21  22
  24
  26
  28
  30
  33  34  35  36
  38  39  40
  42
  44  45  46
  48
  50  51  52
  54  55  56  57  58
  60
		

Crossrefs

See link for prime, composite, squarefree, and nonsquarefree runs/antiruns.
For runs of powers of primes:
- length A174965
- min A373673
- max A373674
- sum A373675
For runs of non-prime-powers:
- length A110969 (firsts A373669, sorted A373670)
- min A373676 (this sequence)
- max A373677
- sum A373678
For antiruns of prime-powers:
- length A373671
- min A120430
- max A006549
- sum A373576
For antiruns of non-prime-powers:
- length A373672
- min A373575
- max A255346
- sum A373679
A000961 lists all powers of primes. A246655 is just prime-powers so lacks 1.
A025528 counts prime-powers up to n.
A057820 gives first differences of consecutive prime-powers, gaps A093555.
A361102 lists all non-prime-powers (A024619 if not including 1).

Programs

  • Mathematica
    Select[Range[100],#==1||!PrimePowerQ[#]&&PrimePowerQ[#-1]&]

A373677 Last element of each maximal run of non-prime-powers.

Original entry on oeis.org

1, 6, 10, 12, 15, 18, 22, 24, 26, 28, 30, 36, 40, 42, 46, 48, 52, 58, 60, 63, 66, 70, 72, 78, 80, 82, 88, 96, 100, 102, 106, 108, 112, 120, 124, 126, 130, 136, 138, 148, 150, 156, 162, 166, 168, 172, 178, 180, 190, 192, 196, 198, 210, 222, 226, 228, 232, 238
Offset: 1

Views

Author

Gus Wiseman, Jun 16 2024

Keywords

Comments

We consider 1 to be a power of a prime and a non-prime-power, but not a prime-power.
A run of a sequence (in this case A000961) is an interval of positions at which consecutive terms differ by one.
The first element of the same run is A373676.
Consists of all non-prime-powers k such that k+1 is a prime-power.

Examples

			The maximal runs of non-prime-powers begin:
   1
   6
  10
  12
  14  15
  18
  20  21  22
  24
  26
  28
  30
  33  34  35  36
  38  39  40
  42
  44  45  46
  48
  50  51  52
  54  55  56  57  58
  60
		

Crossrefs

See link for prime, composite, squarefree, and nonsquarefree runs/antiruns.
For runs of powers of primes:
- length A174965
- min A373673
- max A373674
- sum A373675
For runs of non-prime-powers:
- length A110969 (firsts A373669, sorted A373670)
- min A373676
- max A373677 (this sequence)
- sum A373678
For antiruns of prime-powers:
- length A373671
- min A120430
- max A006549
- sum A373576
For antiruns of non-prime-powers:
- length A373672
- min A373575
- max A255346
- sum A373679
A000961 lists all powers of primes. A246655 is just prime-powers so lacks 1.
A025528 counts prime-powers up to n.
A057820 gives first differences of consecutive prime-powers, gaps A093555.
A361102 lists all non-prime-powers (A024619 if not including 1).

Programs

  • Mathematica
    Select[Range[100],!PrimePowerQ[#]&&PrimePowerQ[#+1]&]

A375707 First differences minus 1 of nonsquarefree numbers.

Original entry on oeis.org

3, 0, 2, 3, 1, 1, 3, 0, 1, 0, 3, 3, 3, 3, 0, 2, 0, 0, 1, 1, 1, 3, 2, 0, 3, 3, 2, 0, 3, 0, 2, 3, 1, 1, 3, 1, 0, 0, 3, 3, 3, 3, 0, 2, 0, 2, 0, 0, 1, 3, 2, 0, 3, 3, 2, 0, 1, 1, 0, 2, 3, 1, 1, 3, 0, 1, 0, 2, 0, 3, 3, 3, 0, 2, 3, 1, 1, 3, 2, 0, 3, 3, 3, 3, 0, 2, 3
Offset: 1

Views

Author

Gus Wiseman, Sep 16 2024

Keywords

Comments

Also the number of squarefree numbers between the nonsquarefree numbers A013929(n) and A013929(n+1).
Delete all 0's to get A120992.
The image is {0,1,2,3}.
Add 1 to all terms for A078147.

Examples

			The runs of squarefree numbers begin:
  (5,6,7)
  ()
  (10,11)
  (13,14,15)
  (17)
  (19)
  (21,22,23)
  ()
  (26)
  ()
  (29,30,31)
  (33,34,35)
		

Crossrefs

Positions of 0, 1, 2, 3 are A375709, A375710, A375711, A375712. This is a set partition of the positive integers into four blocks.
For runs of squarefree numbers:
- length: A120992, anti A373127
- min: A072284, anti A373408
- max: A373415, anti A007674
- sum: A373413, anti A373411
For runs of nonsquarefree numbers:
- length: A053797, anti A373409
- min: A053806, anti A373410
- max: A376164, anti A068781
- sum: A373414, anti A373412
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.
A046933 counts composite numbers between consecutive primes.
A073784 counts primes between consecutive composite numbers.
A093555 counts non-prime-powers between consecutive prime-powers.

Programs

  • Mathematica
    Differences[Select[Range[100],!SquareFreeQ[#]&]]-1
  • PARI
    lista(nmax) = {my(prev = 4); for (n = 5, nmax, if(!issquarefree(n), print1(n - prev - 1, ", "); prev = n));} \\ Amiram Eldar, Sep 17 2024

Formula

Asymptotic mean: lim_{n->oo} (1/n) Sum_{k=1..n} a(k) = 6/(Pi^2-6) = 1.550546... . - Amiram Eldar, Sep 17 2024

A373674 Last element of each maximal run of powers of primes (including 1).

Original entry on oeis.org

5, 9, 11, 13, 17, 19, 23, 25, 27, 29, 32, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 125, 128, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239
Offset: 1

Views

Author

Gus Wiseman, Jun 16 2024

Keywords

Comments

A run of a sequence (in this case A000961) is an interval of positions at which consecutive terms differ by one.
The first element of the same run is A373673.
Consists of all powers of primes k such that k+1 is not a power of primes.

Examples

			The maximal runs of powers of primes begin:
   1   2   3   4   5
   7   8   9
  11
  13
  16  17
  19
  23
  25
  27
  29
  31  32
  37
  41
  43
  47
  49
		

Crossrefs

For prime antiruns we have A001359, min A006512, length A027833.
For composite runs we have A006093, min A008864, length A176246.
For prime runs we have A067774, min A025584, length A251092 or A175632.
For squarefree runs we have A373415, min A072284, length A120992.
For nonsquarefree runs we have min A053806, length A053797.
For runs of prime-powers:
- length A174965
- min A373673
- max A373674 (this sequence)
- sum A373675
For runs of non-prime-powers:
- length A110969 (firsts A373669, sorted A373670)
- min A373676
- max A373677
- sum A373678
For antiruns of prime-powers:
- length A373671
- min A120430
- max A006549
- sum A373576
For antiruns of non-prime-powers:
- length A373672
- min A373575
- max A255346
- sum A373679
A000961 lists all powers of primes (A246655 if not including 1).
A025528 counts prime-powers up to n.
A057820 gives first differences of consecutive prime-powers, gaps A093555.
A361102 lists all non-prime-powers (A024619 if not including 1).

Programs

  • Mathematica
    pripow[n_]:=n==1||PrimePowerQ[n];
    Max/@Split[Select[Range[nn],pripow],#1+1==#2&]//Most

A068140 Smaller of two consecutive numbers each divisible by a cube greater than one.

Original entry on oeis.org

80, 135, 296, 343, 351, 375, 512, 567, 624, 728, 783, 944, 999, 1160, 1215, 1375, 1376, 1431, 1592, 1624, 1647, 1808, 1863, 2024, 2079, 2240, 2295, 2375, 2400, 2456, 2511, 2624, 2672, 2727, 2888, 2943, 3087, 3104, 3159, 3320, 3375, 3429, 3536, 3591
Offset: 1

Views

Author

Amarnath Murthy, Feb 22 2002

Keywords

Comments

Cubeful numbers with cubeful successors. This is to cubes as A068781 is to squares. 1375 is the smallest of three consecutive numbers divisible by a cube, since 1375 = 5^3 * 11 and 1376 = 2^5 * 43 and 1377 = 3^4 * 17. What is the smallest of four consecutive numbers divisible by a cube? Of n consecutive numbers divisible by a cube? - Jonathan Vos Post, Sep 18 2007
22624 is the smallest of four consecutive numbers each divisible by a cube, with factorizations 2^5 * 7 * 101, 5^3 * 181, 2 * 3^3 * 419, and 11^3 * 17. - D. S. McNeil, Dec 10 2010
18035622 is the smallest of five consecutive numbers each divisible by a cube. 4379776620 is the smallest of six consecutive numbers each divisible by a cube. 1204244328624 is the smallest of seven consecutive numbers each divisible by a cube. - Donovan Johnson, Dec 13 2010
The sequence is the union, over all pairs of distinct primes (p,q), of numbers == 0 mod p^3 and == -1 mod q^3 or vice versa. - Robert Israel, Aug 13 2018
The asymptotic density of this sequence is 1 - 2/zeta(3) + Product_{p prime} (1 - 2/p^3) = 1 - 2 * A088453 + A340153 = 0.013077991848467056243... - Amiram Eldar, Feb 16 2021

Examples

			343 is a term as 343 = 7^3 and 344= 2^3 * 43.
		

Crossrefs

Programs

  • Maple
    isA068140 := proc(n)
        isA046099(n) and isA046099(n+1) ;
    end proc:
    for n from 1 to 4000 do
        if isA068140(n) then
            printf("%d,",n) ;
        end if;
    end do: # R. J. Mathar, Dec 08 2015
  • Mathematica
    a = b = 0; Do[b = Max[ Transpose[ FactorInteger[n]] [[2]]]; If[a > 2 && b > 2, Print[n - 1]]; a = b, {n, 2, 5000}]
    Select[Range[2, 6000], Max[Transpose[FactorInteger[ # ]][[2]]] > 2 && Max[Transpose[FactorInteger[ # + 1]][[2]]] > 2 &] (* Jonathan Vos Post, Sep 18 2007 *)
    SequencePosition[Table[If[AnyTrue[Rest[Divisors[n]],IntegerQ[Surd[#,3]]&],1,0],{n,3600}],{1,1}][[All,1]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Apr 18 2020 *)

Formula

{k such that k is in A046099 and k+1 is in A046099}. - Jonathan Vos Post, Sep 18 2007

Extensions

Edited and extended by Robert G. Wilson v, Mar 02 2002
Title edited, cross-references added by Matthew Vandermast, Dec 09 2010
Definition clarified by Harvey P. Dale, Apr 18 2020

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 64 results. Next