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.

Previous Showing 61-70 of 244 results. Next

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]&]

A373678 Sums of maximal runs of non-prime-powers.

Original entry on oeis.org

1, 6, 10, 12, 29, 18, 63, 24, 26, 28, 30, 138, 117, 42, 135, 48, 153, 280, 60, 125, 131, 207, 72, 380, 80, 82, 430, 651, 297, 102, 315, 108, 333, 819, 369, 126, 259, 670, 138, 1296, 150, 770, 800, 495, 168, 513, 880, 180, 1674, 192, 585, 198, 2255, 2387, 675
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 A361102) is an interval of positions at which consecutive terms differ by one.

Examples

			The maximal runs of non-powers of primes 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

A000040 lists the primes, differences A001223.
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).
See link for composite, prime, nonsquarefree, and squarefree runs.
Prime-power runs: A373675, min A373673, max A373674, length A174965.
Non-prime-power runs: A373678, min A373676, max A373677, length A110969.
Prime-power antiruns: A373576, min A120430, max A006549, length A373671.
Non-prime-power antiruns: A373679, min A373575, max A255346, length A373672.

Programs

  • Mathematica
    Total/@Split[Select[Range[100],!PrimePowerQ[#]&],#1+1==#2&]//Most

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

A377286 Numbers k such that there are no prime-powers between prime(k)+1 and prime(k+1)-1.

Original entry on oeis.org

1, 3, 5, 7, 8, 10, 12, 13, 14, 16, 17, 19, 20, 21, 23, 24, 25, 26, 27, 28, 29, 32, 33, 34, 35, 36, 37, 38, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 55, 56, 57, 58, 59, 60, 62, 63, 64, 65, 66, 67, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82
Offset: 1

Views

Author

Gus Wiseman, Oct 25 2024

Keywords

Examples

			Primes 18 and 19 are 61 and 67, and the interval (62, 63, 64, 65, 66) contains the prime-power 64, so 18 is not in the sequence.
		

Crossrefs

The interval from A008864(n) to A006093(n+1) has A046933(n) elements.
For powers of 2 instead of primes see A013597, A014210, A014234, A244508, A304521.
The nearest prime-power before prime(n)-1 is A065514, difference A377289.
These are the positions of 0 in A080101, or 1 in A366833.
The nearest prime-power after prime(n)+1 is A345531, difference A377281.
For at least one prime-power we have A377057.
For one instead of no prime-powers we have A377287.
For two instead of no prime-powers we have A377288.
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.

Programs

  • Mathematica
    Select[Range[100], Length[Select[Range[Prime[#]+1,Prime[#+1]-1],PrimePowerQ]]==0&]
  • Python
    from itertools import count, islice
    from sympy import factorint, nextprime
    def A377286_gen(): # generator of terms
        p, q, k = 2, 3, 1
        for k in count(1):
            if all(len(factorint(i))>1 for i in range(p+1,q)):
                yield k
            p, q = q, nextprime(q)
    A377286_list = list(islice(A377286_gen(),66)) # Chai Wah Wu, Oct 27 2024

A097318 Numbers with more than one prime factor and, in the ordered factorization, the exponent never increases when read from left to right.

Original entry on oeis.org

6, 10, 12, 14, 15, 20, 21, 22, 24, 26, 28, 30, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 51, 52, 55, 56, 57, 58, 60, 62, 63, 65, 66, 68, 69, 70, 72, 74, 76, 77, 78, 80, 82, 84, 85, 86, 87, 88, 91, 92, 93, 94, 95, 96, 99, 100, 102, 104, 105, 106, 110, 111, 112, 114
Offset: 1

Views

Author

Ralf Stephan, Aug 04 2004

Keywords

Comments

If n = Product_{k=1..m} p(k)^e(k), then m > 1, e(1) >= e(2) >= ... >= e(m).
These are numbers whose ordered prime signature is weakly decreasing. Weakly increasing is A304678. Ordered prime signature is A124010. - Gus Wiseman, Nov 10 2019

Examples

			60 is 2^2*3^1*5^1, A001221(60)=3 and 2>=1>=1, so 60 is in sequence.
		

Crossrefs

Programs

  • Maple
    q:= n-> (l-> (t-> t>1 and andmap(i-> l[i, 2]>=l[i+1, 2],
            [$1..t-1]))(nops(l)))(sort(ifactors(n)[2])):
    select(q, [$1..120])[];  # Alois P. Heinz, Nov 11 2019
  • Mathematica
    fQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Length[f] > 1 && Max[Differences[f]] <= 0]; Select[Range[2, 200], fQ] (* T. D. Noe, Nov 04 2013 *)
  • PARI
    for(n=1, 130, F=factor(n); t=0; s=matsize(F)[1]; if(s>1, for(k=1, s-1, if(F[k, 2]
    				

A351201 Numbers whose multiset of prime factors has a permutation without all distinct runs.

Original entry on oeis.org

12, 18, 20, 28, 36, 44, 45, 48, 50, 52, 60, 63, 68, 72, 75, 76, 80, 84, 90, 92, 98, 99, 100, 108, 112, 116, 117, 120, 124, 126, 132, 140, 144, 147, 148, 150, 153, 156, 162, 164, 168, 171, 172, 175, 176, 180, 188, 192, 196, 198, 200, 204, 207, 208, 212, 216
Offset: 1

Views

Author

Gus Wiseman, Feb 12 2022

Keywords

Examples

			The prime factors of 80 are {2,2,2,2,5} and the permutation (2,2,5,2,2) has runs (2,2), (5), and (2,2), which are not all distinct, so 80 is in the sequence. On the other hand, 24 has prime factors {2,2,2,3}, and all four permutations (3,2,2,2), (2,3,2,2), (2,2,3,2), (2,2,2,3) have distinct runs, so 24 is not in the sequence.
The terms and their prime indices begin:
     12: (2,1,1)         76: (8,1,1)        132: (5,2,1,1)
     18: (2,2,1)         80: (3,1,1,1,1)    140: (4,3,1,1)
     20: (3,1,1)         84: (4,2,1,1)      144: (2,2,1,1,1,1)
     28: (4,1,1)         90: (3,2,2,1)      147: (4,4,2)
     36: (2,2,1,1)       92: (9,1,1)        148: (12,1,1)
     44: (5,1,1)         98: (4,4,1)        150: (3,3,2,1)
     45: (3,2,2)         99: (5,2,2)        153: (7,2,2)
     48: (2,1,1,1,1)    100: (3,3,1,1)      156: (6,2,1,1)
     50: (3,3,1)        108: (2,2,2,1,1)    162: (2,2,2,2,1)
     52: (6,1,1)        112: (4,1,1,1,1)    164: (13,1,1)
     60: (3,2,1,1)      116: (10,1,1)       168: (4,2,1,1,1)
     63: (4,2,2)        117: (6,2,2)        171: (8,2,2)
     68: (7,1,1)        120: (3,2,1,1,1)    172: (14,1,1)
     72: (2,2,1,1,1)    124: (11,1,1)       175: (4,3,3)
     75: (3,3,2)        126: (4,2,2,1)      176: (5,1,1,1,1)
		

Crossrefs

The version for run-lengths instead of runs is A024619.
These permutations are counted by A351202.
These rank the partitions counted by A351203, complement A351204.
A005811 counts runs in binary expansion.
A044813 lists numbers whose binary expansion has distinct run-lengths.
A056239 adds up prime indices, row sums of A112798.
A283353 counts normal multisets with a permutation w/o all distinct runs.
A297770 counts distinct runs in binary expansion.
A333489 ranks anti-runs, complement A348612.
A351014 counts distinct runs in standard compositions, firsts A351015.
A351291 ranks compositions without all distinct runs.
Counting words with all distinct runs:
- A351013 = compositions, for run-lengths A329739, ranked by A351290.
- A351016 = binary words, for run-lengths A351017.
- A351018 = binary expansions, for run-lengths A032020, ranked by A175413.
- A351200 = patterns, for run-lengths A351292.

Programs

  • Mathematica
    Select[Range[100],Select[Permutations[Join@@ ConstantArray@@@FactorInteger[#]],!UnsameQ@@Split[#]&]!={}&]

A067871 Number of primes between consecutive terms of A246547 (prime powers p^k, k >= 2).

Original entry on oeis.org

2, 0, 2, 3, 0, 2, 4, 3, 4, 8, 0, 1, 8, 14, 1, 7, 7, 4, 25, 2, 15, 15, 17, 16, 10, 45, 2, 44, 20, 26, 18, 0, 2, 28, 52, 36, 42, 32, 45, 45, 47, 19, 30, 106, 36, 35, 4, 114, 28, 135, 89, 42, 87, 42, 34, 66, 192, 106, 56, 23, 39, 37, 165, 49, 37, 262, 58, 160, 22
Offset: 1

Views

Author

Jon Perry, Mar 07 2002

Keywords

Comments

Does this sequence have any terms appearing infinitely often? In particular, are {2, 5, 11, 32, 77} the only zeros? As an example, {121, 122, 123, 124, 125} is an interval containing no primes, corresponding to a(11) = 0. - Gus Wiseman, Dec 02 2024

Examples

			The first few prime powers A246547 are 4, 8, 9, 16. The first few primes are 2, 3, 5, 7, 11, 13. We have (4), 5, 7, (8), (9), 11, 13, (16) and so the sequence begins with 2, 0, 2.
The initial terms count the following sets of primes: {5,7}, {}, {11,13}, {17,19,23}, {}, {29,31}, {37,41,43,47}, ... - _Gus Wiseman_, Dec 02 2024
		

Crossrefs

For primes between nonsquarefree numbers we have A236575.
For composite instead of prime we have A378456.
A000015 gives the least prime power >= n.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A080101 counts prime powers between primes.
A246547 lists the non prime prime powers, differences A053707.
A246655 lists the prime powers not including 1, complement A361102.

Programs

  • Mathematica
    t = {}; cnt = 0; Do[If[PrimePowerQ[n], If[FactorInteger[n][[1, 2]] == 1, cnt++, AppendTo[t, cnt]; cnt = 0]], {n, 4 + 1, 30000}]; t (* T. D. Noe, May 21 2013 *)
    nn = 2^20; Differences@ Map[PrimePi, Select[Union@ Flatten@ Table[a^2*b^3, {b, nn^(1/3)}, {a, Sqrt[nn/b^3]}], PrimePowerQ] ] (* Michael De Vlieger, Oct 26 2023 *)

Formula

a(n) = A000720(A025475(n+3)) - A000720(A025475(n+2)). - David Wasserman, Dec 20 2002

Extensions

More terms from David Wasserman, Dec 20 2002
Definition clarified by N. J. A. Sloane, Oct 27 2023

A085970 Number of integers ranging from 2 to n that are not prime-powers.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 4, 5, 5, 5, 6, 6, 7, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 13, 14, 15, 16, 17, 17, 18, 19, 20, 20, 21, 21, 22, 23, 24, 24, 25, 25, 26, 27, 28, 28, 29, 30, 31, 32, 33, 33, 34, 34, 35, 36, 36, 37, 38, 38, 39, 40, 41, 41, 42, 42, 43
Offset: 1

Views

Author

Reinhard Zumkeller, Jul 06 2003

Keywords

Comments

For n > 2, a(n) gives the number of duplicate eliminations performed by the Sieve of Eratosthenes when sieving the interval [2, n]. - Felix Fröhlich, Dec 10 2016
Number of terms of A024619 <= n. - Felix Fröhlich, Dec 10 2016
First differs from A082997 at n = 30. - Gus Wiseman, Jul 28 2022

Examples

			The a(30) = 13 numbers: 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30. - _Gus Wiseman_, Jul 28 2022
		

Crossrefs

The complement is counted by A065515, without 1's A025528.
For primes instead of prime-powers we have A065855, with 1's A062298.
Partial sums of A143731.
The version not treating 1 as a prime-power is A356068.
A000688 counts factorizations into prime-powers.
A001222 counts prime-power divisors.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.

Programs

  • Mathematica
    With[{nn = 75}, Table[n - Count[#, k_ /; k < n] - 1, {n, nn}] &@ Join[{1}, Select[Range@ nn, PrimePowerQ]]] (* Michael De Vlieger, Dec 11 2016 *)
  • PARI
    a(n) = my(i=0); forcomposite(c=4, n, if(!isprimepower(c), i++)); i \\ Felix Fröhlich, Dec 10 2016
    
  • Python
    from sympy import primepi, integer_nthroot
    def A085970(n): return n-1-sum(primepi(integer_nthroot(n,k)[0]) for k in range(1,n.bit_length())) # Chai Wah Wu, Aug 20 2024

Formula

a(n) = Max{A024619(k)<=n} k;
a(n) = n - A065515(n) = A085972(n) - A000720(n).

Extensions

Name modified by Gus Wiseman, Jul 28 2022. Normally 1 is not considered a prime-power, cf. A000961, A246655.

A336548 Numbers k such that at least one pair sigma(p_i^e_i), sigma(p_j^e_j) [with i != j] share a prime factor, when k = p_1^e_1 * ... * p_h^e_h, where each p_i^e_i is the maximal power of prime p_i dividing k.

Original entry on oeis.org

10, 15, 21, 22, 30, 33, 34, 35, 39, 40, 42, 46, 51, 52, 55, 57, 58, 60, 65, 66, 69, 70, 77, 78, 82, 84, 85, 87, 88, 90, 91, 93, 94, 95, 98, 102, 105, 106, 110, 111, 114, 115, 118, 119, 120, 123, 129, 130, 132, 133, 135, 136, 138, 140, 141, 142, 143, 145, 152, 154, 155, 156, 159, 160, 161, 164, 165, 166, 168, 170
Offset: 1

Views

Author

Antti Karttunen, Jul 25 2020

Keywords

Comments

Numbers k for which A353802(k) = Product_{p^e||k} A051027(p^e) > A051027(k), i.e. numbers at which points A051027 is not multiplicative. The notation p^e||k means that p^e divides k, but p^(1+e) does not.
If x is present, then also multiples y*x are present for all y for which gcd(x,y) = 1.
Also numbers at which points A062401 and A353750 are not multiplicative. - Antti Karttunen, May 09 2022

Examples

			10 = 2*5 is present as sigma(2) = 3 and sigma(5) = 6, and 3 and 6 share a prime factor (gcd(3,6) = 3). Also we see that sigma(sigma(2))*sigma(sigma(5)) = 4*12 = 48 > sigma(sigma(10)) = 39.
		

Crossrefs

Cf. A336357, A336558, A336560, A336561, A353807 (subsequences).
Positions of nonzero terms in A336562, in A353753 and in A353803.
Positions of terms larger than 1 in A353755, in A353784 and in A353806.
Subsequence of A024619.

Programs

Formula

{k | A336562(k) > 0}. - Antti Karttunen, May 09 2022

Extensions

The old definition moved to comments and replaced with a more generic, but equivalent definition by Antti Karttunen, May 09 2022

A355432 a(n) = number of k < n such that rad(k) = rad(n) and k does not divide n, where rad(k) = A007947(k).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 2, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 4, 0, 2, 0, 2, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Michael De Vlieger, Feb 22 2023

Keywords

Comments

a(n) = 0 for prime powers and squarefree numbers.

Examples

			a(1) = 18, since 18/6 >= 3. We note that rad(12) = rad(18) = 6, yet 12 does not divide 18.
a(2) = 24, since 24/6 >= 3. rad(18) = rad(24) = 6 and 24 mod 18 = 6.
a(3) = 36, since 36/6 >= 3. rad(24) = rad(36) = 6 and 36 mod 24 = 12.
a(6) = 54, since 54/6 >= 3. m in {12, 24, 36, 48} are such that rad(m) = rad(54) = 6, but none divides 54, etc.
		

Crossrefs

Programs

  • Mathematica
    rad[n_] := rad[n] = Times @@ FactorInteger[n][[All, 1]]; Table[Which[PrimePowerQ[n], 0, SquareFreeQ[n], 0, True, r = rad[n]; Count[Select[Range[n], Nor[PrimePowerQ[#], SquareFreeQ[#]] &], _?(And[rad[#] == r, Mod[n, #] != 0] &)]], {n, 120}]
  • PARI
    rad(n) = factorback(factorint(n)[, 1]); \\ A007947
    a(n) = my(rn=rad(n)); sum(k=1, n-1, if (n % k, rad(k)==rn)); \\ Michel Marcus, Feb 23 2023

Formula

a(n) > 0 for n in A360768.
a(n) < A243822(n) < A010846(n).
a(n) = A008479(n) - A005361(n). - Amiram Eldar, Oct 25 2024
Previous Showing 61-70 of 244 results. Next