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.

Showing 1-7 of 7 results.

A057820 First differences of sequence of consecutive prime powers (A000961).

Original entry on oeis.org

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

Views

Author

Labos Elemer, Nov 08 2000

Keywords

Comments

a(n) = 1 iff A000961(n) = A006549(k) for some k. - Reinhard Zumkeller, Aug 25 2002
Also run lengths of distinct terms in A070198. - Reinhard Zumkeller, Mar 01 2012
Does this sequence contain all positive integers? - Gus Wiseman, Oct 09 2024

Examples

			Odd differences arise in pairs in neighborhoods of powers of 2, like {..,2039,2048,2053,..} gives {..,11,5,..}
		

Crossrefs

For perfect-powers (A001597) we have A053289.
For non-perfect-powers (A007916) we have A375706.
Positions of ones are A375734.
Run-compression is A376308.
Run-lengths are A376309.
Sorted positions of first appearances are A376340.
The second (instead of first) differences are A376596, zeros A376597.
Prime-powers:
- terms: A000961 or A246655, complement A024619
- differences: A057820 (this), first appearances A376341
- anti-runs: A373576, A120430, A006549, A373671
Non-prime-powers:
- terms: A361102
- differences: A375708 (ones A375713)
- anti-runs: A373679, A373575, A255346, A373672

Programs

  • Haskell
    a057820_list = zipWith (-) (tail a000961_list) a000961_list
    -- Reinhard Zumkeller, Mar 01 2012
    
  • Maple
    A057820 := proc(n)
            A000961(n+1)-A000961(n) ;
    end proc: # R. J. Mathar, Sep 23 2016
  • Mathematica
    Map[Length, Split[Table[Apply[LCM, Range[n]], {n, 1, 150}]]] (* Geoffrey Critzer, May 29 2015 *)
    Join[{1},Differences[Select[Range[500],PrimePowerQ]]] (* Harvey P. Dale, Apr 21 2022 *)
  • PARI
    isA000961(n) = (omega(n) == 1 || n == 1)
    n_prev=1;for(n=2,500,if(isA000961(n),print(n-n_prev);n_prev=n)) \\ Michael B. Porter, Oct 30 2009
    
  • Python
    from sympy import primepi, integer_nthroot
    def A057820(n):
        def f(x): return int(n+x-1-sum(primepi(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length())))
        m, k = n, f(n)
        while m != k: m, k = k, f(k)
        r, k = m, f(m)+1
        while r != k: r, k = k, f(k)+1
        return r-m # Chai Wah Wu, Sep 12 2024

Formula

a(n) = A000961(n+1) - A000961(n).

Extensions

Offset corrected and b-file adjusted by Reinhard Zumkeller, Mar 03 2012

A376599 Second differences of consecutive non-prime-powers inclusive (A024619). First differences of A375735.

Original entry on oeis.org

-2, 0, -1, 2, -1, -1, 0, 1, 0, 0, 0, 1, -2, 0, 0, 1, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 1, 0, -1, 1, -1, 1, -1, 0, 1, 0, -1, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 1, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 0, 0, 1, -1, 0
Offset: 1

Views

Author

Gus Wiseman, Oct 02 2024

Keywords

Comments

Inclusive means 1 is a prime-power but not a non-prime-power. For the exclusive version, shift left once.

Examples

			The non-prime-powers inclusive (A024619) are:
  6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30, 33, 34, 35, 36, 38, 39, 40, ...
with first differences (A375735):
  4, 2, 2, 1, 3, 2, 1, 1, 2, 2, 2, 2, 3, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, ...
with first differences (A376599):
  -2, 0, -1, 2, -1, -1, 0, 1, 0, 0, 0, 1, -2, 0, 0, 1, -1, 0, 1, 0, -1, 0, 1, 0, ...
		

Crossrefs

The version for A000002 is A376604, first differences of A054354.
For first differences we had A375735, ones A375713(n) - 1.
Positions of zeros are A376600, complement A376601.
A000961 lists prime-powers inclusive, exclusive A246655.
A007916 lists non-perfect-powers.
A057820 gives first differences of prime-powers inclusive, first appearances A376341, sorted A376340.
A321346/A321378 count integer partitions without prime-powers, factorizations A322452.
For non-prime-powers: A024619/A361102 (terms), A375735/A375708 (first differences), A376600 (inflections and undulations), A376601 (nonzero curvature).
For second differences: A036263 (prime), A073445 (composite), A376559 (perfect-power), A376562 (non-perfect-power), A376590 (squarefree), A376593 (nonsquarefree), A376596 (prime-power).

Programs

  • Mathematica
    Differences[Select[Range[100],!(#==1||PrimePowerQ[#])&],2]
  • Python
    from sympy import primepi, integer_nthroot
    def A376599(n):
        def iterfun(f,n=0):
            m, k = n, f(n)
            while m != k: m, k = k, f(k)
            return m
        def f(x): return int(n+1+sum(primepi(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length())))
        return (a:=iterfun(f,n))-((b:=iterfun(lambda x:f(x)+1,a))<<1)+iterfun(lambda x:f(x)+2,b) # Chai Wah Wu, Oct 02 2024

A376311 Position of first appearance of n in the sequence of first differences of squarefree numbers, or the sequence ends if there is none.

Original entry on oeis.org

1, 3, 6, 31, 150, 515, 13391, 131964, 664313, 5392318, 159468672, 134453711, 28728014494, 50131235121, 634347950217, 48136136076258, 1954623227727573, 14433681032814706, 76465679305346797
Offset: 1

Views

Author

Gus Wiseman, Sep 22 2024

Keywords

Examples

			The sequence of squarefree numbers (A005117) is:
  1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 30, ...
The sequence of first differences (A076259) of squarefree numbers is:
  1, 1, 2, 1, 1, 3, 1, 2, 1, 1, 2, 2, 2, 1, 1, 3, 3, 1, 1, 2, 1, 1, 2, 1, ...
The positions of first appearances are a(n).
		

Crossrefs

This is the position of first appearance of n in A076259, ones A375927.
For compression instead of positions of first appearances we have A376305.
For run-lengths instead of first appearances we have A376306.
For run-sums instead of first appearances we have A376307.
For prime-powers instead of squarefree numbers we have A376341.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A003242 counts compressed compositions, ranks A333489.
A005117 lists squarefree numbers, differences A076259.
A013929 lists nonsquarefree numbers, differences A078147.
A116861 counts partitions by compressed sum, by compressed length A116608.

Programs

  • Mathematica
    mnrm[s_]:=If[Min@@s==1,mnrm[DeleteCases[s-1,0]]+1,0];
    q=Differences[Select[Range[10000],SquareFreeQ]];
    Table[Position[q,k][[1,1]],{k,mnrm[q]}]

Extensions

a(11)-a(19) from Amiram Eldar, Sep 24 2024

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

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

A376309 Run-lengths of the sequence of first differences of prime-powers.

Original entry on oeis.org

3, 1, 2, 2, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 2, 1, 1, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Sep 22 2024

Keywords

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, ...
with runs:
  (1,1,1),(2),(1,1),(2,2),(3),(1),(2),(4),(2,2,2,2),(1),(5),(4),(2),(4), ...
with lengths A376309 (this sequence).
		

Crossrefs

For runs of prime-powers increasing by one we have A174965.
For primes instead of prime-powers we have A333254.
For squarefree numbers instead of prime-powers we have A376306.
For compression instead of run-lengths we have A376308.
For run-sums instead of run-lengths 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, 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.
A124767 counts runs in standard compositions, anti-runs A333381.
A238130, A238279, A333755 count compositions by number of runs.
A274174 counts contiguous compositions, ranks A374249.
A373948 encodes compression using compositions in standard order.

Programs

  • Mathematica
    Length/@Split[Differences[Select[Range[100],PrimePowerQ]]]
  • PARI
    up_to = 20000;
    A376309list(up_to) = { my(v=vector(up_to), ppp=2, pd=1, d, rl=0, k=2, i=0); while(i<#v, k++; if(isprimepower(k), d = k-ppp; ppp = k; if(d == pd, rl++, i++; v[i] = rl; rl = 1; pd = d))); (v); };
    v376309 = A376309list(up_to);
    A376309(n) = v376309[n]; \\ Antti Karttunen, Jan 18 2025

Extensions

More terms from Antti Karttunen, Jan 18 2025

A376310 Run-sums of the sequence of first differences of prime-powers.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 22 2024

Keywords

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, ...
with runs:
  (1,1,1),(2),(1,1),(2,2),(3),(1),(2),(4),(2,2,2,2),(1),(5),(4),(2),(4), ...
with sums A376310 (this sequence).
		

Crossrefs

For primes instead of prime-powers we have A373822, halved A373823.
For squarefree numbers instead of prime-powers we have A376307.
For compression instead of run-sums we have A376308.
For run-lengths instead of run-sums we have A376309.
For positions of first appearances we have A376341, sorted A376340.
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.
A124767 counts runs in standard compositions, anti-runs A333381.
A238130, A238279, A333755 count compositions by number of runs.
A274174 counts contiguous compositions, ranks A374249.
A373948 encodes compression using compositions in standard order.

Programs

  • Mathematica
    Total/@Split[Differences[Select[Range[100],PrimePowerQ]]]
Showing 1-7 of 7 results.