A031218 -id:A031218 - cp's OEIS frontend

A377703 First differences of the sequence A345531(k) = least prime-power greater than the k-th prime.

1, 3, 1, 5, 3, 3, 4, 2, 6, 1, 9, 2, 4, 2, 10, 2, 3, 7, 2, 6, 2, 8, 8, 4, 2, 4, 2, 4, 8, 7, 9, 2, 10, 2, 6, 6, 4, 2, 10, 2, 10, 2, 4, 2, 12, 12, 4, 2, 4, 6, 2, 2, 13, 7, 6, 2, 6, 4, 2, 6, 18, 4, 2, 4, 14, 6, 6, 6, 4, 6, 2, 12, 6, 4, 6, 8, 4, 8, 10, 2, 10, 2, 6

Offset: 1

Gus Wiseman, Nov 07 2024

nonn

What is the union of this sequence? In particular, does it contain 17?

First differences of A345531.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect-powers, differences A053289, seconds A376559.
A007916 lists the non-perfect-powers, differences A375706, seconds A376562.
A024619 lists the non-prime-powers, differences A375735, seconds A376599.
A080101 counts prime-powers between primes (exclusive).
A246655 lists the prime-powers, differences A057820 without first term.
A361102 lists the non-powers of primes, differences A375708.
A366833 counts prime-powers between primes, see A053607, A304521, A377057 (positive), A377286 (zero), A377287 (one), A377288 (two).
A377432 counts perfect-powers between primes, see A377434 (one), A377436 (zero), A377466 (multiple).
Cf. A000040, A008864, A031218, A377281, A377282, A377286, A377289, A377468.

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

from sympy import factorint, prime, nextprime
def A377703(n): return -next(filter(lambda m:len(factorint(m))<=1, count((p:=prime(n))+1)))+next(filter(lambda m:len(factorint(m))<=1, count(nextprime(p)+1))) # Chai Wah Wu, Nov 14 2024

A377288 Numbers k such that there are exactly two prime-powers between prime(k)+1 and prime(k+1)-1.

Original entry on oeis.org

4, 9, 30, 327, 3512

Offset: 1

Gus Wiseman, Oct 25 2024

Is this sequence finite? For this conjecture see A053706, A080101, A366833.

Any further terms are > 10^12. - Lucas A. Brown, Nov 08 2024

			Primes 9 and 10 are 23 and 29, and the interval (24, 25, 26, 27, 28) contains the prime-powers 25 and 27, so 9 is in the sequence.

Lucas A. Brown, Python program.

The interval from A008864(n) to A006093(n+1) has A046933 elements.
For powers of 2 instead of primes see A013597, A014210, A014234, A244508, A304521.
The corresponding primes are A053706.
The nearest prime-power before prime(n)-1 is A065514, difference A377289.
The nearest prime-power after prime(n)+1 is A345531, difference A377281.
These are the positions of 2 in A080101, or 3 in A366833.
For at least one prime-power we have A377057, primes A053607.
For no prime-powers we have A377286.
For exactly one prime-power we have A377287.
For squarefree instead of prime-power see A377430, A061398, A377431, A068360.
A000015 gives the least prime-power >= n.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A031218 gives the greatest prime-power <= n.
A246655 lists the prime-powers not including 1, complement A361102.
Cf. A001597, A002808, A024619, A053707, A064113, A065890, A075526, A095195, A224363, A276781, A376596, A376597, A377282.

