A078147 -id:A078147 - cp's OEIS frontend

A376593 Second differences of consecutive nonsquarefree numbers (A013929). First differences of A078147.

-3, 2, 1, -2, 0, 2, -3, 1, -1, 3, 0, 0, 0, -3, 2, -2, 0, 1, 0, 0, 2, -1, -2, 3, 0, -1, -2, 3, -3, 2, 1, -2, 0, 2, -2, -1, 0, 3, 0, 0, 0, -3, 2, -2, 2, -2, 0, 1, 2, -1, -2, 3, 0, -1, -2, 1, 0, -1, 2, 1, -2, 0, 2, -3, 1, -1, 2, -2, 3, 0, 0, -3, 2, 1, -2, 0, 2

Offset: 1

Gus Wiseman, Oct 01 2024

sign

The range is {-3, -2, -1, 0, 1, 2, 3}.

			The nonsquarefree numbers (A013929) are:
  4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 49, 50, 52, ...
with first differences (A078147):
  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, ...
with first differences (A376593):
  -3, 2, 1, -2, 0, 2, -3, 1, -1, 3, 0, 0, 0, -3, 2, -2, 0, 1, 0, 0, 2, -1, -2, ...

The version for A000002 is A376604, first differences of A054354.
The first differences were A078147.
Zeros are A376594, complement A376595.
A000040 lists the prime numbers, differences A001223.
A005117 lists squarefree numbers, differences A076259.
A064113 lists positions of adjacent equal prime gaps.
A114374 counts partitions into nonsquarefree numbers.
A246655 lists prime-powers exclusive, inclusive A000961.
A333254 lists run-lengths of differences between consecutive primes.
For second differences: A036263 (prime), A073445 (composite), A376559 (perfect-power), A376562 (non-perfect-power), A376590 (squarefree), A376596 (prime-power inclusive), A376599 (non-prime-power inclusive).
For nonsquarefree numbers: A013929 (terms), A078147 (first differences), A376594 (inflections and undulations), A376595 (nonzero curvature).
Cf. A061398, A053797, A053806, A120992, A182853, A251092, A373198, A375707, A376305, A376306, A376311, A376312, A376655.

