A131605 -id:A131605 - cp's OEIS frontend

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

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

A378357 Distance from n to the least non perfect power >= n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Nov 24 2024

nonn

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

All terms are <= 2 because the only adjacent perfect powers are 8 and 9.

The version for prime numbers is A007920, subtraction of A159477 or A007918.
The version for perfect powers is A074984, subtraction of A377468.
The version for squarefree numbers is A081221, subtraction of A067535.
Subtracting from n gives A378358, opposite A378363.
The opposite version is A378364.
The version for nonsquarefree numbers is A378369, subtraction of A120327.
The version for prime powers is A378370, subtraction of A000015.
The version for non prime powers is A378371, subtraction of A378372.
The version for composite numbers is A378456, subtraction of A113646.
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.
A377432 counts perfect powers between primes, zeros A377436.
Cf. A013597, A013603, A013632, A031218, A045542, A052410, A081676, A131605, A188951, A216765, A375702, A375736.

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

from sympy import perfect_power
def A378357(n): return 0 if n>1 and perfect_power(n)==False else 1 if perfect_power(n+1)==False else 2 # Chai Wah Wu, Nov 27 2024

a(n) = n - A378358(n).

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

A378358 Least non-perfect-power >= n.

Original entry on oeis.org

2, 2, 3, 5, 5, 6, 7, 10, 10, 10, 11, 12, 13, 14, 15, 17, 17, 18, 19, 20, 21, 22, 23, 24, 26, 26, 28, 28, 29, 30, 31, 33, 33, 34, 35, 37, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 65, 66, 67

Offset: 1

Gus Wiseman, Nov 24 2024

nonn

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

The version for prime-powers is A000015, for non-prime-powers A378372.
The union is A007916, complement A001597.
The version for nonsquarefree numbers is A067535, negative A120327 (subtract A378369).
The version for composite numbers is A113646.
The version for prime numbers is A159477.
The run-lengths are A375706.
Terms appearing only once are A375738, multiple times A375703.
The version for perfect-powers is A377468.
Subtracting from n gives A378357.
The opposite version is A378363, for perfect-powers A081676.
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.
A377432 counts perfect-powers between primes, zeros A377436.
Cf. A007918, A014210, A031218, A045542, A052410, A065514, A070321, A076412, A188951, A216765, A345531, A378033.

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