Select[Range[100], Length[Select[Range[Prime[#]+1,Prime[#+1]-1],PrimePowerQ]]==2&]

prime(a(n)) = A053706(n).

A378035 Greatest perfect power < prime(n).

Original entry on oeis.org

1, 1, 4, 4, 9, 9, 16, 16, 16, 27, 27, 36, 36, 36, 36, 49, 49, 49, 64, 64, 64, 64, 81, 81, 81, 100, 100, 100, 100, 100, 125, 128, 128, 128, 144, 144, 144, 144, 144, 169, 169, 169, 169, 169, 196, 196, 196, 216, 225, 225, 225, 225, 225, 243, 256, 256, 256, 256

Offset: 1

Gus Wiseman, Nov 23 2024

nonn

Perfect powers (A001597) are 1 and numbers with a proper integer root, complement A007916.

			The first number line below shows the perfect powers.
The second shows each positive integer k at position prime(k).
-1-----4-------8-9------------16----------------25--27--------32------36----
===1=2===3===4=======5===6=======7===8=======9==========10==11==========12==

Restriction of A081676 to the primes.
Positions of last appearances are also A377283.
A version for squarefree numbers is A378032.
The opposite is A378249 (run lengths A378251), restriction of A377468 to the primes.
The union is A378253.
Terms appearing exactly once are A378355.
Run lengths are A378356, first differences of A377283, complement A377436.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect powers, differences A053289.
A007916 lists the nonperfect powers, differences A375706.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A080769 counts primes between perfect powers, prime powers A067871.
A131605 lists perfect powers that are not prime powers.
A377432 counts perfect powers between primes, zeros A377436, postpositives A377466.
Cf. A000015, A007918, A031218, A045542, A052410, A065514, A076412, A188951, A216765, A345531, A377434, A378250.

radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
Table[NestWhile[#-1&,Prime[n],radQ[#]&],{n,100}]

a(n) = my(k=prime(n)-1); while (!(ispower(k) || (k==1)), k--); k; \\ Michel Marcus, Nov 25 2024

from sympy import mobius, integer_nthroot, prime
def A378035(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(x-1+sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    m = (p:=prime(n)-1)-f(p)
    return bisection(lambda x:f(x)+m,m,m) # Chai Wah Wu, Nov 25 2024

A378249 Least perfect power > prime(n).

Original entry on oeis.org

4, 4, 8, 8, 16, 16, 25, 25, 25, 32, 32, 49, 49, 49, 49, 64, 64, 64, 81, 81, 81, 81, 100, 100, 100, 121, 121, 121, 121, 121, 128, 144, 144, 144, 169, 169, 169, 169, 169, 196, 196, 196, 196, 196, 216, 216, 216, 225, 243, 243, 243, 243, 243, 256, 289, 289, 289

Offset: 1

Gus Wiseman, Nov 21 2024

nonn

Perfect-powers (A001597) are numbers with a proper integer root, complement A007916.

Which terms appear only once? Just 128, 225, 256, 64009, 1295044?

			The first number line below shows the perfect powers. The second shows each prime.
-1-----4-------8-9------------16----------------25--27--------32------36------------------------49--
===2=3===5===7======11==13======17==19======23==========29==31==========37======41==43======47======

A version for prime powers (but starting with prime(k) + 1) is A345531.
Positions of last appearances are A377283, complement A377436.
Restriction of A377468 to the primes, for prime powers A000015.
The opposite is A378035, restriction of A081676.
The union is A378250.
Run lengths are A378251.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect powers, differences A053289, seconds A376559.
A007916 lists numbers that are not perfect powers, differences A375706, seconds A376562.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A131605 lists perfect powers that are not prime powers.
A377432 counts perfect powers between primes, zeros A377436, postpositives A377466.
Cf. A007918, A023055, A031218, A045542, A052410, A065514, A076412, A188951, A216765, A377434.

radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
Table[NestWhile[#+1&,Prime[n],radQ[#]&],{n,100}]

f(p) = p++; while(!ispower(p), p++); p;
lista(nn) = apply(f, primes(nn)); \\ Michel Marcus, Dec 19 2024

A377781 First differences of A065514(n) = greatest number < prime(n) that is 1 or a prime-power.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 14 2024

nonn

Note 1 is a power of a prime but not a prime-power.

Differences of A065514, which is the restriction of A031218 (differences A377782).
The opposite is A377703 (restriction of A000015), differences of A345531.
The opposite for nonsquarefree is A377784, differences of A377783.
For nonsquarefree we have A378034, differences of A378032 (restriction of A378033).
The opposite for squarefree is A378037, differences of A112926 (restriction of A067535).
For squarefree we have A378038, differences of A112925 (restriction of A070321).
A000040 lists the primes, differences A001223.
A000961 and A246655 list the prime-powers, differences A057820.
A024619 lists the non-prime-powers, differences A375735, seconds A376599.
A361102 lists the non-powers of primes, differences A375708.
Prime-powers between primes:
- A053607 primes
- A080101 count (exclusive)
- A304521 by bits
- A366833 count
- A377057 positive
- A377286 zero
- A377287 one
- A377288 two
Cf. A001597, A007916, A008864, A053289, A120327, A343249, A375706, A377281, A377282, A377289.

Differences[Table[NestWhile[#-1&,Prime[n]-1,#>1&&!PrimePowerQ[#]&],{n,100}]]

A378371 Distance between n and the least non prime power >= n, allowing 1.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 28 2024

nonn

Non prime powers allowing 1 (A361102) are numbers that are not a prime power (A246655), namely 1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, ...

			The least non prime power >= 4 is 6, so a(4) = 2.

Sequences obtained by adding n to each term are placed in parentheses below.
For prime we have A007920 (A151800), strict A013632.
For composite we have A010051 (A113646 except initial terms).
For perfect power we have A074984 (A377468)
For squarefree we have A081221 (A067535).
For nonsquarefree we have (A120327).
For non perfect power we have A378357 (A378358).
The opposite version is A378366 (A378367).
For prime power we have A378370, strict A377282 (A000015).
This sequence is A378371 (A378372).
A000040 lists the primes, differences A001223.
A000961 and A246655 list the prime powers, differences A057820.
A024619 and A361102 list the non prime powers, differences A375708 and A375735.
Prime powers between primes: A053607, A080101, A304521, A366833, A377057.
Cf. A007916, A053707, A065514, A276781 (A031218), A343249, A345531, A376596, A376597, A377051, A377054, A377289, A378457.

Table[NestWhile[#+1&,n,PrimePowerQ[#]&]-n,{n,100}]

a(n) = A378372(n) - n.

A378372 Least non prime power >= n, allowing 1.

Original entry on oeis.org

1, 6, 6, 6, 6, 6, 10, 10, 10, 10, 12, 12, 14, 14, 15, 18, 18, 18, 20, 20, 21, 22, 24, 24, 26, 26, 28, 28, 30, 30, 33, 33, 33, 34, 35, 36, 38, 38, 39, 40, 42, 42, 44, 44, 45, 46, 48, 48, 50, 50, 51, 52, 54, 54, 55, 56, 57, 58, 60, 60, 62, 62, 63, 65, 65, 66, 68

Offset: 1

Gus Wiseman, Nov 29 2024

nonn

Non prime powers allowing 1 (A361102) are numbers that are not a prime power (A246655), namely 1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, ...

			The least non prime power >= 4 is 6, so a(4) = 6.

Sequences obtained by subtracting n from each term are placed in parentheses below.
For prime power we have A000015 (A378370).
For squarefree we have A067535 (A081221).
For composite we have A113646 (A010051).
For nonsquarefree we have A120327.
For prime we have A151800 (A007920), strict (A013632).
Run-lengths are 1 and A375708.
For perfect power we have A377468 (A074984).
For non-perfect power we have A378358 (A378357).
The opposite is A378367, distance A378366.
This sequence is A378372 (A378371).
A000040 lists the primes, differences A001223.
A000961 and A246655 list the prime powers, differences A057820.
A024619 and A361102 list the non prime powers, differences A375708 and A375735.
Prime powers between primes: A053607, A080101, A304521, A366833, A377057.
Cf. A001597, A007916, A031218 (A276781), A053707, A065514, A345531 (A377281), A377051, A377282, A377289, A378457.

Table[NestWhile[#+1&,n,PrimePowerQ[#]&],{n,100}]

a(n) = A378371(n) + n.

A179278 Largest nonprime integer <= n.

Original entry on oeis.org

1, 1, 1, 4, 4, 6, 6, 8, 9, 10, 10, 12, 12, 14, 15, 16, 16, 18, 18, 20, 21, 22, 22, 24, 25, 26, 27, 28, 28, 30, 30, 32, 33, 34, 35, 36, 36, 38, 39, 40, 40, 42, 42, 44, 45, 46, 46, 48, 49, 50, 51, 52, 52, 54, 55, 56, 57, 58, 58, 60, 60, 62, 63, 64, 65, 66, 66, 68, 69, 70, 70, 72

Offset: 1

Reinhard Zumkeller, Jul 08 2010

			From _Gus Wiseman_, Dec 04 2024: (Start)
The nonprime integers <= n:
  1  1  1  4  4  6  6  8  9  10  10  12  12  14  15  16
           1  1  4  4  6  8  9   9   10  10  12  14  15
                 1  1  4  6  8   8   9   9   10  12  14
                       1  4  6   6   8   8   9   10  12
                          1  4   4   6   6   8   9   10
                             1   1   4   4   6   8   9
                                     1   1   4   6   8
                                             1   4   6
                                                 1   4
                                                     1
(End)

Cf. A014684, A113523, A113638, A062298.
For prime we have A007917.
For nonprime we have A179278 (this).
For squarefree we have A070321.
For nonsquarefree we have A378033.
For prime power we have A031218.
For non prime power we have A378367.
For perfect power we have A081676.
For non perfect power we have A378363.
A000040 lists the primes, differences A001223.
A002808 lists the composite numbers, differences A073783.
A018252 lists the nonprimes, differences A065310.
A095195 has row n equal to the k-th differences of the prime numbers.
A113646 gives least nonprime >= n.
A151800 gives the least prime > n, weak version A007918.
A377033 has row n equal to the k-th differences of the composite numbers.

Array[# - Boole[PrimeQ@ #] - Boole[# == 3] &, 72] (* Michael De Vlieger, Oct 13 2018 *)
Table[Max@@Select[Range[n],!PrimeQ[#]&],{n,30}] (* Gus Wiseman, Dec 04 2024 *)

a(n) = if (isprime(n), if (n==3, 1, n-1), n); \\ Michel Marcus, Oct 13 2018

For n > 3: a(n) = A113523(n) = A014684(n);
For n > 0: a(n) = A113638(n). - Georg Fischer, Oct 12 2018
A005171(a(n)) = 1; A010051(a(n)) = 0.
a(n) = A018252(A062298(n)). - Ridouane Oudra, Aug 22 2025

Inequality in the name reversed by Gus Wiseman, Dec 05 2024

A366835 In the pair (A246655(n), A246655(n+1)), how many primes are there?

Original entry on oeis.org

2, 1, 1, 2, 1, 0, 1, 2, 1, 1, 2, 2, 1, 0, 1, 2, 1, 1, 2, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 2, 1, 1, 2, 2, 2, 2, 2, 2, 2, 1, 0, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 2

Offset: 1

Paolo Xausa, Oct 25 2023

nonn

First 0 terms appear at n = 6, 14, 41, 359, 3589, corresponding to consecutive prime powers (8,9), (25,27), (121,125), (2187,2197) and (32761,32768), respectively (cf. A068315 and A068435).

There cannot be primes strictly between consecutive prime powers, so we get the same result considering the whole interval (not just the pair). - Gus Wiseman, Dec 25 2024

			a(1) = 2 because in the first prime power pair (2 and 3) there are two primes.
a(14) = 0 because in the 14th prime power pair (25 and 27) there are no primes.

Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, 1038 X 1038 raster of a(n), n = 1..1077444, read left to right in rows, then top to bottom, showing a(n) = 0 in white, a(n) = 1 in red, and a(n) = 2 in dark blue.

For perfect powers instead of prime powers we have A080769.
Positions of 1 are A379155, indices of A379157.
Positions of 0 are A379156, indices of A068315.
Positions of 2 are A379158, indices of A379541.
A000015 gives the least prime power >= n.
A000040 lists the primes, differences A001223.
A031218 gives the greatest prime power <= n.
A065514 gives the greatest prime power < prime(n), difference A377289.
A080101 and A366833 count prime powers between primes, see A053607, A304521.
A246655 lists the prime powers, differences A057820.
Cf. A000961, A025474, A046933, A067871, A068435, A175106, A178700, A345531, A362965, A377281.

With[{upto=500},Map[Count[#,_?PrimeQ]&,Partition[Select[Range[upto],PrimePowerQ],2,1]]] (* Considers prime powers up to 500 *)

lista(nn) = my(v=[p| p <- [1..nn], isprimepower(p)]); vector(#v-1, k, isprime(v[k]) + isprime(v[k+1])); \\ Michel Marcus, Oct 26 2023

A377467 Number of perfect-powers x in the range 2^n < x < 2^(n+1).

Original entry on oeis.org

0, 0, 0, 1, 2, 2, 4, 6, 7, 10, 15, 23, 31, 41, 60, 81, 117, 165, 230, 321, 452, 634, 891, 1252, 1766, 2486, 3504, 4935, 6958, 9815, 13849, 19537, 27577, 38932, 54971, 77640, 109667, 154921, 218878, 309276, 437046, 617657, 872967, 1233895, 1744152, 2465546, 3485477

Offset: 0

Gus Wiseman, Nov 04 2024

nonn

Perfect-powers (A001597) are numbers with a proper integer root, complement A007916.

Also the number of perfect-powers, except for powers of 2, with n bits.

			The perfect-powers in each prescribed range (rows):
    .
    .
    .
    9
   25   27
   36   49
   81  100  121  125
  144  169  196  216  225  243
  289  324  343  361  400  441  484
  529  576  625  676  729  784  841  900  961 1000
The binary expansions for n >= 3 (columns):
    1001  11001  100100  1010001  10010000  100100001
          11011  110001  1100100  10101001  101000100
                         1111001  11000100  101010111
                         1111101  11011000  101101001
                                  11100001  110010000
                                  11110011  110111001
                                            111100100

The version for squarefree numbers is A077643.
The version for prime-powers is A244508.
For primes instead of powers of 2 we have A377432, zeros A377436.
Including powers of 2 in the range gives A377435.
The version for non-perfect-powers is A377701.
The union of all numbers counted is A377702.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect-powers, differences A053289.
A007916 lists the non-perfect-powers, differences A375706.
A081676 gives the greatest perfect-power <= n.
A131605 lists perfect-powers that are not prime-powers.
A377468 gives the least perfect-power > n.
Cf. A000015, A013597, A014210, A014234, A023055, A031218, A045542, A052410, A065514, A069623, A216765, A246655, A345531.

perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All,2]]>1;
Table[Length[Select[Range[2^n+1,2^(n+1)-1],perpowQ]],{n,0,15}]

from sympy import mobius, integer_nthroot
def A377467(n):
    def f(x): return int(1-sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    return f((1<Chai Wah Wu, Nov 05 2024

For n != 1, a(n) = A377435(n) - 1.

a(26)-a(46) from Chai Wah Wu, Nov 05 2024

cp's OEIS Frontend

A377703 First differences of the sequence A345531(k) = least prime-power greater than the k-th prime.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Python

A377288 Numbers k such that there are exactly two prime-powers between prime(k)+1 and prime(k+1)-1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A378035 Greatest perfect power < prime(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

A378249 Least perfect power > prime(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

A377781 First differences of A065514(n) = greatest number < prime(n) that is 1 or a prime-power.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A378371 Distance between n and the least non prime power >= n, allowing 1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A378372 Least non prime power >= n, allowing 1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A179278 Largest nonprime integer <= n.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

PARI

Formula