A378249 -id:A378249 - cp's OEIS frontend

A378617 First differences of A378249 (next perfect power after prime(n)).

0, 4, 0, 8, 0, 9, 0, 0, 7, 0, 17, 0, 0, 0, 15, 0, 0, 17, 0, 0, 0, 19, 0, 0, 21, 0, 0, 0, 0, 7, 16, 0, 0, 25, 0, 0, 0, 0, 27, 0, 0, 0, 0, 20, 0, 0, 9, 18, 0, 0, 0, 0, 13, 33, 0, 0, 0, 0, 0, 0, 35, 0, 0, 0, 0, 19, 0, 18, 0, 0, 0, 39, 0, 0, 0, 0, 0, 41, 0, 0, 0

Offset: 1

Gus Wiseman, Dec 09 2024

nonn

This is the next perfect power after prime(n+1), minus the next perfect power after prime(n).

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

Positions of positives are A377283.
Positions of zeros are A377436.
The restriction to primes has first differences A377468.
A version for nonsquarefree numbers is A377784, differences of A377783.
The opposite is differences of A378035 (restriction of A081676).
First differences of A378249, run-lengths A378251.
Without zeros we have differences of A378250.
A000040 lists the primes, differences A001223.
A001597 lists the perfect powers, differences A053289.
A007916 lists the non perfect powers, differences A375706.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A377432 counts perfect powers between primes.
A378356 - 1 gives next prime after perfect powers, union A378365 - 1.
Cf. A000015, A007918, A023055, A045542, A052410, A065514, A076412, A188951, A216765, A377431, A377434.

perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All,2]]>1;
Table[NestWhile[#+1&,Prime[n],Not@*perpowQ],{n,100}]//Differences

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

A345531 Smallest prime power greater than the n-th prime.

Original entry on oeis.org

3, 4, 7, 8, 13, 16, 19, 23, 25, 31, 32, 41, 43, 47, 49, 59, 61, 64, 71, 73, 79, 81, 89, 97, 101, 103, 107, 109, 113, 121, 128, 137, 139, 149, 151, 157, 163, 167, 169, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 243, 256, 263, 269, 271, 277

Offset: 1

Dario T. de Castro, Jun 20 2021

nonn

Take the family of correlated prime-indexed conjectures appearing in A343249 - A343253, in which an alternative formula for the p-adic order of positive integers is proposed. There, the general p-indexed conjecture says that v_p(n), the p-adic order of n, is given by the formula: v_p(n) = log_p(n / L_p(k0, n)), where L_p(k0, n) is the lowest common denominator of the elements of the set S_p(k0, n) = {(1/n)*binomial(n, k), with 0 < k <= k0 such that k is not divisible by p}. Evidence suggests that the primality of p is a necessary condition in this general conjecture. So, if a composite number q is used instead of a prime p in the proposed formula for the p-adic (now, q-adic) order of n, the first counterexample (failure) is expected to occur for n = q * a(i), where i is the index of the smallest prime that divides q.

The prime-power a(n) is at most the next prime, so this sequence is strictly increasing. See also A366833. - Gus Wiseman, Nov 06 2024

			a(4) = 8 because the fourth prime number is 7, and the least power of a prime which is greater than 7 is 2^3 = 8.

Robert Israel, Table of n, a(n) for n = 1..10000
Dario T. de Castro, P-adic Order of Positive Integers via Binomial Coefficients, INTEGERS, Electronic J. of Combinatorial Number Theory, Vol. 22, Paper A61, 2022.

Cf. A000040, A007814, A007949, A008864, A343249, A343250, A112765, A343251, A214411, A343252, A343253.
Starting with n instead of prime(n): A000015, A031218, A377468, A377780, A377782.
Opposite (greatest prime-power less than): A065514, A377289, A377781.
For squarefree instead of prime-power: A112926, opposite A112925.
The difference from prime(n) is A377281.
The prime terms have indices A377286(n) - 1.
First differences are A377703.
A version for perfect-powers is A378249.
A000961 and A246655 list the prime-powers, differences A057820.
A024619 and A361102 list the non-prime-powers, differences A375735.
Cf. A366833, A053607, A080101, A304521, A366833, A377057, A377283, A377287, A377288, A377432.

f:= proc(n) local p,x;
  p:= ithprime(n);
  for x from p+1 do
    if nops(numtheory:-factorset(x)) = 1 then return x fi
  od
end proc:
map(f, [$1..100]); # Robert Israel, Aug 25 2024

a[i_]:= Module[{j, k, N = 0, tab={}}, tab = Sort[Drop[DeleteDuplicates[Flatten[Table[ If[Prime[j]^k > Prime[i], Prime[j]^k], {j, 1, i+1}, {k, 1, Floor[Log[Prime[j], Prime[i+1]]]}]]], 1]]; N = Take[tab, 1][[1]]; N];
tabseq = Table[a[i],{i, 1, 100}];
(* second program *)
Table[NestWhile[#+1&,Prime[n]+1, Not@*PrimePowerQ],{n,100}] (* Gus Wiseman, Nov 06 2024 *)

A000015(n) = for(k=n,oo,if((1==k)||isprimepower(k),return(k)));
A345531(n) = A000015(1+prime(n)); \\ Antti Karttunen, Jul 19 2021

from itertools import count
from sympy import prime, factorint
def A345531(n): return next(filter(lambda m:len(factorint(m))<=1, count(prime(n)+1))) # Chai Wah Wu, Oct 25 2024

a(n) = A000015(1+A000040(n)). - Antti Karttunen, Jul 19 2021

a(n) = A000015(A008864(n)). - Omar E. Pol, Oct 27 2021

A065514 Largest power of a prime < prime(n).

Original entry on oeis.org

1, 2, 4, 5, 9, 11, 16, 17, 19, 27, 29, 32, 37, 41, 43, 49, 53, 59, 64, 67, 71, 73, 81, 83, 89, 97, 101, 103, 107, 109, 125, 128, 131, 137, 139, 149, 151, 157, 163, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 243, 256, 257, 263, 269, 271

Offset: 1

Reinhard Zumkeller, Nov 27 2001

nonn

Starting with n instead of prime(n) gives A031218 (A377282, A377782).
The squarefree version is A112925 (A070321, A378038).
The opposite squarefree version is A112926 (A378037, restriction of A067535).
Difference from prime(n) is A377289 (restriction of A276781, opposite A377281).
First differences are A377781.
The nonsquarefree version is A378032 (A377783 (restriction of A378033), A378034, A378040).
The perfect power version is A378035.
A000015 gives the least prime power >= n, differences A377780.
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.
A345531 gives the least prime power > prime(n), differences A377703.
Prime powers between primes: A053607, A080101, A304521, A366833, A377057, A377286.
Cf. A014234, A053707, A065890, A376597, A377051, A378249, A378252.

lpp[n_]:=Module[{k=n-1},While[!PrimePowerQ[k],k--];k]; Join[{1},Table[ lpp[ n],{n,Prime[Range[2,60]]}]] (* Harvey P. Dale, Nov 24 2018 *)

from sympy import factorint, prime
def A065514(n): return next(filter(lambda m:len(factorint(m))<=1, range(prime(n)-1,0,-1))) # Chai Wah Wu, Oct 25 2024

Name edited (1 is technically not a prime power even though it is a power of a prime) by Gus Wiseman, Dec 03 2024.

A377283 Nonnegative integers k such that either k = 0 or there is a perfect power x in the range prime(k) < x < prime(k+1).

Original entry on oeis.org

0, 2, 4, 6, 9, 11, 15, 18, 22, 25, 30, 31, 34, 39, 44, 47, 48, 53, 54, 61, 66, 68, 72, 78, 85, 92, 97, 99, 105, 114, 122, 129, 137, 146, 154, 162, 168, 172, 181, 191, 200, 210, 217, 219, 228, 240, 251, 263, 269, 274, 283, 295, 306, 309, 319, 327, 329, 342, 357

Offset: 1

Gus Wiseman, Nov 21 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==

A version for prime powers is A377057, exclusive A377287.
A version for squarefree numbers is A377431.
Positions of positive terms in A377432 (counts perfect powers between primes).
The case of a unique choice is A377434 (a subset).
The complement (no choices) is A377436.
The case of at least two choices is A377466 (a subset).
Positions of last appearances in A378249.
First-differences are A378251.
This is A378365 - 1, union of A378356 - 1.
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.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A131605 lists perfect powers that are not prime powers.
Cf. A000015, A045542, A065514, A076412, A081676, A188951, A216765, A345531, A377468, A378035, A378250.

perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All,2]]>1;
Select[Range[0,100],#==0||Length[Select[Range[Prime[#]+1,Prime[#+1]-1],perpowQ]]>0&]

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

A378251 Number of primes between consecutive perfect powers, zeros omitted.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 23 2024

First differences of A377283 and A378365. Run-lengths of A378035 and A378249.

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 prime. To get a(n) we count the primes between consecutive perfect powers, skipping the cases where there are none.
-1-----4-------8-9------------16----------------25--27--------32------36----
===2=3===5===7======11==13======17==19======23==========29==31==========37==

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

Same as A080769 with 0's removed (which were at positions A274605).
First differences of A377283 and A378365 (union of A378356).
Run-lengths of A378035 (union A378253) and A378249 (union A378250).
The version for nonprime prime powers is A378373, with zeros A067871.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect powers, differences A053289, run-lengths of A377468.
A007916 lists the non-perfect powers, differences A375706.
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, see A377434, A377436, A377466.
Cf. A007918, A023055, A045542, A052410, A076412, A081676, A188951.
Cf. A216765, A377057, A377431, A378355.

N:= 10^6: # to use perfect powers up to N
PP:= {1,seq(seq(i^j,j=2..ilog[i](N)),i=2..isqrt(N))}:
PP:= sort(convert(PP,list)):
M:= map(numtheory:-pi, PP):
subs(0=NULL, M[2..-1]-M[1..-2]): # Robert Israel, Jan 23 2025

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

A378250 Perfect-powers x > 1 such that it is not possible to choose a prime y and a perfect-power z satisfying x > y > z.

Original entry on oeis.org

4, 8, 16, 25, 32, 49, 64, 81, 100, 121, 128, 144, 169, 196, 216, 225, 243, 256, 289, 324, 343, 361, 400, 441, 484, 512, 529, 576, 625, 676, 729, 784, 841, 900, 961, 1000, 1024, 1089, 1156, 1225, 1296, 1331, 1369, 1444, 1521, 1600, 1681, 1728, 1764, 1849, 1936

Offset: 1

Gus Wiseman, Nov 21 2024

nonn

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

			The first number line below shows the perfect-powers. The second shows the primes. The third is a(n).
-1-----4-------8-9------------16----------------25--27--------32------36----
===2=3===5===7======11==13======17==19======23==========29==31==========37==
       4       8              16                25            32
The terms together with their prime indices begin:
     4: {1,1}
     8: {1,1,1}
    16: {1,1,1,1}
    25: {3,3}
    32: {1,1,1,1,1}
    49: {4,4}
    64: {1,1,1,1,1,1}
    81: {2,2,2,2}
   100: {1,1,3,3}
   121: {5,5}
   128: {1,1,1,1,1,1,1}
   144: {1,1,1,1,2,2}
   169: {6,6}
   196: {1,1,4,4}
   216: {1,1,1,2,2,2}
   225: {2,2,3,3}
   243: {2,2,2,2,2}
   256: {1,1,1,1,1,1,1,1}

A version for prime-powers (but starting with prime(k) + 1) is A345531.
The opposite is union of A378035, restriction of A081676.
Union of A378249, 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 the non-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, positive A377283, postpositive A377466.
Cf. A000015, A007918, A023055, A045542, A052410, A076412, A188951, A216765, A377431, A377434, A377468.

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

A378356 Prime index of the next prime after the n-th perfect power.

Original entry on oeis.org

1, 3, 5, 5, 7, 10, 10, 12, 12, 16, 19, 23, 26, 31, 31, 32, 35, 40, 45, 48, 49, 54, 55, 62, 67, 69, 73, 79, 86, 93, 98, 100, 106, 115, 123, 130, 138, 147, 155, 163, 169, 173, 182, 192, 201, 211, 218, 220, 229, 241, 252, 264, 270, 275, 284, 296, 307, 310, 320

Offset: 1

Gus Wiseman, Dec 05 2024

nonn

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

First differences are A080769.
Union is A378365.
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.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A377432 counts perfect powers between primes, see A377434, A377436, A377466.
A378249 gives the least perfect power > prime(n), restriction of A377468.
Cf. A045542, A052410, A065514, A068361, A076412, A081676, A216765, A345531, A377283, A378035, A378250, A378251, A378253.

Table[PrimePi[NextPrime[n]],{n,Select[Range[1000],perpowQ]}]

a(n) = A000720(A001597(n)) + 1.

A378355 Numbers appearing exactly once in A378035 (greatest perfect power < prime(n)).

Original entry on oeis.org

125, 216, 243, 64000, 1295029, 2535525316, 542939080312

Offset: 1

Gus Wiseman, Nov 26 2024

These are perfect-powers p such that the interval from p to the next perfect power contains a unique prime.

Is this sequence infinite? See A178700.

			We have 125 because 127 is the only prime between 125 and 128.

The next prime is A178700.
Singletons in A378035 (union A378253), restriction of A081676.
The next perfect power is A378374.
Swapping primes and perfect powers gives A379154, unique case of A377283.
A000040 lists the primes, differences A001223.
A001597 lists the perfect powers, differences A053289.
A007916 lists the not perfect powers, differences A375706.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A377432 counts perfect powers between primes, see A377434, A377436, A377466.
A378249 gives least perfect power > prime(n) (run-lengths A378251), restrict of A377468.
Cf. A000961, A031218, A045542, A052410, A067871, A076412, A080769, A131605, A188951, A216765, A378250, A378356.

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

A151800(a(n)) = A178700(n).

cp's OEIS Frontend

A378617 First differences of A378249 (next perfect power after prime(n)).

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A112926 Smallest squarefree integer > the n-th prime.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A345531 Smallest prime power greater than the n-th prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A065514 Largest power of a prime < prime(n).

Views

Author

Keywords

Crossrefs

Programs

Mathematica

Python

Extensions

A377283 Nonnegative integers k such that either k = 0 or there is a perfect power x in the range prime(k) < x < prime(k+1).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A378035 Greatest perfect power < prime(n).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

A378251 Number of primes between consecutive perfect powers, zeros omitted.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A378250 Perfect-powers x > 1 such that it is not possible to choose a prime y and a perfect-power z satisfying x > y > z.

Views

Author

Keywords