from sympy import mobius, integer_nthroot
def A378358(n):
    def f(x): return int(1-sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    a = max(1,n-f(n-1))
    m, k = a, f(a)+a
    while m != k: m, k = k, f(k)+a
    return m # Chai Wah Wu, Nov 26 2024

from sympy import perfect_power
def A378358(n): return n if n>1 and perfect_power(n)==False else n+1 if perfect_power(n+1)==False else n+2 # Chai Wah Wu, Nov 27 2024

a(n) = n - A378357(n).

A378363 Greatest number <= n that is 1 or not a perfect-power.

Original entry on oeis.org

1, 2, 3, 3, 5, 6, 7, 7, 7, 10, 11, 12, 13, 14, 15, 15, 17, 18, 19, 20, 21, 22, 23, 24, 24, 26, 26, 28, 29, 30, 31, 31, 33, 34, 35, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 63, 65, 66, 67

Offset: 1

Gus Wiseman, Nov 24 2024

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

			In the non-perfect-powers ... 5, 6, 7, 10, 11 ... the greatest term <= 8 is 7, so a(8) = 7.

The union is A007916, complement A001597.
The version for prime numbers is A007917 or A151799, opposite A159477.
The version for prime-powers is A031218, opposite A000015.
The version for squarefree numbers is A067535, opposite A070321.
The version for perfect-powers is A081676, opposite A377468.
The version for composite numbers is A179278, opposite A113646.
Terms appearing multiple times are A375704, opposite A375703.
The run-lengths are A375706.
Terms appearing only once are A375739, opposite A375738.
The version for nonsquarefree numbers is A378033, opposite A120327.
The opposite version is A378358.
Subtracting n gives A378364, opposite A378357.
The version for non-prime-powers is A378367 (subtracted A378371), opposite A378372.
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.
A377432 counts perfect-powers between primes, zeros A377436.
Cf. A007918, A014210, A014234, A045542, A052410, A065514, A076412, A162966, A188951, A345531.

perpowQ[n_]:=n==1||GCD@@FactorInteger[n][[All,2]]>1;
Table[NestWhile[#-1&,n,#>1&&perpowQ[#]&],{n,100}]

from sympy import mobius, integer_nthroot
def A378363(n):
    def f(x): return int(1-sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
    a = n-f(n)
    m, k = a, f(a)+a
    while m != k: m, k = k, f(k)+a
    return m # Chai Wah Wu, Nov 26 2024

A368682 Products of primorials that are perfect powers but not prime powers.

Original entry on oeis.org

36, 144, 216, 576, 900, 1296, 1728, 2304, 3600, 5184, 7776, 9216, 13824, 14400, 20736, 27000, 32400, 36864, 44100, 46656, 57600, 82944, 110592, 129600, 147456, 176400, 186624, 216000, 230400, 248832, 279936, 331776, 373248, 518400, 589824, 705600, 746496, 810000

Offset: 1

Michael De Vlieger, Jan 02 2024

nonn

Intersection of A025487 and A131605.
Proper subset of A286708.
Contains A365308 (perfect powers of composite primorials) and A368508 (perfect powers of composite superprimorials).
These numbers are perfect powers of some smaller product of primorials.

			b(n) = A025487(n).
a(1) = b(11) = 36 = 6^2 = b(4)^2,
a(2) = b(19) = 144 = 12^2 = b(6)^2,
a(3) = b(23) = 216 = 6^3 = b(4)^3,
a(4) = b(33) = 576 = 24^2 = b(8)^2,
a(5) = b(38) = 900 = 30^2 = b(9)^2, etc.

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

Cf. A000079, A000961, A001597, A002110, A025487, A131605, A286708, A365308, A368508, A368681.

Select[Range[36, 2^18, 2], And[Union@ Differences@ PrimePi@ #1 == {1}, AllTrue[Union@ Differences@ #2, # <= 0 &], GCD @@ #2 > 1] & @@ Transpose@ FactorInteger[#] &]

This sequence is { A368681 \ A000079 }.

A375397 Numbers divisible by the square of some prime factor other than the least. Non-hooklike numbers.

Original entry on oeis.org

18, 36, 50, 54, 72, 75, 90, 98, 100, 108, 126, 144, 147, 150, 162, 180, 196, 198, 200, 216, 225, 234, 242, 245, 250, 252, 270, 288, 294, 300, 306, 324, 338, 342, 350, 360, 363, 375, 378, 392, 396, 400, 414, 432, 441, 450, 468, 484, 486, 490, 500, 504, 507, 522

Offset: 1

Gus Wiseman, Aug 16 2024

nonn

Contains no squarefree numbers A005117 or prime powers A000961, but some perfect powers A131605.
Also numbers k such that the minima of the maximal anti-runs in the weakly increasing sequence of prime factors of k (with multiplicity) are not identical. Here, an anti-run is a sequence with no adjacent equal parts, and the minima of the maximal anti-runs in a sequence are obtained by splitting it into maximal anti-run subsequences and taking the least term of each. Note the prime factors can alternatively be taken in weakly decreasing order.
Includes all terms of A036785 = non-products of a squarefree number and a prime power.
The asymptotic density of this sequence is 1 - (1/zeta(2)) * (1 + Sum_{p prime} (1/(p^2-p)) / Product_{primes q <= p} (1 + 1/q)) = 0.11514433883... . - Amiram Eldar, Oct 26 2024

			The prime factors of 300 are {2,2,3,5,5}, with maximal anti-runs ((2),(2,3,5),(5)), with minima (2,2,5), so 300 is in the sequence.
The terms together with their prime indices begin:
    18: {1,2,2}
    36: {1,1,2,2}
    50: {1,3,3}
    54: {1,2,2,2}
    72: {1,1,1,2,2}
    75: {2,3,3}
    90: {1,2,2,3}
    98: {1,4,4}
   100: {1,1,3,3}
   108: {1,1,2,2,2}
   126: {1,2,2,4}
   144: {1,1,1,1,2,2}

Amiram Eldar, Table of n, a(n) for n = 1..10000

A superset of A036785.
The complement for maxima is A065200, counted by A034296.
For maxima instead of minima we have A065201, counted by A239955.
A version for compositions is A374520, counted by A374640.
Also positions of non-constant rows in A375128, sums A374706, ranks A375400.
The complement is A375396, counted by A115029.
The complement for distinct minima is A375398, counted by A375134.
For distinct instead of identical minima we have A375399, counts A375404.
Partitions of this type are counted by A375405.
A000041 counts integer partitions, strict A000009.
A003242 counts anti-run compositions, ranks A333489.
A number's prime factors (A027746, reverse A238689) have sum A001414, min A020639, max A006530.
A number's prime indices (A112798, reverse A296150) have sum A056239, min A055396, max A061395.
Both have length A001222, distinct A001221.
Cf. A000005, A013661, A046660, A272919, A319066, A358905, A374686, A374704, A374742, A375133, A375136, A375401.

Select[Range[100],!SameQ@@Min /@ Split[Flatten[ConstantArray@@@FactorInteger[#]],UnsameQ]&]

is(k) = if(k > 1, my(e = factor(k)[, 2]); vecprod(e) > e[1], 0); \\ Amiram Eldar, Oct 26 2024

Name edited by Peter Munn, May 08 2025

A378253 Perfect powers p such that there are no other perfect powers between p and the least prime > p.

Original entry on oeis.org

1, 4, 9, 16, 27, 36, 49, 64, 81, 100, 125, 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

Offset: 1

Gus Wiseman, Nov 26 2024

nonn

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

Each term is the greatest perfect power < prime(k) for some k.

			The first number line below shows the perfect powers. The second shows each prime. To get a(n), we take the last perfect power in each interval between consecutive primes, omitting 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==

Union of A378035, restriction of A081676 to the primes.
The opposite is A378250, union of A378249 (run-lengths A378251).
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.
A080769 counts primes between perfect powers.
A377283 ranks perfect powers between primes, differences A378356.
A377432 counts perfect powers between primes, see A377434, A377436, A377466.
Cf. A000015, A000961, A052410, A067871, A076412, A131605, A188951, A216765, A345531, A377057, A377468, A378355.

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

A379156 Positions in A246655 (prime powers) of terms q such that there is no prime between q and the next prime power.

Original entry on oeis.org

6, 14, 41, 359, 3589

Offset: 1

Gus Wiseman, Dec 22 2024

The powers of primes themselves are 8, 25, 121, 2187, 32761, ... (A068315).

The prime powers themselves are A068315, for just one prime A379157.
For perfect powers instead of prime powers we have A274605.
Positions of 0 in A366835.
For just one prime we have A379155, for perfect powers A378368.
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.
A065514 gives the greatest prime power < prime(n), difference A377289.
A131605 finds perfect powers that are not prime powers.
A246655 lists the prime powers.
A366833 counts prime powers between primes, see A053607, A304521.
Cf. A025474, A046933, A067871, A080101, A080769, A175106, A178700, A345531, A377281, A377287, A377432, A377434.

v=Select[Range[100],PrimePowerQ];
Select[Range[Length[v]-1],FreeQ[Range[v[[#]],v[[#+1]]],_?PrimeQ]&]

A246655(a(n)) = A068315(n).

A377702 Perfect-powers except for powers of 2.

Original entry on oeis.org

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, 1089, 1156, 1225, 1296, 1331, 1369, 1444, 1521, 1600, 1681, 1728, 1764, 1849, 1936, 2025, 2116, 2187, 2197

Offset: 1

Gus Wiseman, Nov 05 2024

nonn

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

			The terms together with their prime indices begin:
     9: {2,2}
    25: {3,3}
    27: {2,2,2}
    36: {1,1,2,2}
    49: {4,4}
    81: {2,2,2,2}
   100: {1,1,3,3}
   121: {5,5}
   125: {3,3,3}
   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}
   289: {7,7}
   324: {1,1,2,2,2,2}

Including the powers of 2 gives A001597, counted by A377435.
For prime-powers we have A061345.
These terms are counted by A377467, for non-perfect-powers A377701.
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.
A081676 gives the greatest perfect-power <= n.
A131605 lists perfect-powers that are not prime-powers.
A188951 counts perfect-powers less than 2^n.
A377468 gives the least perfect-power > n.
Cf. A014210, A014234, A023055, A045542, A052410, A065514, A069623, A246655, A304521, A336416, A345531, A366833.

Select[Range[1000],GCD@@FactorInteger[#][[All,2]]>1&&!IntegerQ[Log[2,#]]&]

from sympy import mobius, integer_nthroot
def A377702(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(n-2+x+(l:=x.bit_length())+sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,l)))
    return bisection(f,n+1,n+1) # Chai Wah Wu, Nov 06 2024

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

Formula

Extensions

A378357 Distance from n to the least non perfect power >= n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Python

Formula

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A378358 Least non-perfect-power >= n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Python

Python

Formula

A378363 Greatest number <= n that is 1 or not a perfect-power.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Python

A368682 Products of primorials that are perfect powers but not prime powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A375397 Numbers divisible by the square of some prime factor other than the least. Non-hooklike numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A378253 Perfect powers p such that there are no other perfect powers between p and the least prime > p.

Views

Author