Differences[Select[Range[100],!SquareFreeQ[#]&],2]

from math import isqrt
from sympy import mobius, factorint
def A376593(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)
    k = next(i for i in range(1,5) if any(d>1 for d in factorint(m+i).values()))
    return next(i for i in range(1-k,5-k) if any(d>1 for d in factorint(m+(k<<1)+i).values())) # Chai Wah Wu, Oct 02 2024

A376312 Run-compression of first differences (A078147) of nonsquarefree numbers (A013929).

Original entry on oeis.org

4, 1, 3, 4, 2, 4, 1, 2, 1, 4, 1, 3, 1, 2, 4, 3, 1, 4, 3, 1, 4, 1, 3, 4, 2, 4, 2, 1, 4, 1, 3, 1, 3, 1, 2, 4, 3, 1, 4, 3, 1, 2, 1, 3, 4, 2, 4, 1, 2, 1, 3, 1, 4, 1, 3, 4, 2, 4, 3, 1, 4, 1, 3, 4, 2, 4, 2, 1, 3, 2, 4, 1, 3, 4, 2, 3, 1, 3, 1, 4, 1, 3, 2, 1, 3, 4, 2

Offset: 1

Gus Wiseman, Sep 24 2024

nonn

We define the run-compression of a sequence to be the anti-run obtained by reducing each run of repeated parts to a single part. Alternatively, we can remove all parts equal to the part immediately to their left. For example, (1,1,2,2,1) has run-compression (1,2,1).

			The sequence of nonsquarefree numbers (A013929) is:
  4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 49, 50, ...
with first differences (A078147):
  4, 1, 3, 4, 2, 2, 4, 1, 2, 1, 4, 4, 4, 4, 1, 3, 1, 1, 2, 2, 2, 4, 3, 1, ...
with runs:
  (4),(1),(3),(4),(2,2),(4),(1),(2),(1),(4,4,4,4),(1),(3),(1,1),(2,2,2), ...
and run-compression (A376312):
  4, 1, 3, 4, 2, 4, 1, 2, 1, 4, 1, 3, 1, 2, 4, 3, 1, 4, 3, 1, 4, 1, 3, 4, ...

For nonprime instead of squarefree numbers we have A037201, halved A373947.
Before compressing we had A078147.
For run-sums instead of compression we have A376264.
For squarefree instead of nonsquarefree we have A376305, ones A376342.
For prime-powers instead of nonsquarefree numbers we have A376308.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A003242 counts compressed compositions, ranks A333489.
A005117 lists squarefree numbers, differences A076259 (ones A375927).
A013929 lists nonsquarefree numbers, differences A078147.
A116861 counts partitions by compressed sum, by compressed length A116608.
Cf. A007674, A053797, A053806, A072284, A112925, A120992, A274174, A373198, A375707, A376306, A376307, A376311.

First/@Split[Differences[Select[Range[100], !SquareFreeQ[#]&]]]

A376264 Run-sums of first differences (A078147) of nonsquarefree numbers (A013929).

Original entry on oeis.org

4, 1, 3, 4, 4, 4, 1, 2, 1, 16, 1, 3, 2, 6, 4, 3, 1, 8, 3, 1, 4, 1, 3, 4, 4, 4, 2, 2, 16, 1, 3, 1, 3, 2, 2, 4, 3, 1, 8, 3, 1, 4, 1, 3, 4, 4, 4, 1, 2, 1, 3, 1, 12, 1, 3, 4, 4, 4, 3, 1, 16, 1, 3, 4, 4, 4, 2, 3, 3, 4, 8, 1, 3, 4, 4, 3, 1, 3, 1, 8, 1, 3, 4, 1, 3, 4

Offset: 1

Gus Wiseman, Sep 26 2024

nonn

Does the image include all positive integers? I have only verified this up to 8.

			The sequence of nonsquarefree numbers (A013929) is:
  4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 49, 50, ...
with first differences (A078147):
  4, 1, 3, 4, 2, 2, 4, 1, 2, 1, 4, 4, 4, 4, 1, 3, 1, 1, 2, 2, 2, 4, 3, 1, ...
with runs:
  (4),(1),(3),(4),(2,2),(4),(1),(2),(1),(4,4,4,4),(1),(3),(1,1),(2,2,2), ...
with sums (A376264):
  4, 1, 3, 4, 4, 4, 1, 2, 1, 16, 1, 3, 2, 6, 4, 3, 1, 8, 3, 1, 4, 1, 3, 4, ...

Before taking run-sums we had A078147.
For nonprime instead of nonsquarefree numbers we have A373822.
Positions of first appearances are A376265, sorted A376266.
For run-lengths instead of run-sums we have A376267.
For squarefree instead of nonsquarefree we have A376307.
For prime-powers instead of nonsquarefree numbers we have A376310.
For compression instead of run-sums we have A376312.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A003242 counts compressed compositions, ranks A333489.
A005117 lists squarefree numbers, differences A076259 (ones A375927).
A013929 lists nonsquarefree numbers, differences A078147.
Cf. A053797, A053806, A061398, A072284, A120992, A373197, A373413, A375707, A376305, A376306, A376311.

Total/@Split[Differences[Select[Range[1000],!SquareFreeQ[#]&]]]//Most

A376267 Run-lengths of first differences (A078147) of nonsquarefree numbers (A013929).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 1, 4, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 4, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 3, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2

Offset: 1

Gus Wiseman, Sep 27 2024

nonn

			The sequence of nonsquarefree numbers (A013929) is:
  4, 8, 9, 12, 16, 18, 20, 24, 25, 27, 28, 32, 36, 40, 44, 45, 48, 49, 50, ...
with first differences (A078147):
  4, 1, 3, 4, 2, 2, 4, 1, 2, 1, 4, 4, 4, 4, 1, 3, 1, 1, 2, 2, 2, 4, 3, 1, ...
with runs:
  (4),(1),(3),(4),(2,2),(4),(1),(2),(1),(4,4,4,4),(1),(3),(1,1),(2,2,2), ...
with lengths (A376267):
  1, 1, 1, 1, 2, 1, 1, 1, 1, 4, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, ...

Robert Israel, Table of n, a(n) for n = 1..10000

For prime instead of nonsquarefree numbers we have A333254.
For run-sums instead of run-lengths we have A376264.
For squarefree instead of nonsquarefree we have A376306.
For prime-powers instead of nonsquarefree numbers we have A376309.
For compression instead of run-lengths we have A376312.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A005117 lists squarefree numbers, differences A076259 (ones A375927).
A013929 lists nonsquarefree numbers, differences A078147.
Cf. A007674, A053797, A053806, A072284, A112925, A120992, A373198, A375707, A376305, A376307, A376311, A380595.

nsf:= remove(numtheory:-issqrfree, [$4..1000]):
S:= nsf[2..-1]-nsf[1..-2]:
R:= NULL: x:= 4: t:= 1:
for i from 2 to nops(S) do
  if S[i] = x then t:= t+1
  else R:= R,t; x:= S[i]; t:= 1
  fi
od:
R; # Robert Israel, Jan 27 2025

Length/@Split[Differences[Select[Range[1000], !SquareFreeQ[#]&]]]//Most

A057820 First differences of sequence of consecutive prime powers (A000961).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 2, 3, 1, 2, 4, 2, 2, 2, 2, 1, 5, 4, 2, 4, 2, 4, 6, 2, 3, 3, 4, 2, 6, 2, 2, 6, 8, 4, 2, 4, 2, 4, 8, 4, 2, 1, 3, 6, 2, 10, 2, 6, 6, 4, 2, 4, 6, 2, 10, 2, 4, 2, 12, 12, 4, 2, 4, 6, 2, 2, 8, 5, 1, 6, 6, 2, 6, 4, 2, 6, 4, 14, 4, 2, 4, 14, 6, 6, 4, 2, 4, 6, 2, 6, 6, 6, 4, 6, 8, 4, 8, 10, 2, 10

Offset: 1

Labos Elemer, Nov 08 2000

nonn

a(n) = 1 iff A000961(n) = A006549(k) for some k. - Reinhard Zumkeller, Aug 25 2002
Also run lengths of distinct terms in A070198. - Reinhard Zumkeller, Mar 01 2012
Does this sequence contain all positive integers? - Gus Wiseman, Oct 09 2024

			Odd differences arise in pairs in neighborhoods of powers of 2, like {..,2039,2048,2053,..} gives {..,11,5,..}

Michael B. Porter, Table of n, a(n) for n = 1..10000

Cf. A000961, A036616, A001223.
For perfect-powers (A001597) we have A053289.
For non-perfect-powers (A007916) we have A375706.
Positions of ones are A375734.
Run-compression is A376308.
Run-lengths are A376309.
Sorted positions of first appearances are A376340.
The second (instead of first) differences are A376596, zeros A376597.
Prime-powers:
- terms: A000961 or A246655, complement A024619
- differences: A057820 (this), first appearances A376341
- runs: A373675, A373673, A373674, A174965
- anti-runs: A373576, A120430, A006549, A373671
Non-prime-powers:
- terms: A361102
- differences: A375708 (ones A375713)
- runs: A373678, A373676, A373677, A110969 (A373669, sorted A373670)
- anti-runs: A373679, A373575, A255346, A373672
Cf. A000015, A014210, A014963, A025475, A025528, A027833, A037201, A046933, A076259, A078147, A093555, A345531, A376310.

a057820_list = zipWith (-) (tail a000961_list) a000961_list
-- Reinhard Zumkeller, Mar 01 2012

A057820 := proc(n)
        A000961(n+1)-A000961(n) ;
end proc: # R. J. Mathar, Sep 23 2016

Map[Length, Split[Table[Apply[LCM, Range[n]], {n, 1, 150}]]] (* Geoffrey Critzer, May 29 2015 *)
Join[{1},Differences[Select[Range[500],PrimePowerQ]]] (* Harvey P. Dale, Apr 21 2022 *)

isA000961(n) = (omega(n) == 1 || n == 1)
n_prev=1;for(n=2,500,if(isA000961(n),print(n-n_prev);n_prev=n)) \\ Michael B. Porter, Oct 30 2009

from sympy import primepi, integer_nthroot
def A057820(n):
    def f(x): return int(n+x-1-sum(primepi(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length())))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    r, k = m, f(m)+1
    while r != k: r, k = k, f(k)+1
    return r-m # Chai Wah Wu, Sep 12 2024

a(n) = A000961(n+1) - A000961(n).

Offset corrected and b-file adjusted by Reinhard Zumkeller, Mar 03 2012

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

N. J. A. Sloane, Apr 07 2000

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)

			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).

Peter Kagey, Table of n, a(n) for n = 1..10000
M. Filaseta and O. Trifonov, On Gaps between Squarefree Numbers. In Analytic Number Theory, Vol 85, 1990, Birkhäuser, Basel, pp. 235-253.
E. Fogels, On the average values of arithmetic functions, Proc. Cambridge Philos. Soc. 1941, 37: 358-372.
L. Marmet, First occurrences of squarefree gaps...
L. Marmet, First occurrences of square-free gaps and an algorithm for their computation, arXiv preprint arXiv:1210.3829 [math.NT], 2012. - From _N. J. A. Sloane_, Jan 01 2013
K. F. Roth, On the gaps between squarefree numbers, J. London Math. Soc. 1951 (2) 26:263-268.
Gus Wiseman, Four statistics for runs and antiruns of prime, nonprime, squarefree, and nonsquarefree numbers

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.
Cf. A007674, A020754, A045882, A053806, A061398, A061399.

SF:= select(numtheory:-issqrfree,[$1..1000]):
map(`-`,select(`>`,SF[2..-1]-SF[1..-2],1),1); # Robert Israel, Sep 22 2015

ReplaceAll[Differences[Select[Range@384, SquareFreeQ]] - 1, 0 -> Nothing] (* Michael De Vlieger, Sep 22 2015 *)

Offset set to 1 by Peter Kagey, Sep 29 2015

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

Leroy Quet, Oct 06 2005

nonn

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

Michael De Vlieger, Table of n, a(n) for n = 1..10000

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).
Cf. A061400, A112928, A112930.
Cf. A007674, A007947, A045882, A053797, A053806, A067535, A072284, A073247, A081217, A175216, A280892, A377783.

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

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 *)

a(n,p=prime(n))=while(!issquarefree(p--),); p \\ Charles R Greathouse IV, Aug 16 2017

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

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

Robert G. Wilson v, Mar 04 2002

nonn

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

			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)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Kaisa Matomäki, Maksym Radziwiłł and Terence Tao, Sign patterns of the Liouville and Möbius functions, Forum of Mathematics, Sigma, Vol. 4. (2016), e14.

Cf. A068780, A068140, A068781, A068782, A068783, A068784, A068785.
Cf. A049535, A077647, A078143, A045882.
Subsequence of A261869.
Cf. A007674, A373409, A373410, A373412, A375707, A376164.
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.
Cf. A049094, A061399, A072284, A120992, A294242, A373573, A375709.

a068781 n = a068781_list !! (n-1)
a068781_list = filter ((== 0) . a261869) [1..]
-- Reinhard Zumkeller, Sep 04 2015

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 *)

isok(m) = !moebius(m) && !moebius(m+1); \\ Michel Marcus, Feb 14 2021

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

A112926 Smallest squarefree integer > the n-th prime.

Original entry on oeis.org

3, 5, 6, 10, 13, 14, 19, 21, 26, 30, 33, 38, 42, 46, 51, 55, 61, 62, 69, 73, 74, 82, 85, 91, 101, 102, 105, 109, 110, 114, 129, 133, 138, 141, 151, 154, 158, 165, 170, 174, 181, 182, 193, 194, 199, 201, 213, 226, 229, 230, 235, 241, 246, 253, 258, 265, 271, 273

Offset: 1

Leroy Quet, Oct 06 2005

nonn

			10 is the smallest squarefree number greater than the 4th prime, 7. So a(4) = 10.
From _Gus Wiseman_, Dec 07 2024: (Start)
The first number line below shows the squarefree numbers. The second shows the primes:
--1--2--3-----5--6--7-------10-11----13-14-15----17----19----21-22-23-------26--
=====2==3=====5=====7==========11====13==========17====19==========23===========
(End)

Michael De Vlieger, Table of n, a(n) for n = 1..10000

Cf. A061400, A112928, A112929, A112930.
Restriction of A067535, differences A378087.
The unrestricted opposite is A070321, differences A378085.
The opposite is A112925, differences A378038.
Subtracting prime(n) from each term gives A240474, opposite A240473.
For nonsquarefree we have A377783, restriction of A120327.
The nonsquarefree differences are A377784, restriction of A378039.
First differences are A378037.
For perfect power we have A378249, A378617, A378250, A378251.
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.
Cf. A007674, A013928, A053797, A053806, A072284, A073247, A175216, A280892, A345531, A378032, A378618.

with(numtheory): a:=proc(n) local p,B,j: p:=ithprime(n): B:={}: for j from p+1 to p+20 do if abs(mobius(j))>0 then B:=B union {j} else B:=B fi od: B[1] end: seq(a(m),m=1..75); # Emeric Deutsch, Oct 10 2005

Do[k = Prime[n] + 1; While[ !SquareFreeQ[k], k++ ]; Print[k], {n, 1, 100}] (* Ryan Propper, Oct 10 2005 *)
With[{k = 120}, Table[SelectFirst[Range[Prime@ n + 1, Prime@ n + k], SquareFreeQ], {n, 58}]] (* Michael De Vlieger, Aug 16 2017 *)

a(n,p=prime(n))=while(!issquarefree(p++),); p \\ Charles R Greathouse IV, Aug 16 2017

a(n) = prime(n) + A240474(n). - Gus Wiseman, Dec 07 2024

More terms from Ryan Propper and Emeric Deutsch, Oct 10 2005

A375706 First differences of non-perfect-powers.

Original entry on oeis.org

1, 2, 1, 1, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 2

Offset: 1

Gus Wiseman, Aug 31 2024

nonn

Non-perfect-powers (A007916) are numbers without a proper integer root.

			The 5th non-perfect-power is 7, and the 6th is 10, so a(5) = 3.

Antti Karttunen, Table of n, a(n) for n = 1..100000

For prime-powers (A000961) we have A057820.
For perfect powers (A001597) we have A053289.
For nonprime numbers (A002808) we have A073783.
For squarefree numbers (A005117) we have A076259.
First differences of A007916.
For nonsquarefree numbers (A013929) we have A078147.
For non-prime-powers (A024619) we have A375708.
Positions of 1s are A375740, complement A375714.
Runs of non-perfect-powers:
- length: A375702 = A053289(n+1) - 1
- first: A375703 (same as A216765 with 2 exceptions)
- last: A375704 (same as A045542 with 8 removed)
- sum: A375705
Cf. A000040, A046933, A052410, A303707, A305630, A305631, A323054, A323088, A323090, A375736.

radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
Differences[Select[Range[100],radQ]]

up_to = 112;
A375706list(up_to) = { my(v=vector(up_to), pk=2, k=2, i=0); while(i<#v, k++; if(!ispower(k), i++; v[i] = k-pk; pk = k)); (v); };
v375706 = A375706list(up_to);
A375706(n) = v375706[n]; \\ Antti Karttunen, Jan 19 2025

from itertools import count
from sympy import mobius, integer_nthroot, perfect_power
def A375706(n):
    def f(x): return int(n+1-sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(2, x.bit_length())))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return next(i for i in count(m+1) if not perfect_power(i))-m # Chai Wah Wu, Sep 09 2024

a(n) = A007916(n+1) - A007916(n).

More terms from Antti Karttunen, Jan 19 2025

cp's OEIS Frontend

A376593 Second differences of consecutive nonsquarefree numbers (A013929). First differences of A078147.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

A376312 Run-compression of first differences (A078147) of nonsquarefree numbers (A013929).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A376264 Run-sums of first differences (A078147) of nonsquarefree numbers (A013929).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A376267 Run-lengths of first differences (A078147) of nonsquarefree numbers (A013929).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A057820 First differences of sequence of consecutive prime powers (A000961).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A053797 Lengths of successive gaps between squarefree numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

A112925 Largest squarefree integer < the n-th prime.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions