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 71-80 of 244 results. Next

A360768 Numbers k that are neither prime powers nor squarefree, such that k/rad(k) >= q, where rad(k) = A007947(k) and prime q = A119288(k).

Original entry on oeis.org

18, 24, 36, 48, 50, 54, 72, 75, 80, 90, 96, 98, 100, 108, 112, 120, 126, 135, 144, 147, 150, 160, 162, 168, 180, 189, 192, 196, 198, 200, 216, 224, 225, 234, 240, 242, 245, 250, 252, 264, 270, 288, 294, 300, 306, 312, 320, 324, 336, 338, 342, 350, 352, 360, 363, 375, 378, 384, 392, 396, 400, 405, 408
Offset: 1

Views

Author

Michael De Vlieger, Feb 22 2023

Keywords

Comments

Proper subsequence of A126706.
Numbers k such that there exists j such that 1 < j < k and rad(j) = rad(k), but j does not divide k.

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. Note: rad(18) = rad(24) = 6 and 24 mod 18 = 6.
a(3) = 36, since 36/6 >= 3. Note: rad(24) = rad(36) = 6 and 36 mod 24 = 12.
a(6) = 54, since 54/6 >= 3. Note: m in {12, 24, 36, 48} are such that rad(m) = rad(54) = 6, but none divides 54, etc.
		

Crossrefs

Programs

  • Mathematica
    Select[Select[Range[120], Nor[SquareFreeQ[#], PrimePowerQ[#]] &], #1/#2 >= #3 & @@ {#1, Times @@ #2, #2[[2]]} & @@ {#, FactorInteger[#][[All, 1]]} &]

Formula

This sequence is { k in A126706 : k/A007947(k) >= A119288(k) }.

A376340 Sorted positions of first appearances in A057820, the sequence of first differences of prime-powers.

Original entry on oeis.org

1, 4, 9, 12, 18, 24, 34, 47, 60, 79, 117, 178, 198, 206, 215, 244, 311, 402, 465, 614, 782, 1078, 1109, 1234, 1890, 1939, 1961, 2256, 2290, 3149, 3377, 3460, 3502, 3722, 3871, 4604, 4694, 6634, 8073, 8131, 8793, 12370, 12661, 14482, 14990, 15912, 17140, 19166
Offset: 1

Views

Author

Gus Wiseman, Sep 22 2024

Keywords

Examples

			The terms together with their prime indices begin:
     1: {}
     4: {1,1}
     9: {2,2}
    12: {1,1,2}
    18: {1,2,2}
    24: {1,1,1,2}
    34: {1,7}
    47: {15}
    60: {1,1,2,3}
    79: {22}
   117: {2,2,6}
   178: {1,24}
   198: {1,2,2,5}
   206: {1,27}
   215: {3,14}
   244: {1,1,18}
		

Crossrefs

For compression instead of sorted firsts we have A376308.
For run-lengths instead of sorted firsts we have A376309.
For run-sums instead of sorted firsts we have A376310.
The version for squarefree numbers is the unsorted version of A376311.
The unsorted version is A376341.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, first differences A057820.
A003242 counts compressed compositions, ranks A333489.
A005117 lists squarefree numbers, differences A076259.
A024619 and A361102 list non-prime-powers, first differences A375708.
A116861 counts partitions by compressed sum, by compressed length A116608.

Programs

  • Mathematica
    q=Differences[Select[Range[100],PrimePowerQ]];
    Select[Range[Length[q]],!MemberQ[Take[q,#-1],q[[#]]]&]

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

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 07 2024

Keywords

Comments

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

Crossrefs

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

Programs

  • Mathematica
    Differences[Table[NestWhile[#+1&, Prime[n]+1,!PrimePowerQ[#]&],{n,100}]]
  • Python
    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

Views

Author

Gus Wiseman, Oct 25 2024

Keywords

Comments

Is this sequence finite? For this conjecture see A053706, A080101, A366833.
Any further terms are > 10^12. - Lucas A. Brown, Nov 08 2024

Examples

			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.
		

Crossrefs

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.

Programs

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

Formula

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

A303606 Powers of composite squarefree numbers that are not squarefree.

Original entry on oeis.org

36, 100, 196, 216, 225, 441, 484, 676, 900, 1000, 1089, 1156, 1225, 1296, 1444, 1521, 1764, 2116, 2601, 2744, 3025, 3249, 3364, 3375, 3844, 4225, 4356, 4761, 4900, 5476, 5929, 6084, 6724, 7225, 7396, 7569, 7776, 8281, 8649, 8836, 9025, 9261, 10000, 10404, 10648, 11025, 11236
Offset: 1

Views

Author

Ilya Gutkovskiy, Apr 26 2018

Keywords

Examples

			196 is in the sequence because 196 = 2^2*7^2.
4900 is in the sequence because 4900 = 2^2*5^2*7^2.
		

Crossrefs

Intersection of A024619 and A072777.
Intersection of A072774 and A126706.
Intersection of A013929 and A182853.

Programs

  • Mathematica
    Select[Range[12000], Length[Union[FactorInteger[#][[All, 2]]]] == 1 && ! SquareFreeQ[#] && ! PrimePowerQ[#] &]
    seq[max_] := Module[{sp = Select[Range[Floor@Sqrt[max]], SquareFreeQ[#] && PrimeNu[#] > 1 &], s = {}}, Do[s = Join[s, sp[[k]]^Range[2, Floor@Log[sp[[k]], max]]], {k, 1, Length[sp]}]; Union@s]; seq[10^4] (* Amiram Eldar, Feb 12 2021 *)
  • Python
    from math import isqrt
    from sympy import mobius, primepi, integer_nthroot
    def A303606(n):
        def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))-primepi(x))
        def f(x): return n-3+x+(y:=x.bit_length())-sum(g(integer_nthroot(x,k)[0]) for k in range(2,y))
        kmin, kmax = 1,2
        while f(kmax) >= kmax:
            kmax <<= 1
        while True:
            kmid = kmax+kmin>>1
            if f(kmid) < kmid:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmax # Chai Wah Wu, Aug 19 2024

Formula

Sum_{n>=1} 1/a(n) = Sum_{n>=1} 1/((A120944(n)-1)*A120944(n)) = Sum_{k>=2} (zeta(k)/zeta(2*k) - P(k) - 1) = 0.07547719891508850482..., where P(k) is the prime zeta function. - Amiram Eldar, Feb 12 2021

A349593 Square array read by downward diagonals: for n >= 0, k >= 1, T(n,k) is the period of {binomial(N,n) mod k: N in Z}.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 4, 1, 1, 4, 3, 4, 1, 1, 5, 8, 9, 8, 1, 1, 6, 5, 8, 9, 8, 1, 1, 7, 12, 5, 16, 9, 8, 1, 1, 8, 7, 36, 5, 16, 9, 8, 1, 1, 9, 16, 7, 72, 25, 16, 9, 16, 1, 1, 10, 9, 16, 7, 72, 25, 16, 9, 16, 1, 1, 11, 20, 27, 32, 7, 72, 25, 32, 27, 16, 1
Offset: 0

Views

Author

Jianing Song, Nov 27 2021

Keywords

Comments

Since binomial(N,n) is defined for all integers N, there is no need to assume that N >= n.
Let Q(N) = 1 if k | binomial(N,n), 0 otherwise. Then T(n,k) is also the period of {Q(N): N in Z}.
By the formula given below, the n-th row is identical to the (n-1)th row if and only if n is not a power of a prime, i.e., n is in A024619. - Jianing Song, Jul 03 2025

Examples

			Rows 0..10:
  1,  1,  1,  1,  1,   1,  1,  1,  1,   1, ...
  1,  2,  3,  4,  5,   6,  7,  8,  9,  10, ...
  1,  4,  3,  8,  5,  12,  7, 16,  9,  20, ...
  1,  4,  9,  8,  5,  36,  7, 16, 27,  20, ...
  1,  8,  9, 16,  5,  72,  7, 32, 27,  40, ...
  1,  8,  9, 16, 25,  72,  7, 32, 27, 200, ...
  1,  8,  9, 16, 25,  72,  7, 32, 27, 200, ...
  1,  8,  9, 16, 25,  72, 49, 32, 27, 200, ...
  1, 16,  9, 32, 25, 144, 49, 64, 27, 400, ...
  1, 16, 27, 32, 25, 432, 49, 64, 81, 400, ...
  1, 16, 27, 32, 25, 432, 49, 64, 81, 400, ...
Example showing that T(4,4) = 16: for N == 0, 1, ..., 15 (mod 16), binomial(N,4) == {0, 0, 0, 0, 1, 1, 3, 3, 2, 2, 2, 2, 3, 3, 1, 1} (mod 4).
Example showing that T(3,10) = 20: for N == 0, 1, ..., 19 (mod 20), binomial(N,3) == {0, 0, 0, 1, 4, 0, 0, 5, 6, 4, 0, 5, 0, 6, 4, 5, 0, 0, 6, 9} (mod 10).
		

Crossrefs

Cf. A022998 (row n = 2), A385555 (row n = 3), A385556 (row n = 4), A385557 (rows n = 5 and 6), A385558 (row n = 7), A385559 (row n = 8), A385560 (rows n = 9 and 10).
Cf. A062383 (2nd column), A064235 (3rd column if offset 0), A385552 (5th column), A385553 (6th column), A385554 (10th column).
Cf. A349221.

Programs

  • Mathematica
    A349593[n_, k_] := If[n == 0 || k == 1, 1, k*Product[p^Floor[Log[p, n]], {p, FactorInteger[k][[All, 1]]}]];
    Table[A349593[k - 1, n - k + 2], {n, 0, 15}, {k, n + 1}] (* Paolo Xausa, Jul 07 2025 *)
  • PARI
    T(n,k) = if(n==0, 1, my(r=1, f=factor(k)); for(j=1, #f[, 1], my(p=f[j, 1], e=f[j, 2]); r *= p^(logint(n,p)+e)); return(r))

Formula

The n-th row is multiplicative with T(n,p^e) = 1 if n = 0, p^(e+floor(log(n)/log(p))) otherwise. In other words, for n > 0, T(n,k) = k * Product_{prime p|k} p^(floor(log(n)/log(p))). See my pdf file for a proof.

A376308 Run-compression of the sequence of first differences of prime-powers.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 20 2024

Keywords

Comments

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, run-compression removes all parts equal to the part immediately to their left. For example, (1,1,2,2,1) has run-compression (1,2,1).

Examples

			The sequence of prime-powers (A246655) is:
  2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, ...
The sequence of first differences (A057820) of prime-powers is:
  1, 1, 1, 2, 1, 1, 2, 2, 3, 1, 2, 4, 2, 2, 2, 2, 1, 5, 4, 2, 4, 2, 4, ...
The run-compression is A376308 (this sequence).
		

Crossrefs

For primes instead of prime-powers we have A037201, halved A373947.
For squarefree numbers instead of prime-powers we have A376305.
For run-lengths instead of compression we have A376309.
For run-sums instead of compression we have A376310.
For positions of first appearances we have A376341, sorted A376340.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A003242 counts compressed compositions, ranks A333489.
A024619 and A361102 list non-prime-powers, differences A375708.
A116861 counts partitions by compressed sum, by compressed length A116608.
A373948 encodes compression using compositions in standard order.

Programs

  • Mathematica
    First/@Split[Differences[Select[Range[100],PrimePowerQ]]]

A373669 Least k such that the k-th maximal run of non-prime-powers has length n. Position of first appearance of n in A110969, and the sequence ends if there is none.

Original entry on oeis.org

1, 5, 7, 12, 18, 190, 28, 109, 40, 28195574, 53
Offset: 1

Views

Author

Gus Wiseman, Jun 14 2024

Keywords

Comments

A run of a sequence (in this case A361102) is an interval of positions at which consecutive terms differ by one.
Are there only 9 terms?
From David A. Corneth, Jun 14 2024: (Start)
No. a(10) exists.
Between the prime 144115188075855859 and 144115188075855872 = 2^57 there are 12 non-prime-powers so a(12) exists. (End)

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

For composite runs we have A073051, sorted A373400, firsts of A176246.
For squarefree runs we have firsts of A120992.
For prime-powers runs we have firsts of A174965.
For prime runs we have firsts of A251092 or A175632.
For squarefree antiruns we have A373128, firsts of A373127.
For nonsquarefree runs we have A373199, firsts of A053797.
The sorted version is A373670.
For antiruns we have firsts of A373672.
For runs of non-prime-powers:
- length A110969
- min A373676
- max A373677
- sum A373678
A000961 lists the powers of primes (including 1).
A005117 lists the squarefree numbers, first differences A076259.
A013929 lists the nonsquarefree numbers, first differences A078147.
A057820 gives first differences of consecutive prime-powers, gaps A093555.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A361102 lists the non-prime-powers, without 1 A024619.

Programs

  • Mathematica
    q=Length/@Split[Select[Range[10000],!PrimePowerQ[#]&],#1+1==#2&]//Most;
    spna[y_]:=Max@@Select[Range[Length[y]],SubsetQ[y,Range[#1]]&];
    Table[Position[q,k][[1,1]],{k,spna[q]}]

A350841 Heinz numbers of integer partitions with a difference < -1 and a conjugate difference < -1.

Original entry on oeis.org

20, 28, 40, 44, 52, 56, 63, 68, 76, 80, 84, 88, 92, 99, 100, 104, 112, 116, 117, 124, 126, 132, 136, 140, 148, 152, 153, 156, 160, 164, 168, 171, 172, 176, 184, 188, 189, 196, 198, 200, 204, 207, 208, 212, 220, 224, 228, 232, 234, 236, 244, 248, 252, 260, 261
Offset: 1

Views

Author

Gus Wiseman, Jan 26 2022

Keywords

Comments

We define a difference of a partition to be a difference of two adjacent parts.

Examples

			The terms together with their prime indices begin:
   20: (3,1,1)
   28: (4,1,1)
   40: (3,1,1,1)
   44: (5,1,1)
   52: (6,1,1)
   56: (4,1,1,1)
   63: (4,2,2)
   68: (7,1,1)
   76: (8,1,1)
   80: (3,1,1,1,1)
   84: (4,2,1,1)
   88: (5,1,1,1)
   92: (9,1,1)
   99: (5,2,2)
		

Crossrefs

Heinz number rankings are in parentheses below.
Taking just one condition gives (A073492) and (A065201), counted by A239955.
These partitions are counted by A350839.
A000041 = integer partitions, strict A000009.
A034296 = partitions with no gaps (A073491), strict A001227 (A073485).
A090858 = partitions with a single gap of size 1 (A325284).
A116931 = partitions with no successions (A319630), strict A003114.
A116932 = partitions with no successions or gaps of size 1, strict A025157.
A350842 = partitions with no gaps of size 1, strict A350844, sets A005314.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Select[Range[100],(Min@@Differences[Reverse[primeMS[#]]]<-1)&&(Min@@Differences[conj[primeMS[#]]]<-1)&]

A080765 Integers m such that m+1 divides lcm(1 through m).

Original entry on oeis.org

5, 9, 11, 13, 14, 17, 19, 20, 21, 23, 25, 27, 29, 32, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 47, 49, 50, 51, 53, 54, 55, 56, 57, 59, 61, 62, 64, 65, 67, 68, 69, 71, 73, 74, 75, 76, 77, 79, 81, 83, 84, 85, 86, 87, 89, 90, 91, 92, 93, 94, 95, 97, 98, 99, 101, 103, 104, 105, 107
Offset: 1

Views

Author

Lekraj Beedassy, Mar 10 2003

Keywords

Comments

Integers m for which A003418(m) = A003418(m+1).
a(n) = A024619(n) - 1. Proof:
If N+1 is a power of a prime (N+1=P^K), then only smaller powers of that prime divide numbers up to N and so lcm(1..N) doesn't have K powers of P; that is, N+1=P^K doesn't divide lcm(1..N).
From Don Reble, Mar 12 2003: (Start)
If N+1 is not a power of a prime, then it has at least two prime factors. Call one of them P, let K be such that P^K divides N+1, but P^(K+1) doesn't, and let N+1=P^K*R. Then
- R is greater than 1 because it is divisible by another prime factor of N+1;
- P^K and R are each less than N+1 because the other is greater than one;
- lcm(P^K,R) divides lcm(1..N) because 1..N includes both numbers;
- lcm(P^K,R)=N+1 because P doesn't divide R;
- N+1 divides lcm(1..N). (End)

Examples

			17 is the sequence because lcm(1,2,...,17)=12252240 and 17+1=18 divides 12252240.
		

Crossrefs

Cf. A003418.

Programs

  • Mathematica
    Select[Range[120], Divisible[LCM @@ Range[#], #+1]&] (* Jean-François Alcover, Jun 21 2018 *)
  • PARI
    a=1;for(n=1,108,a=lcm(a,n);if(a%(n+1)==0,print1(n,","))) \\ Klaus Brockhaus, Jun 11 2004
    
  • PARI
    first(n) = {my(u = max(2*n, 50), charact = vector(u, i, 1), res = List()); forprime(p = 2, 2*n, for(t = 1, logint(u, p), charact[p^t - 1] = 0)); for(i = 1, u, if(charact[i] == 1, listput(res, i); if(#res >= n, return(res)))); res } \\ David A. Corneth, Aug 30 2019
    
  • Sage
    [x - 1 for x in (1..108) if not is_prime_power(n)]  # Peter Luschny, May 23 2013

Formula

a(n) ~ n. - David A. Corneth, Aug 30 2019

Extensions

More terms from Klaus Brockhaus, Jun 11 2004
Previous Showing 71-80 of 244 results. Next